CN106253353A - A kind of dynamic area control deviation allocation strategy improve and optimizate method - Google Patents

Abstract

The invention discloses an improved optimization method for a dynamic regional control deviation allocation strategy, which belongs to the field of electric power system control and comprises the following steps: Aiming at the deficiency of the current dynamic ACE in distributing fault loss power in a fixed proportion, a dynamic ACE considering the regional power grid ACE distribution is proposed The component correction strategy avoids the problem that the ACE in the control area enters the dead zone after the dynamic ACE starts and cannot assist in frequency recovery; a dynamic ACE component correction strategy based on the rotating reserve optimization model is proposed, and the analysis of each control area is carried out by establishing the maximum callable reserve optimization model Actually, the spinning reserve can be called out to provide a decision-making basis for the correction of the dynamic ACE component.

Description

An improved optimization method for dynamic regional control deviation allocation strategy

technical field

The invention relates to an improved optimization method for a dynamic region control deviation allocation strategy, which belongs to the technical field of power systems.

Background technique

When the incoming power outside the high-power area of the receiving end grid is lost or the large unit trips, the power transmission and reception plan in each control area cannot be automatically adjusted according to the grid situation, and the method of manually modifying the tie line exchange plan cannot meet the requirements for rapid recovery of frequency and power. Since the automatic generation control (AGC) of each control area is adjusted according to the regional control deviation (ACE) of the control area, after the fault occurs, the DC drop point or the fault control area will greatly increase the output adjustment due to less transmission or more reception, and the non-DC drop The control area of the non-DC drop point may be the province receiving power from outside the area, and the failure to adjust the contact line plan in time will not be conducive to the frequency and power after the accident. recover.

According to its own characteristics and needs, East China Power Grid proposes dynamic area control deviation (hereinafter referred to as dynamic ACE) technology to improve the power grid’s recovery ability against high power loss and effectively develop the potential of backup and sharing of the entire network. Practical applications show that under high-power loss faults such as DC blocking, dynamic ACE can effectively utilize the backup potential of the entire network to quickly restore power loss due to faults and reduce frequency recovery time. However, during the implementation of dynamic ACE, there are also the following problems:

1) The dynamic ACE does not consider the initial ACE of each control area, so that after the ACE in the control area is superimposed on the fault power allocated by the dynamic ACE, the ACE is positive or in a dead zone, and cannot support power recovery;

2) The dynamic ACE adopts a fixed coefficient to apportion the accident loss power proportionally, without considering the actual adjustable reserve capacity in the control area, and will not be able to provide effective support due to the control area's own adjustable ability or the limitation of the tidal flow section.

The document "Application Analysis of Dynamic Regional Control Errors in East China Power Grid" (Power System Automation, Vol. 28, No. 1, Page 78, 2004) discloses an improved method for dynamic ACE considering the tripping of trans-regional power transmission units and tie-line overrun constraints. , focusing on the current dynamic ACE architecture, the improvement and improvement strategies in the aspects of fault power loss capacity, cascading faults and tie-line safety, in order to enhance the adaptability of dynamic ACE technology. The strategy proposed in the literature is to use the latest sensitivity coefficient and manually adjust the dynamic ACE control tie line overrun strategy, which can effectively avoid the tripping situation caused by the power of the tie line exceeding its safe and stable limit after the dynamic ACE is triggered. This method is manually adjusted and canceled by the dispatcher. When a power grid emergency occurs, it will undoubtedly increase the dispatching work intensity and the probability of miscontrol, and cannot guarantee the timeliness and accuracy of the control.

Contents of the invention

The purpose of the present invention is to overcome the deficiencies in the prior art, provide an improved optimization method for the dynamic area control deviation allocation strategy, and solve the technical problem that the control area ACE enters the dead zone after the dynamic ACE starts and cannot assist frequency recovery.

In order to solve the above-mentioned technical problems, the technical solution adopted in the present invention is: an improved optimization method of a dynamic region control deviation allocation strategy, comprising the following steps:

Step 1: Aiming at the shortage of the current dynamic ACE in distributing fault loss power in a fixed proportion, a dynamic ACE component correction strategy considering the ACE distribution of the regional power grid is proposed;

Step 2: Propose a dynamic ACE component correction strategy based on the spinning reserve optimization model. By establishing the maximum callable reserve optimization model, analyze the actual callable spinning reserve in each control area, and provide decision-making basis for the correction of dynamic ACE components.

The specific method of step 1 is as follows:

1) Calculate the active unbalanced power of the regional power grid after the fault occurs:

${\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}\mathrm{<span\; class="notranslate">\; D.</span>}}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; no</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; \&Delta;P</span>}}_{\mathrm{<span\; class="notranslate">\; the\; s</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; )</span>$

In the formula, ACE _ND is the ACE of the regional power grid; n is the number of control areas to which the regional power grid belongs; ACE _i is the ACE before the failure of the control area i; ΔP _s is the power loss of the regional power grid fault;

2) Allocate the active unbalanced power of the regional power grid according to the rotating reserve coefficient:

${{\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}^{<span\; class="notranslate">\; \&prime;</span><span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; =</span>\frac{{\mathrm{<span\; class="notranslate">\; k</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}{\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; no</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; k</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}\mathrm{<span\; class="notranslate">\; D.</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; )</span>$

In the formula, ACE _i ″ is the control target of the ACE in the control area _i ; ki is the rotation reserve coefficient of the control area i; 3) Calculate the correction component of the dynamic ACE:

ACE _i ″=ACE _i ′+ACE _DI,i ＝(ACE _i +ΔP _s,i )+ACE _DI,i (3)

In the formula, ACE _i ′ is the ACE after superimposing the dynamic ACE component in the control area i; ACE _DI,i is the correction component of the dynamic ACE in the control area i; ΔP _s,i is the dynamic ACE component in the control area i,

Then the correction component of the dynamic ACE is obtained:

ACE _DI,i = ACE _i ″-ACE _i ′ (4);

4) After the correction component ACE _DI,i is superimposed with the dynamic ACE component, it is sent to each control area as the corrected dynamic ACE component and executed. The calculation formula is:

ΔP _s,i '=ΔP _s,i +ACE _DI,i (5)

In the formula, ΔP _s,i ′ is the corrected dynamic ACE component of control area i;

5) The modified dynamic ACE component only replaces the original dynamic ACE component, and the start-up judgment condition, delivery method and recovery time of the dynamic ACE remain unchanged.

The specific method of step two is:

1) Establish the maximum adjustable spinning reserve optimization model:

The optimization objective is to provide the most spinning reserve capacity for each control area and regional power grid in the time period, the spinning reserve release coefficient c _j,t is the decision variable, and its objective function is:

$\mathrm{<span\; class="notranslate">\; min</span>}\mathrm{<span\; class="notranslate">\; f</span>}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 6</span>}<span\; class="notranslate">\; )</span>$

where r _j,t = P _j,max -P _j,t ;

In the formula, r _j,t is the upward spinning reserve capacity of unit j at time t; c _j,t is the spinning reserve release coefficient of unit j at time t; N is the number of units; P _j,max is the upper limit of adjustment of unit j; P _j,min is the adjustment lower limit of unit j; P _j,t is the active output of unit j at time t;

Constraints include line flow constraints and unit output constraints:

${\underset{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; T</span>}}}_{\mathrm{<span\; class="notranslate">\; z</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; z</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; +</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; the\; s</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; z</span>}}<span\; class="notranslate">\; \&le;</span>{\stackrel{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; T</span>}}}_{\mathrm{<span\; class="notranslate">\; z</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 7</span>}<span\; class="notranslate">\; )</span>$

In the formula, f _z,t is the active power flow of transmission line z at time t; N is the number of units; r _j,t is the spinning reserve capacity of unit j at time t; c _j,t is the spinning reserve capacity of unit j at time t release coefficient; s _j,z is the sensitivity coefficient of unit j to transmission line z; T _z,t , are the active lower limit and upper limit of the transmission line z respectively;

$\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; \{</span><span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}^{\mathrm{<span\; class="notranslate">\; r</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; p</span>}}\mathrm{<span\; class="notranslate">\; \&Delta;</span>}\mathrm{<span\; class="notranslate">\; t</span>}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \}</span><span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \&le;</span>\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; \{</span>{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}^{\mathrm{<span\; class="notranslate">\; r</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; p</span>}}\mathrm{<span\; class="notranslate">\; \&Delta;</span>}\mathrm{<span\; class="notranslate">\; t</span>}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \}</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 8</span>}<span\; class="notranslate">\; )</span>$

In the formula, is the climbing rate of unit j; Δt is the statistical period, 0<Δt≤15, rounded to an integer;

The maximum callable spinning reserve is:

${\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 9</span>}<span\; class="notranslate">\; )</span>$

In the formula, R _i,t is the maximum adjustable spinning reserve in the control area i;

2) Using the optimization model of 1) to calculate the maximum adjustable spinning reserve of the regional power grid and each control area, and compare it with the corrected dynamic ACE component;

3) Regional power grid verification:

Let R _t be the maximum adjustable spinning reserve of the regional power grid n is the number of control areas to which the regional power grid belongs;

If ΔP _s ≤ R _t , the regional power grid can cope with the current fault;

If ΔP _s >R _t , the control area under the regional power grid only allocates the maximum transferable spinning reserve capacity, and takes other measures, including: emergency start-up and shutdown, rapid switch-on and withdrawal of pumped storage, or application for outsourced power support. To make up for the excess power;

4) Verification and correction of each control area:

a) If ΔP _s,i ′≤R _i,t , then the spinning reserve in the control area is sufficient to deal with the apportioned fault loss power;

b) If ΔP _s,i ′>R _i,t , then the spinning reserve in the control area is not enough to cope with the shared fault power loss, and other control areas need to be supported;

When the situation b) occurs in a certain control area, based on the principle of rapid fault recovery, the shortfall part (ΔP _s,i ′-R _i,t ) of the control area is allocated to other control areas according to the rotating reserve coefficient; after allocation, if other In the case of b) in the control area, a similar method is adopted to re-allocate in the control area where sufficient spinning reserve can be called out.

Compared with the prior art, the beneficial effects achieved by the present invention are:

1. Under the framework of the existing dynamic ACE control system, the method of apportioning the power loss due to failure in each control area in the dynamic ACE is improved, aiming at the shortage of the current dynamic ACE in distributing the power loss due to failure in a fixed proportion, based on giving full play to the backup sharing advantages of the regional power grid , the principle of safely and reasonably calling out the spinning reserve, and propose a dynamic ACE component correction strategy that considers the ACE distribution of the regional power grid, so as to give full play to the role of the interconnected power grid in the recovery of large power loss disturbances, and avoid the ACE in the control area from entering the dead zone after the dynamic ACE starts. Failure to assist with frequency recovery issues;

2. Propose a dynamic ACE component correction strategy based on the spinning reserve optimization model. By establishing the maximum callable reserve optimization model, analyze the actual callable spinning reserve in each control area, and provide decision-making basis for the correction of dynamic ACE components.

Description of drawings

Figure 1 is a schematic diagram of the relationship between the regional power grid and its control area.

Figure 2 is a schematic diagram of the control area and the power of generator set regulation.

detailed description

The present invention will be further described below in conjunction with the accompanying drawings. The following examples are only used to illustrate the technical solution of the present invention more clearly, but not to limit the protection scope of the present invention.

As shown in Figure 1, it is a schematic diagram of the relationship between the regional power grid and its control area, where: RG is the regional power grid, and CA is the control area that belongs to the regional power grid RG. The CAs under the RG form a backup sharing group to share the power lost due to faults.

The direct adjustment control area and the provincial adjustment control area calculate and adjust the capacity according to ACE and send it to the units under the jurisdiction of the control area. The units under the jurisdiction of the direct regulation control area are distributed in each province. Each province takes the output of the direct regulation unit as its control boundary to form a closed control area; Closed control area. As shown in Figure 2, it is a schematic diagram of the control area and the power of generator set regulation. The direct control area and the provincial control area bear the loss power sharing responsibility according to the rotating reserve ratio.

The present invention is an improved and optimized method for distribution strategy of dynamic regional control deviation, which includes the following steps: aiming at the deficiency of current dynamic ACE in distributing fault loss power in a fixed proportion, a dynamic ACE component correction strategy considering regional power grid ACE distribution is proposed to avoid dynamic ACE start-up The ACE in the rear control area enters the dead zone and cannot assist in frequency recovery; a dynamic ACE component correction strategy based on the spinning reserve optimization model is proposed, and the maximum callable reserve optimization model is established to analyze the actual callable spinning reserve in each control area. Correction of dynamic ACE components provides decision-making basis. Specifically include the following steps:

Step 101: When a large power loss fault such as DC blockage or tripping of a large unit occurs in the regional power grid, the dynamic ACE technology allocates the fault loss power according to the rotating reserve ratio, and superimposes it as the dynamic ACE component of each control area in the reserve sharing group to each control area In ACE, its calculation formula is:

ACE _i '=ACE _i +ΔP _s,i (1)

in,

In the formula, ACE _i ′ is the ACE after superimposing the dynamic ACE component in the control area i; n is the number of control areas to which the regional power grid belongs; ACE _i is the ACE before the failure of the control area i; ΔP _s is the power loss of the regional power grid fault ( Negative value); ΔP _s,i is the dynamic ACE component of control area _i (negative value); ki is the spinning reserve coefficient of control area i (the distribution coefficient agreed by the regional power grid for each control area).

Step 102: Calculate the active unbalanced power of the regional power grid as:

${\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}\mathrm{<span\; class="notranslate">\; D.</span>}}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; no</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}<span\; class="notranslate">\; +</span>{\mathrm{<span\; class="notranslate">\; \&Delta;P</span>}}_{\mathrm{<span\; class="notranslate">\; the\; s</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 2</span>}<span\; class="notranslate">\; )</span>$

In the formula, ACE _ND is the regional grid ACE;

Step 103: Allocate the active unbalanced power of the regional grid according to the rotating reserve coefficient:

${{\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}^{<span\; class="notranslate">\; \&prime;</span><span\; class="notranslate">\; \&prime;</span>}<span\; class="notranslate">\; =</span>\frac{{\mathrm{<span\; class="notranslate">\; k</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}{\underset{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; no</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; k</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}}}{\mathrm{<span\; class="notranslate">\; ACE</span>}}_{\mathrm{<span\; class="notranslate">\; N</span>}\mathrm{<span\; class="notranslate">\; D.</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 3</span>}<span\; class="notranslate">\; )</span>$

In the formula, ACE _i ″ is the control target of the ACE in the control area i;

The distribution method of active and unbalanced power in the regional power grid can ensure that the ACE direction of each control area is consistent, effectively avoid the influence of the initial ACE in the control area, and simultaneously participate in the adjustment of the unbalanced power of the regional power grid.

Step 104: Calculate the corrected component of the dynamic ACE:

ACE _i ″=ACE _i ′+ACE _DI,i ＝(ACE _i +ΔP _s,i )+ACE _DI,i (4)

In the formula, ACE _DI,i is the correction component of the dynamic ACE in the control area i;

Which leads to:

ACE _DI,i = ACE _i "-ACE _i ' (5).

Step 105: After the correction component and the dynamic ACE component are superimposed, the corrected dynamic ACE component is issued to each control area and executed. The calculation formula is:

ΔP _s,i '=ΔP _s,i +ACE _DI,i (6)

In the formula, ΔP _s,i ′ is the corrected dynamic ACE component.

Step 201: Establish a maximally callable spinning reserve optimization model:

The optimization objective is to provide the most spinning reserve capacity for each control area and regional power grid in the time period, the spinning reserve release coefficient c _j,t is the decision variable, and its objective function is:

$\mathrm{<span\; class="notranslate">\; min</span>}\mathrm{<span\; class="notranslate">\; f</span>}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 1</span>}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; )</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 7</span>}<span\; class="notranslate">\; )</span>$

where r _j,t = P _j,max -P _j,t ;

In the formula, r _j,t is the upward spinning reserve capacity of unit j at time t; c _j,t is the spinning reserve release coefficient of unit j at time t; N is the number of units; P _j,max is the upper limit of adjustment of unit j; P _j,min is the adjustment lower limit of unit j; P _j,t is the active output of unit j at time t.

Constraints include line flow constraints and unit output constraints:

${\underset{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; T</span>}}}_{\mathrm{<span\; class="notranslate">\; z</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; f</span>}}_{\mathrm{<span\; class="notranslate">\; z</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; +</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; the\; s</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; z</span>}}<span\; class="notranslate">\; \&le;</span>{\stackrel{<span\; class="notranslate">\; \&OverBar;</span>}{\mathrm{<span\; class="notranslate">\; T</span>}}}_{\mathrm{<span\; class="notranslate">\; z</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 8</span>}<span\; class="notranslate">\; )</span>$

In the formula, f _z,t is the active power flow of transmission line z at time t; N is the number of units; s _j,z is the sensitivity coefficient of unit j to transmission line z; T _z,t , are the active lower limit and upper limit of transmission line z respectively.

$\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; \{</span><span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}^{\mathrm{<span\; class="notranslate">\; r</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; p</span>}}\mathrm{<span\; class="notranslate">\; \&Delta;</span>}\mathrm{<span\; class="notranslate">\; t</span>}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \}</span><span\; class="notranslate">\; \&le;</span>{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \&le;</span>\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; i</span>}\mathrm{<span\; class="notranslate">\; no</span>}<span\; class="notranslate">\; \{</span>{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}}^{\mathrm{<span\; class="notranslate">\; r</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; p</span>}}\mathrm{<span\; class="notranslate">\; \&Delta;</span>}\mathrm{<span\; class="notranslate">\; t</span>}<span\; class="notranslate">\; ,</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; m</span>}\mathrm{<span\; class="notranslate">\; a</span>}\mathrm{<span\; class="notranslate">\; x</span>}}<span\; class="notranslate">\; -</span>{\mathrm{<span\; class="notranslate">\; P</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; \}</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 9</span>}<span\; class="notranslate">\; )</span>$

In the formula, is the climbing rate of unit j; Δt is the statistical period, 0<Δt≤15, rounded to an integer;

The maximum callable spinning reserve is:

${\mathrm{<span\; class="notranslate">\; R</span>}}_{\mathrm{<span\; class="notranslate">\; i</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; =</span>\underset{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; =</span>\mathrm{<span\; class="notranslate">\; 1</span>}}{\overset{\mathrm{<span\; class="notranslate">\; N</span>}}{<span\; class="notranslate">\; \&Sigma;</span>}}{\mathrm{<span\; class="notranslate">\; r</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}{\mathrm{<span\; class="notranslate">\; c</span>}}_{\mathrm{<span\; class="notranslate">\; j</span>}<span\; class="notranslate">\; ,</span>\mathrm{<span\; class="notranslate">\; t</span>}}<span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; -</span><span\; class="notranslate">\; (</span>\mathrm{<span\; class="notranslate">\; 10</span>}<span\; class="notranslate">\; )</span>$

In the formula, R _i,t is the maximum up-adjustable spinning reserve of the control area i.

Step 202: Use the optimization model in step 201 to calculate the maximum adjustable spinning reserve of the regional power grid and each control area, and compare it with the corrected dynamic ACE component.

Step 203: regional grid verification:

Let R _t be the maximum adjustable spinning reserve of the regional power grid n is the number of control areas to which the regional power grid belongs;

If ΔP _s ≤ R _t , the regional power grid can cope with the current fault;

If ΔP _s > R _t , the control area under the regional power grid only allocates the maximum adjustable spinning reserve capacity, and takes other measures, including: emergency start-up and shutdown, rapid switch-on and withdrawal of pumped storage, or application for outsourced power support, to Make up for excess power.

Step 204: Verification and correction of each control area:

a) If ΔP _s,i ′≤R _i,t , then the spinning reserve in the control area is sufficient to deal with the apportioned fault loss power;

b) If ΔP _s,i ′>R _i,t , the spinning reserve in the control area is not enough to cope with the shared fault power loss, and other control areas need to provide support.

When the situation b) occurs in a certain control area, based on the principle of quick recovery from the fault, the shortfall (ΔP _s,i ′-R _i,t ) of the control area is allocated to other control areas according to the rotating reserve coefficient; after allocation, if other In the case of b) in the control area, a similar method is adopted to re-allocate in the control area where sufficient spinning reserve can be called out.

The corrected dynamic ACE component only replaces the original dynamic ACE component. The dynamic ACE component is sent to the provincial control areas within one minute of the fault, and the component starts to decrease in the next full 15 minutes. The connection line plan is gradually adjusted according to the power loss power, and the dynamic ACE decreases to 0 in the next 15 minutes. .

The above is only a preferred embodiment of the present invention, and it should be pointed out that for those of ordinary skill in the art, without departing from the technical principle of the present invention, some improvements and modifications can also be made. It should also be regarded as the protection scope of the present invention.

Claims (3)

1. An improved optimization method for a dynamic regional control deviation allocation strategy is characterized in that: comprises the following steps:

the method comprises the following steps: aiming at the defect that the current dynamic ACE distributes fault loss power in a fixed proportion, a dynamic ACE component correction strategy considering the ACE distribution of a regional power grid is provided;

step two: and providing a dynamic ACE component correction strategy based on a rotary standby optimization model, and analyzing the actually-adjustable rotary standby of each control area by establishing a maximum adjustable standby optimization model so as to provide decision basis for correcting the dynamic ACE component.

2. The improved optimization method of the dynamic zone control deviation allocation strategy according to claim 1, wherein the specific method of the first step is as follows:

1) calculating the active unbalanced power of the regional power grid after the fault occurs:

${\mathrm{ACE}}_{ND}=\underset{i=1}{\overset{n}{\&Sigma;}}{\mathrm{ACE}}_i+{\mathrm{\&Delta;P}}_s---\left(1\right)$

in the formula, ACE_NDIs an area grid ACE; n is the number of control areas to which the regional power grid belongs; ACE_iThe ACE before the control area i fails; delta P_sPower is lost for regional grid faults;

2) distributing the active unbalanced power of the regional power grid according to the rotation standby coefficient:

${\mathrm{ACE}}_i^{\&prime;\&prime;}=\frac{k_i}{\underset{i=1}{\overset{n}{\&Sigma;}}{k}_i}{\mathrm{ACE}}_{ND}---\left(2\right)$

in the formula, ACE_i"is the control target of the ACE of control area i; k is a radical of_iIs the spinning spare coefficient of control area i; 3) calculating the correction component of the dynamic ACE:

ACE_i″＝ACE_i′+ACE_DI,i＝(ACE_i+ΔP_s,i)+ACE_DI,i(3)

in the formula, ACE_i' is ACE after the control area i is superposed with dynamic ACE components; ACE_DI,iA correction component for dynamic ACE for control area i; delta P_s,iTo control the dynamic ACE component of zone i,

and further obtaining a correction component of the dynamic ACE:

ACE_DI,i＝ACE_i″-ACE_i′ (4)；

4) correction component ACE_DI,iAnd the dynamic ACE component is superposed and then is taken as a modified dynamic ACE component, and the modified dynamic ACE component is issued to each control area and executed, and the calculation formula is as follows:

ΔP_s,i′＝ΔP_s,i+ACE_DI,i(5)

in the formula,. DELTA.P_s,i' is the dynamic ACE component after the control area i is corrected;

5) the corrected dynamic ACE component only replaces the original dynamic ACE component, and the starting judgment condition, the issuing mode and the recovery time of the dynamic ACE are kept unchanged.

3. The improved optimization method of the dynamic area control deviation allocation strategy according to claim 2, wherein the specific method in the second step is as follows:

1) establishing a maximum adjustable rotary standby optimization model:

the optimization target provides the maximum spinning reserve capacity and spinning reserve release coefficient c for each control area and regional power grid in the time period_j,tFor decision variables, the objective function is:

$\min F=\underset{j=1}{\overset{N}{\&Sigma;}}{r}_{j,t}(1-c_{j,t})---\left(6\right)$

wherein r is_j,t＝P_j,max-P_j,t；

In the formula, r_j,tThe upper rotation reserve capacity of the unit j at the moment t; c. C_j,tThe rotating standby release coefficient of the unit j at the moment t; n is the number of the units; p_j,maxAdjusting an upper limit for the unit j; p_j,minAdjusting a lower limit for the unit j; p_j,tThe active output of the unit j at the moment t;

the constraint conditions comprise line tide constraint and unit output constraint:

${\underset{\&OverBar;}{T}}_{z,t}\&le;f_{z,t}+\underset{j=1}{\overset{N}{\&Sigma;}}{r}_{j,t}{c}_{j,t}{s}_{j,z}\&le;{\stackrel{\&OverBar;}{T}}_{z,t}---\left(7\right)$

in the formula (f)_z,tThe active power flow of the transmission line z at the moment t; n is the number of the units; r is_j,tThe upper rotation reserve capacity of the unit j at the moment t; c. C_j,tThe rotating standby release coefficient of the unit j at the moment t; s_j,zThe sensitivity coefficient of the unit j to the transmission line z is set;T _z,t、respectively an active lower limit and an active upper limit of the transmission line z;

$min\{-r_j^{ramp}\mathrm{\&Delta;}t,P_{j,\min }-P_{j,t}\}\&le;r_{j,t}{c}_{j,t}\&le;min\{r_j^{ramp}\mathrm{\&Delta;}t,P_{j,max}-P_{j,t}\}---\left(8\right)$

in the formula,the ramp rate of the unit j; Δ t is a statistical period, 0<Delta t is less than or equal to 15, and an integer is taken;

the maximum adjustable rotary standby is as follows:

$R_{i,t}=\underset{j=1}{\overset{N}{\&Sigma;}}{r}_{j,t}{c}_{j,t}---\left(9\right)$

in the formula, R_i,tThe control area i can be adjusted up for standby;

2) calculating the maximum adjustable rotation standby of the regional power grid and each control area by using 1) an optimization model, and comparing the maximum adjustable rotation standby with the corrected dynamic ACE component;

3) checking a regional power grid:

let R_tRotatable reserve for maximum adjustable regional power gridn is the number of control areas to which the regional power grid belongs;

if Δ P_s≤R_tIf so, the regional power grid can deal with the current fault;

if Δ P_s>R_tAnd then, the control area under the regional power grid is only allocated with the maximum adjustable rotating reserve capacity, and other measures are taken, including: emergency starting and stopping, pumped storage fast switching on and off or applying for external power purchase support mode to make up for the excess power;

4) checking and correcting each control area:

a) if Δ P_s,i′≤R_i,tThen the control area spinning reserve is sufficient to handle the apportioned fault loss power;

b) if Δ P_s,i′>R_i,tIf the control area is rotated for standby, the control area is not enough to deal with the distributed fault loss power, and other control areas are required to be supported;

when a condition b) occurs in a certain control area, the shortage part (delta P) of the control area is recovered on the principle of quick fault recovery_s,i′－R_i,t) Distributing the data to other control areas according to the rotating standby coefficient; after allocation, if the other control areas have the condition of b), the other control areas are allocated again by adopting a similar method, and the control areas with enough rotation reserve can be called out.