CN105656024A - Electric power system operation standby safety margin online assessment method - Google Patents

Electric power system operation standby safety margin online assessment method Download PDF

Info

Publication number
CN105656024A
CN105656024A CN201510874414.1A CN201510874414A CN105656024A CN 105656024 A CN105656024 A CN 105656024A CN 201510874414 A CN201510874414 A CN 201510874414A CN 105656024 A CN105656024 A CN 105656024A
Authority
CN
China
Prior art keywords
power
meritorious
forecast failure
node
frequency
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN201510874414.1A
Other languages
Chinese (zh)
Other versions
CN105656024B (en
Inventor
徐泰山
汪际锋
唐卓尧
鲍颜红
黄河
苏寅生
郑亮
刘韶峰
任先成
周海锋
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
China Southern Power Grid Co Ltd
Nari Technology Co Ltd
Nanjing NARI Group Corp
Original Assignee
China Southern Power Grid Co Ltd
Nari Technology Co Ltd
Nanjing NARI Group Corp
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by China Southern Power Grid Co Ltd, Nari Technology Co Ltd, Nanjing NARI Group Corp filed Critical China Southern Power Grid Co Ltd
Priority to CN201510874414.1A priority Critical patent/CN105656024B/en
Publication of CN105656024A publication Critical patent/CN105656024A/en
Application granted granted Critical
Publication of CN105656024B publication Critical patent/CN105656024B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00Circuit arrangements for ac mains or ac distribution networks
    • HELECTRICITY
    • H02GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02JCIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]

Abstract

The invention discloses an electric power system operation standby safety margin online assessment method and belongs to the electric power system operation and control technology field. Through considering time domain simulation of protection and automatic safety device motion, a stable operation state of a faulted system is acquired. Based on a direct current trend equation of the state, active power sensitivity of an operation standby to power transmission equipment and a power transmission channel is calculated. Through solving and considering a controllable operation standby interval, a power transmission equipment overload limit and a power transmission channel safety and stability limit constraint, a maximum operation standby is taken as an optimization function of a target so as to acquire an operation standby capacity maximum value which can be used for controlling a steady-state frequency. Through combining a power frequency static state characteristic coefficient, a maximum controllable amount of the steady-state frequency is calculated and operation standby safety margin online assessment considering the fault is realized.

Description

Electric power system operation standby safety margin online evaluation method
Technical field
The invention belongs to operation and control of electric power system technical field, specifically, the present invention relates to a kind of considering security that is applicable to and stablize constraint to the online quantitative estimation method of frequency control ability of electric power system operation standby after forecast failure.
Background technology
After electric power system operation standby mainly occurs for forecast failure (comprising generator output fluctuation, fluctuation of load, short circuit, tripping operation and direct current locking etc.), safeguards system frequency quality requires to arrange, therefore, need constantly to monitor operation is for subsequent use in Operation of Electric Systems, and assess the virtual rating of its safeguards system frequency quality.
What existing operation in-service monitoring for subsequent use and appraisal procedure adopted is exerts oneself according to the total maximum generation of current running status, the difference of total minimum generated output and total load represents control power system mesomeric state frequency decrease respectively operation reserve capacity and the operation reserve capacity that control power system mesomeric state frequency rises, to run the ratio reflection operation safety margin for subsequent use that reserve capacity accounts for total load or total generated output.
Existing method the size unification that affects of power system power balance is considered with himself by forecast failure, there is the deficiency of two aspects: on the one hand, do not distinguish the impact difference that power system is caused by the identical different forecast failures of power variation, such as, after the different DC line bipolar lockings that transmission power is identical, the power shortage caused after may there is machine of cutting, cutting load because transient state is different with dynamic process is different, even if the power shortage caused is identical, it is also possible to the trend distribution caused is different; On the other hand, do not consider to run yet for subsequent use whether can be unavailable because of the constraint of transmission facility overload margin and passway for transmitting electricity safety and stability limit, there is operation reserve capacity that the method provides can not completely for controlling the operation condition of power system mesomeric state change of frequency. Therefore, adopt the method can not obtain true reflection to consider forecast failure, meet power system safety and stability and run the operation safety margin for subsequent use of constraint requirements.For this reason, need for different forecast failures, by considering the time-domain-simulation of first and second equipment dynamic behaviour and obtain the steady operational status after forecast failure, again based on this running status, consider that the constraint of transmission facility overload margin and passway for transmitting electricity safety and stability limit calculates actual available maximum operation reserve capacity, in conjunction with the upper and lower limit that power-frequency static characteristics coefficient and the systematic steady state frequency quality of this running status require, operation safety margin for subsequent use is carried out quantitative evaluation.
Summary of the invention
It is an object of the invention to: for existing operation in-service monitoring for subsequent use and appraisal procedure, the size unification that affects of power system power balance is considered by forecast failure with it, do not distinguish the impact difference that power system is caused by the identical different forecast failures of power variation, and do not consider the deficiency of the constraint of transmission facility overload margin and passway for transmitting electricity safety and stability limit, after considering forecast failure, power train is unified, the dynamic behaviour of second set, the constraint of transmission facility overload margin and passway for transmitting electricity safety and stability limit in adjustment for subsequent use is run under steady operational status after taking into account forecast failure, a kind of electric power system operation standby safety margin online evaluation method considering forecast failure is proposed, to meet the actual demand of management and running control.
The technology used in the present invention principle is: forecast failure just can be transitioned into stable running status after needing experience transient state, dynamic process after occurring; generator excitation/velocity control system, protection in transient process, peace may action from equipment such as devices; it is therefore necessary to carry out for forecast failure considering that the time-domain-simulation of first and second equipment dynamic behaviour just can obtain after forecast failure closest to actual steady operational status. The speed obtaining steady operational status to improve, in dynamic process, computing node frequency jitter amplitude is less than set(ting)value as criterion premature termination time-domain-simulation. Running the constraint of transmission facility overload margin and passway for transmitting electricity safety and stability limit in adjustment for subsequent use after taking into account forecast failure under steady operational status, the maximum controlled variable for subsequent use that runs calculated is only the operation reserve capacity really can implemented. In conjunction with the upper and lower limit that power-frequency static characteristics coefficient and the systematic steady state frequency quality of system under steady operational status after forecast failure require, can realize the electric power system operation standby safety margin online evaluation considering forecast failure.
Specifically, the present invention adopts following technical scheme to realize, and comprises the following steps:
1) based on the current running status of power system, carry out for forecast failure considering that protection and peace are from the time-domain-simulation of device action, until when emulation is all less than set(ting)value �� to the difference of adjacent maximum value in each computing node dynamic process of frequency and mnm., terminate time-domain-simulation, obtain dynamic process of frequency and each meritorious dynamic process injecting node and transmission facility of each computing node, enter step 2);
Described injection node comprises generator node, load node and equivalent injection node;
Described transmission facility comprises circuit and transformer;
2) for each computing node, respectively using step 1) maximum value at the end of node dynamic process of frequency and the mean value both mnm. as the steady frequency considering each computing node after forecast failure, then using the mean value of the steady frequency of all computing nodes as the steady frequency f considering power system after forecast failurea;
Inject node for each, respectively using step 1) inject the meritorious dynamic process of node at the end of maximum value and the mean value both mnm. as considering that the stable state of this injection node is meritorious after forecast failure;
For each transmission facility, respectively using step 1) maximum value at the end of the meritorious dynamic process of transmission facility and the mean value both mnm. is as considering that the stable state of this transmission facility is meritorious after forecast failure, the passway for transmitting electricity monitored for safety and stability for each again, respectively meritorious for the stable state of transmission facilities all in passway for transmitting electricity sum is gained merit as the stable state of this passway for transmitting electricity after considering forecast failure, enter step 3);
3) according to considering after forecast failure that stable state that each injects node is meritorious and the network topology of power system after forecast failure and model and parameter, form the DC power flow equation of power system, enter step 4);
4) based on the DC power flow equation of power system, for each meritorious controllable node that injection node pre-sets, calculate the sensitivity that each transmission facility in power system is gained merit by meritorious controllable node respectively, enter step 5);
5) if faIt is less than or equal to power system rated frequency frWith the difference of ��/2, then by solving the optimization function that formula (1) represents, after calculating forecast failure, can be used for improving the meritorious injection change total amount �� P of power system mesomeric state frequencyin, enter step 6);
If faIt is more than or equal to frWith ��/2 sum, then by solving the optimization function that formula (2) represents, after calculating forecast failure, can be used for reducing the meritorious injection change total amount �� P of power system mesomeric state frequencyde, enter step 7);
If faIt is greater than frWith the difference of ��/2 and faIt is less than frWith ��/2 sum, then respectively by the optimization function solving formula (1) and formula (2) represents, meritorious injections that can be used for improving power system mesomeric state frequency after calculating forecast failure changes total amount �� PinWith the meritorious injection change total amount �� P that can be used for reducing power system mesomeric state frequencyde, enter step 8);
In formula (1) and (2), I is the sum injecting the meritorious controllable node of node, �� PiFor the meritorious injection variable quantity of i-th node in meritorious controllable node, Pi.0For the meritorious injection rate of stable state of i-th node in controllable node meritorious after consideration forecast failure, Pi.maxFor meritorious injection rate of considering the stable state of i-th node in meritorious controllable node after forecast failure set time long �� t in the upper limit in controlled interval, Pi.minFor meritorious injection rate of considering the stable state of i-th node in meritorious controllable node after forecast failure set time long �� t in the lower limit in controlled interval, J is the passway for transmitting electricity sum monitored for safety and stability, KjFor the sum of transmission facility in the individual passway for transmitting electricity for safety and stability supervision of jth, sj.k.iFor the sensitivity that jth is gained merit by the meritorious injection of i-th node in meritorious controllable node for passway for transmitting electricity kth the transmission facility that safety and stability monitors, �� Pj.kFor the variable quantity that jth is gained merit for kth transmission facility in the passway for transmitting electricity of safety and stability supervision, Pj.k.0For jth after consideration forecast failure is gained merit for the stable state of kth transmission facility in the passway for transmitting electricity of safety and stability supervision, Pj.k.LMTFor the meritorious limit of stable state of kth transmission facility in the individual passway for transmitting electricity for safety and stability supervision of jth, �� PT.jFor the variable quantity that jth is gained merit for the passway for transmitting electricity that safety and stability monitors, PT.j.0For jth after consideration forecast failure is meritorious for the stable state of the passway for transmitting electricity of safety and stability supervision, PT.j.LMTFor the forward security stability limitation of the passway for transmitting electricity that jth after consideration forecast failure monitors for safety and stability, PT'.j.LMTFor the reverse safety and stability limit of the passway for transmitting electricity that jth after consideration forecast failure monitors for safety and stability, the set of the transmission facility except transmission facility in the passway for transmitting electricity being used for safety and stability supervision being designated as B, L is the sum of transmission facility in B, sb.l.iFor the sensitivity that B the l transmission facility is gained merit by the meritorious injection of i-th node in meritorious controllable node, �� Pb.lFor the variable quantity that the l transmission facility in B is meritorious, Pb.l.0For in B after consideration forecast failure, the stable state of l transmission facility is gained merit, Pb.l.LMTFor the meritorious limit of stable state of l transmission facility in B;
First, 6) according to step 5) meritorious injections that can be used for improving power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pin, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (3) and consider �� PinThe variation delta f of power system mesomeric state frequency after implementingin, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (4), method ends;
Δf i n = ΔP i n k p f - - - ( 3 )
η = f a + Δf i n - f d f r - f d × 100 % - - - ( 4 )
In formula, fdFor the lower limit in power system mesomeric state frequency security interval;
First, 7) according to step 5) meritorious injections that can be used for reducing power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pde, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (5) and consider �� PdeThe variation delta f of power system mesomeric state frequency after implementingde, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (6), method ends;
Δf d e = ΔP d e k p f - - - ( 5 )
η = f u - f a - Δf d e f u - f r × 100 % - - - ( 6 )
In formula, fuFor the upper limit in power system mesomeric state frequency security interval;
First, 8) according to step 5) meritorious injections that can be used for improving power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pin, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (7) and consider �� PinThe variation delta f of power system mesomeric state frequency after implementingin, according to step 5) and meritorious injections that can be used for reducing power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pde, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (8) and consider �� PdeThe variation delta f of power system mesomeric state frequency after implementingde, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (9), method ends;
Δf i n = ΔP i n k p f - - - ( 7 )
Δf d e = ΔP d e k p f - - - ( 8 )
η = min ( f a + Δf i n - f d f r - f d × 100 % , f u - f a - Δf d e f u - f r × 100 % ) - - - ( 9 )
In formula, fdFor the lower limit in power system mesomeric state frequency security interval, fuFor the upper limit in power system mesomeric state frequency security interval.
In order to improve operation efficiency, above-mentioned all processes can be processed by PC cluster method, namely each forecast failure concentrated for forecast failure, by the electric power system operation standby safety margin online evaluation using considering single forecast failure as an example, adopt the mode of PC cluster, by above-mentioned all processes, each example is carried out parallel processing, obtain considering respectively the electric power system operation standby safety margin of each forecast failure.
By adopting technique scheme, the present invention achieves following technique effect: by carrying out for forecast failure considering that the time-domain-simulation of first and second equipment dynamic behaviour can obtain after forecast failure closest to actual steady operational status, overcomes existing operation in-service monitoring for subsequent use and ignores transient state with appraisal procedure and dynamic process, foundation stable state analysis technology determine the steady operational status after forecast failure, unrealistic deficiency. In dynamic process, computing node frequency jitter amplitude is less than set(ting)value as criterion premature termination time-domain-simulation, it is to increase obtain the speed of steady operational status after forecast failure. By running the constraint of transmission facility overload margin and passway for transmitting electricity safety and stability limit in adjustment for subsequent use under steady operational status after taking into account forecast failure, the operation calculated maximum controlled variable for subsequent use, reflects the operation reserve capacity really can implemented. In conjunction with the upper and lower limit that power-frequency static characteristics coefficient and the systematic steady state frequency quality of system under steady operational status after forecast failure require, achieve the electric power system operation standby safety margin online evaluation considering forecast failure, meet the actual demand of management and running control.
Accompanying drawing explanation
Fig. 1 is the schema of the inventive method.
Embodiment
Below in conjunction with accompanying drawing, the present invention is described in further detail.
Step 1 in Fig. 1: based on the current running status of power system, carry out for forecast failure considering that protection and peace are from the time-domain-simulation of device action, until when emulation is all less than set(ting)value �� (being usually set to 0.02Hz) to the difference of adjacent maximum value in each computing node dynamic process of frequency and mnm., terminate time-domain-simulation, obtain dynamic process of frequency and each meritorious dynamic process injecting node and transmission facility of each computing node, enter step 2);
Described injection node comprises generator node, load node and equivalent injects node (such as equivalent generator node, etc. duty value node). Described transmission facility comprises circuit and transformer.
Step 2 in Fig. 1: for each computing node, respectively using step 1) maximum value at the end of node dynamic process of frequency and the mean value both mnm. as the steady frequency considering this computing node after forecast failure, then using the mean value of the steady frequency of all computing nodes as the steady frequency f considering power system after forecast failurea��
Inject node for each, respectively using step 1) inject the meritorious dynamic process of node at the end of maximum value and the mean value both mnm. as considering that the stable state of this injection node is meritorious after forecast failure.
For each transmission facility, respectively using step 1) at the end of the meritorious dynamic process of transmission facility maximum value and the mean value both mnm. as considering that the stable state of this transmission facility is meritorious after forecast failure, the passway for transmitting electricity monitored for safety and stability for each again, respectively meritorious for the stable state of transmission facilities all in passway for transmitting electricity sum is gained merit as the stable state of this passway for transmitting electricity after considering forecast failure, enter step 3).
Step 3 in Fig. 1: according to considering after forecast failure that stable state that each injects node is meritorious and the network topology of power system after forecast failure and model and parameter, form the DC power flow equation of power system, enter step 4).
Step 4 in Fig. 1: based on the DC power flow equation of power system, for each meritorious controllable node that injection node pre-sets, calculates the sensitivity that each transmission facility in power system is gained merit by meritorious controllable node respectively, enters step 5).
Step 5 in Fig. 1: if faIt is less than or equal to power system rated frequency frWith the difference of ��/2, then by solving the optimization function that formula (1) represents, after calculating forecast failure, can be used for improving the meritorious injection change total amount �� P of power system mesomeric state frequencyin, enter step 6).
If faIt is more than or equal to frWith ��/2 sum, then by solving the optimization function that formula (2) represents, after calculating forecast failure, can be used for reducing the meritorious injection change total amount �� P of power system mesomeric state frequencyde, enter step 7).
If faIt is greater than frWith the difference of ��/2 and faIt is less than frWith ��/2 sum, then respectively by the optimization function solving formula (1) and formula (2) represents, meritorious injections that can be used for improving power system mesomeric state frequency after calculating forecast failure changes total amount �� PinWith the meritorious injection change total amount �� P that can be used for reducing power system mesomeric state frequencyde, enter step 8);
In formula (1) and (2), I is the sum injecting the meritorious controllable node of node, �� PiFor the meritorious injection variable quantity of i-th node in meritorious controllable node, Pi.0For the meritorious injection rate of stable state of i-th node in controllable node meritorious after consideration forecast failure, Pi.maxFor meritorious injection rate of considering the stable state of i-th node in meritorious controllable node after forecast failure set time long �� t (being usually set to 3 minutes) in the upper limit in controlled interval, Pi.minFor meritorious injection rate of considering the stable state of i-th node in meritorious controllable node after forecast failure set time long �� t in the lower limit in controlled interval, J is the passway for transmitting electricity sum monitored for safety and stability, KjFor the sum of transmission facility in the individual passway for transmitting electricity for safety and stability supervision of jth, sj.k.iFor the sensitivity that jth is gained merit by the meritorious injection of i-th node in meritorious controllable node for passway for transmitting electricity kth the transmission facility that safety and stability monitors, �� Pj.kFor the variable quantity that jth is gained merit for kth transmission facility in the passway for transmitting electricity of safety and stability supervision, Pj.k.0For jth after consideration forecast failure is gained merit for the stable state of kth transmission facility in the passway for transmitting electricity of safety and stability supervision, Pj.k.LMTFor the meritorious limit of stable state of kth transmission facility in the individual passway for transmitting electricity for safety and stability supervision of jth, �� PT.jFor the variable quantity that jth is gained merit for the passway for transmitting electricity that safety and stability monitors, PT.j.0For jth after consideration forecast failure is meritorious for the stable state of the passway for transmitting electricity of safety and stability supervision, PT.j.LMTFor the forward security stability limitation of the passway for transmitting electricity that jth after consideration forecast failure monitors for safety and stability, P'T.j.LMTFor the reverse safety and stability limit of the passway for transmitting electricity that jth after consideration forecast failure monitors for safety and stability, the set of the transmission facility except transmission facility in the passway for transmitting electricity being used for safety and stability supervision being designated as B, L is the sum of transmission facility in B, sb.l.iFor the sensitivity that B the l transmission facility is gained merit by the meritorious injection of i-th node in meritorious controllable node, �� Pb.lFor the variable quantity that the l transmission facility in B is meritorious, Pb.l.0For in B after consideration forecast failure, the stable state of l transmission facility is gained merit, Pb.l.LMTFor the meritorious limit of stable state of l transmission facility in B.
If needing the meritorious constraint considering junctor section between region in the active power adjustment controlled for steady frequency, then junctor section is processed as the passway for transmitting electricity being used for safety and stability supervision, wherein, by PT.j.LMTSubstitute by the meritorious upper limit of this junctor section, by P'T.j.LMTSubstituting with the contrary number of the meritorious lower limit of this junctor section, other parameter constant, can realize.
Wherein, the method of calculation of the meritorious limit of the stable state of transmission facility are as follows: if the ratio of transmission facility electric current and electric current limit under the current running status of power system is greater than set(ting)value (being usually set to 10%), the terminal voltage then remained under the current running status of power system with transmission facility is constant and power factor does not turn into condition, calculate gaining merit when transmission facility electric current reaches electric current limit, the stable state of transmission facility is it can be used as to gain merit limit, otherwise, then constant and under the current running status of power system the terminal voltage of power factor (being usually set to 0.95) of setting is kept not turn into condition with transmission facility, calculate gaining merit when transmission facility electric current reaches electric current limit, the stable state of transmission facility is it can be used as to gain merit limit.
Step 6 in Fig. 1: first, according to step 5) meritorious injections that can be used for improving power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pin, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf(when power system voltage is constant, system always injects meritorious variable quantity and the ratio of systematic steady state frequency variation), is calculated by formula (3) and considers �� PinThe variation delta f of power system mesomeric state frequency after implementingin, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (4), method ends:
Δf i n = ΔP i n k p f - - - ( 3 )
η = f a + Δf i n - f d f r - f d × 100 % - - - ( 4 )
In formula, fdFor the lower limit in power system mesomeric state frequency security interval.
Step 7 in Fig. 1: first, according to step 5) meritorious injections that can be used for reducing power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pde, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (5) and consider �� PdeThe variation delta f of power system mesomeric state frequency after implementingde, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (6), method ends:
Δf d e = ΔP d e k p f - - - ( 5 )
η = f u - f a - Δf d e f u - f r × 100 % - - - ( 6 )
In formula, fuFor the upper limit in power system mesomeric state frequency security interval;
Step 8 in Fig. 1: first, according to step 5) meritorious injections that can be used for improving power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pin, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (7) and consider �� PinThe variation delta f of power system mesomeric state frequency after implementingin, according to step 5) and meritorious injections that can be used for reducing power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pde, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (8) and consider �� PdeThe variation delta f of power system mesomeric state frequency after implementingde, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (9), method ends:
Δf i n = ΔP i n k p f - - - ( 7 )
Δf d e = ΔP d e k p f - - - ( 8 )
η = min ( f a + Δf i n - f d f r - f d × 100 % , f u - f a - Δf d e f u - f r × 100 % ) - - - ( 9 )
In formula, fdFor the lower limit in power system mesomeric state frequency security interval, fuFor the upper limit in power system mesomeric state frequency security interval.
In order to improve operation efficiency, above-mentioned all processes can be processed by PC cluster method, namely each forecast failure concentrated for forecast failure, by the electric power system operation standby safety margin online evaluation using considering single forecast failure as an example, adopt the mode of PC cluster, by above-mentioned all processes, each example is carried out parallel processing, obtain considering respectively the electric power system operation standby safety margin of each forecast failure.
Although the present invention is with better embodiment openly as above, but embodiment is not for limiting the present invention's. Without departing from the spirit and scope of the invention, any equivalence change done or retouching, belong to the protection domain of the present invention equally. Therefore the content that protection scope of the present invention should define taking the claim of the application is as standard.

Claims (3)

1. electric power system operation standby safety margin online evaluation method, it is characterised in that, comprise the steps:
1) based on the current running status of power system, carry out for forecast failure considering that protection and peace are from the time-domain-simulation of device action, until when emulation is all less than set(ting)value �� to the difference of adjacent maximum value in each computing node dynamic process of frequency and mnm., terminate time-domain-simulation, obtain the dynamic process of frequency of each computing node and each meritorious dynamic process injecting node and each transmission facility, enter step 2);
Described injection node comprises generator node, load node and equivalent injection node;
Described transmission facility comprises circuit and transformer.
2) for each computing node, respectively using step 1) maximum value at the end of node dynamic process of frequency and the mean value both mnm. as the steady frequency considering each computing node after forecast failure, then using the mean value of the steady frequency of all computing nodes as the steady frequency f considering power system after forecast failurea;
Inject node for each, respectively using step 1) inject the meritorious dynamic process of node at the end of maximum value and the mean value both mnm. as considering that the stable state of this injection node is meritorious after forecast failure;
For each transmission facility, respectively using step 1) maximum value at the end of the meritorious dynamic process of transmission facility and the mean value both mnm. is as considering that the stable state of this transmission facility is meritorious after forecast failure, the passway for transmitting electricity monitored for safety and stability for each again, respectively meritorious for the stable state of transmission facilities all in passway for transmitting electricity sum is gained merit as the stable state of this passway for transmitting electricity after considering forecast failure, enter step 3);
3) according to considering after forecast failure that stable state that each injects node is meritorious and the network topology of power system after forecast failure and model and parameter, form the DC power flow equation of power system, enter step 4);
4) based on the DC power flow equation of power system, for each meritorious controllable node that injection node pre-sets, calculate the sensitivity that each transmission facility in power system is gained merit by meritorious controllable node respectively, enter step 5);
5) if faIt is less than or equal to power system rated frequency frWith the difference of ��/2, then by solving the optimization function that formula (1) represents, after calculating forecast failure, can be used for improving the meritorious injection change total amount �� P of power system mesomeric state frequencyin, enter step 6);
If faIt is more than or equal to frWith ��/2 sum, then by solving the optimization function that formula (2) represents, after calculating forecast failure, can be used for reducing the meritorious injection change total amount �� P of power system mesomeric state frequencyde, enter step 7);
If faIt is greater than frWith the difference of ��/2 and faIt is less than frWith ��/2 sum, then respectively by the optimization function solving formula (1) and formula (2) represents, meritorious injections that can be used for improving power system mesomeric state frequency after calculating forecast failure changes total amount �� PinWith the meritorious injection change total amount �� P that can be used for reducing power system mesomeric state frequencyde, enter step 8);
In formula (1) and (2), I is the sum injecting the meritorious controllable node of node, �� PiFor the meritorious injection variable quantity of i-th node in meritorious controllable node, Pi.0For the meritorious injection rate of stable state of i-th node in controllable node meritorious after consideration forecast failure, Pi.maxFor meritorious injection rate of considering the stable state of i-th node in meritorious controllable node after forecast failure set time long �� t in the upper limit in controlled interval, Pi.minFor meritorious injection rate of considering the stable state of i-th node in meritorious controllable node after forecast failure set time long �� t in the lower limit in controlled interval, J is the passway for transmitting electricity sum monitored for safety and stability, KjFor the sum of transmission facility in the individual passway for transmitting electricity for safety and stability supervision of jth, sj.k.iFor the sensitivity that jth is gained merit by the meritorious injection of i-th node in meritorious controllable node for passway for transmitting electricity kth the transmission facility that safety and stability monitors, �� Pj.kFor the variable quantity that jth is gained merit for kth transmission facility in the passway for transmitting electricity of safety and stability supervision, Pj.k.0For jth after consideration forecast failure is gained merit for the stable state of kth transmission facility in the passway for transmitting electricity of safety and stability supervision, Pj.k.LMTFor the meritorious limit of stable state of kth transmission facility in the individual passway for transmitting electricity for safety and stability supervision of jth, �� PT.jFor the variable quantity that jth is gained merit for the passway for transmitting electricity that safety and stability monitors, PT.j.0For jth after consideration forecast failure is meritorious for the stable state of the passway for transmitting electricity of safety and stability supervision, PT.j.LMTFor the forward security stability limitation of the passway for transmitting electricity that jth after consideration forecast failure monitors for safety and stability, P 'T.j.LMTFor the reverse safety and stability limit of the passway for transmitting electricity that jth after consideration forecast failure monitors for safety and stability, the set of the transmission facility except transmission facility in the passway for transmitting electricity being used for safety and stability supervision being designated as B, L is the sum of transmission facility in B, sb.l.iFor the sensitivity that B the l transmission facility is gained merit by the meritorious injection of i-th node in meritorious controllable node, �� Pb.lFor the variable quantity that the l transmission facility in B is meritorious, Pb.l.0For in B after consideration forecast failure, the stable state of l transmission facility is gained merit, Pb.l.LMTFor the meritorious limit of stable state of l transmission facility in B;
First, 6) according to step 5) meritorious injections that can be used for improving power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pin, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (3) and consider �� PinThe variation delta f of power system mesomeric state frequency after implementingin, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (4), method ends;
Δf i n = ΔP i n k p f - - - ( 3 )
η = f a + Δf i n - f d f r - f d × 100 % - - - ( 4 )
In formula, fdFor the lower limit in power system mesomeric state frequency security interval;
First, 7) according to step 5) meritorious injections that can be used for reducing power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pde, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (5) and consider �� PdeThe variation delta f of power system mesomeric state frequency after implementingde, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (6), method ends;
Δf d e = ΔP d e k p f - - - ( 5 )
η = f u - f a - Δf d e f u - f r × 100 % - - - ( 6 )
In formula, fuFor the upper limit in power system mesomeric state frequency security interval;
First, 8) according to step 5) meritorious injections that can be used for improving power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pin, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (7) and consider �� PinThe variation delta f of power system mesomeric state frequency after implementingin, according to step 5) and meritorious injections that can be used for reducing power system mesomeric state frequency after the forecast failure that calculates change total amount �� Pde, and the power-frequency static characteristics coefficient k of power system after consideration forecast failurepf, calculated by formula (8) and consider �� PdeThe variation delta f of power system mesomeric state frequency after implementingde, then calculated the electric power system operation standby safety margin �� considering forecast failure by formula (9), method ends;
Δf i n = ΔP i n k p f - - - ( 7 )
Δf d e = ΔP d e k p f - - - ( 8 )
η = min ( f a + Δf i n - f d f r - f d × 100 % , f u - f a - Δf d e f u - f r × 100 % ) - - - ( 9 )
In formula, fdFor the lower limit in power system mesomeric state frequency security interval, fuFor the upper limit in power system mesomeric state frequency security interval.
2. the online cluster appraisal procedure of electric power system operation standby safety margin according to claim 1, it is characterized in that, in step 5) in, if needing the meritorious constraint considering junctor section between region in the active power adjustment controlled for steady frequency, then being processed as the passway for transmitting electricity being used for safety and stability supervision by junctor section, it is by PT.j.LMTSubstitute by the meritorious upper limit of this junctor section, by P 'T.j.LMTSubstitute with the contrary number of the meritorious lower limit of this junctor section.
3. the online cluster appraisal procedure of electric power system operation standby safety margin according to claim 1 and 2, it is characterized in that, for each forecast failure that forecast failure is concentrated, by the electric power system operation standby safety margin online evaluation using considering single forecast failure as an example, adopt the mode of PC cluster, by method according to claim 1, each example is carried out parallel processing, obtain considering respectively the electric power system operation standby safety margin of each forecast failure.
CN201510874414.1A 2015-12-02 2015-12-02 Electric power system operation standby margin of safety online evaluation method Active CN105656024B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201510874414.1A CN105656024B (en) 2015-12-02 2015-12-02 Electric power system operation standby margin of safety online evaluation method

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
CN201510874414.1A CN105656024B (en) 2015-12-02 2015-12-02 Electric power system operation standby margin of safety online evaluation method
JP2018548255A JP6569014B2 (en) 2015-12-02 2016-10-28 Online evaluation method for safety margin of power system operating reserve
PCT/CN2016/103729 WO2017092527A1 (en) 2015-12-02 2016-10-28 Online evaluation method of safety margin of operating reserve in power system

Publications (2)

Publication Number Publication Date
CN105656024A true CN105656024A (en) 2016-06-08
CN105656024B CN105656024B (en) 2018-05-15

Family

ID=56482199

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201510874414.1A Active CN105656024B (en) 2015-12-02 2015-12-02 Electric power system operation standby margin of safety online evaluation method

Country Status (3)

Country Link
JP (1) JP6569014B2 (en)
CN (1) CN105656024B (en)
WO (1) WO2017092527A1 (en)

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017092527A1 (en) * 2015-12-02 2017-06-08 国电南瑞科技股份有限公司 Online evaluation method of safety margin of operating reserve in power system
CN107633271A (en) * 2017-09-30 2018-01-26 电子科技大学 The non-precision method for calculating probability of power system mesomeric state availability
CN108512219A (en) * 2018-03-14 2018-09-07 国家电网公司华中分部 A kind of multichannel receiving-end system under Voltage Stability Constraints is by electric energy power comprehensive estimation method
CN108900368A (en) * 2018-04-12 2018-11-27 广东电网有限责任公司佛山供电局 A kind of regulator control system and peace are from bi-directional data proofreading method between system
CN109696901A (en) * 2018-11-30 2019-04-30 红云红河烟草(集团)有限责任公司 A kind of evaluation of cigarette packet equipment running status and prediction technique
CN109995054A (en) * 2019-03-26 2019-07-09 国网浙江省电力有限公司电力科学研究院 A kind of mains frequency safety opening terminal allowance calculation method and starting-up method
CN111555269A (en) * 2020-04-14 2020-08-18 南方电网科学研究院有限责任公司 Analysis method and device for power grid reserve capacity, terminal equipment and storage medium
CN111555269B (en) * 2020-04-14 2021-12-14 南方电网科学研究院有限责任公司 Analysis method and device for power grid reserve capacity, terminal equipment and storage medium

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120161520A1 (en) * 2010-12-28 2012-06-28 Central Japan Railway Company Double loop structure
CN103401232A (en) * 2013-04-24 2013-11-20 南京南瑞集团公司 Cluster computing method for online emergency control strategy for transient safety stability of power system
CN103514512A (en) * 2013-09-11 2014-01-15 南京南瑞集团公司 Fast auxiliary decision-making method for electric power system on-line isolated island grid connection

Family Cites Families (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPH10290528A (en) * 1997-04-14 1998-10-27 Hitachi Ltd Voltage-reactive power monitoring controlling device
CN102623990B (en) * 2011-08-09 2014-04-30 国电南瑞科技股份有限公司 Low-frequency deloading parameter online checking method of self-adaptation power grid operation mode
CN103715721B (en) * 2013-07-31 2016-02-17 南京南瑞集团公司 The real-time peak regulating method of alternating current-direct current mixing power grid joint
CN103762590B (en) * 2014-01-13 2017-06-30 国家电网公司 Power system UFLS takes turns off-load amount on-line tuning method substantially
CN104578048B (en) * 2014-08-11 2017-01-25 国家电网公司 Method for quickly evaluating transient security and stability of enumerated combined successive fault set on line
CN105656024B (en) * 2015-12-02 2018-05-15 国电南瑞科技股份有限公司 Electric power system operation standby margin of safety online evaluation method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20120161520A1 (en) * 2010-12-28 2012-06-28 Central Japan Railway Company Double loop structure
CN103401232A (en) * 2013-04-24 2013-11-20 南京南瑞集团公司 Cluster computing method for online emergency control strategy for transient safety stability of power system
CN103514512A (en) * 2013-09-11 2014-01-15 南京南瑞集团公司 Fast auxiliary decision-making method for electric power system on-line isolated island grid connection

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2017092527A1 (en) * 2015-12-02 2017-06-08 国电南瑞科技股份有限公司 Online evaluation method of safety margin of operating reserve in power system
CN107633271A (en) * 2017-09-30 2018-01-26 电子科技大学 The non-precision method for calculating probability of power system mesomeric state availability
CN107633271B (en) * 2017-09-30 2020-10-30 电子科技大学 Method for calculating inaccurate probability of steady-state availability of power system
CN108512219A (en) * 2018-03-14 2018-09-07 国家电网公司华中分部 A kind of multichannel receiving-end system under Voltage Stability Constraints is by electric energy power comprehensive estimation method
CN108512219B (en) * 2018-03-14 2020-08-28 国家电网公司华中分部 Method for evaluating power receiving capacity of multi-channel receiving end system under voltage stability constraint
CN108900368B (en) * 2018-04-12 2019-10-29 广东电网有限责任公司佛山供电局 A kind of regulator control system and peace are from bi-directional data proofreading method between system
CN108900368A (en) * 2018-04-12 2018-11-27 广东电网有限责任公司佛山供电局 A kind of regulator control system and peace are from bi-directional data proofreading method between system
CN109696901A (en) * 2018-11-30 2019-04-30 红云红河烟草(集团)有限责任公司 A kind of evaluation of cigarette packet equipment running status and prediction technique
CN109995054A (en) * 2019-03-26 2019-07-09 国网浙江省电力有限公司电力科学研究院 A kind of mains frequency safety opening terminal allowance calculation method and starting-up method
CN111555269A (en) * 2020-04-14 2020-08-18 南方电网科学研究院有限责任公司 Analysis method and device for power grid reserve capacity, terminal equipment and storage medium
CN111555269B (en) * 2020-04-14 2021-12-14 南方电网科学研究院有限责任公司 Analysis method and device for power grid reserve capacity, terminal equipment and storage medium

Also Published As

Publication number Publication date
WO2017092527A1 (en) 2017-06-08
CN105656024B (en) 2018-05-15
JP6569014B2 (en) 2019-08-28
JP2018536382A (en) 2018-12-06

Similar Documents

Publication Publication Date Title
CN105656024A (en) Electric power system operation standby safety margin online assessment method
Bonfiglio et al. Design and implementation of a variable synthetic inertia controller for wind turbine generators
CN106099906B (en) A kind of power system security stablizes online integrated campaign method
CN103248055B (en) A kind of solve regional power grid off-the-line after correct choosing cut the method for generator and load
Suvire et al. Wind farm: Dynamic model and impact on a weak power system
Teymouri et al. Cyber security risk assessment of solar PV units with reactive power capability
CN106154084A (en) Network load exception and operation risk real-time monitoring and early warning method
Caldognetto et al. Improving microgrid performance by cooperative control of distributed energy sources
CN102412575B (en) Simulation evaluation method for configuration scheme of safe automatic device of transmitting-end power grid high-frequency cutting machine
CN104934971A (en) Dynamic section control method based on power flow transfer ratio
CN102623990A (en) Low-frequency deloading parameter online checking method of self-adaptation power grid operation mode
CN106469908A (en) The static security synthesis aid decision-making method that a kind of multiple faults is coordinated
Parniani et al. SCADA based under frequency load shedding integrated with rate of frequency decline
Berizzi et al. Architecture and functionalities of a smart Distribution Management System
CN108400597B (en) Direct-current fault static safety analysis method considering security control and primary frequency modulation characteristics
CN105900308B (en) The method for keeping electrical network stable by load break
van der Meer et al. Stability assessment of VSC-HVDC connected large-scale offshore wind power: A North-Sea region case study
CN104505855A (en) Online evaluation system for grid-connected characteristic of wind power plant
Kaur et al. Review and analysis of voltage collapse in power system
CN111525598A (en) Optimal generator tripping capacity obtaining method for inhibiting high-frequency problem of transmission-end power grid
Hassan et al. Under-frequency load shedding: Towards a smarter smart house with a consumer level controller
Gunadin et al. Analysis of Transmission Line Stability for Sulselbar Interconnection System with the Penetration of Renewable Energy to Prevent Voltage Collapse
CN104253433B (en) A kind of method reducing massive blackout accident
Fan et al. Steady frequency prediction algorithm for power system with governor deadband
Silva et al. Control architectures to perform voltage regulation on low voltage networks using DG

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant