CN103065060A - Computing method for transient stability limit of transmission section of power system - Google Patents

Computing method for transient stability limit of transmission section of power system Download PDF

Info

Publication number
CN103065060A
CN103065060A CN2013100346363A CN201310034636A CN103065060A CN 103065060 A CN103065060 A CN 103065060A CN 2013100346363 A CN2013100346363 A CN 2013100346363A CN 201310034636 A CN201310034636 A CN 201310034636A CN 103065060 A CN103065060 A CN 103065060A
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.)
Granted
Application number
CN2013100346363A
Other languages
Chinese (zh)
Other versions
CN103065060B (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

Images

Landscapes

  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention provides a computing method for a transient stability limit of a transmission section of a power system. The method comprises the following steps of: S1, dividing a power grid into a power transmitting region and a power receiving region, increasing the output of the power transmitting region, reducing the output of the power receiving region, and increasing the transmission power of the section; S2, obtaining a transient stability limit, a limit operation mode and a control fault according to a set output changing scheme; S3, recognizing an instable cluster and a stable cluster of the control fault in the limit mode; S4, generating a new output changing scheme, increasing the output of a machine set in the instable cluster or starting a machine set which is not started, and reducing the output of the machine set in the stable cluster until the machine set stops; S5, obtaining a transient stability limit, a limit operation mode and a control fault according to the new output changing scheme; and S6, ending if the control fault in the step S5 is same as that in the step S2 and taking the transient stability limit obtained in the step S5 as a result, or otherwise, executing the step S3 to regenerate an output changing scheme and search a transient stability limit. The computing method has the advantages of high accuracy and good stability.

Description

A kind of computing method of the electric system transmission cross-section transient stability limit
Technical field
The present invention relates to the Power System Analysis technical field, particularly a kind of computing method of the electric system transmission cross-section transient stability limit.
Background technology
Transient stability is the important content of power system safety and stability.Operation of power networks department need to guarantee that electrical network can keep transient stability under forecast failure.In operation of power networks, generally guarantee the transient stability of system by control transmission cross-section power.Transmission cross-section is the set of some branch road in the electrical network, and these branch roads are divided into two disconnected zones with electrical network after disconnecting, and the power on the transmission cross-section represents the power of two interregional transmission.When the power that transmits on the transmission cross-section is excessive, system may be under certain fault transient state unstable, therefore, there is a transient stability limit in transmission cross-section, when section power during less than the limit, for all concentrated faults of forecast failure, system is transient stability, when the section increased power to greater than the limit time, some faults of concentrating for forecast failure, system at first transient state is unstable, and this fault is called the control fault.In operation of power networks, control section power is no more than the transient stability limit, can guarantee the transient stability of system.
The Conventional Calculation Method of the section transient stability limit is: take a certain method of operation as the basis, keep load and other regional generator outputs constant, this section of corresponding increase and decrease send the electricity district and is subjected to electricity district unit, and the fault that forecast failure is concentrated carried out stability Calculation, be met the section limit of power system transient stability.How increasing and decreasing unit (be generator output change scheme) has a significant impact the result of calculation tool of the final transient stability limit.When generator output changes scheme not simultaneously, the stability limit that obtains is likely different, this is so that the result of calculation of stability limit depends on artificial setting in advance, can not guarantee to obtain the stability limit in the most serious situation (generator output of stability limit minimum shifts gears), i.e. the most conservative stability limit.This is disadvantageous for the security that ensures electrical network.Therefore wish to find the most dangerous generator output to change scheme, the transient stability limit that obtains under this scheme is minimum, and is 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 provides at least a kind of useful commerce to select.For this reason, the object of the invention is to propose a kind of computing method of the electric system transmission cross-section transient stability limit, the method can guarantee the conservative property of transient stability limit result of calculation, ensures power grid security.
To achieve these goals, computing method according to the electric system transmission cross-section transient stability limit of the embodiment of the invention comprise: S1: at the transmission cross-section place that analyzes electrical network is divided into and send the electricity district and be subjected to the electricity district, set a kind of initial generator output and change scheme, increase the described electricity district generated output that send, equivalent reduces the described electricity district generated output that is subjected to, to increase the section through-put power; S2: change scheme for the generator output of setting, the transient stability limit of search section obtains the transient stability limit, the limit method of operation of section and controls accordingly fault; S3: under the described limit method of operation, for the control fault, carry out the identification of a transient stability unstability group of planes and a stable group of planes, obtain a unstability group of planes and a stable group of planes; S4: generate new generator output and change scheme, the unit increase in the described unstability group of planes is exerted oneself or the not start unit in its region is started shooting successively, and the unit minimizing that will stablize accordingly in the group of planes is exerted oneself until shut down; S5: change scheme for new generator output, the described transient stability limit of search section obtains the transient stability limit, the limit method of operation of section and controls accordingly fault; S6: if the control fault phase that the control fault that step S5 obtains and step S2 obtain together, then calculate and finish, the transient stability limit that obtains take step S5 is as the final section limit, otherwise turns step S3, regenerates generator output and changes scheme and search for the transient stability limit.
In an embodiment of the present invention, the method of search section stability limit comprises among the described step S2: S201: to the 1st concentrated fault of forecast failure, adopt dichotomy to carry out complete stability limit search, obtain section stability limit corresponding to this fault, the control fault is this fault; S202: to the corresponding limit method of operation of the described stability limit of step S201, scan successively consequent malfunction i, wherein i is that fault sequence number and i are the positive integer greater than 1; S203: if system stability under the 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 the fault i, then this fault is re-used dichotomy and carry out the stability limit search, obtain new section stability limit, the limit method of operation, the control fault becomes fault i, then to the new limit method of operation, continue scanning consequent malfunction i+1; S205: when scanning last fault, just obtain resolving face stability limit and corresponding control fault.
In an embodiment of the present invention, described employing dichotomy carries out the search of complete stability limit and comprises: S2011: system is stable to the setting fault under the basic method of operation, is designated as stationary mode S a, this moment, the section through-put power was P aSend the electricity district and be subjected to electricity district generator output with larger step-size change, so that system is unstable under the setting fault, be designated as unstable mode S b, this moment, the section through-put power was designated as P bS2012: if P b-P aLess than predefined threshold value, then turn step S2014, to send the electricity district and to be subjected to electricity district generator output variable quantity be original half otherwise change, and regenerates the new method of operation so that the section through-put power is (P a+ P b)/2, this mode is designated as S c, under this mode, carry out time-domain-simulation 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, unstable mode S then 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 aBe final section stability limit.
In an embodiment of the present invention, the method of an identification unstability group of planes and a stable group of planes comprises among the described step S3: S301: under the limit method of operation that step S2 obtains, the control fault is carried out time-domain-simulation, obtain the neutrality track, then slightly increase the fault clearance time, so that system transient modelling is unstable, emulation obtains critical unstable track; S302: estimate leading unstable equilibrium point according to described neutrality track and described critical unstable track, concrete grammar is for from fault clearance constantly, search for respectively first kinetic energy minimum point on described neutrality track and the described critical unstable track, then with 2 mean value as leading unstable equilibrium point; S303: calculate that energy margin is to the sensitivity of each generator output under this fault, 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, iThe expression energy margin is to the sensitivity of the mechanical output of generator i, δ Cl, iBe illustrated in generator i under the center of inertia coordinate at fault clearance merit angle constantly, δ Uep, iBe illustrated under the center of inertia coordinate generator i at the merit angle at described leading unstable equilibrium point place, ω Cl, iBe illustrated in generator i under the center of inertia coordinate at fault clearance rotating speed constantly, t ClExpression begins time to fault clearance, M from fault iThe moment of inertia of expression generator i, P M, iThe mechanical output of expression generator i, P Ecl, iRepresent the generator i electromagnetic power of moment after fault clearance.S304: the result according to sensitivity, unit is hived off, a part of unit of sensitivity minimum is a unstability group of planes, and a part of unit of sensitivity maximum is for to stablize a group of planes, and the quantity of unit is determined according to concrete electrical network in the group of planes.
According to electric system transmission cross-section transient stability Limit Calculation Method of the present invention, become the most dangerous generator output to change scheme with the stabilizing machine all living creatures based on a unstability group of planes, guarantee the conservative property of transient stability limit result of calculation, ensure power grid security.Described method has improved analysis and the operation control technology level of electrical network, and the operation of power networks personnel can accurately grasp stability limit and the stability margin of each transmission cross-section, carry out can guaranteeing when the method for operation is adjusted the transient stability of system.
Additional aspect of the present invention and advantage in the following description part provide, and part will become obviously from the following description, or recognize by practice of the present invention.
Description of drawings
Of the present invention and/or additional aspect and advantage are from obviously and easily understanding becoming the description of embodiment below in conjunction with accompanying drawing, wherein:
Fig. 1 is the process flow diagram according to the computing method of the electric system transmission cross-section transient stability limit of the embodiment of the invention;
Fig. 2 is the process flow diagram according to the method for the search section stability limit of the computing method of the electric system transmission cross-section transient stability limit of the embodiment of the invention;
Fig. 3 is the process flow diagram according to the dichotomy of the search stability limit of the computing method of the electric system transmission cross-section transient stability limit of the embodiment of the invention; And
Fig. 4 is the process flow diagram according to the method for the identification unstability group of planes of the computing method of the electric system transmission cross-section transient stability limit of the embodiment of the invention and a stable group of planes.
Embodiment
The below describes embodiments of the invention in detail, and the example of described embodiment is shown in the drawings, and wherein identical or similar label represents identical or similar element or the element with identical or similar functions from start to finish.Be exemplary below by the embodiment that is described with reference to the drawings, only be used for explaining the present invention, and can not be interpreted as limitation of the present invention.On the contrary, embodiments of the invention comprise spirit and interior all changes, modification and the equivalent of intension scope that falls into additional claims.
In description of the invention, it will be appreciated that term " first ", " second " etc. only are used for describing purpose, and can not be interpreted as indication or hint relative importance.In description of the invention, need to prove that unless clear and definite regulation and restriction are arranged in addition, term " links to each other ", " connection " should do broad understanding, for example, can be to be fixedly connected with, and also can be to removably connect, or connect integratedly; Can be mechanical connection, also can be to be electrically connected; Can be directly to link to each other, also can indirectly link to each other by intermediary.For the ordinary skill in the art, can concrete condition understand above-mentioned term concrete meaning in the present invention.In addition, in description of the invention, except as otherwise noted, the implication of " a plurality of " is two or more.
According to the described computing method for the electric system transmission cross-section transient stability limit of the embodiment of the invention, at 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 and change scheme, an increase unstability group of planes is exerted oneself and is reduced to stablize a group of planes and exerts oneself, because increasing a unstability group of planes exerts oneself and reduces to stablize the exert oneself transient stability of reduction system that can be the fastest of a group of planes, therefore compare with initial generator output change scheme, new scheme is a kind of more dangerous scheme, the stability limit of section is less under this scheme, by generator output being changed the iteration adjustment of scheme, guarantee the conservative property of section stability limit.
As shown in Figure 1, the embodiment of the invention may further comprise the steps:
S1: at the transmission cross-section place that analyzes electrical network is divided into two parts, is respectively and send the electricity district and be subjected to the electricity district, set a kind of initial generator output and change scheme, increase and send electricity district generated output, reduce being subjected to electricity district generated output, to increase the section through-put power.
S2: change scheme for the generator output of setting, the transient stability limit of search section obtains the transient stability limit, the limit method of operation of section and controls accordingly fault.
Particularly, as shown in Figure 2, the method for search section stability limit further may further comprise the steps:
S201: to the 1st fault that forecast failure is concentrated, adopt dichotomy to carry out complete stability limit search, obtain section stability limit corresponding to this fault, the control fault is this fault.
S202: to the corresponding limit method of operation of step S201 stability limit, scan successively consequent malfunction i, wherein i is that fault sequence number and i are the positive integer greater than 1.
S203: if system stability under the 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 the fault i, then this fault is re-used dichotomy and carry out the stability limit search, obtain new section stability limit, the limit method of operation, the control fault becomes fault i, then to the new limit method of operation, continue scanning consequent malfunction i+1.
S205: when scanning last fault, just obtain resolving face stability limit and corresponding control fault.
Wherein, the employing dichotomy of mentioning among the step S201 carries out complete stability limit search as shown in Figure 3, further may further comprise the steps:
S2011: system is stable to setting fault under the basic method of operation, is designated as stationary mode S a, this moment, the section through-put power was P aSend the electricity district and be subjected to electricity district generator output with larger step-size change, so that system is unstable under the setting fault, be designated as unstable mode S b, this moment, the section through-put power was designated as P b
S2012: if P b-P aLess than predefined threshold value, wherein threshold value generally can be taken as 10MW, then turns step S2015, and to send the electricity district and to be subjected to electricity district generator output variable quantity be original half otherwise change, and regenerates the new method of operation so that the section through-put power is (P a+ P b)/2, this mode is designated as S c, under this mode, calculate and set fault, 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, unstable mode S then 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 aBe final section stability limit.
S3: under the limit method of operation, for the control fault, carry out the identification of a transient stability unstability group of planes and a stable group of planes, obtain a unstability group of planes and a stable group of planes.
Particularly, as shown in Figure 4, the method for an identification unstability group of planes and a stable group of planes is as follows:
S301: under the limit method of operation that step S2 obtains, the control fault is carried out time-domain-simulation, obtain the neutrality track, then slightly increase the fault clearance time, generally increase 0.01s, so that system transient modelling is unstable, obtain critical unstable track.
S302: estimate leading unstable equilibrium point according to neutrality track and critical unstable track, concrete grammar is for from fault clearance constantly, search for respectively first kinetic energy minimum point on neutrality track and the critical unstable track, then with 2 mean value as leading 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 gExpression generator bus number, ω iThe rotating speed of expression generator.
S303: calculate that energy margin is to the sensitivity of each generator output under this fault, 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, iThe expression energy margin is to the sensitivity of the mechanical output of generator i, δ Cl, iBe illustrated in generator i under the center of inertia coordinate at fault clearance merit angle constantly, δ Uep, iBe illustrated in generator i under the center of inertia coordinate at the merit angle at leading unstable equilibrium point place, ω Cl, iBe illustrated in generator i under the center of inertia coordinate at fault clearance rotating speed constantly, t ClExpression begins time to fault clearance, M from fault iThe moment of inertia of expression generator i, P M, iThe mechanical output of expression generator i, P Ecl, iRepresent the generator i electromagnetic power of moment after fault clearance.
S304: the result according to 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 be to stablize a group of planes.In this step, minimum sensitivity is for negative, and maximum sensitivity is for just.Be that the generator of negative value is as the alternative generator of a unstability group of planes with all sensitivity.When specifically choosing, can the sensitivity that all alternative generators are corresponding sort, the generator that the numerical value of corresponding sensitivity is less more may be critical machine.Generally speaking, for easy to operate consideration, do not need all alternative generators are identified as a unstability group of planes, like this can be according to the ranking results of sensitivity, the generator of the front predetermined quantity that the numerical value of corresponding sensitivity is less is identified as a unstability group of planes.Predetermined quantity determines that according to concrete electrical network those skilled in the art know its numerical value of how to confirm, do not repeat them here.The choosing method of stablizing a group of planes is similar.
S4: generate new generator output and change scheme, the unit increase in the unstability group of planes is exerted oneself or the not start unit in its region is started shooting successively, and the unit minimizing that will stablize accordingly in the group of planes is exerted oneself until shut down.
S5: change scheme for new generator output, the transient stability limit of search section obtains the transient stability limit, the limit method of operation of section and controls accordingly fault.
S6: if the control fault phase that the control fault that S5 obtains and S2 obtain with calculate to finish, the transient stability limit that obtains take S5 is as the final section limit, otherwise turns S3, regenerates generator output and changes scheme and search for the transient stability limit.
Electric system transmission cross-section transient stability Limit Calculation Method of the present invention, become the less generator output of stability limit to change scheme with the stabilizing machine all living creatures based on a unstability group of planes, and by iteration adjustment generator output change scheme, guarantee the final limit method of operation corresponding be that the most dangerous generator output changes scheme, stable to all faults and stable to the control fault critical simultaneously, be the stability limit that a while the most conservative can guarantee grid stability.The method has improved analysis and the operation control technology level of electrical network, and the operation of power networks personnel can accurately grasp stability limit and the stability margin of each transmission cross-section, carry out can guaranteeing when the method for operation is adjusted the transient stability of system.
Describe and to be understood in the process flow diagram or in this any process of otherwise describing or method, expression comprises module, fragment or the part of code of the executable instruction of the step that one or more is used to realize specific logical function or process, and the scope of preferred implementation of the present invention comprises other realization, wherein can be not according to order shown or that discuss, comprise according to related function by the mode of basic while or by opposite order, carry out function, this should be understood by the embodiments of the invention person of ordinary skill in the field.
In process flow diagram the expression or in this logic of otherwise describing and/or step, for example, can be considered to the sequencing tabulation for the executable instruction that realizes logic function, may be embodied in any computer-readable medium, use for instruction execution system, device or equipment (such as the computer based system, comprise that the system of processor or other can and carry out the system of instruction from instruction execution system, device or equipment instruction fetch), or use in conjunction with these instruction execution systems, device or equipment.With regard to this instructions, " computer-readable medium " can be anyly can comprise, storage, communication, propagation or transmission procedure be for instruction execution system, device or equipment or the device that uses in conjunction with these instruction execution systems, device or equipment.The more specifically example of computer-readable medium (non-exhaustive list) comprises following: the electrical connection section (electronic installation) with one or more wirings, portable computer diskette box (magnetic device), random-access memory (ram), ROM (read-only memory) (ROM), the erasable ROM (read-only memory) (EPROM or flash memory) of editing, fiber device, and portable optic disk ROM (read-only memory) (CDROM).In addition, computer-readable medium even can be paper or other the suitable media that to print described program thereon, because can be for example by paper or other media be carried out optical scanning, then edit, decipher or process to obtain described program in the electronics mode with other suitable methods in case of necessity, then it is stored in the 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, a plurality of steps or method can realize with being stored in the storer and by software or firmware that suitable instruction execution system is carried out.For example, if realize with hardware, the same in another embodiment, can realize with the combination of each or they in the following technology well known in the art: have for the discrete logic of data-signal being realized the logic gates of logic function, special IC with suitable combinational logic gate circuit, programmable gate array (PGA), field programmable gate array (FPGA) etc.
Those skilled in the art are appreciated that and realize that all or part of step that above-described embodiment method is carried is to come the relevant hardware of instruction to finish by program, described program can be stored in a kind of computer-readable recording medium, this program comprises step of embodiment of the method one or a combination set of when carrying out.
In addition, each functional unit in each embodiment of the present invention can be integrated in the processing module, also can be that the independent physics of unit exists, and also can be integrated in the module two or more unit.Above-mentioned integrated module both can adopt the form of hardware to realize, also can adopt the form of software function module to realize.If described integrated module realizes with the form of software function module and during as independently production marketing or use, also can be stored in the computer read/write memory medium.
In the description of this instructions, the description of reference term " embodiment ", " some embodiment ", " example ", " concrete example " or " some examples " etc. means to be contained at least one embodiment of the present invention or the example in conjunction with specific features, structure, material or the characteristics of this embodiment or example description.In this manual, the schematic statement of above-mentioned term not necessarily referred to identical embodiment or example.And the specific features of description, structure, material or characteristics can be with suitable mode combinations in any one or more embodiment or example.
Although illustrated and described embodiments of the invention, those having ordinary skill in the art will appreciate that: can carry out multiple variation, modification, replacement and modification to these embodiment in the situation that does not break away from principle of the present invention and aim, scope of the present invention is limited by claim and equivalent thereof.

Claims (4)

1. the computing method of an electric system transmission cross-section transient stability limit is characterized in that, may further comprise the steps:
S1: at the transmission cross-section place that analyzes electrical network is divided into and send the electricity district and to be subjected to the electricity district, set a kind of initial generator output and change scheme, increase the described electricity district generated output that send, equivalent reduces the described electricity district generated output that is subjected to, with increase section through-put power;
S2: change scheme for the generator output of setting, the transient stability limit of search section obtains the transient stability limit, the limit method of operation of section and controls accordingly fault;
S3: under the described limit method of operation, for the control fault, carry out the identification of a transient stability unstability group of planes and a stable group of planes, obtain a unstability group of planes and a stable group of planes;
S4: generate new generator output and change scheme, the unit increase in the described unstability group of planes is exerted oneself or the not start unit in its region is started shooting successively, and the unit minimizing that will stablize accordingly in the group of planes is exerted oneself until shut down;
S5: change scheme for new generator output, the described transient stability limit of search section obtains the transient stability limit, the limit method of operation of section and controls accordingly fault; And
S6: if the control fault phase that the control fault that step S5 obtains and step S2 obtain together, then calculate and finish, the transient stability limit that obtains take step S5 is as the final section limit, otherwise turns step S3, regenerates generator output and changes scheme and search for the transient stability limit.
2. the method for claim 1 is characterized in that, the method for search section stability limit among the described step S2 may further comprise the steps:
S201: to the 1st fault that forecast failure is concentrated, adopt dichotomy to carry out complete stability limit search, obtain section stability limit corresponding to this fault, the control fault is this fault;
S202: to the corresponding limit method of operation of the described stability limit of step S201, scan successively consequent malfunction i, wherein i is that fault sequence number and i are the positive integer greater than 1;
S203: if system stability under the 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 the fault i, then this fault is re-used dichotomy and carry out the stability limit search, obtain new section stability limit, the limit method of operation, the control fault becomes fault i, then to the new limit method of operation, continue scanning consequent malfunction i+1; And
S205: when scanning last fault, obtained section stability limit and corresponding control fault.
3. method as claimed in claim 2 is characterized in that, adopts dichotomy to carry out complete stability limit search described in the described step S201, may further comprise the steps:
S2011: system is stable to setting fault under the basic method of operation, is designated as stationary mode S a, this moment, the section through-put power was P aSend the electricity district and be subjected to electricity district generator output with larger step-size change, so that system is unstable under the setting fault, be designated as unstable mode S b, this moment, the section through-put power was designated as P b
S2012: if P b-P aLess than predefined threshold value, then turn step S2015, to send the electricity district and to be subjected to electricity district generator output variable quantity be original half otherwise change, and regenerates the new method of operation so that the section through-put power is (P a+ P b)/2, this mode is designated as S c, under this mode, carry out time-domain-simulation 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, unstable mode S then 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 aBe final section stability limit.
4. the method for claim 1 is characterized in that, the method for an identification unstability group of planes and a stable group of planes among the described step S3 may further comprise the steps:
S301: under the limit method of operation that step S2 obtains, the control fault is carried out time-domain-simulation, obtain the neutrality track, then slightly increase the fault clearance time, so that system transient modelling is unstable, emulation obtains critical unstable track;
S302: estimate leading unstable equilibrium point according to described neutrality track and described critical unstable track, concrete grammar is for from fault clearance constantly, search for respectively first kinetic energy minimum point on described neutrality track and the described critical unstable track, then with 2 mean value as leading unstable equilibrium point;
S303: calculate that energy margin is to the sensitivity of each generator output under this fault, 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, iThe expression energy margin is to the sensitivity of the mechanical output of generator i, δ Cl, iBe illustrated in generator i under the center of inertia coordinate at fault clearance merit angle constantly, δ Uep, iBe illustrated under the center of inertia coordinate generator i at the merit angle at described leading unstable equilibrium point place, ω Cl, iBe illustrated in generator i under the center of inertia coordinate at fault clearance rotating speed constantly, t ClExpression begins time to fault clearance, M from fault iThe moment of inertia of expression generator i, P M, iThe mechanical output of expression generator i, P Ecl, iRepresent the generator i electromagnetic power of moment after fault clearance; And
S304: the result according to sensitivity, unit is hived off, a part of unit of sensitivity minimum is a unstability group of planes, and a part of unit of sensitivity maximum is for to stablize a group of planes, and the quantity of unit is determined according to concrete electrical network in the group of planes.
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 true CN103065060A (en) 2013-04-24
CN103065060B 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)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413043A (en) * 2013-08-09 2013-11-27 中国南方电网有限责任公司 Solving method for power system multi-dimensional space cross section thermal stability limit boundary
CN103823998A (en) * 2014-03-19 2014-05-28 国家电网公司 Weak cross section determination method taking influence of network topology changes on power transmission capacity into consideration
CN107767028A (en) * 2017-09-26 2018-03-06 广东电网有限责任公司电力调度控制中心 A kind of transmission cross-section limit automatic identification based on xml and out-of-limit computational methods
CN110943453A (en) * 2019-12-23 2020-03-31 北京交通大学 Power system fault sample generation and model construction method facing transfer learning
CN115454024A (en) * 2022-08-29 2022-12-09 华北电力科学研究院有限责任公司 Security and stability control system strategy verification method and device

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
LINGYU REN,ET AL.: "The Online Application and Fast Solving Method for Critical Clearing Time of Three-Phase Short Circuit in Power System", 《2012 (2ND) INTERNATIONAL CONFERENCE ON SMART GRID AND CLEAN ENERGY TECHNOLOGIES(ICSGCE 2012)》 *
XIA Y,ET AL.: "Direct nonlinear primal-dual interior-point method for transient stability constrained optimal power flow", 《IEEE PROCEEDINGS: GENERATION TRANSMISSION AND DISTRIBUTION》 *
YORINO N, ET AL.: "A new method for obtaining critical clearing time for transient stability", 《POWER SYSTEMS, IEEE TRANSACTIONS ON》 *
吴政球等: "电力系统静态电压稳定极限及裕度计算综述", 《电力系统及其自动化学报》 *
张晓东等: "基于三道防线的故障集自动批量生成方法", 《电力系统及其自动化》 *

Cited By (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103413043A (en) * 2013-08-09 2013-11-27 中国南方电网有限责任公司 Solving method for power system multi-dimensional space cross section thermal stability limit boundary
CN103413043B (en) * 2013-08-09 2016-04-20 中国南方电网有限责任公司 A kind of method for solving of electric system hyperspace section thermal stability limit boundary
CN103823998A (en) * 2014-03-19 2014-05-28 国家电网公司 Weak cross section determination method taking influence of network topology changes on power transmission capacity into consideration
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
CN107767028A (en) * 2017-09-26 2018-03-06 广东电网有限责任公司电力调度控制中心 A kind of transmission cross-section limit automatic identification based on xml and out-of-limit computational methods
CN107767028B (en) * 2017-09-26 2021-09-07 广东电网有限责任公司电力调度控制中心 Transmission section limit automatic identification and out-of-limit calculation method based on xml
CN110943453A (en) * 2019-12-23 2020-03-31 北京交通大学 Power system fault sample generation and model construction method facing transfer learning
CN115454024A (en) * 2022-08-29 2022-12-09 华北电力科学研究院有限责任公司 Security and stability control system strategy verification method and device

Also Published As

Publication number Publication date
CN103065060B (en) 2015-10-28

Similar Documents

Publication Publication Date Title
CN103065060A (en) Computing method for transient stability limit of transmission section of power system
Jadhav et al. Performance verification of PID controller in an interconnected power system using particle swarm optimization
Moghaddam et al. Optimal sizing of hybrid energy storage systems to mitigate wind power fluctuations
CN103296688A (en) Method for implementing power system stabilizer
Gao et al. Improving frequency stability in low inertia power systems using synthetic inertia from wind turbines
Setiadi et al. Enabling resilient wide‐area POD at BESS in Java, Indonesia 500 kV power grid
Jaber et al. A new parameters identification of single area power system based LFC using Segmentation Particle Swarm Optimization (SePSO) algorithm
CN104378033A (en) Method for tracing maximum power point of switch reluctance motor wind power generation system
CN106875293B (en) A kind of wind power plant booster stations main transformer failure generated energy loss acquisition methods
CN103344881A (en) Grid fault diagnosing method based on combined type cross entropy algorithm
JP2015532578A (en) Multiple input single output combined power generation system
CN104868490A (en) Transient stability index-based direct current modulation controller parameter optimization method
CN113359468A (en) Wind turbine generator fault-tolerant control method based on robust self-adaption and sliding mode variable structure control
Arioua et al. Control of high voltage direct current links with overall large‐scale grid objectives
CN112801533A (en) Power system operation reliability assessment method considering decision dependence uncertainty
JP2016067101A (en) Change mitigation storage battery system for renewable energy
CN104238624B (en) A kind of tracking of maximum power point and device
JP2020137164A (en) Power system stabilization system
Murat et al. Validation of hydroelectric power plant model for speed governor development studies
Gerin-Lajoie et al. Hydro-Quebec experience with PSS tuning
Feng et al. Sensor fault tolerant control of a wind turbine via Takagi-Sugeno fuzzy observer and model predictive control
CN106911323A (en) Triangular-wave generator and triangular wave generation system
Chamorro et al. Dynamic measurements of the wind power impact on power system inertia and stability
Hashimoto et al. A new control method of multiple sets of battery energy storage system for load frequency control in power system with massive integration of renewable energy resources
Wang et al. DFIG WT electromechanical transient behaviour influenced by PLL: Modelling and analysis

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