CN111832180A - Method and system for terminating numerical integration in advance in transient power angle stability analysis - Google Patents

Method and system for terminating numerical integration in advance in transient power angle stability analysis Download PDF

Info

Publication number
CN111832180A
CN111832180A CN202010691505.2A CN202010691505A CN111832180A CN 111832180 A CN111832180 A CN 111832180A CN 202010691505 A CN202010691505 A CN 202010691505A CN 111832180 A CN111832180 A CN 111832180A
Authority
CN
China
Prior art keywords
criterion
generator
quantitative analysis
now
leading
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
CN202010691505.2A
Other languages
Chinese (zh)
Other versions
CN111832180B (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
Original Assignee
State Grid Corp of China SGCC
State Grid Hubei Electric Power Co Ltd
NARI Group Corp
Nari Technology Co Ltd
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 filed Critical State Grid Corp of China SGCC
Priority to CN202010691505.2A priority Critical patent/CN111832180B/en
Publication of CN111832180A publication Critical patent/CN111832180A/en
Application granted granted Critical
Publication of CN111832180B publication Critical patent/CN111832180B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00Details relating to the application field
    • G06F2113/04Power grid distribution networks

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention discloses a method for terminating a transient power angle stability analysis numerical integration in advance, which comprises the steps of judging whether the current simulation moment reaches the system fault clearing moment; responding to the current simulation time reaching the system fault clearing time, and obtaining a leading image according to the maximum rotor angle clearance at the current simulation time; acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set; judging whether the stable state/unstable state of the leading image is mature or not; responding to the maturity of the stable state/unstable state of the dominant image, and carrying out quantitative analysis on the stability/instability of the power angle; judging whether the quantitative analysis conclusion is credible; and responding to the credibility of the quantitative analysis conclusion, outputting the quantitative analysis conclusion, and ending the transient power angle stability analysis. A responsive system is also disclosed. The invention realizes the termination of numerical integration in advance on the premise of ensuring the accuracy of quantitative analysis of the transient power angle stability, and the required calculated amount can be reduced to at least 50% of the original calculated amount.

Description

Method and system for terminating numerical integration in advance in transient power angle stability analysis
Technical Field
The invention relates to a method and a system for early terminating a transient power angle stability analysis numerical integration, belonging to the technical field of power systems and automation thereof.
Background
Energy conversion, while driving a "clean alternative," would significantly increase system operational complexity and uncertainty. The synchronous rising of the capacity of the source load side and the frequent interaction of the internal and external information enable the transmission capacity of the key connecting line to continuously approach the limit, increase the disturbed frequency and increase the spread degree.
And the system stability characteristics caused by disturbance are quickly analyzed, and the premise that risk examples are screened in time and the control strategy is effectively matched is provided. The method aims to filter out credible risk-free or high-risk examples, and obtains wide attention of academic and engineering fields.
Most of the existing researches are carried out by following data-driven and model-driven research models and constructing classifiers to realize example screening. Following the data-driven research paradigm, the feature quantities on individual time sections are taken as input parameters of the classifier: in the off-line training stage, the mapping relation between the input parameters and the actual transient stability conclusion is learned by self; and in the online application stage, a qualitative conclusion of transient stability or instability is directly deduced from the characteristic quantity obtained in the online analysis stage. Following the model-driven research paradigm, the multiple bases on causal relationships perform a stability approximation: and the qualitative analysis of most of the examples is realized by simplifying the element model, refining the algorithm flow and constructing a simplified algorithm. Following a research paradigm of a fusion model and data driving, on the basis of deep mining of causal elements affecting analysis accuracy, the causal elements are introduced into a classifier so as to further improve screening robustness.
On the premise of trying to ensure the detailed analysis of uncertain risk calculation examples, most of the existing screening technologies achieve better effects. However, based on the current example screening technology, even if the stability characteristics of each example are determined, the work of strategy formulation, plan adjustment and the like required for further improving the transient stability characteristics of the system cannot be executed. For the example judged as transient instability, it is not clear which emergency control strategy (such as cutting machine, load cutting and the like) and optimization technology thereof need to be adopted, so that the system recovers transient stability; for the example determined as transient stability, since the margin of the current situation from the critical situation cannot be obtained, how to deal with the stability analysis under the coupling of potential multiple uncertain factors is a difficult problem.
This requires the goal of improving the analysis utility, not only to utilize or merge various analysis-driven ways to accurately classify risk-free or high-risk examples, but also to further obtain transient details such as (nearly accurate) quantitative margins and dominant patterns (the latter being a necessary condition for influencing stability control strategy formulation and even optimization). For this reason, it is necessary to study the early termination method of numerical integration for the purpose of obtaining accurate (or more accurate) dominant mode and quantization margin.
The existing early termination methods (such as an algorithm for automatically terminating numerical integration and engineering application thereof, estimation of subsequent stability of a stable track, a numerical integration early termination method based on dynamic characteristics of a disturbed track and the like) mostly start from the essence that time-varying factors in the essence of a power system influence the stability of a power angle, and start from motion state quantities (such as omega, and the like) and energy (such as E)kEtc.) view angle reveals the degree of system time variation, and further design criteria are supplemented by integral termination. However, the transient power angle stability characteristic is a global characteristic, which is reflected by the phase between the motion state quantities of the generatorsMutual containment; particularly for the study calculation examples with complex stabilization mechanism and remarkable time-varying characteristics, the power angle stabilization situation based on the motion state quantity has a remarkable time-lag effect; in other words, the relative motion situation of each machine rotor angle is a mature manifestation of the current situation of power angle stability, and is not a potential representation of the evolution trend thereof. The time-series evolution view based on the motion state quantity is not favorable for studying and judging the stable situation of the power angle in advance, so that a method for terminating numerical integration in advance on the premise of ensuring the quantitative analysis precision of the transient power angle is urgently needed.
Disclosure of Invention
The invention provides a method and a system for early terminating a transient power angle stability analysis numerical integration, which solve the problems in the background technology.
In order to solve the technical problems, the technical scheme adopted by the invention is as follows:
a method for terminating numerical integration in advance for transient power angle stability analysis includes,
in response to the current simulation time reaching the system fault clearing time, identifying a leading group and a remaining group of the generator set according to the maximum rotor angle gap at the current simulation time to obtain a leading image;
acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set;
judging whether the leading image stable state/unstable state is mature or not according to a preset maturation criterion, the bus voltage amplitude of each generator and the corresponding generator set;
responding to the maturity of the stable state/unstable state of the dominant image, and carrying out quantitative analysis on the stability/instability of the power angle;
judging whether the quantitative analysis conclusion is credible or not according to a preset credible criterion;
and responding to the credibility of the quantitative analysis conclusion, outputting the quantitative analysis conclusion, and ending the transient power angle stability analysis.
The maturation criterion comprises a criterion for leading the maturation of the stable situation of the image and a criterion for leading the maturation of the unstable situation;
the criterion for leading the map to be mature in stable situation comprises the following steps:
criterion 1: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front collar group has the same trend along with the change of time, and the current simulation time t isnowThe highest value of bus voltage of the group generator set before the leader is not less than the lowest value of bus voltage of the rest group generator sets, wherein delta t is a numerical integration time step;
criterion 2: current simulation time tnowThe lowest voltage value of each generator bus is greater than 0.75;
criteria governing maturation of destabilization potentials include:
criterion 3: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generator buses in the front leading group is consistent with the change trend of time, the highest value of the bus voltage of the front leading group generator is not more than the lowest value of the bus voltage of the rest group generator, and the generator connected with the bus corresponding to the lowest value of the bus is the generator set in the front leading group;
criterion 4: current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.20; or the current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.40, and the time zone t isnow-Δt,tnow]And the voltage amplitude of the bus of any leading group generator set is reduced with time, and the gradient is less than-3.00.
In response to the satisfaction of the criteria 1 and 2, judging that the stable situation of the leading image is mature; and in response to the satisfaction of the criteria 3 and 4, judging that the instability situation of the dominant image is mature.
The credible criterion comprises a credible criterion of a quantitative analysis conclusion of power angle stability and a credible criterion of a quantitative analysis conclusion of power angle instability;
the credible criterion of the quantitative analysis conclusion of the power angle stability is as follows:
current simulation time tnowThe minimum value of the bus voltage of each generator is greater than a first threshold value, and the quantitative analysis conclusion is greater than a second threshold value;
credible criterion for quantitative analysis conclusion of power angle instability:
current simulation time tnowThe minimum value of the bus voltage of each generator is less than a third thresholdAnd the quantitative analysis conclusion is less than the fourth threshold.
A system for early terminating a transient power angle stability analysis numerical integration comprises,
a leading image identification module: in response to the current simulation time reaching the system fault clearing time, identifying a leading group and a remaining group of the generator set according to the maximum rotor angle gap at the current simulation time to obtain a leading image;
a voltage information acquisition module: acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set;
a stable situation/unstable situation maturity judging module: judging whether the leading image stable state/unstable state is mature or not according to a preset maturation criterion, the bus voltage amplitude of each generator and the corresponding generator set;
a quantitative analysis module: responding to the maturity of the stable state/unstable state of the dominant image, and carrying out quantitative analysis on the stability/instability of the power angle;
a quantitative analysis conclusion judgment module: judging whether the quantitative analysis conclusion is credible or not according to a preset credible criterion;
a result module: and responding to the credibility of the quantitative analysis conclusion, outputting the quantitative analysis conclusion, and ending the transient power angle stability analysis.
The maturation criterion comprises a criterion for leading the maturation of the stable situation of the image and a criterion for leading the maturation of the unstable situation;
the criterion for leading the map to be mature in stable situation comprises the following steps:
criterion 1: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front collar group has the same trend along with the change of time, and the current simulation time t isnowThe highest value of bus voltage of the group generator set before the leader is not less than the lowest value of bus voltage of the rest group generator sets, wherein delta t is a numerical integration time step;
criterion 2: current simulation time tnowThe lowest voltage value of each generator bus is greater than 0.75;
criteria governing maturation of destabilization potentials include:
criterion 3: in the time zoneSegment [ t ]now-Δt,tnow]The voltage amplitude of any two generator buses in the front leading group is consistent with the change trend of time, the highest value of the bus voltage of the front leading group generator is not more than the lowest value of the bus voltage of the rest group generator, and the generator connected with the bus corresponding to the lowest value of the bus is the generator set in the front leading group;
criterion 4: current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.20; or the current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.40, and the time zone t isnow-Δt,tnow]And the voltage amplitude of the bus of any leading group generator set is reduced with time, and the gradient is less than-3.00.
The stable situation/unstable situation maturity judging module comprises a dominant image stable situation maturity judging module and a dominant image unstable situation maturity judging module; a dominant mapping stable situation maturity judging module: in response to the satisfaction of the criteria 1 and 2, judging that the stable situation of the leading image is mature; a dominant image instability situation maturity judgment module: and in response to the satisfaction of the criteria 3 and 4, judging that the instability situation of the dominant image is mature.
The credible criterion comprises a credible criterion of a quantitative analysis conclusion of power angle stability and a credible criterion of a quantitative analysis conclusion of power angle instability;
the credible criterion of the quantitative analysis conclusion of the power angle stability is as follows:
current simulation time tnowThe minimum value of the bus voltage of each generator is greater than a first threshold value, and the quantitative analysis conclusion is greater than a second threshold value;
credible criterion for quantitative analysis conclusion of power angle instability:
current simulation time tnowAnd the minimum value of the bus voltage of each generator is smaller than a third threshold, and the quantitative analysis conclusion is smaller than a fourth threshold.
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 a transient power angle stability analysis numerical integration early termination method.
A computing device comprising one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing a transient power angle stability analysis numerical integration early termination method.
The invention achieves the following beneficial effects: the method and the device realize that numerical integration is terminated in advance on the premise of ensuring the accuracy of quantitative analysis of the transient power angle stability, the required calculated amount can be reduced to at least 50% of the original calculated amount, a foundation is laid for realizing rapid analysis of the transient power angle stability under the influence of multiple coupling uncertain factors, and certain theoretical and engineering significance is achieved.
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.
Transient voltage stability mostly presents local characteristics, the bus voltage of each generator is jointly determined by the reactive power support obtained by the bus voltage and the reactive power demand required to be consumed, and weak interaction influence situation is presented among all generator groups. Therefore, based on the incidence relation between the transient power angle stable situation and the transient voltage stable situation, the potential motion situation information of the generator rotor angle can be revealed in advance based on the time sequence evolution information of the generator terminal bus voltage. If the time sequence evolution information of the generator terminal bus voltage keeps consistent with the current relative motion situation information of the rotor angle of each generator, the current relative motion situation of the rotor angle of each generator is proved to have evolved and mature (the subsequent relative motion situation of each generator is approximately kept unchanged), thereby guiding the early termination of numerical integration. The specific process is shown in fig. 1, and a method for terminating a transient power angle stability analysis numerical integration early comprises the following steps:
step 1, judging whether the current simulation time reaches the system fault clearing time.
In the transient power angle stability analysis, an element model (such as a generator classical model or a detailed model, an induction motor load model or a constant impedance load model or a comprehensive load model and the like) and calculation parameters (such as an analysis step length and the like) required by simulation are set, and a numerical integration method is adopted to simulate the system.
Step 2, if the current simulation time reaches the system fault clearing time, identifying a leading group and a remaining group of the generator set according to the maximum rotor angle gap of the current simulation time to obtain a leading image, and turning to step 3; otherwise, simulating one step and turning to the step 1.
And 3, acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set.
And 4, judging whether the dominant image stable state/unstable state is mature or not according to a preset maturation criterion, the bus voltage amplitude of each generator and the corresponding generator set.
The maturation criterion comprises a criterion for leading the maturation of the stable situation of the image and a criterion for leading the maturation of the unstable situation.
The criterion for leading the map to be mature in stable situation comprises the following steps:
criterion 1: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front collar group has the same trend along with the change of time, and the current simulation time t isnowThe highest value of bus voltage of the group generator set before the lead is not less than the lowest value of bus voltage of the rest group generator sets, wherein delta t is a numerical integration time step from the current simulation time tnowThe time section formed by tracing back an integral time step is tnow-Δt,tnow];
Criterion 2: current simulation time tnowThe lowest voltage value of each generator bus is greater than 0.75 (per unit value);
criteria governing maturation of destabilization potentials include:
criterion 3: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front-leading group has the same trend along with the change of time, and the maximum value of the bus voltage of the front-leading group generator is not more than that of the rest groups for power generationThe generator set connected with the bus with the lowest voltage value is a generator set in the front group;
criterion 4: current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.20 (per unit value); or the current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.40 (per unit value) and is in a time section [ tnow-Δt,tnow]And the voltage amplitude of the bus of any leading group generator set is reduced with time, and the gradient is less than-3.00.
In response to the satisfaction of the criteria 1 and 2, judging that the stable situation of the leading image is mature; and in response to the satisfaction of the criteria 3 and 4, judging that the instability situation of the dominant image is mature.
Step 5, if the dominant map stable state/instability state is mature, carrying out quantitative analysis on power angle stability/instability, and turning to step 6; otherwise, simulating one step and turning to the step 1.
And 6, judging whether the quantitative analysis conclusion is credible or not according to a preset credibility criterion.
The credible criterion comprises a credible criterion for quantitative analysis conclusion of power angle stability and a credible criterion for quantitative analysis conclusion of power angle instability.
The credible criterion of the quantitative analysis conclusion of the power angle stability is as follows:
current simulation time tnowThe minimum value of the bus voltage of each generator is larger than a first threshold value1And the quantitative analysis conclusion is greater than the second threshold value2
Credible criterion for quantitative analysis conclusion of power angle instability:
current simulation time tnowThe minimum value of the bus voltage of each generator is less than a third threshold value3And the quantitative analysis conclusion is less than the fourth threshold value4
Step 7, outputting a quantitative analysis conclusion if the quantitative analysis conclusion (the quantitative analysis conclusion of the power angle stability is credible and the quantitative analysis conclusion of the power angle instability is credible) is credible, and ending the transient power angle stability analysis; otherwise, simulating one step and turning to the step 1.
The method realizes that numerical integration is terminated in advance (reflected in step 7) on the premise of ensuring the accuracy of quantitative analysis of the transient power angle stability, the required calculated amount can be reduced to at least 50% of the original calculated amount, a foundation is laid for realizing rapid analysis of the transient power angle stability under the influence of multiple coupling uncertain factors, and the method has certain theoretical and engineering significance.
A system for early terminating a transient power angle stability analysis numerical integration comprises,
a time judgment module: and judging whether the current simulation time reaches the system fault clearing time.
A leading image identification module: and identifying a leading group and a remaining group of the generator set according to the maximum rotor angle clearance at the current simulation moment to obtain a leading image in response to the fact that the current simulation moment reaches the system fault clearing moment.
A voltage information acquisition module: and acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set.
A stable situation/unstable situation maturity judging module: and judging whether the leading image stable state/unstable state is mature or not according to a preset maturation criterion, the bus voltage amplitude of each generator and the corresponding generator set.
The maturation criterion comprises a criterion for leading the maturation of the stable situation of the image and a criterion for leading the maturation of the unstable situation;
the criterion for leading the map to be mature in stable situation comprises the following steps:
criterion 1: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front collar group has the same trend along with the change of time, and the current simulation time t isnowThe highest value of bus voltage of the group generator set before the leader is not less than the lowest value of bus voltage of the rest group generator sets, wherein delta t is a numerical integration time step;
criterion 2: current simulation time tnowThe lowest voltage value of each generator bus is greater than 0.75 (per unit value);
criteria governing maturation of destabilization potentials include:
criterion 3: in the time zone [ tnow-Δt,tnow]Bus voltage amplitude of any two generators in inner and leading groups along with timeThe variation trends are consistent, the highest value of the bus voltage of the cluster generator set before the lead is not larger than the lowest value of the bus voltage of the rest cluster generator sets, and the generator sets connected with the buses corresponding to the lowest value of the bus voltage are the generator sets in the cluster before the lead;
criterion 4: current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.20 (per unit value); or the current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.40 (per unit value) and is in a time section [ tnow-Δt,tnow]And the voltage amplitude of the bus of any leading group generator set is reduced with time, and the gradient is less than-3.00.
The stable situation/unstable situation maturity judging module comprises a dominant image stable situation maturity judging module and a dominant image unstable situation maturity judging module; a dominant mapping stable situation maturity judging module: in response to the satisfaction of the criteria 1 and 2, judging that the stable situation of the leading image is mature; a dominant image instability situation maturity judgment module: and in response to the satisfaction of the criteria 3 and 4, judging that the instability situation of the dominant image is mature.
A quantitative analysis module: and responding to the maturity of the stable state/the unstable state of the dominant image, and carrying out quantitative analysis on the stability/the instability of the power angle.
A quantitative analysis conclusion judgment module: and judging whether the quantitative analysis conclusion is credible or not according to a preset credible criterion.
The credible criterion comprises a credible criterion of a quantitative analysis conclusion of power angle stability and a credible criterion of a quantitative analysis conclusion of power angle instability;
the credible criterion of the quantitative analysis conclusion of the power angle stability is as follows:
current simulation time tnowThe minimum value of the bus voltage of each generator is greater than a first threshold value, and the quantitative analysis conclusion is greater than a second threshold value;
credible criterion for quantitative analysis conclusion of power angle instability:
current simulation time tnowAnd the minimum value of the bus voltage of each generator is smaller than a third threshold, and the quantitative analysis conclusion is smaller than a fourth threshold.
A result module: and responding to the credibility of the quantitative analysis conclusion, outputting the quantitative analysis conclusion, and ending the transient power angle stability analysis.
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 a transient power angle stability analysis numerical integration early termination method.
A computing device comprising one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing a transient power angle stability analysis numerical integration early termination 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 terminating numerical integration in advance in transient power angle stability analysis is characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
in response to the current simulation time reaching the system fault clearing time, identifying a leading group and a remaining group of the generator set according to the maximum rotor angle gap at the current simulation time to obtain a leading image;
acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set;
judging whether the leading image stable state/unstable state is mature or not according to a preset maturation criterion, the bus voltage amplitude of each generator and the corresponding generator set;
responding to the maturity of the stable state/unstable state of the dominant image, and carrying out quantitative analysis on the stability/instability of the power angle;
judging whether the quantitative analysis conclusion is credible or not according to a preset credible criterion;
and responding to the credibility of the quantitative analysis conclusion, outputting the quantitative analysis conclusion, and ending the transient power angle stability analysis.
2. The method according to claim 1, wherein the method comprises: the maturation criterion comprises a criterion for leading the maturation of the stable situation of the image and a criterion for leading the maturation of the unstable situation;
the criterion for leading the map to be mature in stable situation comprises the following steps:
criterion 1: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front collar group has the same trend along with the change of time, and the current simulation time t isnowThe highest value of bus voltage of the group generator set before the leader is not less than the lowest value of bus voltage of the rest group generator sets, wherein delta t is a numerical integration time step;
criterion 2: current simulation time tnowThe lowest voltage value of each generator bus is greater than 0.75;
criteria governing maturation of destabilization potentials include:
criterion 3: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generator buses in the front leading group is consistent with the change trend of time, the highest value of the bus voltage of the front leading group generator is not more than the lowest value of the bus voltage of the rest group generator, and the generator connected with the bus corresponding to the lowest value of the bus is the generator set in the front leading group;
criterion 4: current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.20; or the current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.40, and the time zone t isnow-Δt,tnow]And the voltage amplitude of the bus of any leading group generator set is reduced with time, and the gradient is less than-3.00.
3. The method according to claim 2, wherein the method comprises: in response to the satisfaction of the criteria 1 and 2, judging that the stable situation of the leading image is mature; and in response to the satisfaction of the criteria 3 and 4, judging that the instability situation of the dominant image is mature.
4. The method according to claim 1, wherein the method comprises: the credible criterion comprises a credible criterion of a quantitative analysis conclusion of power angle stability and a credible criterion of a quantitative analysis conclusion of power angle instability;
the credible criterion of the quantitative analysis conclusion of the power angle stability is as follows:
current simulation time tnowThe minimum value of the bus voltage of each generator is greater than a first threshold value, and the quantitative analysis conclusion is greater than a second threshold value;
credible criterion for quantitative analysis conclusion of power angle instability:
current simulation time tnowAnd the minimum value of the bus voltage of each generator is smaller than a third threshold, and the quantitative analysis conclusion is smaller than a fourth threshold.
5. A transient state power angle stability analysis numerical integration early termination system is characterized in that: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
a leading image identification module: in response to the current simulation time reaching the system fault clearing time, identifying a leading group and a remaining group of the generator set according to the maximum rotor angle gap at the current simulation time to obtain a leading image;
a voltage information acquisition module: acquiring the bus voltage amplitude of each generator at the current simulation moment and the corresponding generator set;
a stable situation/unstable situation maturity judging module: judging whether the leading image stable state/unstable state is mature or not according to a preset maturation criterion, the bus voltage amplitude of each generator and the corresponding generator set;
a quantitative analysis module: responding to the maturity of the stable state/unstable state of the dominant image, and carrying out quantitative analysis on the stability/instability of the power angle;
a quantitative analysis conclusion judgment module: judging whether the quantitative analysis conclusion is credible or not according to a preset credible criterion;
a result module: and responding to the credibility of the quantitative analysis conclusion, outputting the quantitative analysis conclusion, and ending the transient power angle stability analysis.
6. The system according to claim 5, wherein the system comprises: the maturation criterion comprises a criterion for leading the maturation of the stable situation of the image and a criterion for leading the maturation of the unstable situation;
the criterion for leading the map to be mature in stable situation comprises the following steps:
criterion 1: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generators in the front collar group has the same trend along with the change of time, and the current simulation time t isnowThe highest value of bus voltage of the group generator set before the leader is not less than the lowest value of bus voltage of the rest group generator sets, wherein delta t is a numerical integration time step;
criterion 2: current simulation time tnowThe lowest voltage value of each generator bus is greater than 0.75;
criteria governing maturation of destabilization potentials include:
criterion 3: in the time zone [ tnow-Δt,tnow]The voltage amplitude of any two generator buses in the front leading group is consistent with the change trend of time, the highest value of the bus voltage of the front leading group generator is not more than the lowest value of the bus voltage of the rest group generator, and the generator connected with the bus corresponding to the lowest value of the bus is the generator set in the front leading group;
criterion 4: current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.20; or the current simulation time tnowThe lowest value of the bus voltage of the generator set is less than 0.40, and the time zone t isnow-Δt,tnow]And the voltage amplitude of the bus of any leading group generator set is reduced with time, and the gradient is less than-3.00.
7. The system of claim 6, wherein the system comprises: the stable situation/unstable situation maturity judging module comprises a dominant image stable situation maturity judging module and a dominant image unstable situation maturity judging module; a dominant mapping stable situation maturity judging module: in response to the satisfaction of the criteria 1 and 2, judging that the stable situation of the leading image is mature; a dominant image instability situation maturity judgment module: and in response to the satisfaction of the criteria 3 and 4, judging that the instability situation of the dominant image is mature.
8. The system according to claim 5, wherein the system comprises: the credible criterion comprises a credible criterion of a quantitative analysis conclusion of power angle stability and a credible criterion of a quantitative analysis conclusion of power angle instability;
the credible criterion of the quantitative analysis conclusion of the power angle stability is as follows:
current simulation time tnowThe minimum value of the bus voltage of each generator is greater than a first threshold value, and the quantitative analysis conclusion is greater than a second threshold value;
credible criterion for quantitative analysis conclusion of power angle instability:
current simulation time tnowAnd the minimum value of the bus voltage of each generator is smaller than a third threshold, and the quantitative analysis conclusion is smaller than a fourth threshold.
9. A computer readable storage medium storing one or more programs, characterized in that: the one or more programs include instructions that, when executed by a computing device, cause the computing device to perform any of the methods of claims 1-4.
10. A computing device, characterized by: comprises the steps of (a) preparing a mixture of a plurality of raw materials,
one or more processors, memory, and one or more programs stored in the memory and configured to be executed by the one or more processors, the one or more programs including instructions for performing any of the methods of claims 1-4.
CN202010691505.2A 2020-07-17 2020-07-17 Method and system for terminating numerical integration in advance in transient power angle stability analysis Active CN111832180B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010691505.2A CN111832180B (en) 2020-07-17 2020-07-17 Method and system for terminating numerical integration in advance in transient power angle stability analysis

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010691505.2A CN111832180B (en) 2020-07-17 2020-07-17 Method and system for terminating numerical integration in advance in transient power angle stability analysis

Publications (2)

Publication Number Publication Date
CN111832180A true CN111832180A (en) 2020-10-27
CN111832180B CN111832180B (en) 2022-07-01

Family

ID=72924313

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010691505.2A Active CN111832180B (en) 2020-07-17 2020-07-17 Method and system for terminating numerical integration in advance in transient power angle stability analysis

Country Status (1)

Country Link
CN (1) CN111832180B (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
CN111832180B (en) 2022-07-01

Similar Documents

Publication Publication Date Title
US7827122B1 (en) Data mining of unfiltered controller data
CN115097920B (en) Method and device for adjusting temperature of server case, server and medium
CN107838917A (en) Robot simulation apparatus
CN111369226A (en) BIM-based project progress supervision system
CN102509152A (en) Switched reluctance motor on-line modeling method based RBF neural network
Wessel et al. Ontology-based traceability system for interoperable data acquisition in battery cell manufacturing
CN110348540B (en) Clustering-based method and device for screening transient power angle stability faults of power system
CN110928267B (en) Method and device for evaluating stability of automobile electric control unit
CN111832180B (en) Method and system for terminating numerical integration in advance in transient power angle stability analysis
CN104008046A (en) Method and device for testing programs
CN107122285A (en) A kind of multithread process control Dynamic Announce and intelligence change System and method for
CN114064012B (en) Dynamic and static combined interface code generation method and system and electronic equipment
CN115829122A (en) Electric power customer service work order complaint early warning method and system
US20230124678A1 (en) Storage medium to store transmission data setting support program, gateway device, and transmission data setting supporting method
CN111796513B (en) Data processing method and device
CN111291464A (en) Dynamic equivalence method and device for power system
Zou et al. Analysis and evaluation of control action end-state impact on manufacturing system future production
CN114002948B (en) Method and device for accessing third-party numerical control machine tool to service platform
CN110286948A (en) Process based on the read-write of fictitious host computer storage device hangs up method and relevant device
CN117922538B (en) Measurement and control method and system for hybrid electric vehicle based on digital twin technology
US20240169209A1 (en) Computer-implemented method for producing a software image, suitable for a numerical simulation, of at least part of a real controller
CN111126620B (en) Characteristic fingerprint generation method for time sequence and application
JP2023002962A (en) Information processing apparatus, and program and method for generating model
CN116148610A (en) Method, system, equipment and storage medium for generating partial discharge signal of switch cabinet
CN113900784A (en) Method and device for determining task baseline time, electronic equipment and storage medium

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