The following relates to the technical field of electric power systems control. It is applicable to establishment and implementation of adaptive control strategy in case of equipment overload in electric power systems.
With enhanced grid structure and increase of grid transmission capacity, equipment overload has gradually become a main issue that constraints safe, reliable, and economic operation of grid. In the stage of fast grid development, complicated grid frame structure, and varying operation modes, for equipment (line, transformer) overload, it is still difficult to use offline analysis to establish reasonable control strategy. To solve the problems of insufficient dispatching operation control execution time and control quantity, over high control cost, over centralized load control measures, etc. it is urgent to study equipment overload adaptive control method based on centralized real-time decision-making.
Patent “Centralized decision-making and real-time emergency control method for eliminating overload in large power grid” (ZL200710135096.2) proposes the following method: according to information measured by security control device (SCD), obtain grid operation state data through adjustment of EMS state estimation data, and based on this, according to sensitivity of eliminating equipment overload by candidate control measures (generator and load etc.), taking into consideration cost of control measures and extent of equipment overload, calculate comprehensive indexes of control measures; next, according to sequence of absolute values of comprehensive indexes, select generator units of positive index and load nodes of negative index for participation in emergency control, and bisection method is used to calculate optimized emergency control scheme. This method can quickly calculate control strategy, but the search strategy of control measures according to sequence of comprehensive indexes cannot ensure optimization, and does not involve judgment of equipment actual overload extent and permissible continuous operation time, nor coordination between emergency control and dispatching operation control.
Different from prior art, the present invention discloses the method of real-time estimation of equipment long-term allowable current and continuous operation time according to measured current and temperature information, and determination of control modes according to continuous operation time of equipment instead of the long-term allowable current. The decision optimization models of dispatching operation control and SCD emergency control aimed to minimize control costs have been established, taking into consideration discrete and decentralized requirements on load control. The control strategy of successive approximation according to equipment state is adopted, introducing mixed integer nonlinear programming algorithm, thus solving the problems of insufficient dispatching operation control execution time and control quantity, over high control cost, over centralized load control measures, etc. while realizing comprehensive application of and coordination between emergency control measures and dispatching operation control measures for eliminating equipment overload.
The purpose of the present invention is to propose an equipment overload adaptive emergency control method based on centralized real-time decision-making, that can automatically judge equipment overload state according to present and history operation state of the equipment, comprehensively apply emergency control and dispatching operation control means to eliminate equipment overload, and solve the problems of insufficient dispatching operation control execution time and control quantity, over high control cost, over centralized load control measures, etc., and ensure reliability and optimization of equipment overload control.
The present invention is realized by the following technical scheme, comprising the following steps:
1) Based on equipment current acquired in real-time by security control device (SCD), and temperature information acquired in real-time by equipment temperature monitoring device, estimate in real-time long-term allowable current and continuous operation time of the equipment under present operating environment;
2) If minimum value of continuous operation time of all equipment is greater than the set dispatching operation control time limit, neither control strategy calculation nor control of the electric system will be performed. If this minimum value is less than or equal to the dispatching operation control time limit but greater than the emergency control time limit, go to step 3) for optimal dispatching operation control strategy calculation and implementation. If this minimum value is less than or equal to the emergency control time limit, go to step 4) for optimal emergency control strategy calculation and implementation;
3) First, based on grid beyond the dispatching management (abbreviated as External Network) state estimation data or typical operation mode data, carry out static equivalence of the External Network according to the criterion that electric distance between External Network side buses of tie-line between External Network and Internal Network (grid under dispatching management) shall exceed the set value. Next, based on Internal Network state estimation data and External Network equivalence data, carry out the operation profile data integration according to real-time power flow data of equipment of continuous operation time less than or equal to dispatching operation control time limit (such equipment is referred to as overload equipment in alarm state) and power flow data acquired in real-time by SCD. Later, for the grid after equivalence, using the objective function aimed to minimize total control cost and comprehensive indexes of proportions of load control quantity in different regions, taking into consideration discreteness and cost of dispatching operation control measures, and under the restraint of electric system power flow, and based on mixed integer nonlinear programming algorithm, centralized decision-making optimization and successive approximation control strategy are adopted for calculation of real-time optimal dispatching operation control strategy for equipment overload in alarm state, which will be implemented by dispatching operator; after implementation of control measures, return to step 1);
4) First, based on External Network state estimation data or typical operation mode data, carry out static equivalence of the External Network according to the criterion that electric distance between External Network side buses of tie-line between External Network and Internal Network (grid under dispatching management) shall exceed the set value. Next, based on Internal Network state estimation data and External Network equivalence data, carry out the operation profile data integration according to real-time power flow data of equipment of continuous operation time less than or equal to emergency control time limit (such equipment is referred to as overload equipment in emergent state) and power flow data acquired in real-time by SCD. Later, for the grid after equivalence, using the objective function aimed to minimize total control cost and comprehensive indexes of proportions of load control quantity in different regions, taking into consideration discreteness and cost of emergency control measures, and under the restraint of electric system power flow, and based on mixed integer nonlinear programming algorithm, centralized decision-making optimization and successive approximation control strategy are adopted for calculation of real-time optimal emergency control strategy for equipment overload in emergent state, which will be implemented by SCD; after implementation of control measures, return to step 1).
Since equipment long-term allowable current and continuous operation time under certain set current condition are related to equipment operating environment (e.g. ambient temperature, sunshine, and wind velocity etc.), the present invention has solved the problem by real-time estimation of equipment long-term allowable current and continuous operation time according to present measured information and history current and temperature information. In prior art, equipment overload control mode is determined according to equipment current exceeding different preset values. As a contrast, the present invention decides control mode according to equipment continuous operation time instead of current, thus meeting actual requirements on selection of control mode. The present invention has established decision-making optimization models of dispatching operation control and SCD emergency control aimed to minimize control costs respectively, taking into consideration discrete and decentralized requirements on load control, adopting successive approximation control strategy according to equipment state, and introducing mixed integer nonlinear programming algorithm, thus solving the problems of insufficient dispatching operation control execution time and control quantity, over high control cost, over centralized load control measures, etc. while realizing comprehensive application of and coordination between emergency control and dispatching operation control means that eliminate equipment overload, to ensure reliable and optimized equipment overload control.
The following describes method of the present invention in greater details in combination with
In
In
For equipment for which the temperature can be actually measured, real-time estimation of Ir and Δt is only carried out if the measured temperature rises and the ratio of measured current to preset rated current exceeds a set threshold. For such equipment that does not satisfy these two conditions, Ir is taken as its rated current and Δt is set to long-term. In collected equipment current and temperature history information (I(t), T(t)), take history data of two periods starting from the most recent measurement time point (earlier than this point). Assume that the first period provides data of m time points in chronological order, namely [(I1.i(t1.i), T1.i(t1.i)), i=1, 2, . . . , m], and that the second period provides data of n time points in chronological order, namely [(I2.j(t2.j), T2.j(t2.j)), j=1, 2, . . . , n]. Then, equipment Ir is estimated using formula (1), where
and k1 is a correction factor.
According to equipment current I(trt) and corresponding temperature T(trt) of equipment acquired in real-time at the most recent time point, in history information (I(t), T(t)) of current and temperature of this equipment, starting from the most recent measurement time point backward (to earlier time), find measurement time points of T(t) less than T(trt) in sequence. If the period between two time points (trt−t) exceeds a set value and the equipment current difference between them is less than a set value, then formula (2) is used to estimate Δt of this equipment. In this formula, Tcr is the highest permissible operation temperature of the equipment under present environment, and k2 is a correction factor.
For equipment for which temperature is not actually measured, if the function (Δt=f(I)) of equipment Δt to current in present operating environment is available, equipment Ir is taken as the value of current corresponding to dispatching operation control time limit td (e.g. 15 min) timed by a coefficient less than 1. If only equipment Δt to current correspondence table (Δtk, Ik) under present operating environment is available, curve fitting will be carried out according to such time to current correspondence points, to obtain function (Δt=f(I)) of Δt to current of this equipment, and then equipment Ir is taken as the value of current corresponding to td timed by a coefficient less than 1.
If the present measured current of equipment exceeds Ir, in equipment current history information collected, take history data of a period starting from the most recent measurement time point (to earlier time). Assume that this period contains data of m time points in chronological order, namely [Ii(ti), i=1, 2, . . . , m], it is required that current at the first time point of this period I1(t1) is less than or equal to Ir, and that current at the second time point of this period I2(t2) exceeds Ir. Formula (3) is used to estimate Δt of this equipment. If the present measured current of the equipment is less than or equal to Ir, Δt of this equipment is set to long-term.
In
In
In
In
Sub-step 1: External Network static equivalence: to be carried out based on External Network state estimation data or typical operation mode data, and according to the criterion that electric distance between External Network side buses of tie-line between External Network and Internal Network (grid under dispatching management) shall exceed the set value.
Sub-step 2: Internal Network and External Network operation profile data integration: based on Internal Network state estimation data and External Network equivalence data, and according to real-time power flow data of equipment of Δt less than td (abbreviated as overload equipment in alarm state) and power flow data acquired in real-time by SCD.
Sub-step 3: optimal dispatching operation control strategy calculation: for the grid after equivalence, using the objective function aimed to minimize total control cost and comprehensive indexes of proportions of load control quantity in different regions, taking into consideration discreteness and cost of dispatching operation control measures, and under the restraint of electric power system power flow, and based on mixed integer nonlinear programming algorithm, centralized decision-making optimization and successive approximation control strategy are adopted to calculate real-time optimal dispatching operation control strategy for equipment overload in alarm state. Its objective function is formula (4). In this formula, term 1 is the cost of generator active power adjustment, PGi and P′Gi are active power output before and after adjustment of generator i in this round of control strategy optimization calculation, Ci is unit active power adjustment cost of generator i, and G is total number of generators that can be adjusted. Term 2 is the load control cost: if load j is shed, Lj is 1; otherwise it is 0; CLj is load shedding cost of load j, and L is total number of loads that can be controlled. Term 3 reflects decentralized load control requirements: N is total number of regions for which power outage effect index is examined, PLj is active power of the load at the most recent time point after occurrence of this overload event and before taking the control measures, Zk is the kth examined region, x is a set coefficient (greater than 1), and k3 is the factor of converting power outage effect to control cost.
Corresponding restraint conditions are given in formula (5), including power flow equation restraint (including bus voltage restraint), equipment overload restraint, and dispatching operation control measures space restraint to be selected. In this formula, Irj is long-term allowable current of equipment j estimated in step 2, Ij0 is current in equipment j at the most recent time point before this round of control strategy calculation, λd is set dispatching operation control successive approximation coefficient, and Md is number of equipment for which Δr estimated in real-time in sub-step 2 is less than k′dtd (k′d is greater than 1). This optimization calculation of dispatching operation control strategy for overload event has considered the restraint that generator active power output adjustment direction cannot be reversed and that after load control, such load cannot be restored.
In
In
In
Sub-step 1: External Network static equivalence: to be carried out based on External Network state estimation data or typical operation mode data, and according to the criterion that electric distance between External Network side buses of tie-line between External Network and Internal Network (grid under dispatching management) shall exceed the set value.
Sub-step 2: Internal Network and External Network operation profile data integration: based on Internal Network state estimation data and External Network equivalence data, and according to real-time power flow data of equipment of Δt less than te (abbreviated as overload equipment in emergent state) and power flow data acquired in real-time by SCD.
Sub-step 3: optimal emergency control strategy calculation: for the grid after equivalence, using the objective function aimed to minimize total control cost and comprehensive indexes of proportions of load control quantity in different regions, taking into consideration of discreteness and cost of emergency control measures, and under the restraint of electric power system power flow, and based on mixed integer nonlinear programming algorithm, centralized decision-making optimization and successive approximation control strategy are adopted to calculate optimal emergency control strategy for equipment overload in emergent state. Its objective function is formula (6). In this formula, term 1 is generator emergency control cost, comprising two parts: power adjustment cost and shutdown cost. If generator i is tripped, then Gi is 1; otherwise it is 0. CGi is shutdown cost of generator i. Ge is total number of generators controlled by the SCD. Term 2 is the load control cost: Le is total number of loads controlled by the SCD. Term 3 reflects decentralized load control requirements: y is a set coefficient (greater than 1 and less than or equal to x). Meanings of other variables are the same as described earlier.
Corresponding restraint conditions are given in formula (7), including power flow equation restraint (including bus voltage restraint), equipment overload restraint, and emergency control measures space restraint to be selected. In this formula, Irj is long-term allowable current of equipment j estimated in step 2, Ij0 is current in equipment j at the most recent time point before this round of control strategy calculation, λe is set emergency control successive approximation coefficient, and Me is number of equipment for which Δt estimated in real-time in sub-step 2 is less than k′ete (k′e is greater than 1).
In
In
In
Number | Date | Country | Kind |
---|---|---|---|
201110445359.6 | Dec 2011 | CN | national |
This application claims priority to PCT Application No. PCT/CN2012/082158, having a filing date of Sep. 27, 2012, based off of CN Application No. 201110445359.6, having a filing date of Dec. 28, 2011, the entire contents of which are hereby incorporated by reference.
Filing Document | Filing Date | Country | Kind |
---|---|---|---|
PCT/CN2012/082158 | 9/27/2012 | WO | 00 |