CN106655202A - Interruptible-load-based multiple-target overloading assistant decision making method - Google Patents

Abstract

The invention discloses an interruptible-load-based multiple-target overloading assistant decision making method. Contribution indexes of an interruptible load for reducing two target functions, namely an interruptible load cutting total amount and an interruptible load cutting total number, are respectively computed; a comprehensive contribution index of the interruptible load to the two target functions is computed by utilizing linear weighting factors; the linear weighting factors are divided into different grades; an overloading assistant decision making measure corresponding to each grade of the linear weighting factor is computed; and the final measure is selected via a scheduling operator. By applying the interruptible-load-based multiple-target overloading assistant decision making method, multiple-target optimization can be performed according to the contribution indexes of the interruptible load to the target functions, an effective solution set meeting the computation precision is given out, and coordinate optimization of reduction of the two target functions, namely the interruptible load cutting total amount and the interruptible load cutting total number, is implemented.

Description

A kind of multiple target based on interruptible load transships aid decision-making method

Technical field

The present invention relates to a kind of multiple target based on interruptible load transships aid decision-making method, belong to power system automatic Change technical field.

Background technology

Receiving end grid generator it is standby deficiency in the case of, need by spatial load forecasting to eliminate electrical network in equipment mistake Carry.Ratio of the extra-high voltage direct-current transmission capacity in receiving-end system scale constantly increases, and direct current locking failure causes high-power punching Hit, in the presence of system frequency modulation characteristic, generator output increase may cause Section Tidal Current of Power Transmission to get over stability limitation, to load The scale and response speed of control is put forward higher requirement.

Traditional spatial load forecasting does not account for the practical situation of user side, and the concentration cutting load of extensive indiscriminate property is not Huge economic losses can be only caused, the acceptance level of user is relatively low under Power Market.Scope is modeled by EMS system Limit, provincial scheduling institution is difficult to directly perceive the controllable burden capacity in electrical network, can only be according to the ratio for pre-setting each Share the cutting load amount for eliminating overload between districts and cities, draw road to carry out spatial load forecasting by accident by regional dispatch mechanism.When certain When regional dispatch mechanism controls amount is not enough, provincial scheduling institution needs the system running state of the execution end that road is drawn according to accident Aid decision calculating is carried out again, it is impossible to meet rapidity requirement of the complex fault scene to accident handling.

The content of the invention

In order to solve above-mentioned technical problem, the invention provides a kind of multiple target overload auxiliary based on interruptible load is determined Plan method.

In order to achieve the above object, the technical solution adopted in the present invention is：

A kind of multiple target based on interruptible load transships aid decision-making method, comprises the following steps,

Step 1, obtains the controllable capacity data of interruptible load and sets up interruptible load and energy at control centre station The corresponding relation of the medium duty value of management system；

Step 2, if control centre station detects there is overload prevention device under current state, into step 3；Otherwise, terminate This method；

Step 3, detects that overload prevention device associates the adjustable space of electromotor, if all of overload prevention device cannot be by increasing Itself association electromotor it is active reducing overload degree, then into step 4；Otherwise, method ends；

Step 4, calculates interruptible load and shifts distribution factor to the power of overload prevention device, and power transfer distribution factor is big Decision space is added in the interruptible load of setting threshold value；

Step 5, shifts distribution factor and calculates interruptible load to reducing interruptible load excision total amount target according to power Function f₁Contribution degree index, total object function f is cut off to interruptible load according to controllable calculation of capacity interruptible load₂'s Contribution degree index, interruptible load is calculated to object function f using the linear weighted function factor₁And f₂Comprehensive contribution degree index；

Step 6, C gear is equally divided into according to the check figure of parallel computing platform by the excursion of the linear weighted function factor； Wherein, the corresponding linear weighted function factors of gear m are κ_m=m/ (C-1), m=0 ..., C-1；

Step 7, is ranked up according to comprehensive contribution degree index to interruptible load, calculate respectively each gear linear weighted function because Son corresponding overload aid decision measure；Overload aid decision measure include needing the interruptible load of excision and with can interrupt The total flux matched generated power Adjusted Option of load excision；

Step 8, calculates interruptible load excision total amount E of each gear_mSum W is cut off with interruptible load_mIf there is phase The difference of adjacent gear interruptible load excision total amount is more than computational accuracy ξ, then execution step 9；Otherwise, execution step 10；

Step 9, cuts off total amount difference true more than the linear weighted function factor region of search of computational accuracy ξ according to interruptible load Fixed newly-increased gear, return to step 7；

Step 10, with interruptible load excision total amount E that gear is 0₀Sum W is cut off with interruptible load₀As reference, Variation delta E that interruptible load cuts off total amount is calculated successively for other effective gear n_n=E₀-E_nCut with interruptible load Except variation delta W of sum_n=W₀-W_n, E_nAnd W_nThe interruptible load excision total amount of respectively gear n and interruptible load excision Sum, by the overload aid decision measure that management and running personnel selection is final.

Determine that overload prevention device associates the process of the adjustable space of electromotor and is,

3-1) electromotor being calculated respectively, distribution factor is shifted to the power of overload prevention device；

3-2) filter out electromotor of the power transfer distribution factor absolute value more than setting threshold value, using the electromotor as Overload prevention device associates electromotor；

3-3) meet overload prevention device association electromotor of the power transfer distribution factor less than 0 exerts oneself adjustment direction to increase Power, calculates and increases the adjustable space that overload prevention device of exerting oneself associates electromotor；

Wherein, Δ P_g,jTo increase the adjustable space that overload prevention device of exerting oneself associates electromotor g,WithRespectively increase and exert oneself Overload prevention device association electromotor g's is maximum active and current active；

If 3-4) all adjustable spaces for increasing overload prevention device association electromotor of exerting oneself are 0, cannot be by increasing Load equipment associates generator output to reduce the overload degree of overload prevention device.

Calculate interruptible load is to the power transfer distribution factor process of overload prevention device：Calculate corresponding with interruptible load Duty value distribution factor is shifted to the power of overload prevention device, the power is shifted distribution factor as interruptible load to mistake The power transfer distribution factor of load equipment.

Interruptible load is calculated to object function f₁And f₂The process of comprehensive contribution degree index be,

5-1) distribution factor is shifted to the power of all overload prevention devices according to interruptible load, calculate the control of interruptible load Performance indications processed；

Wherein, γ_iFor interruptible load L_iControl performance standard, N for overload prevention device sum, λ_jFor overload prevention device Q_j Load factor, γ_i,jFor interruptible load L_iTo overload prevention device Q_jPower transfer distribution factor；

5-2) choose in interruptible load and be worth on the basis of the maximum of Control performance standard, by the controlling of interruptible load Energy index is normalized, and calculates interruptible load to reducing interruptible load excision total amount object function f₁Contribution degree Index；

Wherein,For interruptible load L_iTo reducing interruptible load excision total amount object function f₁Contribution degree index, γ_max=max (γ_i) for the maximum of Control performance standard in interruptible load；

5-3) choose in interruptible load and be worth on the basis of the maximum of controllable capacity, the controllable capacity of interruptible load is entered Row normalized, calculates interruptible load to reducing the total object function f of interruptible load excision₂Contribution degree index；

Wherein,For interruptible load L_iTo reducing the total object function f of interruptible load excision₂Contribution degree index, P_iFor interruptible load L_iControllable capacity, P_max=max (P_i) it is interruptible load L_iIn controllable capacity the maximum；

5-4) interruptible load is calculated to object function f using the linear weighted function factor₁And f₂Comprehensive contribution degree index；

Wherein,For the linear weighted function factor be κ when, interruptible load L_iTo object function f₁And f₂Comprehensive contribution Degree index.

The process of corresponding overload aid decision measure is when the calculating linear weighted function factor is κ,

7-1) AC power flow check number of times A is initialized as into 0, it is ε to arrange AC power flow and check step-length₂；

7-2) determine that overload prevention device overload is completely eliminated corresponding interruptible load resection successively, meet：

Wherein, E_j(κ) it is overload prevention device Q_jOverload is completely eliminated corresponding interruptible load resection, M_jCan for what is cut off Interruptible load number,For overload prevention device Q_jPower more limit the quantity, P_jFor overload prevention device Q_jActive power, For alternating current circuit,It is according to overload prevention device Q_jThe power limit that rated current and current voltage are converted into, for transformation Device,It is according to overload prevention device Q_jThe power limit that rated capacity and current voltage are converted into；For transmission cross-section, For power allowances of the transmission cross-section under the current method of operation, γ_i+1,jFor interruptible load L_i+1To overload prevention device Q_jPower turn Move distribution factor；

7-3) maximum of all overload prevention devices correspondence interruptible load resection is examined as interruptible load resection Consider frequency modulation frequency modulation during frequency modulation frequency modulation association electromotor and interruptible load response condition, to electromotor and do not cut off can in The active of disconnected load is adjusted, and control measure are checked by AC power flow；

7-4) AC power flow check number of times A is updated to into A+1, the power that all overload prevention devices are calculated respectively is more limited the quantity, really The corresponding overload prevention device Q of absolute value of the minima that fixed all power are more limited the quantity_d；

Process according to the following steps：

(1) ifThen by C_AAs final interruptible load control program；

Overload prevention device Q when number of times is A is checked for AC power flow_dPower more limit the quantity, C_ACheck for AC power flow secondary The maximum of number all overload prevention devices correspondence interruptible load resection when being A；

(2) if A>1、AndThen by C_AAs final interruptible load control program；

Overload prevention device Q when number of times is A-1 is checked for AC power flow_dPower more limit the quantity；

(3) if A>1、AndThen by C_A-1As final interruptible load control program；

C_A-1For the maximum that AC power flow checks all overload prevention device correspondence interruptible load resection when number of times is A-1；

(4) ifThen by overload prevention device Q_dPower get over limitation and be updated to Δ P_d=Δ P_d+ε₂, return to step 7-1；

(5) ifThen by overload prevention device Q_dPower get over limitation and be updated to Δ P_d=Δ P_d-ε₂, return to step 7- 1。

In step 7-3, it is to electromotor and the active process that is adjusted for not cutting off interruptible load,

7-3-1) Control performance standard of the computation overload equipment association electromotor to all overload prevention devices；

Wherein, τ_gControl performance standards of the electromotor g to all overload prevention devices, τ are associated for overload prevention device_g,jSet for overload Standby association electromotor g is to overload prevention device Q_jPower transfer distribution factor；

7-3-2) according to system frequency and the frequency departure of rated frequency computing system；

Δ f=f-f₀

Wherein, Δ f is system frequency f and rated frequency f₀The frequency departure of computing system；

7-3-3) compare the corresponding frequency departure threshold value Δ f in AGC operating dead zones_AGCWith the size of Δ f, if Δ f >=-Δ f_AGC, then select overload prevention device association electromotor drop to exert oneself by Control performance standard is descending, adjust total amount and C_AUnanimously, and Terminate adjustment process；Otherwise, simulate frequency modulation frequency modulation course of action frequency modulation frequency modulation is associated electromotor and not cutting off can interrupt it is negative Lotus is adjusted, and step is as follows：

(1) the active variable quantity that frequency modulation frequency modulation associates electromotor is participated at the end of calculating controlling cycle T；

Wherein, Δ P_g′The active variable quantity that frequency modulation frequency modulation associates electromotor g ', V are participated at the end of for controlling cycle T_g′For Frequency modulation frequency modulation associates the climbing rate of electromotor g ',WithRespectively frequency modulation frequency modulation association electromotor g ' it is maximum active with It is current active；

(2) calculate controlling cycle T at the end of system frequency variable quantity；

Wherein, Δ f_TThe variable quantity of system frequency, K at the end of for controlling cycle T_sNot cut off total tune of load in system Frequency coefficient, N_GTo participate in frequency modulation frequency modulation association electromotor quantity；

(3) duty value B after frequency change is calculated_kActive variable quantity；

Wherein,For duty value B after frequency change_kActive variable quantity,For duty value B_kFrequency change Active power coefficient.

Total amount difference is cut off according to interruptible load newly-increased more than the linear weighted function factor region of search determination of computational accuracy ξ The process of gear be,

It is t and t+1 9-1) to define adjacent gear positions, then the corresponding linear weighted function factor is respectively κ_tAnd κ_t+1；

The gear number for 9-2) newly increasingE_t、E_t+1Be divided into gear for t and t+1 can in Disconnected load cuts off total amount；

9-3) by the region of search [κ_t, κ_t+1] press M_tRespectively obtain the linear weighted function factor of newly-increased gear.

The beneficial effect that the present invention is reached：It is total to reducing interruptible load excision that the present invention calculates respectively interruptible load The contribution degree index of amount, total two object functions of interruptible load excision, using the linear weighted function factor interruptible load is calculated Comprehensive contribution degree index to two object functions；The linear weighted function factor is divided into into different gears, each gear line is calculated respectively Property weighter factor it is corresponding overload aid decision measure, by the measure that management and running personnel selection is final；It is given linear calculating During spatial load forecasting scheme under weighter factor, distribution factor estimation interruptible load excision is shifted according to the power of interruptible load Total amount, it is considered to the response condition of frequency modulation frequency modulation association electromotor and interruptible load during frequency modulation frequency modulation, to electromotor and not Excision the active of interruptible load is adjusted, and gets over limitation to the power of overload prevention device using AC power flow check and is modified, Approach by inchmeal meets computational accuracy and the aid decision measure of apparatus overload in electrical network can be completely eliminated；Using present invention side Method, can carry out multiple-objection optimization according to interruptible load to the contribution degree index of object function, be given and meet computational accuracy Efficient solution set, realizes reducing interruptible load excision total amount and the coordination of total two object functions of interruptible load excision is excellent Change.

Description of the drawings

Fig. 1 is the flow chart of the present invention.

Specific embodiment

Below in conjunction with the accompanying drawings the invention will be further described.Following examples are only used for clearly illustrating the present invention Technical scheme, and can not be limited the scope of the invention with this.

As shown in figure 1, a kind of multiple target based on interruptible load transships aid decision-making method, its principle is：Can interrupt Load shifts distribution factor to the power of overload prevention device, reflects the excision interruptible load total to reducing interruptible load excision The contribution of amount, the capacity of interruptible load reflects the tribute for cutting off the interruptible load to reducing interruptible load excision sum Offer；Interruptible load can be calculated to reducing interruptible load excision total amount and interruptible load excision using the linear weighted function factor The comprehensive contribution degree index of sum, the change of the linear weighted function factor reflects different contribution degree indexs in comprehensive contribution degree index Weight, interruptible load is ranked up by comprehensive contribution degree index, then can be answered by solving different linear weighted function factor pairs Interruptible load control program can obtain effective solution of Model for Multi-Objective Optimization；Should given linear weighted function factor pair is calculated Interruptible load control program when, can according to power shift distribution factor estimate interruptible load control problem, by hand over Stream trend check carrys out Approach by inchmeal and meets computational accuracy and the aid decision measure of apparatus overload in electrical network can be completely eliminated； Before AC power flow is checked, it is considered to the response condition of electromotor and interruptible load is associated during frequency modulation frequency modulation, to electromotor and Do not cut off the active of interruptible load to be adjusted, the trend that can simulate system after spatial load forecasting is distributed and guarantees in system Power-balance；The corresponding interruptible load controlling party of multiple gear linear weighted function factors can simultaneously be calculated by parallel computation Case, improves the calculating speed of overload aid decision.

Comprise the following steps that：

Step 1, the collection of user side installation data and interruptible load control device in interruptible load, in control centre Stand and obtain the controllable capacity data of interruptible load and to set up interruptible load corresponding with the medium duty value of EMS Relation.

, if control centre station detects has overload prevention device under current state, that is, there is alternating current circuit, transformator in step 2 Overload or transmission cross-section get over the situation of stability limitation, then into step 3；Otherwise, method ends.

Step 3, detects that overload prevention device associates the adjustable space of electromotor, if all of overload prevention device cannot be by increasing Itself association electromotor it is active reducing overload degree, then into step 4；Otherwise, method ends.

Determine that overload prevention device associates the process of the adjustable space of electromotor and is：

3-1) electromotor being calculated respectively, distribution factor is shifted to the power of overload prevention device；

3-2) filter out electromotor of the power transfer distribution factor absolute value more than setting threshold value, using the electromotor as Overload prevention device associates electromotor；

3-3) meet overload prevention device association electromotor of the power transfer distribution factor less than 0 exerts oneself adjustment direction to increase Power, calculates and increases the adjustable space that overload prevention device of exerting oneself associates electromotor；

Wherein, Δ P_g,jTo increase the adjustable space that overload prevention device of exerting oneself associates electromotor g,WithRespectively increase and exert oneself Overload prevention device association electromotor g's is maximum active and current active；

If 3-4) all adjustable spaces for increasing overload prevention device association electromotor of exerting oneself are 0, cannot be by increasing Load equipment associates generator output to reduce the overload degree of overload prevention device.

Step 4, calculates interruptible load and shifts distribution factor to the power of overload prevention device, and power transfer distribution factor is big Decision space is added in the interruptible load of setting threshold value.

Calculate interruptible load is to the power transfer distribution factor process of overload prevention device：Calculate corresponding with interruptible load Duty value distribution factor is shifted to the power of overload prevention device, the power is shifted distribution factor as interruptible load to mistake The power transfer distribution factor of load equipment.

Step 5, shifts distribution factor and calculates interruptible load to reducing interruptible load excision total amount target according to power Function f₁Contribution degree index, total object function f is cut off to interruptible load according to controllable calculation of capacity interruptible load₂'s Contribution degree index, interruptible load is calculated to object function f using the linear weighted function factor₁And f₂Comprehensive contribution degree index.

Interruptible load is calculated to object function f₁And f₂The process of comprehensive contribution degree index be：

5-1) distribution factor is shifted to the power of all overload prevention devices according to interruptible load, calculate the control of interruptible load Performance indications processed；

Wherein, γ_iFor interruptible load L_iControl performance standard, N for overload prevention device sum, λ_jFor overload prevention device Q_j Load factor, γ_i,jFor interruptible load L_iTo overload prevention device Q_jPower transfer distribution factor；

5-2) choose in interruptible load and be worth on the basis of the maximum of Control performance standard, by the controlling of interruptible load Energy index is normalized, and calculates interruptible load to reducing interruptible load excision total amount object function f₁Contribution degree Index；

Wherein,For interruptible load L_iTo reducing interruptible load excision total amount object function f₁Contribution degree index, γ_max=max (γ_i) for the maximum of Control performance standard in interruptible load；

5-3) choose in interruptible load and be worth on the basis of the maximum of controllable capacity, the controllable capacity of interruptible load is entered Row normalized, calculates interruptible load to reducing the total object function f of interruptible load excision₂Contribution degree index；

Wherein,For interruptible load L_iTo reducing the total object function f of interruptible load excision₂Contribution degree index, P_iFor interruptible load L_iControllable capacity, P_max=max (P_i) it is interruptible load L_iIn controllable capacity the maximum；

5-4) interruptible load is calculated to object function f using the linear weighted function factor₁And f₂Comprehensive contribution degree index；

Wherein,For the linear weighted function factor be κ when, interruptible load L_iTo object function f₁And f₂Comprehensive contribution Degree index, the excursion of κ is [0,1].

Step 6, C gear is equally divided into according to the check figure of parallel computing platform by the excursion of the linear weighted function factor； Wherein, the corresponding linear weighted function factors of gear m are κ_m=m/ (C-1), m=0 ..., C-1.

Step 7, is ranked up according to comprehensive contribution degree index to interruptible load, calculate respectively each gear linear weighted function because Son corresponding overload aid decision measure；Overload aid decision measure include needing the interruptible load of excision and with can interrupt The total flux matched generated power Adjusted Option of load excision.

The process of corresponding overload aid decision measure is when the calculating linear weighted function factor is κ：

7-1) AC power flow check number of times A is initialized as into 0, it is ε to arrange AC power flow and check step-length₂；

7-2) determine that overload prevention device overload is completely eliminated corresponding interruptible load resection successively, meet：

Wherein, E_j(κ) it is overload prevention device Q_jOverload is completely eliminated corresponding interruptible load resection, M_jCan for what is cut off Interruptible load number,For overload prevention device Q_jPower more limit the quantity, P_jFor overload prevention device Q_jActive power, For alternating current circuit,It is according to overload prevention device Q_jThe power limit that rated current and current voltage are converted into, for transformation Device,It is according to overload prevention device Q_jThe power limit that rated capacity and current voltage are converted into；For transmission cross-section, For power allowances of the transmission cross-section under the current method of operation, γ_i+1,jFor interruptible load L_i+1To overload prevention device Q_jPower turn Move distribution factor；

7-3) maximum of all overload prevention devices correspondence interruptible load resection is examined as interruptible load resection Consider frequency modulation frequency modulation during frequency modulation frequency modulation association electromotor and interruptible load response condition, to electromotor and do not cut off can in The active of disconnected load is adjusted, and control measure are checked by AC power flow；

It is to electromotor and the active process that is adjusted for not cutting off interruptible load,

7-3-1) Control performance standard of the computation overload equipment association electromotor to all overload prevention devices；

Wherein, τ_gControl performance standards of the electromotor g to all overload prevention devices, τ are associated for overload prevention device_g,jSet for overload Standby association electromotor g is to overload prevention device Q_jPower transfer distribution factor；

7-3-2) according to system frequency and the frequency departure of rated frequency computing system；

Δ f=f-f₀

Wherein, Δ f is system frequency f and rated frequency f₀The frequency departure of computing system；

7-3-3) compare the corresponding frequency departure threshold value Δ f in AGC operating dead zones_AGCWith the size of Δ f, if Δ f >=-Δ f_AGC, then select overload prevention device association electromotor drop to exert oneself by Control performance standard is descending, adjust total amount and C_AUnanimously, and Terminate adjustment process；Otherwise, simulate frequency modulation frequency modulation course of action frequency modulation frequency modulation is associated electromotor and not cutting off can interrupt it is negative Lotus is adjusted, and step is as follows：

(1) the active variable quantity that frequency modulation frequency modulation associates electromotor is participated at the end of calculating controlling cycle T；

Wherein, Δ P_g′The active variable quantity that frequency modulation frequency modulation associates electromotor g ', V are participated at the end of for controlling cycle T_g′For Frequency modulation frequency modulation associates the climbing rate of electromotor g ',WithRespectively frequency modulation frequency modulation association electromotor g ' it is maximum active with It is current active；

(2) calculate controlling cycle T at the end of system frequency variable quantity；

Wherein, Δ f_TThe variable quantity of system frequency, K at the end of for controlling cycle T_sNot cut off total tune of load in system Frequency coefficient, N_GTo participate in frequency modulation frequency modulation association electromotor quantity；

(3) duty value B after frequency change is calculated_kActive variable quantity；

Wherein,For duty value B after frequency change_kActive variable quantity,For duty value B_kFrequency change Active power coefficient；

7-4) AC power flow check number of times A is updated to into A+1, the power that all overload prevention devices are calculated respectively is more limited the quantity, really The corresponding overload prevention device Q of absolute value of the minima that fixed all power are more limited the quantity_d；

Process according to the following steps：

(1) ifThen by C_AAs final interruptible load control program；

Overload prevention device Q when number of times is A is checked for AC power flow_dPower more limit the quantity, C_ACheck for AC power flow secondary The maximum of number all overload prevention devices correspondence interruptible load resection when being A；

(2) if A>1、AndThen by C_AAs final interruptible load control program；

Overload prevention device Q when number of times is A-1 is checked for AC power flow_dPower more limit the quantity；

(3) if A>1、AndThen by C_A-1As final interruptible load control program；

C_A-1For the maximum that AC power flow checks all overload prevention device correspondence interruptible load resection when number of times is A-1；

(4) ifThen by overload prevention device Q_dPower get over limitation and be updated to Δ P_d=Δ P_d+ε₂, return to step 7- 1；

(5) ifThen by overload prevention device Q_dPower get over limitation and be updated to Δ P_d=Δ P_d-ε₂, return to step 7- 1。

Step 8, calculates interruptible load excision total amount E of each gear_mSum W is cut off with interruptible load_mIf there is phase The difference of adjacent gear interruptible load excision total amount is more than computational accuracy ξ, then execution step 9；Otherwise, execution step 10.

Step 9, cuts off total amount difference true more than the linear weighted function factor region of search of computational accuracy ξ according to interruptible load Fixed newly-increased gear, return to step 7.

Total amount difference is cut off according to interruptible load newly-increased more than the linear weighted function factor region of search determination of computational accuracy ξ The process of gear be：

It is t and t+1 9-1) to define adjacent gear positions, then the corresponding linear weighted function factor is respectively κ_tAnd κ_t+1；

The gear number for 9-2) newly increasingE_t、E_t+1Be divided into gear for t and t+1 can in Disconnected load cuts off total amount；

9-3) by the region of search [κ_t, κ_t+1] press M_tRespectively obtain the linear weighted function factor of newly-increased gear.

Step 10, with interruptible load excision total amount E that gear is 0₀Sum W is cut off with interruptible load₀As reference, Variation delta E that interruptible load cuts off total amount is calculated successively for other effective gear n_n=E₀-E_nCut with interruptible load Except variation delta W of sum_n=W₀-W_n, E_nAnd W_nThe interruptible load excision total amount of respectively gear n and interruptible load excision Sum, by the overload aid decision measure that management and running personnel selection is final.

Said method calculates respectively interruptible load to reducing interruptible load excision total amount, interruptible load excision sum The contribution degree index of two object functions, using the linear weighted function factor comprehensive tribute of the interruptible load to two object functions is calculated Degree of offering index；The linear weighted function factor is divided into into different gears, the corresponding overload of each gear linear weighted function factor is calculated respectively auxiliary Policy making steps are helped, by the measure that management and running personnel selection is final；Spatial load forecasting side in the case where the given linear weighted function factor is calculated During case, distribution factor estimation interruptible load excision total amount is shifted according to the power of interruptible load, it is considered to frequency modulation frequency modulation process The response condition of middle frequency modulation frequency modulation association electromotor and interruptible load, to electromotor and do not cut off interruptible load it is active enter Row adjustment, gets over limitation and is modified using AC power flow check to the power of overload prevention device, and Approach by inchmeal meets computational accuracy simultaneously The aid decision measure of apparatus overload in electrical network can be completely eliminated；Using the inventive method, can be according to interruptible load pair The contribution degree index of object function carries out multiple-objection optimization, provides the efficient solution set for meeting computational accuracy, realizes that reduction can be interrupted Load cuts off the coordination optimization of total amount and total two object functions of interruptible load excision.

The above is only the preferred embodiment of the present invention, it is noted that for the ordinary skill people of the art For member, on the premise of without departing from the technology of the present invention principle, some improvement and deformation can also be made, these improve and deform Also should be regarded as protection scope of the present invention.

Claims (7)

1. a kind of multiple target based on interruptible load transships aid decision-making method, it is characterised in that：Comprise the following steps,

Step 1, obtains the controllable capacity data of interruptible load and sets up interruptible load and energy management at control centre station The corresponding relation of the medium duty value of system；

Step 2, if control centre station detects there is overload prevention device under current state, into step 3；Otherwise, we are terminated Method；

Step 3, detects that overload prevention device associates the adjustable space of electromotor, if all of overload prevention device cannot be by increasing itself Association electromotor it is active reducing overload degree, then into step 4；Otherwise, method ends；

Step 4, calculates interruptible load and shifts distribution factor to the power of overload prevention device, power transfer distribution factor is more than and is set The interruptible load of set threshold value adds decision space；

Step 5, shifts distribution factor and calculates interruptible load to reducing interruptible load excision total amount object function f according to power₁ Contribution degree index, total object function f is cut off to interruptible load according to controllable calculation of capacity interruptible load₂Contribution degree Index, interruptible load is calculated to object function f using the linear weighted function factor₁And f₂Comprehensive contribution degree index；

Step 6, C gear is equally divided into according to the check figure of parallel computing platform by the excursion of the linear weighted function factor；Wherein, The corresponding linear weighted function factors of gear m are κ_m=m/ (C-1), m=0 ..., C-1；

Step 7, is ranked up according to comprehensive contribution degree index to interruptible load, and each gear linear weighted function factor pair is calculated respectively The overload aid decision measure answered；The measure of overload aid decision includes needing the interruptible load and and interruptible load of excision The total flux matched generated power Adjusted Option of excision；

Step 8, calculates interruptible load excision total amount E of each gear_mSum W is cut off with interruptible load_mIf there are adjacent shelves The difference of position interruptible load excision total amount is more than computational accuracy ξ, then execution step 9；Otherwise, execution step 10；

Step 9, cuts off total amount difference and determines newly more than the linear weighted function factor region of search of computational accuracy ξ according to interruptible load The gear of increasing, return to step 7；

Step 10, with interruptible load excision total amount E that gear is 0₀Sum W is cut off with interruptible load₀As reference, for Other effective gear n calculate successively variation delta E that interruptible load cuts off total amount_n=E₀-E_nIt is total with interruptible load excision Several variation deltas W_n=W₀-W_n, E_nAnd W_nThe interruptible load excision total amount of respectively gear n and interruptible load excision are total Number, by the overload aid decision measure that management and running personnel selection is final.

2. a kind of multiple target based on interruptible load according to claim 1 transships aid decision-making method, and its feature exists In：Determine that overload prevention device associates the process of the adjustable space of electromotor and is,

3-1) electromotor being calculated respectively, distribution factor is shifted to the power of overload prevention device；

Electromotor of the power transfer distribution factor absolute value more than setting threshold value 3-2) is filtered out, using the electromotor as overload Equipment associates electromotor；

The adjustment direction of exerting oneself for 3-3) meeting overload prevention device association electromotor of the power transfer distribution factor less than 0 is exerted oneself for increasing, Calculate and increase the adjustable space that overload prevention device of exerting oneself associates electromotor；

${\mathrm{\ΔP}}_{g,j}=P_g^{max}-P_g^0$

Wherein, Δ P_g,jTo increase the adjustable space that overload prevention device of exerting oneself associates electromotor g,WithRespectively increase overload of exerting oneself Equipment association electromotor g's is maximum active and current active；

If 3-4) all adjustable spaces for increasing overload prevention device association electromotor of exerting oneself are 0, cannot be set by increasing overload It is standby to associate generator output to reduce the overload degree of overload prevention device.

3. a kind of multiple target based on interruptible load according to claim 1 transships aid decision-making method, and its feature exists In：Calculate interruptible load is to the power transfer distribution factor process of overload prevention device：Calculate corresponding with interruptible load etc. Duty value shifts distribution factor to the power of overload prevention device, the power is shifted into distribution factor overload is set as interruptible load Standby power transfer distribution factor.

4. a kind of multiple target based on interruptible load according to claim 1 transships aid decision-making method, and its feature exists In：Interruptible load is calculated to object function f₁And f₂The process of comprehensive contribution degree index be,

5-1) distribution factor is shifted to the power of all overload prevention devices according to interruptible load, calculate the controlling of interruptible load Can index；

${\mathrm{\γ}}_i=\underset{j=1}{\overset{N}{\Σ}}{\mathrm{\γ}}_{i,j}{\mathrm{\λ}}_j$

Wherein, γ_iFor interruptible load L_iControl performance standard, N for overload prevention device sum, λ_jFor overload prevention device Q_jLoad Rate, γ_i,jFor interruptible load L_iTo overload prevention device Q_jPower transfer distribution factor；

5-2) choose in interruptible load and be worth on the basis of the maximum of Control performance standard, the control performance of interruptible load is referred to Mark is normalized, and calculates interruptible load to reducing interruptible load excision total amount object function f₁Contribution degree refer to Mark；

$S_i^E={\mathrm{\γ}}_i/{\mathrm{\γ}}_{max}$

Wherein,For interruptible load L_iTo reducing interruptible load excision total amount object function f₁Contribution degree index, γ_max =max (γ_i) for the maximum of Control performance standard in interruptible load；

5-3) choose in interruptible load and be worth on the basis of the maximum of controllable capacity, the controllable capacity of interruptible load is returned One change is processed, and calculates interruptible load to reducing the total object function f of interruptible load excision₂Contribution degree index；

$S_i^W=P_i/P_{max}$

Wherein,For interruptible load L_iTo reducing the total object function f of interruptible load excision₂Contribution degree index, P_iFor Interruptible load L_iControllable capacity, P_max=max (P_i) it is interruptible load L_iIn controllable capacity the maximum；

5-4) interruptible load is calculated to object function f using the linear weighted function factor₁And f₂Comprehensive contribution degree index；

$S_i^C\left(\mathrm{\κ}\right)={\mathrm{\κS}}_i^E+(1-\mathrm{\κ})S_i^W$

Wherein,For the linear weighted function factor be κ when, interruptible load L_iTo object function f₁And f₂Comprehensive contribution degree refer to Mark.

5. a kind of multiple target based on interruptible load according to claim 1 transships aid decision-making method, and its feature exists In：The process of corresponding overload aid decision measure is when the calculating linear weighted function factor is κ,

7-1) AC power flow check number of times A is initialized as into 0, it is ε to arrange AC power flow and check step-length₂；

7-2) determine that overload prevention device overload is completely eliminated corresponding interruptible load resection successively, meet：

$\left\{\begin{array}{c}{E}_j\left(\mathrm{\&kappa;}\right)=\underset{i=1}{\overset{M_j}{\&Sigma;}}{P}_i\\ \underset{i=1}{\overset{M_j}{\&Sigma;}}{\mathrm{\&gamma;}}_{i,j}{P}_i>{\mathrm{\&Delta;P}}_j\\ \underset{i=1}{\overset{M_j-1}{\&Sigma;}}{\mathrm{\&gamma;}}_{i,j}{P}_i<{\mathrm{\&Delta;P}}_j\\ {\mathrm{\&gamma;}}_{i,j}\&GreaterEqual;{\mathrm{\&gamma;}}_{i+1,j}\end{array}\right.$

Wherein, E_j(κ) it is overload prevention device Q_jOverload is completely eliminated corresponding interruptible load resection, M_jFor interrupting for excision Load number,For overload prevention device Q_jPower more limit the quantity, P_jFor overload prevention device Q_jActive power, for Alternating current circuit,It is according to overload prevention device Q_jThe power limit that rated current and current voltage are converted into, for transformator,It is according to overload prevention device Q_jThe power limit that rated capacity and current voltage are converted into；For transmission cross-section,For defeated Power allowances of the electric section under the current method of operation, γ_i+1,jFor interruptible load L_i+1To overload prevention device Q_jPower transfer point The cloth factor；

7-3) using the maximum of all overload prevention devices correspondence interruptible load resection as interruptible load resection, it is considered to two The response condition of frequency modulation frequency modulation association electromotor and interruptible load, can interrupt negative to electromotor and not cutting off in secondary frequency-modulating process The active of lotus is adjusted, and control measure are checked by AC power flow；

7-4) AC power flow check number of times A is updated to into A+1, the power that all overload prevention devices are calculated respectively is more limited the quantity, and determines institute There is the corresponding overload prevention device Q of absolute value of the minima of the out-of-limit amount of power_d；

Process according to the following steps：

(1) ifThen by C_AAs final interruptible load control program；

Overload prevention device Q when number of times is A is checked for AC power flow_dPower more limit the quantity, C_AIt is A to check number of times for AC power flow The maximum of Shi Suoyou overload prevention devices correspondence interruptible load resection；

(2) if A>1、AndThen by C_AAs final interruptible load control program；

Overload prevention device Q when number of times is A-1 is checked for AC power flow_dPower more limit the quantity；

(3) if A>1、AndThen by C_A-1As final interruptible load control program；

C_A-1For the maximum that AC power flow checks all overload prevention device correspondence interruptible load resection when number of times is A-1；

(4) ifThen by overload prevention device Q_dPower get over limitation and be updated to Δ P_d=Δ P_d+ε₂, return to step 7-1；

(5) ifThen by overload prevention device Q_dPower get over limitation and be updated to Δ P_d=Δ P_d-ε₂, return to step 7-1.

6. a kind of multiple target based on interruptible load according to claim 5 transships aid decision-making method, and its feature exists In：In step 7-3, it is to electromotor and the active process that is adjusted for not cutting off interruptible load,

7-3-1) Control performance standard of the computation overload equipment association electromotor to all overload prevention devices；

${\mathrm{\&tau;}}_g=\underset{j=1}{\overset{N}{\&Sigma;}}{\mathrm{\&tau;}}_{g,j}{\mathrm{\&lambda;}}_j$

Wherein, τ_gControl performance standards of the electromotor g to all overload prevention devices, τ are associated for overload prevention device_g,jFor overload prevention device association Electromotor g is to overload prevention device Q_jPower transfer distribution factor；

7-3-2) according to system frequency and the frequency departure of rated frequency computing system；

Δ f=f-f₀

Wherein, Δ f is system frequency f and rated frequency f₀The frequency departure of computing system；

7-3-3) compare the corresponding frequency departure threshold value Δ f in AGC operating dead zones_AGCWith the size of Δ f, if Δ f >=-Δ f_AGC, Then select overload prevention device association electromotor drop to exert oneself by Control performance standard is descending, adjust total amount and C_AUnanimously, and terminate Adjustment process；Otherwise, the course of action for simulating frequency modulation frequency modulation associates electromotor to frequency modulation frequency modulation and does not cut off interruptible load Row adjustment, step is as follows：

(1) the active variable quantity that frequency modulation frequency modulation associates electromotor is participated at the end of calculating controlling cycle T；

${\mathrm{\&Delta;P}}_{g^{\&prime;}}=min(V_{g^{\&prime;}}T,P_{g^{\&prime;}}^{\max }-P_{g^{\&prime;}}^0)$

Wherein, Δ P_g′The active variable quantity that frequency modulation frequency modulation associates electromotor g ', V are participated at the end of for controlling cycle T_g′For secondary tune Frequency associates the climbing rate of electromotor g ',WithRespectively frequency modulation frequency modulation associates the maximum active of electromotor g ' and currently has Work(；

(2) calculate controlling cycle T at the end of system frequency variable quantity；

${\mathrm{\&Delta;f}}_T=(C_A+\underset{g^{\&prime;}=1}{\overset{N_G}{\&Sigma;}}{\mathrm{\&Delta;P}}_{g^{\&prime;}})/K_s$

Wherein, Δ f_TThe variable quantity of system frequency, K at the end of for controlling cycle T_sNot cut off total frequency modulation system of load in system Number, N_GTo participate in frequency modulation frequency modulation association electromotor quantity；

(3) duty value B after frequency change is calculated_kActive variable quantity；

${\mathrm{\&Delta;P}}_{B_k}=K_{B_k}{\mathrm{\&Delta;f}}_T$

Wherein,For duty value B after frequency change_kActive variable quantity,For duty value B_kFrequency change it is active Coefficient.

7. a kind of multiple target based on interruptible load according to claim 1 transships aid decision-making method, and its feature exists In：Total amount difference is cut off according to interruptible load and determines newly-increased shelves more than the linear weighted function factor region of search of computational accuracy ξ Position process be,

It is t and t+1 9-1) to define adjacent gear positions, then the corresponding linear weighted function factor is respectively κ_tAnd κ_t+1；

The gear number for 9-2) newly increasingE_t、E_t+1It is that interrupting for t and t+1 is negative to be divided into gear Lotus cuts off total amount；

9-3) by the region of search [κ_t, κ_t+1] press M_tRespectively obtain the linear weighted function factor of newly-increased gear.