CN103065060B - A kind of computing method of transient stability limit of transmission section of power system - Google Patents

A kind of computing method of transient stability limit of transmission section of power system Download PDF

Info

Publication number
CN103065060B
CN103065060B CN201310034636.3A CN201310034636A CN103065060B CN 103065060 B CN103065060 B CN 103065060B CN 201310034636 A CN201310034636 A CN 201310034636A CN 103065060 B CN103065060 B CN 103065060B
Authority
CN
China
Prior art keywords
fault
limit
section
stability limit
planes
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.)
Active
Application number
CN201310034636.3A
Other languages
Chinese (zh)
Other versions
CN103065060A (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.)
NORTHEAST GRID CO Ltd
Tsinghua University
Original Assignee
NORTHEAST GRID CO Ltd
Tsinghua University
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 NORTHEAST GRID CO Ltd, Tsinghua University filed Critical NORTHEAST GRID CO Ltd
Priority to CN201310034636.3A priority Critical patent/CN103065060B/en
Publication of CN103065060A publication Critical patent/CN103065060A/en
Application granted granted Critical
Publication of CN103065060B publication Critical patent/CN103065060B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The computing method that the present invention proposes transient stability limit of transmission section of power system comprise: electrical network is divided into power transmission district and powered district by S1, increase power transmission district and exert oneself, reduce powered district and exert oneself, increase section through-put power; S2, for the change scheme of exerting oneself of setting, obtains transient stability limit, the limit method of operation and controls fault; S3, under limit mode, to control fault, carries out the identification of a unstability group of planes and a stable group of planes; S4 generates new change scheme of exerting oneself, and the unit increase in a unstability group of planes is exerted oneself or will not start shooting unit start, is reduced by the unit in a stable group of planes and exerts oneself to shutdown; S5, for new change scheme of exerting oneself, obtains transient stability limit, the limit method of operation and controls fault; If the control fault of S6 S5 is identical with S2, terminate, with the transient stability limit of S5 for result, otherwise turn S3, regenerate and exert oneself change scheme and search for transient stability limit.The present invention has that accuracy is high, the advantage of good stability.

Description

A kind of computing method of transient stability limit of transmission section of power system
Technical field
The present invention relates to Power System Analysis technical field, particularly a kind of computing method of transient stability limit of transmission section of power system.
Background technology
Transient stability is the important content of power system safety and stability.Operation of power networks department needs to ensure that electrical network can keep transient stability under forecast failure.In operation of power networks, generally ensure the transient stability of system by controlling transmission cross-section power.Transmission cross-section is the set of some branch road in electrical network, and electrical network is divided into two disconnected regions after disconnecting by these branch roads, and the power on transmission cross-section represents the power of two interregional transmission.When the power that transmission cross-section transmits is excessive, system may under certain fault transient state unstable, therefore, there is a transient stability limit in transmission cross-section, when section power is less than the limit, for all faults that forecast failure is concentrated, system is transient stability, when section power be increased to be greater than the limit time, for forecast failure concentrate some faults, system first transient state is unstable, and this fault is called control fault.In operation of power networks, control section power is no more than transient stability limit, can ensure the transient stability of system.
The conventional Calculation Method of section transient stability limit is: based on a certain method of operation, keep the generator output in load and other regions constant, this section power transmission district of corresponding increase and decrease and powered district unit, and stability Calculation is carried out to the fault that forecast failure is concentrated, be met the section limit of power system transient stability.How to increase and decrease unit (namely generator output changes scheme) to have a significant impact the result of calculation tool of final transient stability limit.When generator output changes scheme difference, the stability limit obtained is likely different, this makes the result of calculation of stability limit depend on artificial setting in advance, can not ensure to obtain the stability limit under most serious conditions (stability limit minimum generator output shift gears), namely the most conservative stability limit.This security for guarantee electrical network is disadvantageous.Therefore wish to find the most dangerous generator output to change scheme, the transient stability limit obtained under this scenario is minimum, namely the most conservative, to ensure the security of electrical network.
Summary of the invention
The present invention one of is intended to solve the problems of the technologies described above at least to a certain extent or at least provides a kind of useful business to select.For this reason, the object of the invention is to the computing method proposing a kind of transient stability limit of transmission section of power system, the method can ensure the conservative property of transient stability limit result of calculation, ensures power grid security.
To achieve these goals, computing method according to the transient stability limit of transmission section of power system of the embodiment of the present invention comprise: S1: at analyzed transmission cross-section place, electrical network is divided into power transmission district and powered district, set a kind of initial generator output and change scheme, increase described power transmission district generated output, equivalent reduces described powered district generated output, to increase section through-put power; S2: the generator output for setting changes scheme, the transient stability limit of search section, obtains the transient stability limit of section, the limit method of operation and controls fault accordingly; S3: under the described limit method of operation, for control fault, carries out the identification of a transient stability unstability group of planes and a stable group of planes, obtains a unstability group of planes and a stable group of planes; S4: generate new generator output and change scheme, the unit increase in a described unstability group of planes is exerted oneself or started shooting successively by the unit of not starting shooting in its region, reduces the unit in a stable group of planes accordingly and exerts oneself until shut down; S5: scheme is changed for new generator output, the described transient stability limit of search section, obtains the transient stability limit of section, the limit method of operation and controls fault accordingly; S6: if the control fault phase that obtains of the control fault that obtains of step S5 and step S2 with, then calculate end, the transient stability limit obtained with step S5 is the final section limit, otherwise goes to step S3, regenerates generator output and changes scheme and search for transient stability limit.
In an embodiment of the present invention, the method of searching for section stability limit in described step S2 comprises: S201: the 1st fault concentrated forecast failure, adopt dichotomy to carry out complete stability limit search, obtain the section stability limit that this fault is corresponding, controlling fault is this fault; S202: to the limit method of operation corresponding to stability limit described in step S201, successively scan consequent malfunction i, wherein i be fault sequence number and i be greater than 1 positive integer; S203: if system stability under fault i, then section stability limit, the limit method of operation and control fault are all constant, continue next fault i+1 of scanning; S204: if system is unstable under fault i, then dichotomy is re-used to this fault and carry out stability limit search, obtain new section stability limit, the limit method of operation, control fault and become fault i, then to the new limit method of operation, scanning consequent malfunction i+1 is continued; S205: when scanning last fault, just obtains the control fault resolving face stability limit and correspondence.
In an embodiment of the present invention, described employing dichotomy carries out the search of complete stability limit and comprises: S2011: under the basic method of operation, system is stable to setting fault, is designated as stationary mode S a, now section through-put power is P a; With larger step-size change power transmission district and powered district generator output, make system unstable under setting fault, be designated as unstable mode S b, now section through-put power is designated as P b; S2012: if P b-P abe less than the threshold value preset, then go to step S2014, otherwise change power transmission district and powered district generator output variable quantity are original half, regenerate the new method of operation and make section through-put power for (P a+ P b)/2, which is designated as S c, in this mode time-domain-simulation is carried out to setting fault; S2013: if system stability, then stationary mode S abecome mode S c, P abecome (P a+ P b)/2, unstable mode is constant, returns step S2012 and continues to calculate; S2014: if system is unstable, then unstable mode S bbecome mode S c, P bbecome (P a+ P b)/2, stationary mode is constant, returns step S2012 and continues to calculate; S2015: real stability limit is positioned at P aand P bbetween, note P afor final section stability limit.
In an embodiment of the present invention, identify in described step S3 that the method for a unstability group of planes and a stable group of planes comprises: S301: under the limit method of operation that step S2 obtains, time-domain-simulation is carried out to control fault, obtain neutrality track, then slightly the fault clearance time is increased, make system transient modelling unstable, emulation obtains critical unstable track; S302: estimate controlling unstable equilibrium point according to described neutrality track and described critical unstable track, concrete grammar is for from the fault clearance moment, search for first kinetic energy minimum point on described neutrality track and described critical unstable track respectively, then using the mean value of 2 as controlling unstable equilibrium point; S303: under calculating this fault, energy margin is to the sensitivity of each generator output, and computing formula is as follows:
S g , i = ( δ cl , i - δ uep , i ) - ω cl , i t cl + 1 2 M i ( P m , i - P ecl , i ) t cl 2
Wherein, S g, irepresent that energy margin is to the sensitivity of the mechanical output of generator i, δ cl, irepresent that generator i is at the merit angle in fault clearance moment under center of inertia coordinate, δ uep, irepresent that generator i is at the merit angle at described controlling unstable equilibrium point place under center of inertia coordinate, ω cl, irepresent that generator i is at the rotating speed in fault clearance moment under center of inertia coordinate, t clrepresent to the time of fault clearance from fault, M irepresent the moment of inertia of generator i, P m,irepresent the mechanical output of generator i, P ecl, ithe electromagnetic power of expression generator i moment after fault clearance.S304: according to the result of sensitivity, hive off to unit, the minimum a part of unit of sensitivity is a unstability group of planes, and the maximum a part of unit of sensitivity is for stablize a group of planes, and in a group of planes, the quantity of unit is determined according to concrete electrical network.
According to transient stability limit of transmission section of power system computing method of the present invention, become the most dangerous generator output to change scheme based on a unstability group of planes with stabilizing machine all living creatures, ensure the conservative property of transient stability limit result of calculation, ensure power grid security.Described method improves the analysis of electrical network and runs control technology level, and operation of power networks personnel accurately can grasp stability limit and the stability margin of each transmission cross-section, can ensure the transient stability of system when carrying out method of operation adjustment.
Additional aspect of the present invention and advantage will part provide in the following description, and part will become obvious from the following description, or be recognized by practice of the present invention.
Accompanying drawing explanation
Of the present invention and/or additional aspect and advantage will become obvious and easy understand from the following description of the accompanying drawings of embodiments, wherein:
Fig. 1 is the process flow diagram of the computing method of transient stability limit of transmission section of power system according to the embodiment of the present invention;
Fig. 2 is the process flow diagram of the method for the search section stability limit of the computing method of transient stability limit of transmission section of power system according to the embodiment of the present invention;
Fig. 3 is the process flow diagram of the dichotomy of the search stability limit of the computing method of transient stability limit of transmission section of power system according to the embodiment of the present invention; And
Fig. 4 is an identification unstability group of planes for computing method for transient stability limit of transmission section of power system according to the embodiment of the present invention and the process flow diagram of the method for a stable group of planes.
Embodiment
Be described below in detail embodiments of the invention, the example of described embodiment is shown in the drawings, and wherein same or similar label represents same or similar element or has element that is identical or similar functions from start to finish.Being exemplary below by the embodiment be described with reference to the drawings, only for explaining the present invention, and can not limitation of the present invention being interpreted as.On the contrary, embodiments of the invention comprise fall into attached claims spirit and intension within the scope of all changes, amendment and equivalent.
In describing the invention, it is to be appreciated that term " first ", " second " etc. are only for describing object, and instruction or hint relative importance can not be interpreted as.In describing the invention, it should be noted that, unless otherwise clearly defined and limited, term " is connected ", " connection " should be interpreted broadly, such as, can be fixedly connected with, also can be removably connect, or connect integratedly; Can be mechanical connection, also can be electrical connection; Can be directly be connected, also indirectly can be connected by intermediary.For the ordinary skill in the art, concrete condition above-mentioned term concrete meaning in the present invention can be understood.In addition, in describing the invention, except as otherwise noted, the implication of " multiple " is two or more.
For the computing method of transient stability limit of transmission section of power system according to the embodiment of the present invention, first set a kind of initial generator output and change scheme, then under the limit method of operation, carry out the identification of a unstability group of planes and a stable group of planes, and generate new generator output change scheme, an increase unstability group of planes is exerted oneself and is reduced to stablize a group of planes and exerts oneself, exert oneself and reduce to stablize a group of planes to exert oneself the transient stability of reduction system that can be the fastest owing to increasing a unstability group of planes, therefore change scheme with initial generator output to compare, new scheme is a kind of more dangerous scheme, under the program, the stability limit of section is less, by changing the iteration adjustment of scheme to generator output, ensure the conservative property of section stability limit.
As shown in Figure 1, the embodiment of the present invention comprises the following steps:
S1: electrical network is divided into two parts at analyzed transmission cross-section place, is respectively power transmission district and powered district, sets a kind of initial generator output and changes scheme, increase power transmission district generated output, reduce powered district generated output, to increase section through-put power.
S2: the generator output for setting changes scheme, the transient stability limit of search section, obtains the transient stability limit of section, the limit method of operation and controls fault accordingly.
Particularly, as shown in Figure 2, the method for searching for section stability limit is further comprising the steps:
S201: the 1st fault concentrated forecast failure, adopt dichotomy to carry out complete stability limit search, obtain the section stability limit that this fault is corresponding, controlling fault is this fault.
S202: to the limit method of operation corresponding to step S201 stability limit, successively scan consequent malfunction i, wherein i be fault sequence number and i be greater than 1 positive integer.
S203: if system stability under fault i, then section stability limit, the limit method of operation and control fault are all constant, continue next fault i+1 of scanning.
S204: if system is unstable under fault i, then dichotomy is re-used to this fault and carry out stability limit search, obtain new section stability limit, the limit method of operation, control fault and become fault i, then to the new limit method of operation, scanning consequent malfunction i+1 is continued.
S205: when scanning last fault, just obtains the control fault resolving face stability limit and correspondence.
Wherein, the employing dichotomy mentioned in step S201 carries out complete stability limit search as shown in Figure 3, further comprising the steps:
S2011: under the basic method of operation, system is stable to setting fault, is designated as stationary mode S a, now section through-put power is P a; With larger step-size change power transmission district and powered district generator output, make system unstable under setting fault, be designated as unstable mode S b, now section through-put power is designated as P b.
S2012: if P b-P abe less than the threshold value preset, wherein threshold value generally can be taken as 10MW, then go to step S2015, otherwise change power transmission district and powered district generator output variable quantity are original half, regenerates the new method of operation and makes section through-put power for (P a+ P b)/2, which is designated as S c, calculate setting fault in this mode, carry out time-domain-simulation.
S2013: if system stability, then stationary mode S abecome mode S c, P abecome (P a+ P b)/2, unstable mode is constant, returns step S2012 and continues to calculate.
S2014: if system is unstable, then unstable mode S bbecome mode S c, P bbecome (P a+ P b)/2, stationary mode is constant, returns step S2012 and continues to calculate.
S2015: real stability limit is positioned at P aand P bbetween, note P afor final section stability limit.
S3: under the limit method of operation, for control fault, carries out the identification of a transient stability unstability group of planes and a stable group of planes, obtains a unstability group of planes and a stable group of planes.
Particularly, as shown in Figure 4, identify that the method for a unstability group of planes and a stable group of planes is as follows:
S301: under the limit method of operation that step S2 obtains, carries out time-domain-simulation to control fault, obtains neutrality track, then slightly increase the fault clearance time, generally increase 0.01s, makes system transient modelling unstable, obtains critical unstable track.
S302: estimate controlling unstable equilibrium point according to neutrality track and critical unstable track, concrete grammar is for from the fault clearance moment, search for first kinetic energy minimum point on neutrality track and critical unstable track respectively, then using the mean value of 2 as controlling unstable equilibrium point.
The computing formula of electric system kinetic energy is as follows:
E = Σ i = 1 n g 1 2 M i ω i 2 ;
Wherein, E represents electric system kinetic energy, n grepresent generator bus number, ω irepresent the rotating speed of generator.
S303: under calculating this fault, energy margin is to the sensitivity of each generator output, and computing formula is as follows:
S g , i = ( δ cl , i - δ uep , i ) - ω cl , i t cl + 1 2 M i ( P m , i - P ecl , i ) t cl 2
Wherein, S g, irepresent that energy margin is to the sensitivity of the mechanical output of generator i, δ cl, irepresent that generator i is at the merit angle in fault clearance moment under center of inertia coordinate, δ uep, irepresent that generator i is at the merit angle at controlling unstable equilibrium point place under center of inertia coordinate, ω cl, irepresent that generator i is at the rotating speed in fault clearance moment under center of inertia coordinate, t clrepresent to the time of fault clearance from fault, M irepresent the moment of inertia of generator i, P m, irepresent the mechanical output of generator i, P ecl, ithe electromagnetic power of expression generator i moment after fault clearance.
S304: according to the result of sensitivity, unit is hived off, the minimum a part of unit of sensitivity be a unstability group of planes, the maximum a part of unit of sensitivity for stablize a group of planes.In this step, minimum sensitivity is negative, maximum sensitivity is just.Be the alternative generator of generator as a unstability group of planes of negative value using all sensitivity.When specifically choosing, sensitivity corresponding for all alternative generators can be sorted, the generator that the numerical value of corresponding sensitivity is less may be more critical machine.Generally, for easy to operate consideration, do not need all alternative generators to be identified as a unstability group of planes, like this can according to the ranking results of sensitivity, the generator of predetermined quantity before less for the numerical value of corresponding sensitivity is identified as a unstability group of planes.Predetermined quantity is determined according to concrete electrical network, and those skilled in the art know how to determine its numerical value, do not repeat them here.The choosing method stablizing a group of planes is similar.
S4: generate new generator output and change scheme, the unit increase in a unstability group of planes is exerted oneself or started shooting successively by the unit of not starting shooting in its region, reduces the unit in a stable group of planes accordingly and exerts oneself until shut down.
S5: change scheme for new generator output, the transient stability limit of search section, obtains the transient stability limit of section, the limit method of operation and controls fault accordingly.
S6: if the control fault phase that obtains of the control fault that obtains of S5 and S2 with; calculate end, the transient stability limit obtained with S5 is the final section limit, otherwise turns S3, regenerates generator output and changes scheme and search for transient stability limit.
Transient stability limit of transmission section of power system computing method of the present invention, scheme is changed with the generator output that stabilizing machine all living creatures becomes stability limit less based on a unstability group of planes, and change scheme by iteration adjustment generator output, ensure the final limit method of operation corresponding be that the most dangerous generator output changes scheme, simultaneously stable to all faults and control fault critical is stablized, be one the most conservative while can ensure the stability limit of grid stability.The method increase the analysis of electrical network and run control technology level, operation of power networks personnel accurately can grasp stability limit and the stability margin of each transmission cross-section, can ensure the transient stability of system when carrying out method of operation adjustment.
Describe and can be understood in process flow diagram or in this any process otherwise described or method, represent and comprise one or more for realizing the module of the code of the executable instruction of the step of specific logical function or process, fragment or part, and the scope of the preferred embodiment of the present invention comprises other realization, wherein can not according to order that is shown or that discuss, comprise according to involved function by the mode while of basic or by contrary order, carry out n-back test, this should understand by embodiments of the invention person of ordinary skill in the field.
In flow charts represent or in this logic otherwise described and/or step, such as, the sequencing list of the executable instruction for realizing logic function can be considered to, may be embodied in any computer-readable medium, for instruction execution system, device or equipment (as computer based system, comprise the system of processor or other can from instruction execution system, device or equipment instruction fetch and perform the system of instruction) use, or to use in conjunction with these instruction execution systems, device or equipment.With regard to this instructions, " computer-readable medium " can be anyly can to comprise, store, communicate, propagate or transmission procedure for instruction execution system, device or equipment or the device that uses in conjunction with these instruction execution systems, device or equipment.The example more specifically (non-exhaustive list) of computer-readable medium comprises following: the electrical connection section (electronic installation) with one or more wiring, portable computer diskette box (magnetic device), random-access memory (ram), ROM (read-only memory) (ROM), erasablely edit ROM (read-only memory) (EPROM or flash memory), fiber device, and portable optic disk ROM (read-only memory) (CDROM).In addition, computer-readable medium can be even paper or other suitable media that can print described program thereon, because can such as by carrying out optical scanning to paper or other media, then carry out editing, decipher or carry out process with other suitable methods if desired and electronically obtain described program, be then stored in computer memory.
Should be appreciated that each several part of the present invention can realize with hardware, software, firmware or their combination.In the above-described embodiment, multiple step or method can with to store in memory and the software performed by suitable instruction execution system or firmware realize.Such as, if realized with hardware, the same in another embodiment, can realize by any one in following technology well known in the art or their combination: the discrete logic with the logic gates for realizing logic function to data-signal, there is the special IC of suitable combinational logic gate circuit, programmable gate array (PGA), field programmable gate array (FPGA) etc.
Those skilled in the art are appreciated that realizing all or part of step that above-described embodiment method carries is that the hardware that can carry out instruction relevant by program completes, described program can be stored in a kind of computer-readable recording medium, this program perform time, step comprising embodiment of the method one or a combination set of.
In addition, each functional unit in each embodiment of the present invention can be integrated in a processing module, also can be that the independent physics of unit exists, also can be integrated in a module by two or more unit.Above-mentioned integrated module both can adopt the form of hardware to realize, and the form of software function module also can be adopted to realize.If described integrated module using the form of software function module realize and as independently production marketing or use time, also can be stored in a computer read/write memory medium.
In the description of this instructions, specific features, structure, material or feature that the description of reference term " embodiment ", " some embodiments ", " example ", " concrete example " or " some examples " etc. means to describe in conjunction with this embodiment or example are contained at least one embodiment of the present invention or example.In this manual, identical embodiment or example are not necessarily referred to the schematic representation of above-mentioned term.And the specific features of description, structure, material or feature can combine in an appropriate manner in any one or more embodiment or example.
Although illustrate and describe embodiments of the invention, those having ordinary skill in the art will appreciate that: can carry out multiple change, amendment, replacement and modification to these embodiments when not departing from principle of the present invention and aim, scope of the present invention is by claim and equivalents thereof.

Claims (3)

1. computing method for transient stability limit of transmission section of power system, is characterized in that, comprise the following steps:
S1: electrical network is divided into power transmission district and powered district at analyzed transmission cross-section place, set a kind of initial generator output and change scheme, increase described power transmission district generated output, equivalent reduces described powered district generated output, to increase section through-put power;
S2: the generator output for setting changes scheme, the transient stability limit of search section, obtains the transient stability limit of section, the limit method of operation and controls fault accordingly;
S3: under the described limit method of operation, for control fault, carries out the identification of a transient stability unstability group of planes and a stable group of planes, obtains a unstability group of planes and a stable group of planes;
S4: generate new generator output and change scheme, the unit increase in a described unstability group of planes is exerted oneself or started shooting successively by the unit of not starting shooting in its region, reduces the unit in a stable group of planes accordingly and exerts oneself until shut down;
S5: scheme is changed for new generator output, the described transient stability limit of search section, obtains the transient stability limit of section, the limit method of operation and controls fault accordingly; And
S6: if the control fault phase that obtains of the control fault that obtains of step S5 and step S2 with, then calculate end, the transient stability limit obtained with step S5 is the final section limit, otherwise goes to step S3, regenerate generator output change scheme and search for transient stability limit
Wherein, search for the method for the transient stability limit of section in described step S2, comprise the following steps:
S201: the 1st fault that forecast failure is concentrated, adopt dichotomy to carry out complete stability limit search, obtain the section stability limit that this fault is corresponding, controlling fault is this fault;
S202: to the limit method of operation corresponding to stability limit described in step S201, successively scan consequent malfunction i, wherein i be fault sequence number and i be greater than 1 positive integer;
S203: if system stability under fault i, then section stability limit, the limit method of operation and control fault are all constant, continue next fault i+1 of scanning;
S204: if system is unstable under fault i, then dichotomy is re-used to this fault and carry out stability limit search, obtain new section stability limit, the limit method of operation, control fault and become fault i, then to the new limit method of operation, scanning consequent malfunction i+1 is continued; And
S205: when scanning last fault, obtains the control fault of section stability limit and correspondence.
2. the method for claim 1, is characterized in that, adopts dichotomy to carry out complete stability limit search, comprise the following steps described in described step S201:
S2011: under the basic method of operation, system is stable to setting fault, is designated as stationary mode S a, now section through-put power is P a; Change power transmission district and powered district generator output, make system unstable under setting fault, be designated as unstable mode S b, now section through-put power is designated as P b;
S2012: if P b-P abe less than the threshold value preset, then go to step S2015, otherwise change power transmission district and powered district generator output variable quantity are original half, regenerate the new method of operation and make section through-put power for (P a+ P b)/2, which is designated as S c, in this mode time-domain-simulation is carried out to setting fault;
S2013: if system stability, then stationary mode S abecome mode S c, P abecome (P a+ P b)/2, unstable mode is constant, returns step S2012 and continues to calculate;
S2014: if system is unstable, then unstable mode S bbecome mode S c, P bbecome (P a+ P b)/2, stationary mode is constant, returns step S2012 and continues to calculate;
S2015: real stability limit is positioned at P aand P bbetween, note P afor final section stability limit.
3. the method for claim 1, is characterized in that, identifies the method for a unstability group of planes and a stable group of planes, comprise the following steps in described step S3:
S301: under the limit method of operation that step S2 obtains, carries out time-domain-simulation to control fault, obtains neutrality track, then increase the fault clearance time, and make system transient modelling unstable, emulation obtains critical unstable track;
S302: estimate controlling unstable equilibrium point according to described neutrality track and described critical unstable track, concrete grammar is for from the fault clearance moment, search for first kinetic energy minimum point on described neutrality track and described critical unstable track respectively, then using the mean value of 2 as controlling unstable equilibrium point;
S303: under calculating this fault, energy margin is to the sensitivity of each generator output, and computing formula is as follows:
S g , i = ( δ cl , i - δ uep , i ) - ω cl , i t cl + 1 2 M i ( P m , i - P ecl , i ) t cl 2
Wherein, S g, irepresent that energy margin is to the sensitivity of the mechanical output of generator i, δ cl, irepresent that generator i is at the merit angle in fault clearance moment under center of inertia coordinate, δ uep, irepresent that generator i is at the merit angle at described controlling unstable equilibrium point place under center of inertia coordinate, ω cl, irepresent that generator i is at the rotating speed in fault clearance moment under center of inertia coordinate, t clrepresent to the time of fault clearance from fault, M irepresent the moment of inertia of generator i, p m, irepresent the mechanical output of generator i, P ecl, ithe electromagnetic power of expression generator i moment after fault clearance; And
S304: according to the result of sensitivity, hive off to unit, the minimum a part of unit of sensitivity is a unstability group of planes, and the maximum a part of unit of sensitivity is for stablize a group of planes, and in a group of planes, the quantity of unit is determined according to concrete electrical network.
CN201310034636.3A 2013-01-29 2013-01-29 A kind of computing method of transient stability limit of transmission section of power system Active CN103065060B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201310034636.3A CN103065060B (en) 2013-01-29 2013-01-29 A kind of computing method of transient stability limit of transmission section of power system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201310034636.3A CN103065060B (en) 2013-01-29 2013-01-29 A kind of computing method of transient stability limit of transmission section of power system

Publications (2)

Publication Number Publication Date
CN103065060A CN103065060A (en) 2013-04-24
CN103065060B true CN103065060B (en) 2015-10-28

Family

ID=48107688

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201310034636.3A Active CN103065060B (en) 2013-01-29 2013-01-29 A kind of computing method of transient stability limit of transmission section of power system

Country Status (1)

Country Link
CN (1) CN103065060B (en)

Families Citing this family (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413043B (en) * 2013-08-09 2016-04-20 中国南方电网有限责任公司 A kind of method for solving of electric system hyperspace section thermal stability limit boundary
CN103823998B (en) * 2014-03-19 2017-05-10 国家电网公司 Weak cross section determination method taking influence of network topology changes on power transmission capacity into consideration
CN107767028B (en) * 2017-09-26 2021-09-07 广东电网有限责任公司电力调度控制中心 Transmission section limit automatic identification and out-of-limit calculation method based on xml
CN110943453B (en) * 2019-12-23 2022-12-02 北京交通大学 Power system fault sample generation and model construction method facing transfer learning

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN101969202A (en) * 2010-11-23 2011-02-09 国网电力科学研究院 Method for recognizing transmission section power limit interval of power system transient stability
JP2011050165A (en) * 2009-08-26 2011-03-10 Chugoku Electric Power Co Inc:The Method and device for calculation of transient stability limit value, and program
CN102035203A (en) * 2010-10-20 2011-04-27 中国电力科学研究院 Fast calculation method of tie-line transient state stable limiting transmission power

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011050165A (en) * 2009-08-26 2011-03-10 Chugoku Electric Power Co Inc:The Method and device for calculation of transient stability limit value, and program
CN102035203A (en) * 2010-10-20 2011-04-27 中国电力科学研究院 Fast calculation method of tie-line transient state stable limiting transmission power
CN101969202A (en) * 2010-11-23 2011-02-09 国网电力科学研究院 Method for recognizing transmission section power limit interval of power system transient stability

Non-Patent Citations (5)

* Cited by examiner, † Cited by third party
Title
A new method for obtaining critical clearing time for transient stability;Yorino N, et al.;《Power Systems, IEEE transactions on》;20101231;第25卷(第3期);1620-1626 *
Direct nonlinear primal-dual interior-point method for transient stability constrained optimal power flow;Xia Y,et al.;《IEEE Proceedings: Generation transmission and distribution》;20051231;第152卷(第1期);11-16 *
The Online Application and Fast Solving Method for Critical Clearing Time of Three-Phase Short Circuit in Power System;Lingyu Ren,et al.;《2012 (2nd) International Conference on Smart Grid and Clean Energy Technologies(ICSGCE 2012)》;20121231;94-99 *
基于三道防线的故障集自动批量生成方法;张晓东等;《电力系统及其自动化》;20121231;第34卷(第5期);53-55 *
电力系统静态电压稳定极限及裕度计算综述;吴政球等;《电力系统及其自动化学报》;20100228;第22卷(第1期);126-131 *

Also Published As

Publication number Publication date
CN103065060A (en) 2013-04-24

Similar Documents

Publication Publication Date Title
CN105425782B (en) Real-time fault simulation platform of high-speed train traction transmission control system
CN106650072B (en) Semi-physical simulation-based virtual grid-connected detection system and method for same generator
CN103065060B (en) A kind of computing method of transient stability limit of transmission section of power system
CN108736465B (en) Future-state power grid model construction method and device
Baldwin et al. Distribution load flow methods for shipboard power systems
CN107546769B (en) Method for obtaining transient stability of grid-connected inverter type distributed power supply
Feng et al. A two-level forced oscillations source location method based on phasor and energy analysis
CN103197556A (en) Half period repetitive control method based on attractive rule
CN102346442A (en) Method for carrying out parallel transient simulation computing on microgrid suitable for nodal analysis frameworks
CN103700036A (en) Transient stability projection integral method suitable for multi-time scale of electrical power system
Dalla Vedova et al. Evaluation of the correlation coefficient as a prognostic indicator for electromechanical servomechanism failures
KR101980685B1 (en) Appratus and method of analyzing power system network
Oerter et al. Experience with first smart, autonomous LV-grids in Germany
JP2015532578A (en) Multiple input single output combined power generation system
Bélanger et al. A modern and open real-time digital simulator of contemporary power systems
Mallem et al. Dynamic equivalent of neighbor power system for day-ahead stability studies
Gerin-Lajoie et al. Hydro-Quebec experience with PSS tuning
Alashter et al. Fault Diagnosis of Broken Rotor Bar in AC Induction Motor based on A Qualitative Simulation Approach
Li et al. Automatic test generation of large boolean expressions in computer based interlocking system
Mannah et al. Hardware in the loop simulation for optimal management of electrical power converters
Chen et al. A Convex–Concave Procedure-Based Method for Optimal Power Flow of Offshore Wind Farms
Nikiforov The through simulation of devices on the basis of the structural linguistic method
Taleires Filho et al. Motor test bench control for load profile emulation
Senyuk et al. A PMU-based algorithm of synchronous generator stability prediction during a disturbance
Yang et al. A unified scheme for power system transient stability simulation and transient stability assessment

Legal Events

Date Code Title Description
C06 Publication
PB01 Publication
C10 Entry into substantive examination
SE01 Entry into force of request for substantive examination
C14 Grant of patent or utility model
GR01 Patent grant