CN112614015A - Method and system for judging transient stability of uncertain class arithmetic example - Google Patents

Method and system for judging transient stability of uncertain class arithmetic example Download PDF

Info

Publication number
CN112614015A
CN112614015A CN202011526534.XA CN202011526534A CN112614015A CN 112614015 A CN112614015 A CN 112614015A CN 202011526534 A CN202011526534 A CN 202011526534A CN 112614015 A CN112614015 A CN 112614015A
Authority
CN
China
Prior art keywords
track
disturbed
uncertain
fep
dsp
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Granted
Application number
CN202011526534.XA
Other languages
Chinese (zh)
Other versions
CN112614015B (en
Inventor
黄天罡
薛禹胜
彭慧敏
刘庆龙
宋晓芳
赖业宁
李威
吕睿
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Hubei Electric Power Co Ltd
NARI Group Corp
Nari Technology Co Ltd
State Grid Electric Power Research Institute
Original Assignee
State Grid Corp of China SGCC
State Grid Hubei Electric Power Co Ltd
NARI Group Corp
Nari Technology Co Ltd
State Grid Electric Power Research Institute
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Hubei Electric Power Co Ltd, NARI Group Corp, Nari Technology Co Ltd, State Grid Electric Power Research Institute filed Critical State Grid Corp of China SGCC
Priority to CN202011526534.XA priority Critical patent/CN112614015B/en
Publication of CN112614015A publication Critical patent/CN112614015A/en
Application granted granted Critical
Publication of CN112614015B publication Critical patent/CN112614015B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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
    • 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
    • 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
    • G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
    • Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation

Landscapes

  • Business, Economics & Management (AREA)
  • Human Resources & Organizations (AREA)
  • Engineering & Computer Science (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Theoretical Computer Science (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Educational Administration (AREA)
  • Marketing (AREA)
  • Development Economics (AREA)
  • Health & Medical Sciences (AREA)
  • Tourism & Hospitality (AREA)
  • Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • General Physics & Mathematics (AREA)
  • Public Health (AREA)
  • Primary Health Care (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Game Theory and Decision Science (AREA)
  • Operations Research (AREA)
  • Quality & Reliability (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a method and a system for judging the transient stability of an uncertain algorithm, which are used for implementing numerical integration on the uncertain algorithm to obtain a disturbed track, and judging the transient stability by mining track quantization information or track form information of the disturbed track to realize the rapid judgment of the transient stability of the uncertain algorithm, and can reduce the calculated amount by at least 60 percent compared with the traditional stability judgment mode based on the detailed model small-step numerical integration.

Description

Method and system for judging transient stability of uncertain class arithmetic example
Technical Field
The invention relates to a method and a system for judging transient stability of an uncertain classification example, belonging to the technical field of power systems and automation thereof.
Background
The patent ZL201310132812.7 discloses a method for rapidly screening an expected fault set for transient stability evaluation of an electric power system (patent number: ZL201310132812.7), which is based on an extended equal-area criterion (EEAC), uses transient stability analysis algorithms with different simplification degrees, namely Static EEAC (SEAAC) and dynamic EEAC (DEEAC), reflects transient stability margins of research examples and further obtained time-varying degrees of the examples as indexes, combines fault information of the research examples, matches and combines different screening criteria, hierarchically screens expected fault subsets meeting the screening criteria from an expected fault complete set, greatly reduces the number of expected faults needing detailed transient stability analysis, and improves the rapidity of transient stability evaluation on the premise of ensuring the accuracy of the transient stability evaluation, thereby coordinating the precision and the benefit speed of the transient stability analysis properly.
Furthermore, a patent of 'a rapid and robust classification method for expected faults in power system transient stability assessment' (patent number ZL201410271454.2) optimizes and enriches identification rules, a complete set of the examples is divided into five types, namely, stable, quasi-stable, critical, quasi-instability and instability, and the screening efficiency is improved due to the fact that identification of the instability examples is supplemented.
The arithmetic cases identified by the existing rapid screening method (such as the two patents) are uncertain arithmetic cases, which are called uncertain arithmetic cases, and for the uncertain arithmetic cases which are not screened out, the stable characteristics of the uncertain arithmetic cases still need to be obtained based on detailed numerical integration, and still large calculation cost is consumed, and rapid judgment cannot be carried out.
Disclosure of Invention
The invention provides a method and a system for judging transient stability of an uncertain algorithm, which solve the problem that the transient stability of the uncertain algorithm cannot be quickly judged.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a method for determining transient stability of uncertain arithmetic examples includes,
performing numerical integration on the uncertain arithmetic embodiment to obtain a disturbed track taking the fault occurrence moment as a starting point;
if the disturbed track meets the preset track quantization information mining requirement, mining track quantization information according to the disturbed track;
if the track quantization information is credible, directly judging the transient stability of the uncertain analogy example;
if the track quantization information is not credible or the disturbed track meets the preset track form information mining requirement, mining track form information according to the disturbed track;
and if the track form information is mature at the current moment, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain class examples based on the overall process disturbed track.
The track quantization information mining requirements are that the integral step number of the disturbed track is greater than or equal to a step number threshold value, and the integral time is greater than or equal to a time threshold value; the track form information mining requirements are that the integral step number of the disturbed track is greater than or equal to a step number threshold value, and the integral time is smaller than a time threshold value.
If the track quantization information is credible, the transient stability of the uncertain classification example is directly judged, and the specific process is,
according to the track quantization information, obtaining a track swing quantization margin index, and judging whether the corresponding swing FEP/DSP is ill-conditioned;
and judging the transient stability of the uncertain calculation examples according to the track swing quantization margin index and the corresponding swing FEP/DSP ill-condition judgment result.
The process of judging whether the swing FEP/DSP is ill-conditioned is as follows,
calculating the ratio of the kinetic energy in the group at the pendulum FEP/DSP to the kinetic energy between the groups;
if the ratio is smaller than the ratio threshold, the FEP/DSP is in a non-sick state;
if the ratio is larger than or equal to the ratio threshold value, the FEP/DSP ill condition is determined.
Calculating the ratio of the kinetic energy in the group at the pendulum FEP/DSP to the kinetic energy between the groups, wherein the formula is as follows,
Figure BDA0002850767660000031
wherein σFEP(DSP)The ratio of the intra-cluster kinetic energy to the inter-cluster kinetic energy at the pendulum FEP/DSP, MiInertia of unit i, omegai(tFEP(DSP)) To reach the swing FEP/DSP time tFEP(DSP)Angular speed of rotor of processing unit i, MeqIs equivalent unit inertia, omegaeq(tFEP(DSP)) To reach the swing FEP/DSP time tFEP(DSP)And (4) processing the rotor angular speed of the equivalent unit.
If the current moment of the track form information is judged to be mature according to the difference of the characteristic values at the adjacent time sections and the criterion, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain calculation examples based on the whole process disturbed track; the characteristic value is a characteristic value representing the track form information.
The characteristic values include, for example,
rotor angle envelope distance:
δminmax(ti)＝δmax(ti)-δmin(ti)
wherein, deltamax(ti) Is a time section tiAt maximum value of rotor angle, delta, of each machinemin(ti) Is a time section tiAt minimum value of rotor angle of each machine, deltaminmax(ti) Is a time section tiThe rotor angle envelope distance;
sub-maximum gap to maximum gap ratio:
Figure BDA0002850767660000032
wherein, deltafir-max(ti) Is a time section tiMaximum value of angular clearance of rotor, deltasec-max(ti) Is a time section tiAt the sub-maximum value of the rotor angular gap, UMV (t)i) Is a time section tiThe ratio of the next maximum gap to the maximum gap;
first s large gap mean square error:
Figure BDA0002850767660000041
wherein, deltaj(ti) Is a time section tiAt the jth gap value of rotor angle, j ∈ [1, s ]]，
Figure BDA0002850767660000042
Is deltaj(ti) Mean value of Ss(ti) The mean square error of the previous s large gap values.
The criteria may include,
criterion one is as follows: at 3 continuous adjacent time sections after a short time of fault clearing, the rotor angle enveloping distances in the characteristic values tend to be consistent;
criterion two: at 2 continuous adjacent time sections after a short time of fault clearing, the ratio of the sub-maximum clearance to the maximum clearance in the characteristic values tends to be consistent;
criterion three: at 5 continuous adjacent time sections after a short time of fault clearing, the mean square deviations of the first s large gap values in the characteristic values tend to be consistent
And if any criterion is met, judging that the track form information is mature at the current moment.
Based on the whole process disturbed track, the transient stability of the uncertain analogy example is judged, and the specific process is as follows:
if disturbed trajectory convergence obtained by implementing numerical integration, overall process disturbed trajectory convergence, maximum value of single machine kinetic energy in the current time period is smaller than a kinetic energy threshold value, and quantitative stability margin obtained by the overall process disturbed trajectory is larger than or equal to a margin threshold value, determining that the transient stability of the uncertain class example is stable; the current time interval is a plurality of integration step lengths traced back from the current time section;
and if the disturbed trajectory acquired by implementing numerical integration is unstable, the disturbed trajectory in the whole process is unstable, the maximum value of the single-machine kinetic energy in the current time period is greater than the kinetic energy threshold value, and the quantitative stability margin acquired by the disturbed trajectory in the whole process is smaller than the margin threshold value, determining the transient state instability of the uncertain similar algorithm.
A transient stability determination system for uncertain arithmetic examples comprises,
a disturbed track module: performing numerical integration on the uncertain arithmetic embodiment to obtain a disturbed track taking the fault occurrence moment as a starting point;
the track quantization information mining module: if the disturbed track meets the preset track quantization information mining requirement, mining track quantization information according to the disturbed track;
a first determination module: if the track quantization information is credible, directly judging the transient stability of the uncertain analogy example;
the track form information mining module: if the track quantization information is not credible or the disturbed track meets the preset track form information mining requirement, mining track form information according to the disturbed track;
a second determination module: and if the track form information is mature at the current moment, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain class examples based on the overall process disturbed track.
The invention achieves the following beneficial effects: the method implements numerical integration on the uncertain arithmetic embodiments to obtain disturbed tracks, and carries out transient stability judgment by mining track quantization information or track form information of the disturbed tracks, thereby realizing the rapid judgment of the transient stability of the uncertain arithmetic embodiments, and reducing the calculated amount by at least 60 percent compared with the traditional stability judgment mode based on the detailed model small-step numerical integration.
Drawings
FIG. 1 is a flow chart of the method of the present invention.
Detailed Description
The invention is further described below with reference to the accompanying drawings. The following examples are only for illustrating the technical solutions of the present invention more clearly, and the protection scope of the present invention is not limited thereby.
As shown in fig. 1, a method for determining transient stability of an uncertain classification example includes the following steps:
step 1, implementing proper numerical integration on the uncertain calculation example based on a detailed model and a small integration step length to obtain a disturbed track taking the fault occurrence moment as a starting point.
Step 2, if the disturbed track meets the preset track quantization information mining requirement, mining track quantization information according to the disturbed track, and turning to step 3; and if the disturbed track meets the preset track form information mining requirement, turning to the step 4.
Track quantization information mining requirements: the integral step number of the disturbed track is greater than or equal to the step number threshold value, and the integral time is greater than or equal to the time threshold value; track form information mining requirements: the integral step number of the disturbed track is more than or equal to the step number threshold value, and the integral time is less than the time threshold value; wherein, the threshold value of the number of steps is 5 steps, and the time threshold value is 0.5 second.
Step 3, judging whether the track quantization information is credible, and if so, directly judging the transient stability of the uncertain algorithms; if not, go to step 4.
The transient stability of the uncertain classification example is directly judged, and the specific process comprises the following steps:
31) according to the track quantization information, obtaining a track swing quantization margin index, and judging whether a corresponding swing FEP/DSP (farthest point/dynamic saddle point) is ill-conditioned;
and (3) judging whether the FEP/DSP is ill or not:
1) calculating the ratio of the kinetic energy in the group at the pendulum FEP/DSP to the kinetic energy between the groups;
Figure BDA0002850767660000061
wherein σFEP(DSP)The ratio of the intra-cluster kinetic energy to the inter-cluster kinetic energy at the pendulum FEP/DSP, MiInertia of unit i, omegai(tFEP(DSP)) To reach the swing FEP/DSP time tFEP(DSP)Angular speed of rotor of processing unit i, MeqIs equivalent unit inertia, omegaeq(tFEP(DSP)) To reach the swing FEP/DSP time tFEP(DSP)Processing the rotor angular speed of the equivalent unit;
2) if the ratio is smaller than the ratio threshold, the FEP/DSP is in a non-sick state;
3) if the ratio is larger than or equal to the ratio threshold value, the FEP/DSP ill condition is determined.
32) And judging the transient stability of the uncertain calculation examples according to the track swing quantization margin index and the corresponding swing FEP/DSP ill-condition judgment result.
And 4, mining track form information according to the disturbed track.
Step 5, if the current moment of the track form information is judged to be mature according to the difference of the characteristic values at the adjacent time sections and the criterion, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain analogy example based on the whole process disturbed track; if the track is not mature, increasing the step number, and repeating the steps, wherein the characteristic value is the characteristic value representing the track form information.
Specific characteristic values include:
1) rotor angle envelope distance:
δminmax(ti)＝δmax(ti)-δmin(ti)
wherein, deltamax(ti) Is a time section tiAt maximum value of rotor angle, delta, of each machinemin(ti) Is a time section tiAt minimum value of rotor angle of each machine, deltaminmax(ti) Is a time section tiThe rotor angle envelope distance;
2) sub-maximum gap to maximum gap ratio:
Figure BDA0002850767660000071
wherein, deltafir-max(ti) Is a time section tiMaximum value of angular clearance of rotor, deltasec-max(ti) Is a time section tiAt the sub-maximum value of the rotor angular gap, UMV (t)i) Is a time section tiThe ratio of the next maximum gap to the maximum gap;
3) first s large gap mean square error:
Figure BDA0002850767660000072
wherein, deltaj(ti) Is a time section tiAt the jth gap value of rotor angle, j ∈ [1, s ]]，
Figure BDA0002850767660000081
Is deltaj(ti) Mean value of Ss(ti) The former s is the mean square error of the large gap value, and the value of s is 5.
The criteria for the maturation of the trajectory morphology information include,
criterion one is as follows: at 3 continuous adjacent time sections after a short time of fault clearing, the rotor angle enveloping distances in the characteristic values tend to be consistent;
the expression is as follows:
Figure BDA0002850767660000082
Figure BDA0002850767660000083
Figure BDA0002850767660000084
tk＞τ+0.05
criterion two: at 2 continuous adjacent time sections after a short time of fault clearing, the ratio of the sub-maximum clearance to the maximum clearance in the characteristic values tends to be consistent;
the expression is as follows:
Figure BDA0002850767660000085
Figure BDA0002850767660000086
tk＞τ+0.05
criterion three: at 5 continuous adjacent time sections after a short time of fault clearing, the mean square deviations of the first s large gap values in the characteristic values tend to be consistent;
the expression is as follows:
Figure BDA0002850767660000091
Figure BDA0002850767660000092
Figure BDA0002850767660000093
Figure BDA0002850767660000094
Figure BDA0002850767660000095
tk＞τ+0.05
wherein, tk、tk-1、tk-2、tk-3、tk-4、tk-5All the time sections are time sections, tau is fault clearing time, tau +0.05 is a small section of time for fault clearing, and if any criterion is met, the current time of the track form information is judged to be mature.
Based on the whole process disturbed track, the transient stability of the uncertain analogy example is judged, and the specific process is as follows:
1) if disturbed trajectory convergence obtained by implementing numerical integration, overall process disturbed trajectory convergence, maximum value of single machine kinetic energy in the current time period is smaller than a kinetic energy threshold value, and quantitative stability margin obtained by the overall process disturbed trajectory is larger than or equal to a margin threshold value, determining that the transient stability of the uncertain class example is stable; the current time interval is a plurality of integral step lengths traced back from the current time section, and 4 step lengths are generally traced back;
2) and if the disturbed trajectory acquired by implementing numerical integration is unstable, the disturbed trajectory in the whole process is unstable, the maximum value of the single-machine kinetic energy (backtracking 4 steps from the current time section forward) in the current time period is greater than the kinetic energy threshold value, and the quantitative stability margin acquired by the disturbed trajectory in the whole process is smaller than the margin threshold value, determining the transient instability of the uncertain analogy example.
The software system corresponding to the method comprises,
a disturbed track module: performing numerical integration on the uncertain arithmetic embodiment to obtain a disturbed track taking the fault occurrence moment as a starting point;
the track quantization information mining module: if the disturbed track meets the preset track quantization information mining requirement, mining track quantization information according to the disturbed track;
a first determination module: if the track quantization information is credible, directly judging the transient stability of the uncertain analogy example;
the track form information mining module: if the track quantization information is not credible or the disturbed track meets the preset track form information mining requirement, mining track form information according to the disturbed track;
a second determination module: and if the track form information is mature at the current moment, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain class examples based on the overall process disturbed track.
The method implements numerical integration on the uncertain algorithm to obtain the disturbed track, carries out transient stability judgment by mining track quantization information or track form information of the disturbed track, flexibly applies numerical integration and Taylor expansion to realize quick judgment of the transient stability of the uncertain algorithm, and can reduce the calculated amount by at least 60 percent compared with the traditional stability judgment mode based on the detailed model small-step numerical integration.
The above method may be a supplement to an existing transient stability determination method based only on taylor expansion, that is, a transient stability determination method, and specifically includes the following steps:
A1) for all the example sets, screening out the examples of which the stability property can be identified only based on Taylor expansion, and remaining uncertain algorithms;
for all the example sets, the method for screening the expected fault set for the transient stability assessment of the power system and the method for rapidly and robustly classifying the expected faults for the transient stability assessment of the power system (patent numbers ZL201310132812.7 and ZL201410271454.2) are used for screening the examples, and the filtering stability property is clear;
A2) and traversing all uncertain algorithms, and judging the transient stability of the uncertain algorithms by adopting the transient stability judging method of the uncertain algorithms.
A computer readable storage medium storing one or more programs, the one or more programs comprising instructions, which when executed by a computing device, cause the computing device to perform an uncertain-class-example transient stability determination method or a transient stability determination method.
A computing device comprising one or more processors, one or more memories, and one or more programs stored in the one or more memories and configured to be executed by the one or more processors, the one or more programs including instructions for performing an uncertain-class-example transient stability determination method or a transient stability determination method.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
The present invention is not limited to the above embodiments, and any modifications, equivalent replacements, improvements, etc. made within the spirit and principle of the present invention are included in the scope of the claims of the present invention which are filed as the application.

Claims (10)

1. A method for judging transient stability of uncertain analogy examples is characterized by comprising the following steps: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
performing numerical integration on the uncertain arithmetic embodiment to obtain a disturbed track taking the fault occurrence moment as a starting point;
if the disturbed track meets the preset track quantization information mining requirement, mining track quantization information according to the disturbed track;
if the track quantization information is credible, judging the transient stability of the uncertain analogy example;
if the track quantization information is not credible or the disturbed track meets the preset track form information mining requirement, mining track form information according to the disturbed track;
and if the track form information is mature at the current moment, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain class examples based on the overall process disturbed track.
2. The method of claim 1, wherein the method comprises: the track quantization information mining requirements are that the integral step number of the disturbed track is greater than or equal to a step number threshold value, and the integral time is greater than or equal to a time threshold value; the track form information mining requirements are that the integral step number of the disturbed track is greater than or equal to a step number threshold value, and the integral time is smaller than a time threshold value.
3. The method of claim 1, wherein the method comprises: if the track quantization information is credible, the transient stability of the uncertain classification example is determined, and the specific process is,
according to the track quantization information, obtaining a track swing quantization margin index, and judging whether the corresponding swing FEP/DSP is ill-conditioned;
and judging the transient stability of the uncertain calculation examples according to the track swing quantization margin index and the corresponding swing FEP/DSP ill-condition judgment result.
4. The method of claim 3, wherein the method comprises: the process of judging whether the swing FEP/DSP is ill-conditioned is as follows,
calculating the ratio of the kinetic energy in the group at the pendulum FEP/DSP to the kinetic energy between the groups;
if the ratio is smaller than the ratio threshold, the FEP/DSP is in a non-sick state;
if the ratio is larger than or equal to the ratio threshold value, the FEP/DSP ill condition is determined.
5. The method of claim 4, wherein the method comprises: calculating the ratio of the kinetic energy in the group at the pendulum FEP/DSP to the kinetic energy between the groups, wherein the formula is as follows,
Figure FDA0002850767650000021
wherein σFEP(DSP)The ratio of the intra-cluster kinetic energy to the inter-cluster kinetic energy at the pendulum FEP/DSP, MiInertia of unit i, omegai(tFEP(DSP)) To reach the swing FEP/DSP time tFEP(DSP)Angular speed of rotor of processing unit i, MeqIs equivalent unit inertia, omegaeq(tFEP(DSP)) To reach the swing FEP/DSP time tFEP(DSP)And (4) processing the rotor angular speed of the equivalent unit.
6. The method of claim 1, wherein the method comprises: if the current moment of the track form information is judged to be mature according to the difference of the characteristic values at the adjacent time sections and the criterion, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain calculation examples based on the whole process disturbed track; the characteristic value is a characteristic value representing the track form information.
7. The method of claim 6, wherein the method comprises: the characteristic values include, for example,
rotor angle envelope distance:
δmin max(ti)＝δmax(ti)-δmin(ti)
wherein, deltamax(ti) Is a time section tiAt maximum value of rotor angle, delta, of each machinemin(ti) Is a time section tiAt minimum value of rotor angle of each machine, deltamin max(ti) Is a time section tiThe rotor angle envelope distance;
sub-maximum gap to maximum gap ratio:
Figure FDA0002850767650000031
wherein, deltafir-max(ti) Is a time section tiMaximum value of angular clearance of rotor, deltasec-max(ti) Is a time section tiAt the sub-maximum value of the rotor angular gap, UMV (t)i) Is a time section tiThe ratio of the next maximum gap to the maximum gap;
first s large gap mean square error:
Figure FDA0002850767650000032
wherein, deltaj(ti) Is a time section tiAt the jth gap value of rotor angle, j ∈ [1, s ]]，
Figure FDA0002850767650000033
Is deltaj(ti) Mean value of Ss(ti) The mean square error of the previous s large gap values.
8. The uncertain-class-example transient stability determination method according to claim 6 or 7, wherein: the criteria may include,
criterion one is as follows: at 3 continuous adjacent time sections after a short time of fault clearing, the rotor angle enveloping distances in the characteristic values tend to be consistent;
criterion two: at 2 continuous adjacent time sections after a short time of fault clearing, the ratio of the sub-maximum clearance to the maximum clearance in the characteristic values tends to be consistent;
criterion three: at 5 continuous adjacent time sections after a short time of fault clearing, the mean square deviations of the first s large gap values in the characteristic values tend to be consistent
And if any criterion is met, judging that the track form information is mature at the current moment.
9. The method of claim 1, wherein the method comprises: based on the whole process disturbed track, the transient stability of the uncertain analogy example is judged, and the specific process is as follows:
if disturbed trajectory convergence obtained by implementing numerical integration, overall process disturbed trajectory convergence, maximum value of single machine kinetic energy in the current time period is smaller than a kinetic energy threshold value, and quantitative stability margin obtained by the overall process disturbed trajectory is larger than or equal to a margin threshold value, determining that the transient stability of the uncertain class example is stable; the current time interval is a plurality of integration step lengths traced back from the current time section;
and if the disturbed trajectory acquired by implementing numerical integration is unstable, the disturbed trajectory in the whole process is unstable, the maximum value of the single-machine kinetic energy in the current time period is greater than the kinetic energy threshold value, and the quantitative stability margin acquired by the disturbed trajectory in the whole process is smaller than the margin threshold value, determining the transient state instability of the uncertain similar algorithm.
10. An uncertain classification example transient stability determination system is characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
a disturbed track module: performing numerical integration on the uncertain arithmetic embodiment to obtain a disturbed track taking the fault occurrence moment as a starting point;
the track quantization information mining module: if the disturbed track meets the preset track quantization information mining requirement, mining track quantization information according to the disturbed track;
a first determination module: if the track quantization information is credible, directly judging the transient stability of the uncertain analogy example;
the track form information mining module: if the track quantization information is not credible or the disturbed track meets the preset track form information mining requirement, mining track form information according to the disturbed track;
a second determination module: and if the track form information is mature at the current moment, acquiring a subsequent disturbed track by adopting Taylor expansion, and judging the transient stability of the uncertain class examples based on the overall process disturbed track.
CN202011526534.XA 2020-12-22 2020-12-22 Method and system for judging transient stability of uncertain analogy example Active CN112614015B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202011526534.XA CN112614015B (en) 2020-12-22 2020-12-22 Method and system for judging transient stability of uncertain analogy example

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202011526534.XA CN112614015B (en) 2020-12-22 2020-12-22 Method and system for judging transient stability of uncertain analogy example

Publications (2)

Publication Number Publication Date
CN112614015A true CN112614015A (en) 2021-04-06
CN112614015B CN112614015B (en) 2022-07-08

Family

ID=75244039

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202011526534.XA Active CN112614015B (en) 2020-12-22 2020-12-22 Method and system for judging transient stability of uncertain analogy example

Country Status (1)

Country Link
CN (1) CN112614015B (en)

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108269017A (en) * 2018-01-19 2018-07-10 国电南瑞科技股份有限公司 A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number
CN109492286A (en) * 2018-10-30 2019-03-19 南瑞集团有限公司 Numerical integration based on disturbed track dynamic characteristic shifts to an earlier date terminating method

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN108269017A (en) * 2018-01-19 2018-07-10 国电南瑞科技股份有限公司 A kind of fast transient Method of Stability Analysis based on Adaptive Integral step number
CN109492286A (en) * 2018-10-30 2019-03-19 南瑞集团有限公司 Numerical integration based on disturbed track dynamic characteristic shifts to an earlier date terminating method

Also Published As

Publication number Publication date
CN112614015B (en) 2022-07-08

Similar Documents

Publication Publication Date Title
US20200257543A1 (en) Aggregate Features For Machine Learning
WO2021151292A1 (en) Corpus monitoring method based on mask language model, corpus monitoring apparatus, device, and medium
CN109933778B (en) Word segmentation method, word segmentation device and computer readable storage medium
US10002142B2 (en) Method and apparatus for generating schema of non-relational database
CN105706092B (en) The method and system of four values simulation
CN108880915B (en) Electric power information network safety alarm information false alarm determination method and system
US20160170961A1 (en) Perceptual associative memory for a neuro-linguistic behavior recognition system
WO2021109874A1 (en) Method for generating topology diagram, anomaly detection method, device, apparatus, and storage medium
CN113626400A (en) Log event extraction method and system based on log tree and analytic tree
CN112614015B (en) Method and system for judging transient stability of uncertain analogy example
CN113890821A (en) Log association method and device and electronic equipment
Zhang et al. MolFeSCue: enhancing molecular property prediction in data-limited and imbalanced contexts using few-shot and contrastive learning
CN112100400A (en) Node recommendation method and device based on knowledge graph
CN104731889B (en) A kind of method for estimating query result size
US20180349322A1 (en) Generating models for real time embedded systems that approximate non-embedded models while reducing complexity associated with the non-embedded models
CN114281950B (en) Data retrieval method and system based on multi-graph weighted fusion
CN113495901B (en) Quick retrieval method for variable-length data blocks
CN109902178A (en) A kind of multistage file classification method and system
CN113485878B (en) Multi-data center fault detection method
CN114969761A (en) Log anomaly detection method based on LDA theme characteristics
Zhang et al. Small files storing and computing optimization in Hadoop parallel rendering
Zhang et al. Log Anomaly Detection Based on Bi-LSTM Feature Extraction
WO2021047576A1 (en) Log record processing method and apparatus, and device and machine-readable storage medium
Shang et al. A Multi-view Clustering Model for Event Detection in Twitter
Yang et al. A Parallel Gene Expression Clustering Algorithm Based on Producer-Consumer Model

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