CN108269017A - A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number - Google Patents
A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number Download PDFInfo
- Publication number
- CN108269017A CN108269017A CN201810055727.8A CN201810055727A CN108269017A CN 108269017 A CN108269017 A CN 108269017A CN 201810055727 A CN201810055727 A CN 201810055727A CN 108269017 A CN108269017 A CN 108269017A
- Authority
- CN
- China
- Prior art keywords
- algorithms
- deeac
- analysis
- step number
- stability analysis
- 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.)
- Pending
Links
- 238000004458 analytical method Methods 0.000 title claims abstract description 94
- 238000000034 method Methods 0.000 title claims abstract description 60
- 230000001052 transient effect Effects 0.000 title claims abstract description 36
- 230000003044 adaptive effect Effects 0.000 title claims abstract description 14
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 108
- 230000010354 integration Effects 0.000 claims abstract description 32
- 238000004364 calculation method Methods 0.000 claims abstract description 15
- 238000012360 testing method Methods 0.000 claims description 21
- 238000009833 condensation Methods 0.000 claims description 8
- 230000005494 condensation Effects 0.000 claims description 8
- 230000008569 process Effects 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 claims description 7
- 238000005516 engineering process Methods 0.000 claims description 5
- 230000008859 change Effects 0.000 claims description 2
- 230000000694 effects Effects 0.000 description 11
- 238000011160 research Methods 0.000 description 7
- 238000013461 design Methods 0.000 description 4
- 230000007246 mechanism Effects 0.000 description 4
- 241001269238 Data Species 0.000 description 3
- 238000004445 quantitative analysis Methods 0.000 description 3
- 238000012216 screening Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 230000004888 barrier function Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000000605 extraction Methods 0.000 description 2
- 230000006872 improvement Effects 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 238000005457 optimization Methods 0.000 description 2
- 238000004451 qualitative analysis Methods 0.000 description 2
- 238000013139 quantization Methods 0.000 description 2
- 230000003068 static effect Effects 0.000 description 2
- 230000006978 adaptation Effects 0.000 description 1
- 230000002776 aggregation Effects 0.000 description 1
- 238000004220 aggregation Methods 0.000 description 1
- 230000003466 anti-cipated effect Effects 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000295 complement effect Effects 0.000 description 1
- 238000000205 computational method Methods 0.000 description 1
- 238000000354 decomposition reaction Methods 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 235000013399 edible fruits Nutrition 0.000 description 1
- 230000005611 electricity Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 230000002045 lasting effect Effects 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 230000006641 stabilisation Effects 0.000 description 1
- 238000011105 stabilization Methods 0.000 description 1
- 239000002699 waste material Substances 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/06—Energy or water supply
Landscapes
- Business, Economics & Management (AREA)
- Human Resources & Organizations (AREA)
- Engineering & Computer Science (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- General Physics & Mathematics (AREA)
- Development Economics (AREA)
- Health & Medical Sciences (AREA)
- Educational Administration (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Theoretical Computer Science (AREA)
- Tourism & Hospitality (AREA)
- Physics & Mathematics (AREA)
- General Business, Economics & Management (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Game Theory and Decision Science (AREA)
- Public Health (AREA)
- Water Supply & Treatment (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a kind of fast transient Method of Stability Analysis based on Adaptive Integral step number, belong to Power System and its Automation technical field.The present invention is based on extended equal area criterions, the different nargin information included between the different each fast transient stability analysis algorithm of step number and its comparison result are integrated, come image study example time-varying degree by deep excavate, according to its power further determine whether that few integration step number need to be increased to improve analysis precision.The present invention can be that each example appropriately integrates step number according to the strong and weak Auto-matching of its time-varying degree, former base is compared in the fixed fast transient Method of Stability Analysis for integrating step number less, using smaller calculating increment as cost, under the premise of maintaining the weaker sample calculation analysis precision of time-varying degree constant, the analysis precision of the stronger example of time-varying degree is improved, to further coordinating the transient stability analysis problem of on-line transient stability analysis precision and rate, solution definitely and under uncertain factor with great theory and engineering significance.
Description
Technical field
The invention belongs to Power System and its Automation technical fields more particularly to a kind of based on Adaptive Integral step number
Fast transient Method of Stability Analysis.
Background technology
Electric power system transient stability is (referred to as:Temporarily surely) analysis quantization algorithm experienced from based on the assumption that the hair based on track
Exhibition course puts forth effort on the strong adaptation of various complex models and the rapid extraction of the temporary steady information of high-precision.Extra-high voltage and intelligence electricity
Net the iterative method built so that force device type and number continue to increase, this proposes temporary steady analytical performance more tight
The requirement of lattice.
Numerical computations are the bases of quantitative analysis, and common technology is hiding-trapezium integral method, and analysis precision and rate are always
It is conflict.Existing research is from the angles of optimization numerical computations to accelerating computation rate:Can increase integration step, it is common from
Moment point, the expansion of automatic changing step, decomposition aggregation, fast higher order Taylor series etc.;Also it can reduce integrating range, including base
Integral Technology is respectively terminated in advance in stable mechanism;The combination of novel integration method or distinct methods can be also introduced, such as implicit essence
Thin integration, Duhamel numerical integrations, projecting integral's method.The studies above has coordinated analysis precision and rate to a certain extent, so right
It is objective reality and can not evades in specific temporarily steady algorithm, the negative correlation between integration step and precision.
The technology that screens out for stablizing example is to coordinate the effective way of analysis precision and rate, it by preferred feature parameter,
Dependency rule is extracted, highly stable example is identified, so as to eliminate the calculation amount performed to it needed for detailed analysis.It is existing
Most of example screens out method under the premise of making every effort to retain catastrophe failure emulation in detail, reduces calculating to a certain extent
Cost.But, in order under different test systems and operating mode, ensure that higher example screens out performance, the design for screening out rule should
With relatively stringent theoretical foundation, i.e., need to be typically based on statistical analysis screen out incorporate in problem solving way it is certain because
Fruit is analyzed.
As three steps complementary in EEAC algorithm frames, static EEAC (SEEAC), dynamic EEAC (DEEAC)
Information excavating is carried out to gradually accurate disturbed track with three kinds of algorithms of integrated EEAC (IEEAC), asks for being drawn during disturbed track
The assumed condition entered is by by force to weak until disappearing.Wherein, deviation angle between each unit of same group is freezed by segmentation relaxation SEEAC algorithms
Approximating assumption, classical DEEAC is (referred to as:DEEAC (2)) algorithm established in failure, each 2 step Taylor series expansion after failure
On the basis of, the analysis precision of most of example is improved to increase a small amount of computation burden as cost, compared with SEEAC algorithms, embodies it
Coordinate the bridge effect of SEEAC algorithms rapidity and IEEAC algorithm accuracies, thus preferably coordinate the analysis of algorithm entirety
Precision and rate are the embodiments of EEAC algorithm marrow.Utilize its area with SEEAC algorithms in terms of time-varying factor ability is reflected
Not, research example time-varying factor can be roughly characterized, and further according to time-varying factor near using smaller calculating cost
Positively related causality analysis design screens out rule between introducing error like algorithm.
Accordingly, the scheme that example screens out is performed using quantitatively designing approximate data and being aided with certain mechanism deduction, specially
Sharp " transient stability evaluation in power system forecast failure collection rapid screening method " (grant number:ZL201310132812.7 it) devises
Stablize 3 criterions of example necessary condition, and form layering and screen out the strong efficiently identification devoid of risk example of frame, test result
Show the high efficiency and universality for screening out frame.On its basis, " transient stability evaluation in power system forecast failure is fast for patent
The strong sorting technique of speed " (grant number:ZL201410271454.2) further by example each in forecast failure complete or collected works be divided into stablize,
Doubtful stabilization, critical, doubtful unstability, 5 class of unstability, and apparent line temporarily surely analyzes the severe to calculating time requirement, flexibly
The example classification that can be screened out is selected, the example of fairly constant and fairly unstability is reliably screened out simultaneously so that need to perform in detail
The forecast failure number temporarily surely analyzed further reduces.
Approximate data is designed by dexterously selecting to simplify element, and deduced based on mechanism, excavate different approximate datas
Analytical conclusions information is to estimate its error range, is that the core of said program is thought so as to estimate the confidence level of approximate analysis conclusion
Road.Obviously, influence screens out and the key of classification performance is the design of approximate data.Ensureing the qualitative premise correctly classified
Under, for further reduction need to perform the example number of detailed analysis, dependent on the bridge effect for giving full play to DEEAC algorithms, into one
Step ground larger analysis precision lifting capacity is replaced with smaller calculating cost increment.
Invention content
The purpose of the present invention is:It screens out and taxonomy model performance, need to fill to further improve the example based on EEAC algorithms
The bridge effects of DEEAC algorithms is waved in distribution, research on utilization example time-varying degree and to meet the integration needed for certain integral accuracy
Positively related causality analysis between step number on the basis of the classical DEEAC algorithms based on the fixed step number of integration less, proposes a kind of base
In the fast transient Method of Stability Analysis of Adaptive Integral step number.
Specifically, the present invention adopts the following technical solutions realize, include the following steps:
1) after transient stability quick analysis system starts, some example in test example complete or collected works is taken out, is calculated using SEEAC
Method carries out nargin calculating to the example;
2) application classics DEEAC algorithms carry out nargin calculating to the example;
3) the example time-varying degree is reflected according to the difference of SEEAC algorithms and classics DEEAC algorithm nargin result of calculations;
4) if threshold epsilon is less than or equal to by the time-varying degree of step 3) reflection1, then it is assumed that apply classics DEEAC algorithms
It carries out nargin calculating acquired results and has met quick analysis precision, terminate the quick analysis process to the example, perform step 6),
Otherwise next step is performed;
5) application enhancements DEEAC algorithms carry out nargin calculating to the example, and think that acquired results have met quick analysis
Precision terminates the quick analysis process to the example, performs step 6);
6) if testing each example in example complete or collected works is complete quick analysis, terminates quick analysis process, otherwise take
Next example performs step 1).
Furthermore, application SEEAC algorithms are to the method for example progress nargin calculating in the step 1):
Assuming that two groups of inside of example complementation are the preferable people having the same aspiration and interest, leading Infinite bus power system is obtained based on model condensation technique
Image system is (referred to as:Leading OMIB image systems), using fault clearance time τ as step-length, asked based on single step Taylor series expansion
Rotor angle and acceleration of the leading OMIB image systems at the τ moment are taken, the kinetic energy that further parsing is acquired under SEEAC algorithms accelerates
AreaWith kinetic energy retardation areaStability margin is asked for by formula (1):
In above-mentioned formula, ηSE(τ) is the counted stability margin of example application SEEAC algorithms.
Furthermore, application classics DEEAC algorithms are to the method for example progress nargin calculating in the step 2):
Stage in failure and after failure is based on two step Taylor series expansions and asks for each machine rotor angle, and utilize track
Condensation technique obtains transient state parameter of the leading OMIB image systems at each integration moment, and the transient state parameter includes rotor angle with adding
Speed, the kinetic energy that further piecewise analytic is acquired under classical DEEAC algorithms accelerate areaWith kinetic energy retardation areaStability margin is asked for by formula (2):
In above-mentioned formula, ηDE(2)(τ) applies the counted stability margin of classics DEEAC algorithms for the example.
Furthermore, by comparing SEEAC algorithms and classics DEEAC algorithm nargin result of calculations in the step 3)
Difference is to reflect the method for the example time-varying degree:
The difference value σ of SEEAC algorithms and classics DEEAC algorithm nargin result of calculations is asked for by formula (3)1(τ) reflects
The example time-varying degree:
Furthermore, in the step 4), threshold epsilon1Value be:ε1=0.75.
Furthermore, application enhancements DEEAC algorithms are to the method for example progress nargin calculating in the step 5):
It in the stage in failure and after failure, is based respectively on two steps and five step Taylor series expansions asks for each machine rotor angle, and
Transient state parameter of the leading OMIB image systems at each integration moment is obtained using track condensation technique, the transient state parameter includes turning
Sub- angle and acceleration, further piecewise analytic acquire the kinetic energy improved under DEEAC algorithms and accelerate areaSubtract with kinetic energy
Fast areaStability margin is asked for by formula (4):
In above-mentioned formula, ηDE(5)(τ) is the counted stability margin of example application enhancements DEEAC algorithms.
Furthermore, when the example in test example complete or collected works is taken out in step 1), according to example arrangement sequence successively
It takes out.
Compared with the prior art, beneficial effects of the present invention are embodied in:It is calculated the present invention is based on EEAC algorithms, and according to research
Example time-varying degree and to meet positively related causality analysis between integration step number needed for certain integral accuracy, is integrated based on single step
SEEAC algorithms and DEEAC (2) algorithm based on 4 steps integration on the basis of, deeply excavate based on each fast of different integration step numbers
The different information that surely parser is included fast temporarily and its comparison result carry out image study example time-varying degree, are calculated for each test
Example integrates the different classical DEEAC algorithms (DEEAC (2) algorithm) of step number according to the different Auto-matchings of its time-varying degree power, changes
Into DEEAC algorithms (DEEAC (5) algorithm).It is theoretical unanimously to show with testing for any test example, if its time-varying degree compared with
It is weak, on the basis of DEEAC (2) algorithm, integration step number is further increased almost without helping being obviously improved for analysis precision, it is believed that
Quick analysis precision requirement is met at this time;If its time-varying degree is stronger, have on the basis of DEEAC (2) algorithm, by increasing
The integration step number (i.e. DEEAC (5) algorithm) of limit can promote former analysis precision, to meet quick analysis precision requirement.According to this hair
When quick temporary surely analysis is performed to any example, the quick temporary of step number is integrated using based on fixed less for the bright technical solution
Steady parser (SEEAC and DEEAC (2) algorithm) a priori judges whether to further increase limited integration step number with full
The quick analysis precision requirement of foot, so as to using minimum calculating increment as cost, maintain the weaker sample calculation analysis precision of time-varying degree
Under the premise of constant, the analysis precision of the stronger example of time-varying degree is improved.So as to be conducive to further be promoted based on EEAC algorithms
Example screen out and taxonomy model performance, to further coordinate online temporarily steady analysis precision and rate, solution definitely and it is uncertain because
Temporary steady problem analysis under element has great theory and engineering significance.
Description of the drawings
Fig. 1 is the flow chart of the method for the present invention.
Specific embodiment
The present invention is described in detail with specific embodiment below in conjunction with the accompanying drawings.
Embodiment 1:
One embodiment of the present of invention discloses a kind of fast transient stability analysis side based on Adaptive Integral step number
Method, basic principle are as follows:
Increase integration step number is more to count and study example time-varying factor and then improve the effective way of analysis precision,
It is also DEEAC (2) algorithm performance, further plays the only way of its bridge effect.If on the basis of DEEAC (2) algorithm
On blindly increase integration step number, then probably due to lacking the support of stable mechanism and leading to the waste of part calculation amount;Just
When method be the regularity of distribution based on time-varying factor in whole transient processes, the specific temporarily steady stage is targetedly increased
Add appropriate integration step number.
Further research shows that, no matter for three-phase symmetrical or unbalanced fault, after failure time-varying degree be better than therefore
In barrier.Therefore further to play the bridge effect of DEEAC algorithms, should put forth effort on to the time-varying degree stronger stage (stage after failure)
Thus the temporarily careful extraction of steady information promotes the analysis precision of the stronger example of time-varying degree.
A large amount of emulation statistics show that, for the stronger example of time-varying degree, analysis precision can be with integration step after failure
Several increases up to " saturation " state;For the weaker example of time-varying degree, analysis precision is little with integration step after failure
Several increase and become.And for any example, perform calculatings cost needed for analysis then with the increase of integration step number after failure and
Lasting increase.It on the basis of DEEAC (2) algorithm, keeps integration step number in failure constant for 2, step number will be integrated after failure by 2
When increasing to 5, the required precision of most of stronger examples of time-varying degree can have been met, and step number is integrated after further increasing failure
For precision improvement effect then " close to saturation ", and can cause calculate cost increase.It is by DEEAC (2) algorithm improvement
During DEEAC (5) algorithm, the analysis precision of the stronger example of time-varying degree can be met significantly.
The studies above on the basis of magnanimity emulation is performed, return by angle, fusion cause and effect and statistical analysis from posterior analysis
Receive to obtain the related conclusions for further playing DEEAC algorithm bridge effects.If however, can be a priori different each of time-varying degree
Each quick temporarily steady algorithm of the example matching based on appropriate integration step number, to further improving, example screens out and taxonomy model performance is anticipated
Justice is great.
The present embodiment be according to research example time-varying degree and between the integration step number needed for meeting certain integral accuracy just
Relevant causality analysis designs the quick temporarily steady analysis method based on Adaptive Integral step number, to contain time-varying factor degree of strength
Each quick temporarily steady algorithm of the different each example matchings based on different integration step numbers, to give full play to the bridge of DEEAC algorithms
Effect.Specifically, the step of the present embodiment, is as shown in Figure 1:
After step 1 describes the startup of transient stability quick analysis system in Fig. 1, certain in test example complete or collected works is taken out successively
A example carries out the example using SEEAC algorithms temporarily steady nargin and calculates, obtains stability margin ηSE(τ)。
Particularly, it is assumed that two groups of inside of example complementation are the preferable people having the same aspiration and interest, are dominated based on model condensation technique
Infinite bus power system (OMIB) image system is asked for dominating with the fault clearance time (τ) for step-length, based on single step Taylor series expansion
Rotor angle and acceleration of the OMIB image systems at the τ moment, the kinetic energy that further parsing is acquired under SEEAC algorithms accelerate areaWith kinetic energy retardation areaStability margin is asked for by formula (1):
Step 2 is described using classics DEEAC algorithms (referred to as in Fig. 1:DEEAC (2) algorithm) it is abundant to example progress
Degree calculates, and obtains stability margin ηDE(2)(τ)。
Particularly, in the stage in failure and after failure, two step Taylor series expansions is based on and ask for each machine rotor angle, and
Transient state parameter (such as rotor angle and acceleration) of the leading OMIB image systems at each integration moment is obtained using track condensation technique,
The kinetic energy that further piecewise analytic is acquired under DEEAC (2) algorithm accelerates areaWith kinetic energy retardation areaStability margin is asked for by formula (2):
Step 3 discloses a kind of computational methods of image study example time-varying degree in Fig. 1:By comparing ηSE(τ) and
ηDE(2)Difference between (τ) reflects the example time-varying degree, by formula (3) Suo Shi:
Quick temporarily steady algorithms (DEEAC (5) calculations being whether to call based on more integration step numbers of step 4 description in Fig. 1
Method) decision rule, if the research example time-varying degree σ acquired by step 31(τ) is less than or equal to threshold epsilon1, then it is assumed that application
DEEAC (2) algorithm carries out nargin and calculates gained ηDE(2)(τ) has met quick analysis precision, can terminate quick point to the example
Flow is analysed, performs step 6, otherwise performs step 5.
ε in the step1For static threshold, it is a large amount of typical examples according to different real systems, headed by reliability
Principle optimization is wanted to acquire, there is robustness for different system, different operating modes and different faults, different system, model and
It is all constant under failure.Threshold epsilon1Value be 0.75.
Step 5 describes application enhancements DEEAC algorithms (referred to as in Fig. 1:DEEAC (5) algorithm) it is abundant to example progress
Degree calculates, and obtains stability margin ηDE(5)(τ), and think that it has met quick analysis precision, terminate the quick analysis stream to the example
Journey performs step 6.
Particularly, in the stage in failure and after failure, two steps is based respectively on and five step Taylor series expansions ask for each machine
Rotor angle, and obtain transient state parameter (such as rotor angle of the leading OMIB image systems at each integration moment using track condensation technique
With acceleration), the kinetic energy that further piecewise analytic is acquired under DEEAC (5) algorithm accelerates areaWith kinetic energy deceleration face
ProductStability margin is asked for by formula (4):
Step 6 describes in Fig. 1, if each example is complete quick analysis process in test example complete or collected works, ties
Beam is quickly analyzed, and is otherwise removed an example and is performed step 1.
As the specific calculating of the present embodiment, with Hainan (2009 annual data), Shandong (2004 and 2012 annual datas, difference
Be denoted as Shandong A and Shandong B), Jiangxi (2011 annual data), (2012 and 2013 annual datas, are denoted as Zhejiang A and Zhejiang respectively in Zhejiang
B), Henan (2011 annual data), Xinjiang (2012 annual data) and southern net (2012 annual data) 9 systems are in original operating mode, modification
Test example of the circuit three-phase permanent short failure as three-phase symmetrical failure under operating mode, tests the event of example by totally 1652
Barrier place is mostly system core node, the random sampling between 0.08~0.50s of fault clearance time.
With circuit asymmetrical three-phase permanent short failure of above-mentioned 9 systems under original operating mode, modification operating mode (including list
Mutually three kinds of ground connection, two phase ground, two-phase phase fault situations) test example as asymmetrical three-phase failure, totally 1620, survey
The position of fault of tentative calculation example and the selection principle of checkout time are same as above.
Using above-mentioned three-phase symmetrical and unbalanced fault as test example complete or collected works, it is abbreviated as SU。
For SUIn each example, with the stability margin η acquired based on IEEAC algorithmsIE(τ) is benchmark, picks out that
A little stability margin η acquired based on DEEAC (2) algorithmDE(2)(τ) same to ηIEThe opposite example of (τ) qualitative conclusions is opposite for three
Title and unbalanced fault, respectively have 103 and 92 examples to be picked, these examples are considered due to containing strong time-varying factor, causing to answer
Qualitative erroneous judgement is generated during temporarily steady analysis quick with DEEAC (2) algorithm performs.These examples can be classified as containing respectively by fault type
The symmetrical and unbalanced fault test example subset of strong time-varying factor, is referred to as the test example subset containing strong time-varying factor, letter
It is denoted as SSTVD.Remaining example is then referred to as the test example subset containing weak time-varying factor, is abbreviated as SWTVD.Obviously, SSTVD∪SWTVD
=SU,
Through the present embodiment method to SUIn each example when carrying out quick temporarily steady analysis, each example is applied based on different integration steps
The statistics situation of several DEEAC algorithm performs analyses is as shown in table 1:
Table 1 is directed to SUIn the statistical result of each example application the method for the present invention when performing quick temporarily steady analysis
As shown in Table 1, do not occurred belonging to SSTVDIn example can only apply DEEAC (2) algorithm performs temporarily surely analysis be
The situation of quick analysis precision can be met (in fact, belonging to SSTVDIn whole examples be both needed to using DEEAC (5) algorithm performs
Temporarily steady analysis can just meet quick analysis precision), and there is sub-fraction to belong to SWTVDExample be mistaken for still needing to using DEEAC
(5) algorithm performs, which are temporarily surely analyzed, can just meet quick analysis precision.The former ensures example complete or collected works' analysis precision, and the latter then slightly increases
Calculating cost is added.
Further statistical analysis is it is found that using the present embodiment method to SUIn each example carry out it is required during quick temporarily steady analysis
Calculation amount, compare only apply DEEAC (2) algorithm performs when averagely increase about 3.98% calculating cost, calculate cost
Incrementss are smaller.
It has addressed above, for SWTVDIn each example, it is quickly temporary based on DEEAC (2) algorithms and DEEAC (5) algorithm performs
During steady analysis, it will not almost cause the change of analysis precision;When performing quick temporarily steady analysis via the present embodiment method, compare only
Using DEEAC (2) algorithm, it is mainly used for improving SSTVDIn each sample calculation analysis precision, promoted effect it is as shown in table 2.
The effectiveness of 2 the method for the present invention of table
Table 2 is from qualitative and quantitative two angle changing rates effectiveness of the DEEAC (2) with DEEAC (5) algorithm.Wherein, N0It represents
SSTVDIn, the stability margin η that is acquired by DEEAC (2) (DEEAC (5)) algorithmDE(2)(τ)(ηDE(5)(τ)) it is asked with by IEEAC algorithms
The stability margin η obtainedIE(τ) is qualitative to sentence the opposite example number of steady conclusion;In SSTVDIn, for ηDE(2)(τ)(ηDE(5)(τ)) same to ηIE
(τ) is qualitative to sentence those identical examples of steady conclusion, and each example η is characterized with DDE(2)(τ)(ηDE(5)(τ)) value is the same as corresponding ηIE(τ) value
Between whole difference, shown in definition such as formula (5):
In formula (5) N be in symmetrical (asymmetry) the fault test example subset containing strong time-varying factor, the η of each exampleDE(2)
(τ)(ηDE(5)(τ)) with corresponding ηIE(τ) is qualitative to sentence the identical example number of steady conclusion, characterized respectively using D DEEAC (2),
During DEEAC (5) algorithm effectiveness, k values are respectively 2,5 in formula (5).
As shown in Table 2, it for the three-phase symmetrical failure containing strong time-varying factor, in terms of qualitative analysis, compares and only applies
The quick temporarily steady analysis of DEEAC (2) algorithm performs can at least reduce by 50% qualitative False Rate using the present embodiment method;For
The quick temporarily steady analysis gained correct example of qualitative conclusions is performed using DEEAC (2) algorithms and the present embodiment method, is compared only
Using quantization analytical conclusions obtained by DEEAC (2) algorithm, using the present embodiment method conclusion closer to IEEAC algorithms
Quantitative analysis conclusion.For asymmetrical three-phase failure, in terms of Qualitative and quantitative analysis, performed using the present embodiment method quick
Temporarily utility of income is superior to only apply DEEAC (2) algorithm during steady analysis.
To sum up, using the present embodiment method to SUIn each example when performing quick temporarily steady analysis, with minimum calculating increment
For cost, other (time-varying degree is weaker) is being maintained to improve time-varying degree and calculate more by force under the premise of sample calculation analysis precision is constant
The analysis precision of example.DEEAC algorithm bridge effects have further been played, example has thus been further improved and screens out and taxonomy model
Can, for coordinating the transient stability analysis problem of on-line transient stability analysis precision and rate, solution definitely and under uncertain factor
With great theory and engineering significance.
Although the present invention has been described by way of example and in terms of the preferred embodiments, embodiment is not for limiting the present invention's.Not
It is detached from the spirit and scope of the present invention, any equivalence changes done or retouching, also belongs to the protection domain of the present invention.Cause
This protection scope of the present invention should be using the content that claims hereof is defined as standard.
Claims (7)
1. a kind of fast transient Method of Stability Analysis based on Adaptive Integral step number, which is characterized in that include the following steps:
1) after transient stability quick analysis system starts, some example in test example complete or collected works is taken out, using SEEAC algorithms pair
The example carries out nargin calculating;
2) application classics DEEAC algorithms carry out nargin calculating to the example;
3) the example time-varying degree is reflected according to the difference of SEEAC algorithms and classics DEEAC algorithm nargin result of calculations;
4) if threshold epsilon is less than or equal to by the time-varying degree of step 3) reflection1, then it is assumed that application classics DEEAC algorithms carry out abundant
Degree calculates acquired results and has met quick analysis precision, terminates the quick analysis process to the example, performs step 6), otherwise holds
Row next step;
5) application enhancements DEEAC algorithms carry out nargin calculating to the example, and think that acquired results have met quick analysis precision,
Terminate the quick analysis process to the example, perform step 6);
6) if testing each example in example complete or collected works is complete quick analysis, terminates quick analysis process, otherwise remove one
Example performs step 1).
2. the fast transient Method of Stability Analysis according to claim 1 based on Adaptive Integral step number, which is characterized in that
Application SEEAC algorithms are to the method for example progress nargin calculating in the step 1):
Assuming that two groups of inside of example complementation are the preferable people having the same aspiration and interest, leading Infinite bus power system image is obtained based on model condensation technique
System, referred to as leading OMIB image systems, using fault clearance time τ as step-length, asks for leading based on single step Taylor series expansion
Rotor angle and acceleration of the OMIB image systems at the τ moment are led, the kinetic energy that further parsing is acquired under SEEAC algorithms accelerates areaWith kinetic energy retardation areaStability margin is asked for by formula (1):
In above-mentioned formula, ηSE(τ) is the counted stability margin of example application SEEAC algorithms.
3. the fast transient Method of Stability Analysis according to claim 1 based on Adaptive Integral step number, which is characterized in that
Application classics DEEAC algorithms are to the method for example progress nargin calculating in the step 2):
Stage in failure and after failure is based on two step Taylor series expansions and asks for each machine rotor angle, and agglomerated using track
Technology obtains transient state parameter of the leading OMIB image systems at each integration moment, and the transient state parameter includes rotor angle with accelerating
Degree, the kinetic energy that further piecewise analytic is acquired under classical DEEAC algorithms accelerate areaWith kinetic energy retardation areaStability margin is asked for by formula (2):
In above-mentioned formula, ηDE(2)(τ) applies the counted stability margin of classics DEEAC algorithms for the example.
4. the fast transient Method of Stability Analysis according to claim 1 based on Adaptive Integral step number, which is characterized in that
By comparing the difference of SEEAC algorithms and classics DEEAC algorithm nargin result of calculations come when reflecting the example in the step 3)
The method of change degree into:
The difference value σ of SEEAC algorithms and classics DEEAC algorithm nargin result of calculations is asked for by formula (3)1(τ) reflects the example
Time-varying degree:
5. the fast transient Method of Stability Analysis according to claim 1 based on Adaptive Integral step number, which is characterized in that
In the step 4), threshold epsilon1Value be:ε1=0.75.
6. a kind of fast transient Method of Stability Analysis based on Adaptive Integral step number according to claim 1, feature
It is, the method that application enhancements DEEAC algorithms carry out nargin calculating to the example in the step 5) is:
It in the stage in failure and after failure, is based respectively on two steps and five step Taylor series expansions asks for each machine rotor angle, and utilize
Track condensation technique obtains transient state parameter of the leading OMIB image systems at each integration moment, and the transient state parameter includes rotor angle
With acceleration, further piecewise analytic acquires the kinetic energy improved under DEEAC algorithms and accelerates areaWith kinetic energy deceleration face
ProductStability margin is asked for by formula (4):
In above-mentioned formula, ηDE(5)(τ) is the counted stability margin of example application enhancements DEEAC algorithms.
It is 7. special according to any fast transient Method of Stability Analysis based on Adaptive Integral step number of claim 1~7
Sign is:When the example in test example complete or collected works is taken out in step 1), the sequence arranged according to example is taken out successively.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810055727.8A CN108269017A (en) | 2018-01-19 | 2018-01-19 | A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810055727.8A CN108269017A (en) | 2018-01-19 | 2018-01-19 | A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number |
Publications (1)
Publication Number | Publication Date |
---|---|
CN108269017A true CN108269017A (en) | 2018-07-10 |
Family
ID=62776199
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810055727.8A Pending CN108269017A (en) | 2018-01-19 | 2018-01-19 | A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108269017A (en) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492286A (en) * | 2018-10-30 | 2019-03-19 | 南瑞集团有限公司 | Numerical integration based on disturbed track dynamic characteristic shifts to an earlier date terminating method |
CN111832180A (en) * | 2020-07-17 | 2020-10-27 | 南瑞集团有限公司 | Method and system for terminating numerical integration in advance in transient power angle stability analysis |
CN112287284A (en) * | 2020-10-28 | 2021-01-29 | 山东电力研究院 | Transient stability fault screening method and system considering N-m fault time interval |
CN112614015A (en) * | 2020-12-22 | 2021-04-06 | 南瑞集团有限公司 | Method and system for judging transient stability of uncertain class arithmetic example |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336994A (en) * | 2013-04-16 | 2013-10-02 | 南京南瑞集团公司 | Rapid screening method of expected failure sets in electric power system transient stability assessment |
CN104252685A (en) * | 2014-06-17 | 2014-12-31 | 国网电力科学研究院 | Quick robust classification method of preconceived faults for power system transient stability assessment |
CN104331837A (en) * | 2014-08-13 | 2015-02-04 | 国网电力科学研究院 | Simplification method for optimal cutting machine control strategy search of electric system transient stability |
CN105160586A (en) * | 2015-02-17 | 2015-12-16 | 国网电力科学研究院 | Quick assessment method for changeableness of instability mode of electric power system |
-
2018
- 2018-01-19 CN CN201810055727.8A patent/CN108269017A/en active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103336994A (en) * | 2013-04-16 | 2013-10-02 | 南京南瑞集团公司 | Rapid screening method of expected failure sets in electric power system transient stability assessment |
CN104252685A (en) * | 2014-06-17 | 2014-12-31 | 国网电力科学研究院 | Quick robust classification method of preconceived faults for power system transient stability assessment |
CN104331837A (en) * | 2014-08-13 | 2015-02-04 | 国网电力科学研究院 | Simplification method for optimal cutting machine control strategy search of electric system transient stability |
CN105160586A (en) * | 2015-02-17 | 2015-12-16 | 国网电力科学研究院 | Quick assessment method for changeableness of instability mode of electric power system |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109492286A (en) * | 2018-10-30 | 2019-03-19 | 南瑞集团有限公司 | Numerical integration based on disturbed track dynamic characteristic shifts to an earlier date terminating method |
CN109492286B (en) * | 2018-10-30 | 2022-07-01 | 南瑞集团有限公司 | Numerical integration early termination method based on disturbed trajectory dynamic characteristics |
CN111832180A (en) * | 2020-07-17 | 2020-10-27 | 南瑞集团有限公司 | Method and system for terminating numerical integration in advance in transient power angle stability analysis |
CN112287284A (en) * | 2020-10-28 | 2021-01-29 | 山东电力研究院 | Transient stability fault screening method and system considering N-m fault time interval |
CN112614015A (en) * | 2020-12-22 | 2021-04-06 | 南瑞集团有限公司 | Method and system for judging transient stability of uncertain class arithmetic example |
CN112614015B (en) * | 2020-12-22 | 2022-07-08 | 南瑞集团有限公司 | Method and system for judging transient stability of uncertain analogy example |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108269017A (en) | A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number | |
Fioravanti et al. | A study on fault-proneness detection of object-oriented systems | |
CN104331837B (en) | The optimal method for simplifying for cutting the search of machine control strategy of electric power system transient stability | |
Reger et al. | A pattern-based approach to parametric specification mining | |
CN105868116A (en) | Semantic mutation operator based test case generation and optimization method | |
Padmanaban et al. | Efficient identification of (critical) testable path delay faults using decision diagrams | |
CN112866292B (en) | Attack behavior prediction method and device for multi-sample combination attack | |
CN116401532B (en) | Method and system for recognizing frequency instability of power system after disturbance | |
CN111782532B (en) | Software fault positioning method and system based on network abnormal node analysis | |
CN112560269A (en) | Rhapbody state machine-based high fault tolerance electronic system task reliability simulation analysis method | |
CN106021116A (en) | Inaccessible function call path detection method in complex system | |
CN112882899B (en) | Log abnormality detection method and device | |
CN104252685B (en) | The quick strong sorting technique of transient stability evaluation in power system forecast failure | |
Sayre et al. | Partition testing with usage models | |
CN115185920B (en) | Method, device and equipment for detecting log type | |
CN105987823A (en) | Acceleration test method for communication equipment based on vibration stress transformation | |
CN113886832A (en) | Intelligent contract vulnerability detection method, system, computer equipment and storage medium | |
CN111091194B (en) | Operation system identification method based on CAVWBB _ KL algorithm | |
CN105989095B (en) | Take the correlation rule significance test method and device of data uncertainty into account | |
CN107644145A (en) | A kind of failure behavior emulation mode based on Monte Carlo and decision logic | |
CN115935795B (en) | Power system transient voltage stability evaluation model construction method and evaluation method | |
CN107203469A (en) | Complier test accelerated method based on machine learning | |
CN115239971A (en) | GIS partial discharge type recognition model training method, recognition method and system | |
CN117112847B (en) | Data generation method and device of graph database based on community model | |
Zhang et al. | A Software Multi-Fault Clustering Ensemble Technology |
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 | ||
RJ01 | Rejection of invention patent application after publication | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20180710 |