CN106981872B - A kind of Operation of Electric Systems reliability determination method and device - Google Patents
A kind of Operation of Electric Systems reliability determination method and device Download PDFInfo
- Publication number
- CN106981872B CN106981872B CN201710237221.4A CN201710237221A CN106981872B CN 106981872 B CN106981872 B CN 106981872B CN 201710237221 A CN201710237221 A CN 201710237221A CN 106981872 B CN106981872 B CN 106981872B
- Authority
- CN
- China
- Prior art keywords
- power grid
- probability
- malfunction
- load
- power
- 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
Links
- 238000000034 method Methods 0.000 title claims abstract description 30
- 230000007257 malfunction Effects 0.000 claims abstract description 83
- 238000005520 cutting process Methods 0.000 claims abstract description 37
- 230000005611 electricity Effects 0.000 claims description 11
- 230000000694 effects Effects 0.000 claims description 10
- 230000007812 deficiency Effects 0.000 claims description 8
- 230000005540 biological transmission Effects 0.000 claims description 7
- 240000002853 Nelumbo nucifera Species 0.000 claims description 5
- 235000006508 Nelumbo nucifera Nutrition 0.000 claims description 5
- 235000006510 Nelumbo pentapetala Nutrition 0.000 claims description 5
- 238000013459 approach Methods 0.000 claims description 5
- 238000005303 weighing Methods 0.000 claims 2
- 230000008859 change Effects 0.000 abstract description 3
- 238000012546 transfer Methods 0.000 abstract description 2
- 238000010586 diagram Methods 0.000 description 13
- 238000004590 computer program Methods 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 238000012545 processing Methods 0.000 description 5
- 238000005516 engineering process Methods 0.000 description 4
- 230000006870 function Effects 0.000 description 4
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000003860 storage Methods 0.000 description 2
- FNMKZDDKPDBYJM-UHFFFAOYSA-N 3-(1,3-benzodioxol-5-yl)-7-(3-methylbut-2-enoxy)chromen-4-one Chemical compound C1=C2OCOC2=CC(C2=COC=3C(C2=O)=CC=C(C=3)OCC=C(C)C)=C1 FNMKZDDKPDBYJM-UHFFFAOYSA-N 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 238000004422 calculation algorithm Methods 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000012502 risk assessment Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Classifications
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J3/00—Circuit arrangements for ac mains or ac distribution networks
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02H—EMERGENCY PROTECTIVE CIRCUIT ARRANGEMENTS
- H02H7/00—Emergency protective circuit arrangements specially adapted for specific types of electric machines or apparatus or for sectionalised protection of cable or line systems, and effecting automatic switching in the event of an undesired change from normal working conditions
- H02H7/26—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured
- H02H7/261—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations
- H02H7/262—Sectionalised protection of cable or line systems, e.g. for disconnecting a section on which a short-circuit, earth fault, or arc discharge has occured involving signal transmission between at least two stations involving transmissions of switching or blocking orders
-
- H—ELECTRICITY
- H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
- H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
- H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
- H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
Landscapes
- Engineering & Computer Science (AREA)
- Power Engineering (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The present invention relates to a kind of Operation of Electric Systems reliability determination method and devices, which comprises determines the probability of malfunction of each element in power grid；The probability of malfunction of power grid, loss load and cutting load coefficient when according to the probability of malfunction of each element and determining each element fault using dynamic power flow；Operation of Electric Systems reliability is determined according to the probability of malfunction of power grid, loss load and cutting load coefficient when each element fault；Technical solution provided by the invention, it is contemplated that system sectionalizing or it is grid-connected after frequency change procedure and failure occur after power flow transfer situation, the element operating status being calculated is more accurate compared to conventional Load Flow；Simultaneously, this method utilizes the operational reliability model of power system component, probability of malfunction of the element under different operating statuses is considered, more accurate multiple failure probability can be obtained, the operational reliability result being calculated has highly important directive significance to Operation of Electric Systems scheduling.
Description
Technical field
The present invention relates to Operation of Electric Systems analysis fields, and in particular to a kind of determination side of Operation of Electric Systems reliability
Method and device.
Background technique
With the rapid growth of Electricity Demand, transferring electricity from the west to the east, on national network implementation, the interconnection of China's power grid
Scale is growing, complexity is continuously improved.It improves large-scale interconnected power system operational reliability, guarantee that power grid security has become China
The key and urgency problem that electric system faces.
Electric power system operation risk assessment is one of modern interconnected network planning key technology, the technology by quantitative analysis not
Carry out influence of the uncertain factors such as electric power demand growth, element fault to power grid pattern.Conventional reliability algorithm is based on element
Stable state statistical parameter, state is carried out to typical system running pattern and enumerates or simulates, to assess its state consequence, but this
When kind method is analyzed in large-scale electrical power system, the abundant intensity of system is very big, it is contemplated that all spare capacities and extension wire
Put into operation, however Operation of Electric Systems analyze during, specific run mode of the people more concerned with electric system, it is desirable to certain
A little methods of operation are assessed and the assessment of operational reliability.
Summary of the invention
The present invention provides a kind of Operation of Electric Systems reliability determination method and device, and the purpose is to more accurately point
The operating status for analysing element after failure occurs then obtains the probability of malfunction after failure occurs.
The purpose of the present invention is adopt the following technical solutions realization:
A kind of Operation of Electric Systems reliability determination method, it is improved in that including:
Determine the probability of malfunction of each element in power grid；
The probability of malfunction of power grid when according to the probability of malfunction of each element and determining each element fault using dynamic power flow,
Lose load and cutting load coefficient；
Electric system is determined according to the probability of malfunction of power grid, loss load and cutting load coefficient when each element fault
Operational reliability.
Preferably, in the determining power grid each element probability of malfunction, comprising:
If element is route in power grid, the probability of malfunction P (L) of route is determined as the following formula:
In above formula,For the annual probability of malfunction of route, L is line transmission power,For under route nominal situation most
Big transimission power, LmaxFor the maximum transmission power that can be born under route extreme operating condition；
If element is generator in power grid, the probability of malfunction P (V) of generator is determined as the following formula:
In above formula,For the annual probability of malfunction of generator, VgFor generator voltage,For the normal work of generator
Minimum under condition, Vg,minFor the minimum that can be born under generator extreme operating condition,It is normal for generator
Maximum working voltage under operating condition, Vg,maxFor the maximum working voltage that can be born under generator extreme operating condition.
Preferably, according to the probability of malfunction of each element and when determining each element fault using dynamic power flow power grid therefore
Hinder probability, loss load and cutting load coefficient, comprising:
The probability of malfunction P of power grid when a. initializing k-th of element faultk=FkAnd the loss load L of power gridk=0,
FkFor the probability of malfunction of k-th of element in power grid；
B. electric network data when k-th of element fault is read in；
C. electric network swim is calculated using Newton-Raphson approach, obtains the active-power P of equalization pointswing, and judge equalization point
Active PswingWhether preset condition is met, if so, step e is gone to, if it is not, then enabling accelerating power Pacc=Pswing-P0, P0
For the setting power of balance nodes；
D. the power output of the generator of renewal frequency deviation delta f and power grid interior joint and load power output, and return step c；Its
In, frequency deviation f initial value is 0；
E. output electric network swim is as a result, and update system frequency f=f0Δ f, wherein f0For the reference frequency of system；
If f. the voltage of each node and system frequency f are out-of-limit in power grid, cutting load is carried out to power grid, and update k-th
The loss load L of power grid when element faultk=Lk+Lx, wherein LxFor minimum tangential load amount, if in power grid each node electricity
Pressure and system frequency f be not out-of-limit, then does not carry out cutting load to power grid；
G. the out-of-limit element of power grid is searched for, and if it exists, then disconnect the out-of-limit element x selected at random, and update the failure of power grid
Probability Pk=Pk*Fx, go to step b, wherein FxIt is then gone to if it does not exist for the probability of malfunction of the out-of-limit element x selected at random
Step h；
H. according to the loss load L of power grid when k-th of element faultkWhen determining k-th of element fault power grid cut it is negative
Lotus coefficient Ik；
I. judge whether k is less than or equal to m, if so, k=k+1 and return step a are enabled, if it is not, then end operation, output
Final calculation result.
Further, the preset condition is | Pswing-Plaxt|≤εmax, temporary variable Plast=P0, εmaxFor preset mistake
Poor upper limit value.
Further, renewal frequency deviation delta f as the following formula:
In above formula, KGiCorresponding revolving speed offset, KL when changing for the unit load of i-th of nodeiFor having for i-th node
The frequency mediating effect+6 coefficient of workload, N are grid nodes sum.
Further, in the step d, the power output and load power output for updating the generator of power grid interior joint include:
The power output PG of the generator of i-th of node in power grid is updated as the following formulai:
PGi=PGi0+KGi·Δf
In above formula, PGi0For the original power output of the generator of i-th of node in power grid, Δ f is frequency departure, KGiIt is i-th
Corresponding revolving speed offset when changing of the unit load of a node；
The load power output PL of i-th of node in power grid is updated as the following formulai:
PLi=PLi0+KLi·Δf
In above formula, PLi0For the original power output of the load of i-th of node in power grid, KLiFor the burden with power of i-th of node
Frequency mediating effect+6 coefficient.
Preferably, described according to the probability of malfunction of power grid, loss load and cutting load coefficient when each element fault
Determine Operation of Electric Systems reliability, comprising:
Determine that the mistake load of electric system is general according to the probability of malfunction of power grid when each element fault and cutting load coefficient
Rate；
Determine that the electric energy of electric system is insufficient according to the probability of malfunction of power grid when each element fault and loss load
Desired value.
Further, the load-loss probability LOLP of electric system is determined as the following formula:
In above formula, PkThe loss load of power grid when for k-th of element fault, IkPower grid when for k-th of element fault
Cutting load coefficient, m are component population in power grid；
The electric energy deficiency desired value EENS of electric system is determined as the following formula:
In above formula, LkThe loss load of power grid when for k-th of element fault.
A kind of determining device of Operation of Electric Systems reliability, it is improved in that described device includes:
First determining module, for determining the probability of malfunction of each element in power grid；
Second determining module for the probability of malfunction according to each element and determines each element fault using dynamic power flow
When the probability of malfunction of power grid, loss load and cutting load coefficient；
Third determining module, for according to the probability of malfunction of power grid when each element fault, lose load and cut negative
Lotus coefficient determines Operation of Electric Systems reliability.
Preferably, the third determining module, comprising:
First determination unit, for determining electricity according to the probability of malfunction and cutting load coefficient of power grid when each element fault
The load-loss probability of Force system；
Second determination unit, for determining electricity according to the probability of malfunction and loss load of power grid when each element fault
The electric energy deficiency desired value of Force system.
Beneficial effects of the present invention:
Technical solution provided by the invention is based on dynamic power flow technology, when calculating Operation of Electric Systems reliability, uses
Dynamic power flow analyzes operating status of the power grid under different faults in real time, then according to operation of power networks state, in conjunction with electric system
The probability of malfunction of element obtains the probability that electric network fault occurs, the operational reliability of electric system is calculated.This method considers
System sectionalizing or it is grid-connected after frequency change procedure and failure occur after power flow transfer situation, the element that is calculated
Operating status is more accurate compared to conventional Load Flow；Meanwhile this method utilizes the probability of malfunction of power system component, Neng Gougeng
Add the probability of malfunction for accurately determining power grid, the operational reliability result being calculated has very Operation of Electric Systems scheduling
Important directive significance, further, technical solution provided by the invention can make up existing Operation of Electric Systems reliability meter
The deficiency of calculation provides effective calculating instrument for operational reliability calculating, operational reliability computational accuracy is improved, to electric system
Safe and stable operation has highly important use value and meaning.
Detailed description of the invention
Fig. 1 is a kind of flow chart of Operation of Electric Systems reliability determination method of the present invention；
Fig. 2 is the operational reliability model characteristics curve synoptic diagram of route in the embodiment of the present invention；
Fig. 3 is the operational reliability model characteristics curve synoptic diagram of generator in the embodiment of the present invention；
Fig. 4 is that Operation of Electric Systems reliability determines flow diagram in the embodiment of the present invention；
Fig. 5 is corresponding revolving speed offset characteristic curve synoptic diagram when unit load changes in the embodiment of the present invention；
Fig. 6 is the frequency mediating effect+6 coefficient feature curve synoptic diagram of burden with power in the embodiment of the present invention；
Fig. 7 is a kind of schematic device of Operation of Electric Systems reliability determination method of the present invention.
Specific embodiment
It elaborates with reference to the accompanying drawing to a specific embodiment of the invention.
In order to make the object, technical scheme and advantages of the embodiment of the invention clearer, below in conjunction with the embodiment of the present invention
In attached drawing, technical scheme in the embodiment of the invention is clearly and completely described, it is clear that described embodiment is
A part of the embodiment of the present invention, instead of all the embodiments.Based on the embodiments of the present invention, those of ordinary skill in the art
All other embodiment obtained without making creative work, shall fall within the protection scope of the present invention.
A kind of Operation of Electric Systems reliability determination method provided by the invention is based on dynamic power flow technology, is calculating
When Operation of Electric Systems reliability, operating status of the power grid under different faults is analyzed in real time using dynamic power flow, then basis
Operation of power networks state obtains the probability of malfunction of different elements in conjunction with the operational reliability model of power system component, thus to obtain
The probability that multiple failure occurs, is calculated the operational reliability of electric system, as shown in Figure 1, comprising:
101. determining the probability of malfunction of each element in power grid；
102. the failure of power grid when according to the probability of malfunction of each element and determining each element fault using dynamic power flow
Probability, loss load and cutting load coefficient；
103. determining electric power according to the probability of malfunction of power grid, loss load and cutting load coefficient when each element fault
System operation reliability.
Specifically, in the step 101, can use route operational reliability model and generator it is reliable for operation
Property model determines the probability of malfunction of each element in power grid, wherein operational reliability model linearity curve schematic diagram such as Fig. 2 institute of route
Show, the operational reliability model linearity curve schematic diagram of generator is as shown in figure 3, specifically include:
If element is route in power grid, the probability of malfunction P (L) of route is determined as the following formula:
In above formula,For the annual probability of malfunction of route, L is line transmission power,For under route nominal situation
Maximum transmission power, LmaxFor the maximum transmission power that can be born under route extreme operating condition；
If element is generator in power grid, the probability of malfunction P (V) of generator is determined as the following formula:
In above formula,For the annual probability of malfunction of generator, VgFor generator voltage,For the normal work of generator
Minimum under condition, Vg,minFor the minimum that can be born under generator extreme operating condition,It is normal for generator
Maximum working voltage under operating condition, Vg,maxFor the maximum working voltage that can be born under generator extreme operating condition.
It determines in power grid after the probability of malfunction of each element, the step 102, as shown in Figure 4, comprising:
The probability of malfunction P of power grid when a. initializing k-th of element faultk=FkAnd the loss load L of power gridk=0,
FkFor the probability of malfunction of k-th of element in power grid；
Wherein, electric network data includes title, the voltage class of power grid median generatrix node；Line length, impedance parameter；Transformation
Device impedance parameter；Generated power is idle maximin, active setting value, difference coefficient KG；Active and idle, the frequency of load
Rate mediating effect+6 COEFFICIENT KL；
KG is that can quantitatively show corresponding revolving speed offset when unit load changes, characteristic curve as shown in figure 5, its
In, PGNFor generating set active power, fNFor generating set frequency, f0=50Hz is the reference frequency of system, for example, working as KG=
5%, if active power changes 1%, frequency shift (FS) 0.05%；Change 20% if any power, then frequency shift (FS) 1% (0.5Hz).
Steam turbine KG=3%~5%, hydraulic turbine KG=2%~4%；
KL is the frequency mediating effect+6 coefficient of burden with power or the frequency mediating effect+6 of abbreviation load, and characteristic curve is such as
Shown in Fig. 6, wherein PLNFor load active power, general electric system KL=1~3, when it indicates frequency variation 1%, load has
It is 1%~3% that function power changes accordingly；
B. electric network data when k-th of element fault is read in；
C. electric network swim is calculated using Newton-Raphson approach, obtains the active-power P of equalization pointswing, and judge equalization point
Active PswingWhether preset condition is met, if so, step e is gone to, if it is not, then enabling accelerating power Pacc=Pswing-P0, P0
For the setting power of balance nodes；
Wherein, the preset condition is | Pswing-Plaxt|≤εmax, temporary variable Plast=P0, εmaxFor in preset error
Limit value.
D. the power output of the generator of renewal frequency deviation delta f and power grid interior joint and load power output, and return step c；Its
In, frequency deviation f initial value is 0；
E. output electric network swim is as a result, and update system frequency f=f0Δ f, wherein f0For the reference frequency of system；
If f. the voltage of each node and system frequency f are out-of-limit in power grid, cutting load is carried out to power grid, and update k-th
The loss load L of power grid when element faultk=Lk+Lx, wherein LxFor minimum tangential load amount, if in power grid each node electricity
Pressure and system frequency f be not out-of-limit, then does not carry out cutting load to power grid；
G. the out-of-limit element of power grid is searched for, and if it exists, then disconnect the out-of-limit element x selected at random, and update the failure of power grid
Probability Pk=Pk*Fx, go to step b, wherein FxIt is then gone to if it does not exist for the probability of malfunction of the out-of-limit element x selected at random
Step h；
Wherein, out-of-limit element includes: the out-of-limit route of load and voltage or frequency out-of-limit generator；
H. according to the loss load L of power grid when k-th of element faultkWhen determining k-th of element fault power grid cut it is negative
Lotus coefficient Ik；
I. judge whether k is less than or equal to m, if so, k=k+1 and return step a are enabled, if it is not, then end operation, output
Final calculation result.
Further, the Load flow calculation in the step is conventional Load Flow calculating, and Newton-Raphson approach, PQ solution can be used
The methods of coupling method is calculated, and the acquisition of the Load Flow Calculation Softwares such as BPA, PSASP, PSS/E also can be used；
In the step d, renewal frequency deviation delta f as the following formula:
In above formula, KGiCorresponding revolving speed offset, KL when changing for the unit load of i-th of nodeiFor having for i-th node
The frequency mediating effect+6 coefficient of workload, N are grid nodes sum.
In the step d, the power output and load power output for updating the generator of power grid interior joint include:
The power output PG of the generator of i-th of node in power grid is updated as the following formulai:
PGi=PGi0+KGi·Δf
In above formula, PGi0For the original power output of the generator of i-th of node in power grid, Δ f is frequency departure, KGiIt is i-th
Corresponding revolving speed offset when changing of the unit load of a node；
The load power output PL of i-th of node in power grid is updated as the following formulai:
PLi=PLi0+KLi·Δf
In above formula, PLi0For the original power output of the load of i-th of node in power grid, KLiFor the burden with power of i-th of node
Frequency mediating effect+6 coefficient.
When obtaining each element fault after the probability of malfunction of power grid, loss load and cutting load coefficient, the step 103,
Include:
Determine that the mistake load of electric system is general according to the probability of malfunction of power grid when each element fault and cutting load coefficient
Rate；
Determine that the electric energy of electric system is insufficient according to the probability of malfunction of power grid when each element fault and loss load
Desired value.
Wherein, the load-loss probability LOLP of electric system is determined as the following formula:
In above formula, PkThe loss load of power grid when for k-th of element fault, IkPower grid when for k-th of element fault
Cutting load coefficient, m are component population in power grid；
The electric energy deficiency desired value EENS of electric system is determined as the following formula:
In above formula, LkThe loss load of power grid when for k-th of element fault.
The present invention also provides a kind of determining devices of Operation of Electric Systems reliability, as shown in fig. 7, described device includes:
First determining module, for determining the probability of malfunction of each element in power grid；
Second determining module for the probability of malfunction according to each element and determines each element fault using dynamic power flow
When the probability of malfunction of power grid, loss load and cutting load coefficient；
Third determining module, for according to the probability of malfunction of power grid when each element fault, lose load and cut negative
Lotus coefficient determines Operation of Electric Systems reliability.
Wherein, the third determining module, comprising:
First determination unit, for determining electricity according to the probability of malfunction and cutting load coefficient of power grid when each element fault
The load-loss probability of Force system；
Second determination unit, for determining electricity according to the probability of malfunction and loss load of power grid when each element fault
The electric energy deficiency desired value of Force system.
It should be understood by those skilled in the art that, embodiments herein can provide as method, system or computer program
Product.Therefore, complete hardware embodiment, complete software embodiment or reality combining software and hardware aspects can be used in the application
Apply the form of example.Moreover, it wherein includes the computer of computer usable program code that the application, which can be used in one or more,
The computer program implemented in usable storage medium (including but not limited to magnetic disk storage, CD-ROM, optical memory etc.) produces
The form of product.
The application is referring to method, the process of equipment (system) and computer program product according to the embodiment of the present application
Figure and/or block diagram describe.It should be understood that every one stream in flowchart and/or the block diagram can be realized by computer program instructions
The combination of process and/or box in journey and/or box and flowchart and/or the block diagram.It can provide these computer programs
Instruct the processor of general purpose computer, special purpose computer, Embedded Processor or other programmable data processing devices to produce
A raw machine, so that being generated by the instruction that computer or the processor of other programmable data processing devices execute for real
The device for the function of being specified in present one or more flows of the flowchart and/or one or more blocks of the block diagram.
These computer program instructions, which may also be stored in, is able to guide computer or other programmable data processing devices with spy
Determine in the computer-readable memory that mode works, so that it includes referring to that instruction stored in the computer readable memory, which generates,
Enable the manufacture of device, the command device realize in one box of one or more flows of the flowchart and/or block diagram or
The function of being specified in multiple boxes.
These computer program instructions also can be loaded onto a computer or other programmable data processing device, so that counting
Series of operation steps are executed on calculation machine or other programmable devices to generate computer implemented processing, thus in computer or
The instruction executed on other programmable devices is provided for realizing in one or more flows of the flowchart and/or block diagram one
The step of function of being specified in a box or multiple boxes.
Finally it should be noted that: the above embodiments are merely illustrative of the technical scheme of the present invention and are not intended to be limiting thereof, to the greatest extent
Invention is explained in detail referring to above-described embodiment for pipe, it should be understood by those ordinary skilled in the art that: still
It can be with modifications or equivalent substitutions are made to specific embodiments of the invention, and without departing from any of spirit and scope of the invention
Modification or equivalent replacement, should all cover within the scope of the claims of the present invention.
Claims (9)
1. a kind of Operation of Electric Systems reliability determination method, which is characterized in that the described method includes:
Determine the probability of malfunction of each element in power grid；
The probability of malfunction of power grid, loss when according to the probability of malfunction of each element and determining each element fault using dynamic power flow
Load and cutting load coefficient；
Operation of Electric Systems is determined according to the probability of malfunction of power grid, loss load and cutting load coefficient when each element fault
Reliability；
The probability of malfunction according to each element and the probability of malfunction of power grid when determining each element fault using dynamic power flow,
Lose load and cutting load coefficient, comprising:
The probability of malfunction P of power grid when a. initializing k-th of element faultk=FkAnd the loss load L of power gridk=0, FkFor electricity
The probability of malfunction of k-th of element in net；
B. electric network data when k-th of element fault is read in；
C. electric network swim is calculated using Newton-Raphson approach, obtains the active-power P of equalization pointswing, and judge having for equalization point
Function PswingWhether preset condition is met, if so, step e is gone to, if it is not, then enabling accelerating power Pacc=Pswing-P0, P0It is flat
The setting power of weighing apparatus node；
D. the power output of the generator of renewal frequency deviation delta f and power grid interior joint and load power output, and return step c；Wherein, frequency
Rate deviation delta f initial value is 0；
E. output electric network swim is as a result, and update system frequency f=f0Δ f, wherein f0For the reference frequency of system；
If f. the voltage of each node and system frequency f are out-of-limit in power grid, cutting load is carried out to power grid, and update k-th of element
The loss load L of power grid when failurek=Lk+Lx, wherein LxFor minimum tangential load amount, if in power grid the voltage of each node and
System frequency f is not out-of-limit, then does not carry out cutting load to power grid；
G. the out-of-limit element of power grid is searched for, and if it exists, then disconnect the out-of-limit element x selected at random, and update the probability of malfunction P of power gridk
=Pk*Fx, go to step b, wherein FxStep h is then gone to if it does not exist for the probability of malfunction of the out-of-limit element x selected at random；
H. according to the loss load L of power grid when k-th of element faultkThe cutting load coefficient of power grid when determining k-th of element fault
Ik；
I. judge whether k is less than or equal to m, if so, enabling k=k+1 and return step a, if it is not, then end operation, output is final
Calculated result.
2. the method as described in claim 1, which is characterized in that the probability of malfunction of each element in the determining power grid, comprising:
If element is route in power grid, the probability of malfunction P (L) of route is determined as the following formula:
In above formula,For the annual probability of malfunction of route, L is line transmission power,It is maximum under route nominal situation
Transimission power, LmaxFor the maximum transmission power that can be born under route extreme operating condition；
If element is generator in power grid, the probability of malfunction P (V) of generator is determined as the following formula:
In above formula,For the annual probability of malfunction of generator, VgFor generator voltage,For under generator nominal situation
Minimum, Vg,minFor the minimum that can be born under generator extreme operating condition,For generator nominal situation
Lower maximum working voltage, Vg,maxFor the maximum working voltage that can be born under generator extreme operating condition.
3. the method as described in claim 1, which is characterized in that the preset condition is | Pswing-Plaxt|≤εmax, interim to become
Measure Plast=P0, εmaxFor preset error upper limit value.
4. the method as described in claim 1, which is characterized in that renewal frequency deviation delta f as the following formula:
In above formula, KGiCorresponding revolving speed offset, KL when changing for the unit load of i-th of nodeiIt is born for i-th of the active of node
The frequency mediating effect+6 coefficient of lotus, N are grid nodes sum.
5. the method as described in claim 1, which is characterized in that in the step d, update going out for the generator of power grid interior joint
Power and load power output include:
The power output PG of the generator of i-th of node in power grid is updated as the following formulai:
PGi=PGi0+KGi·Δf
In above formula, PGi0For the original power output of the generator of i-th of node in power grid, Δ f is frequency departure, KGiFor i-th of section
Corresponding revolving speed offset when changing of the unit load of point；
The load power output PL of i-th of node in power grid is updated as the following formulai:
PLi=PLi0+KLi·Δf
In above formula, PLi0For the original power output of the load of i-th of node in power grid, KLiFor the frequency of the burden with power of i-th of node
Rate mediating effect+6 coefficient.
6. the method as described in claim 1, which is characterized in that the failure according to power grid when each element fault is general
Rate, loss load and cutting load coefficient determine Operation of Electric Systems reliability, comprising:
The load-loss probability of electric system is determined according to the probability of malfunction of power grid when each element fault and cutting load coefficient；
The electric energy deficiency expectation of electric system is determined according to the probability of malfunction of power grid when each element fault and loss load
Value.
7. method as claimed in claim 6, which is characterized in that determine the load-loss probability LOLP of electric system as the following formula:
In above formula, PkThe probability of malfunction of power grid, I when for k-th of element faultkThe cutting load system of power grid when for k-th of element fault
Number, m are component population in power grid；
The electric energy deficiency desired value EENS of electric system is determined as the following formula:
In above formula, LkThe loss load of power grid when for k-th of element fault.
8. a kind of determining device of Operation of Electric Systems reliability, which is characterized in that described device includes:
First determining module, for determining the probability of malfunction of each element in power grid；
Second determining module, electricity when for according to the probability of malfunction of each element and determining each element fault using dynamic power flow
Probability of malfunction, loss load and the cutting load coefficient of net；
Third determining module, for according to the probability of malfunction of power grid, loss load and cutting load system when each element fault
Number determines Operation of Electric Systems reliability；
The probability of malfunction according to each element and the probability of malfunction of power grid when determining each element fault using dynamic power flow,
Lose load and cutting load coefficient, comprising:
The probability of malfunction P of power grid when a. initializing k-th of element faultk=FkAnd the loss load L of power gridk=0, FkFor electricity
The probability of malfunction of k-th of element in net；
B. electric network data when k-th of element fault is read in；
C. electric network swim is calculated using Newton-Raphson approach, obtains the active-power P of equalization pointswing, and judge having for equalization point
Function PswingWhether preset condition is met, if so, step e is gone to, if it is not, then enabling accelerating power Pacc=Pswing-P0, P0It is flat
The setting power of weighing apparatus node；
D. the power output of the generator of renewal frequency deviation delta f and power grid interior joint and load power output, and return step c；Wherein, frequency
Rate deviation delta f initial value is 0；
E. output electric network swim is as a result, and update system frequency f=f0Δ f, wherein f0For the reference frequency of system；
If f. the voltage of each node and system frequency f are out-of-limit in power grid, cutting load is carried out to power grid, and update k-th of element
The loss load L of power grid when failurek=Lk+Lx, wherein LxFor minimum tangential load amount, if in power grid the voltage of each node and
System frequency f is not out-of-limit, then does not carry out cutting load to power grid；
G. the out-of-limit element of power grid is searched for, and if it exists, then disconnect the out-of-limit element x selected at random, and update the probability of malfunction P of power gridk
=Pk*Fx, go to step b, wherein FxStep h is then gone to if it does not exist for the probability of malfunction of the out-of-limit element x selected at random；
H. according to the loss load L of power grid when k-th of element faultkThe cutting load coefficient of power grid when determining k-th of element fault
Ik；
I. judge whether k is less than or equal to m, if so, enabling k=k+1 and return step a, if it is not, then end operation, output is final
Calculated result.
9. device as claimed in claim 8, which is characterized in that the third determining module, comprising:
First determination unit, for determining power train according to the probability of malfunction and cutting load coefficient of power grid when each element fault
The load-loss probability of system；
Second determination unit, for determining power train according to the probability of malfunction and loss load of power grid when each element fault
The electric energy deficiency desired value of system.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710237221.4A CN106981872B (en) | 2017-04-12 | 2017-04-12 | A kind of Operation of Electric Systems reliability determination method and device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201710237221.4A CN106981872B (en) | 2017-04-12 | 2017-04-12 | A kind of Operation of Electric Systems reliability determination method and device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN106981872A CN106981872A (en) | 2017-07-25 |
| CN106981872B true CN106981872B (en) | 2019-11-01 |
Family
ID=59344274
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201710237221.4A Active CN106981872B (en) | 2017-04-12 | 2017-04-12 | A kind of Operation of Electric Systems reliability determination method and device |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN106981872B (en) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108376984B (en) * | 2018-01-25 | 2021-03-23 | 云南电网有限责任公司 | Automatic adjustment method for node voltage of power system |
| CN108446816B (en) * | 2018-02-01 | 2024-07-02 | 中国电力科学研究院有限公司 | Reliability assessment method and system for energy hub model |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2520800A1 (en) * | 2011-05-02 | 2012-11-07 | Kabushiki Kaisha Toshiba | Wind power generation system and control method for the same |
| CN103426056A (en) * | 2013-07-18 | 2013-12-04 | 清华大学 | Power system weak link identification method based on risk evaluation |
| CN103440400A (en) * | 2013-07-29 | 2013-12-11 | 清华大学 | Power system short-term risk determination method taking disaster factors into account |
| CN106327034A (en) * | 2015-06-18 | 2017-01-11 | 中国电力科学研究院 | Cascading failure search and weak link analysis method based on operation reliability model |
-
2017
- 2017-04-12 CN CN201710237221.4A patent/CN106981872B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| EP2520800A1 (en) * | 2011-05-02 | 2012-11-07 | Kabushiki Kaisha Toshiba | Wind power generation system and control method for the same |
| CN103426056A (en) * | 2013-07-18 | 2013-12-04 | 清华大学 | Power system weak link identification method based on risk evaluation |
| CN103440400A (en) * | 2013-07-29 | 2013-12-11 | 清华大学 | Power system short-term risk determination method taking disaster factors into account |
| CN106327034A (en) * | 2015-06-18 | 2017-01-11 | 中国电力科学研究院 | Cascading failure search and weak link analysis method based on operation reliability model |
Non-Patent Citations (3)
| Title |
|---|
| Reliability assessment of distribution system with innovative Smart Grid technology implementation;Ahmed S. Al-Abdulwahab;《2011 IEEE PES Conference on Innovative Smart Grid Technologies - Middle East》;20111220;正文1-4页 * |
| 含分布式电源配电网可靠性评估的点估计法;芦晶晶，赵渊;《电网技术》;20130831;第37卷(第8期);2250-2257 * |
| 电力***随机潮流及其安全评估应用研究综述;朱星阳，刘文霞;《电工技术学报》;20131031;第28卷(第10期);257-269 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN106981872A (en) | 2017-07-25 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Sun et al. | An online dynamic security assessment scheme using phasor measurements and decision trees | |
| Wang et al. | Decision tree based online stability assessment scheme for power systems with renewable generations | |
| CN104333005B (en) | Based on frequency dynamic Forecasting Methodology after the Power System Disturbances of support vector regression | |
| CN105610192B (en) | Consider the online methods of risk assessment of large-scale wind power access | |
| CN109102132A (en) | A kind of micro-capacitance sensor vulnerability assessment method and device based on probability theory | |
| CN108075470A (en) | A kind of power grid transient frequency feature Forecasting Methodology based on SFR and ELM | |
| CN109752629B (en) | Intelligent diagnosis method and system for power grid measurement problems | |
| CN106355308B (en) | A method of wind power integration system core equipment is recognized based on decision tree | |
| CN111969658B (en) | Defensive-conventional coordination planning method for power generation and transmission system considering wind power | |
| CN102510060B (en) | Computing method of frequency characteristic coefficient of electric system | |
| CN110492524A (en) | Active-control method for frequency and system without the communication soft direct join net of marine wind electric field | |
| CN102868224A (en) | Secondary network measurement and multimode decision-making method and device for intelligent substation | |
| CN106981872B (en) | A kind of Operation of Electric Systems reliability determination method and device | |
| CN103887795A (en) | Electrical power system real-time probabilistic load flow online computing method | |
| CN103344881A (en) | Grid fault diagnosing method based on combined type cross entropy algorithm | |
| CN202872464U (en) | Secondary network measurement and multimode decision apparatus of intelligent substation | |
| Chen et al. | Distribution system state estimation: A survey of some relevant work | |
| CN107346889A (en) | Consider that the load of a frequency modulation frequency modulation and minimum frequency deviation cuts down Optimized model construction method | |
| CN110445171A (en) | A kind of active power and frequency control method and system based on the soft lineal system of marine wind electric field | |
| Kamisan et al. | Forecasting Wind Speed Data by Using a Combination of ARIMA Model with Single Exponential Smoothing. | |
| CN105701265A (en) | Double-fed wind generator modeling method and apparatus | |
| CN117424210A (en) | Method, system and equipment for evaluating vulnerability of power CPS under typhoon disaster | |
| CN108183481B (en) | Method and system for rapidly judging stability of power grid based on deep learning | |
| CN108092272A (en) | A kind of voltage stabilization on-line monitoring method based on the Kalman filtering that fades | |
| CN114707789A (en) | Low-voltage distribution network high-frequency measurement data processing and topology identification method and system |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |