Transient time domain simulation switching method and system of power system and storage medium
By adaptively switching electromechanical transient simulation and electromagnetic transient simulation in the power system, the credibility of the simulation results is used to identify the credibility of the simulation results, the problem of simulation accuracy and speed bottleneck in the existing technology is solved, and more accurate and efficient Transient time domain simulation of power system.
Patent Information
- Application Number
- CN202411539384.4
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- Filing Date
- 2024-10-31
- Publication Date
- 2025-02-14
AI Technical Summary
The electromechanical transient simulation algorithms of existing power systems have natural limitations in terms of model, integral step size and electrical quantity, making it difficult to accurately describe the multi-time scale characteristics, sequential switching characteristics and transient stability characteristics of the power system, resulting in online safety and stability analysis. Conclusion There is a possibility of scenario mismatch.
Provide a transient time domain simulation switching method for power system adaptive analysis and computing scenarios. By constructing an electromagnetic transient simulation system consistent with the electromechanical transient simulation system, electromagnetic transient simulation and electromechanical transient simulation are carried out simultaneously, and classifiers are used. Identify the credibility of the simulation results and adaptive switching simulation method.
It achieves the improvement of simulation calculation speed and adaptability while ensuring simulation accuracy, avoid scene mismatch, improves the level of grid characteristics cognition, and has engineering practicality.
Smart Images

Figure CN119442665A_ABST
Abstract
Description
Technical Field
[0001] The present invention relates to power system simulation technology, and in particular to a method, system and storage medium for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario. Background Art
[0002] The high proportion of renewable energy and power electronic equipment connected to the grid has aggravated the strong time-varying and uncertain factors of the power system. This is reflected in the research of online safety and stability time-domain simulation technology, which is manifested in the improvement of the accuracy of the required power grid data model and the expansion of the scale of the expected fault set. At present, researchers have carried out a lot of research on improving the accuracy of calculation data, improving the simulation calculation speed and reducing the number of simulation examples, so as to improve the computing power of online safety and stability analysis to adapt to the development of the power grid. However, the electromechanical transient simulation algorithm commonly used in the existing technology has natural limitations in terms of model, integration step and electrical quantity, and it is difficult to accurately describe the multi-time scale characteristics, sequential switching characteristics and transient stability characteristics of the power system. The online safety and stability analysis conclusion may have scenario mismatch. Electromagnetic transient simulation uses a three-phase instantaneous value model to simulate the change process of the electric field, magnetic field and corresponding voltage and current in each component of the power system. The nonlinear and frequency-related characteristics of the components can be considered in detail, and the simulation accuracy is high. The use of electromagnetic transient refinement model for grid connection analysis of power electronic equipment such as new energy and DC transmission can improve the level of understanding of power grid characteristics. However, electromagnetic transient simulation of large power grids has not yet broken through the bottleneck of simulation scale and computing speed, and it is difficult to completely replace traditional electromechanical transient simulation. Therefore, it is necessary to comprehensively consider electromechanical transient simulation and electromagnetic transient simulation and switch the simulation according to the situation. Summary of the invention
[0003] In view of the problems existing in the prior art, the purpose of the present invention is to provide a method, system, device and storage medium for switching transient time-domain simulation of a power system that can adaptively analyze and calculate scenarios.
[0004] In order to achieve the above-mentioned object of the invention, the present invention provides the following technical solutions:
[0005] A method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario, comprising:
[0006] According to the calculation data of the electromechanical transient simulation system, an electromagnetic transient simulation system consistent with the electromechanical transient simulation system is constructed;
[0007] The electromagnetic transient simulation system and the electromechanical transient simulation system are respectively used to synchronously perform the electromagnetic transient simulation and the electromechanical transient simulation considering the electromagnetic transient process of the simulation object until the fault is restored and a preset time period has passed;
[0008] According to the simulation results of electromagnetic transient simulation and electromechanical transient simulation carried out simultaneously, the power transmitted by the simulation object to the grid through the grid-connected port is calculated, and the characteristic quantity for identifying whether the electromechanical transient simulation result is credible is obtained according to the power calculation;
[0009] The feature quantity is input into a classifier generated offline in advance. If the output result of the classifier is credible, the electromechanical transient simulation is continued and the electromagnetic transient simulation is terminated. Otherwise, the electromagnetic transient simulation is switched to perform subsequent simulation processes and the electromechanical transient simulation is terminated. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the feature quantity.
[0010] Furthermore, when the simulation is switched, if a DC commutation failure or a new energy grid disconnection event occurs in the previous electromagnetic transient simulation, the electromagnetic transient simulation is always performed.
[0011] Furthermore, the construction of an electromagnetic transient simulation system consistent with the electromechanical transient simulation system based on the calculation data of the electromechanical transient simulation system specifically includes:
[0012] Obtain the calculation data of the electromechanical transient simulation system, including the power flow file, stability file and fault file of the electromechanical transient simulation system;
[0013] According to the parameters of nodes, lines, transformers in the power flow file and the parameters of dynamic components in the stability file, the parameters of the corresponding equipment in the electromagnetic transient simulation system are determined, and the topology of the electromagnetic transient simulation system is determined according to the network parameters in the stability file, thereby establishing the electromagnetic transient simulation system;
[0014] According to the fault information in the fault file, add corresponding timing switches in the electromagnetic transient simulation system;
[0015] The power flow calculation of the electromagnetic transient simulation system is performed based on the power flow file, and the relevant equipment in the electromagnetic transient simulation is initialized according to the power flow calculation results.
[0016] Furthermore, the electromagnetic transient simulation system and the electromechanical transient simulation system are respectively used to simultaneously perform electromagnetic transient simulation and electromechanical transient simulation considering the electromagnetic transient process of the simulation object, specifically including
[0017] Set the electromagnetic transient simulation step size to Δt, the electromechanical transient simulation step size to ΔT, and start the electromagnetic transient simulation system and the electromechanical transient simulation system at the same time; where ΔT = k·Δt, k is a positive integer;
[0018] Starting from the first electromechanical transient simulation step after the simulation starts, to the electromechanical transient simulation step before the preset time period after the fault is restored, every time an electromechanical transient simulation step ΔT is completed, according to the simulation results of the electromagnetic transient simulation system and the electromechanical transient simulation system in the previous ΔT time period, the computer electrical transient simulation will reach the compensation injection current of each simulation object when it is consistent with the electromagnetic transient simulation results, and the calculated compensation injection current will be injected into the electromechanical transient simulation system, so as to perform electromechanical transient simulation considering the electromagnetic transient process of the simulation object in the next ΔT time period.
[0019] Furthermore, the compensation injection current of each simulation object when the computer electrical transient simulation is to be consistent with the electromagnetic transient simulation result according to the simulation results of the electromagnetic transient simulation system and the electromechanical transient simulation system in the previous ΔT time period specifically includes:
[0020] According to the simulation results of the electromagnetic transient simulation system in the previous ΔT time period, the active power and reactive power transmitted to the grid through the grid-connected port of each simulation object in the electromagnetic transient simulation are calculated;
[0021] According to the simulation results of the electromechanical transient simulation system in the previous ΔT time period, the active power and reactive power transmitted to the grid through the grid-connected port by each simulation object in the computer electrical transient simulation;
[0022] The active power difference and reactive power difference of each simulation object in the electromagnetic transient simulation and the electromechanical transient simulation are calculated, and the compensation injection current of each simulation object is calculated in combination with the grid connection point voltage of each simulation object in the electromechanical transient simulation.
[0023] Further, the calculation of the active power and reactive power transmitted to the power grid through the grid-connected port of each simulation object in the electromagnetic transient simulation according to the simulation result of the electromagnetic transient simulation system in the previous ΔT time period specifically includes:
[0024] From the simulation results of the electromagnetic transient simulation system in the previous ΔT time period, the grid connection point voltage and injection current of each simulation object are obtained;
[0025] Perform fast Fourier transformation on the single-phase instantaneous values of the grid connection point voltage and injection current of each simulation object to obtain the ABC three-phase fundamental frequency phasor;
[0026] The positive-sequence voltage and positive-sequence current are extracted from the fundamental-frequency phasors of the three-phase ABC based on the symmetrical component method, and the active power and reactive power are calculated based on the positive-sequence voltage and positive-sequence current.
[0027] Furthermore, the calculation formula of the compensation injection current is:
[0028]
[0029] In the formula, IR i,j ,II i,j They represent the real and imaginary parts of the compensation injection current of the simulation object i at the jth ΔT moment, respectively. represents the active power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, It represents the reactive power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT moment in the electromechanical transient simulation. They represent the real part and imaginary part of the grid connection point voltage of simulation object i at the jth ΔT moment in the electromechanical transient simulation, and n is the total number of simulation objects.
[0030] Furthermore, the characteristic quantity for identifying whether the electromechanical transient simulation result is credible is obtained according to the power calculation, specifically including:
[0031] According to the active power and reactive power in the power, the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment are calculated according to the following formula:
[0032]
[0033]
[0034] Where ΔP i,j , ΔQ i,j They represent the active power difference ratio and reactive power difference ratio of the simulation object i at the jth ΔT moment, respectively. They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT in the electromechanical transient simulation, n is the total number of simulation objects, and m is the total number of ΔTs from the start to the end of the simulation;
[0035] According to the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment, the active power intermediate characteristic quantity and reactive power intermediate characteristic quantity of each simulation object are calculated respectively according to the following formula;
[0036]
[0037] In the formula, D(ΔP i)、D(ΔQ i ) represent the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of the simulation object i, respectively. They represent the average active power difference ratio and average reactive power difference ratio of simulation object i respectively;
[0038] According to the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of each simulation object, the active power characteristic quantity and the reactive power characteristic quantity are calculated according to the following formula, which are used as characteristic quantities to identify whether the electromechanical transient simulation results are credible:
[0039]
[0040] In the formula, They represent the active power characteristic quantity and the reactive power characteristic quantity respectively.
[0041] Furthermore, the classifier is generated in the following manner:
[0042] Obtain calculation data of several electromechanical transient simulation systems with anticipated faults;
[0043] According to the calculation data of each electromechanical transient simulation system, a corresponding electromagnetic transient simulation system consistent with the electromechanical transient simulation system is constructed;
[0044] For each electromechanical transient simulation system, two electromechanical transient simulation processes Process1 and Process2 are started, and an electromagnetic transient simulation process Process3 is started in the corresponding electromagnetic transient simulation system until the fault is restored and a preset time period has passed, wherein Process1 performs the whole process electromechanical transient simulation, and Process2 performs the electromechanical transient simulation considering the electromagnetic transient simulation process Process3 of the simulation object;
[0045] According to the simulation results, characteristic quantities for identifying whether the electromechanical transient simulation results are credible are calculated;
[0046] After the fault is recovered and a preset time period has passed, Process3 is terminated, and Process2 is modified to no longer consider the electromagnetic transient simulation process Process3 of the simulation object, and the electromechanical transient simulation is continued until the simulation termination time;
[0047] According to the simulation results of Process2, the first transient power angle stability margin η under fault is calculated hybrid , according to the simulation results of Process1, calculate the second transient power angle stability margin η under fault origin ;
[0048] According to η hybrid and η origin The margin difference of the simulation object is calculated as Δη=ηhybrid -η origin |, if Δη>Δη th , it means that the electromechanical transient simulation results are unreliable, otherwise it means that the electromechanical transient simulation results are reliable, where Δη th is the margin difference threshold;
[0049] Through machine learning methods, a nonlinear mapping relationship between the characteristic quantity and whether the electromechanical transient simulation results are credible is established as a classifier for identifying whether the electromechanical transient simulation results are credible.
[0050] A power system transient time domain simulation switching system for adaptive analysis and calculation scenarios, comprising:
[0051] A mirror simulation system building module is used to build an electromagnetic transient simulation system consistent with the electromechanical transient simulation system based on the calculation data of the electromechanical transient simulation system;
[0052] The electromechanical and electromagnetic synchronous simulation module is used to respectively use the electromagnetic transient simulation system and the electromechanical transient simulation system to synchronously perform electromagnetic transient simulation and electromechanical transient simulation considering the electromagnetic transient process of the simulation object until a preset time period has passed after the fault is restored;
[0053] A feature quantity information extraction module is used to calculate the power transmitted by the simulation object to the power grid through the grid-connected port according to the simulation results of the synchronous electromagnetic transient simulation and the electromechanical transient simulation, and obtain a feature quantity for identifying whether the electromechanical transient simulation result is credible according to the power calculation;
[0054] The simulation method discrimination module is used to input the feature quantity into a classifier pre-generated offline. If the output result of the classifier is credible, the electromechanical transient simulation continues to be used and the electromagnetic transient simulation is terminated. Otherwise, the electromagnetic transient simulation is switched to perform subsequent simulation processes and the electromechanical transient simulation is terminated. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the feature quantity.
[0055] A power system transient time domain simulation switching system for adaptive analysis and calculation scenarios, characterized by comprising:
[0056] A mirror simulation system building module is used to build an electromagnetic transient simulation system consistent with the electromechanical transient simulation system based on the calculation data of the electromechanical transient simulation system;
[0057] The electromechanical and electromagnetic synchronous simulation module is used to respectively use the electromagnetic transient simulation system and the electromechanical transient simulation system to synchronously perform electromagnetic transient simulation and electromechanical transient simulation considering the electromagnetic transient process of the simulation object until a preset time period has passed after the fault is restored;
[0058] A feature quantity information extraction module is used to calculate the power transmitted by the simulation object to the power grid through the grid-connected port according to the simulation results of the synchronous electromagnetic transient simulation and the electromechanical transient simulation, and obtain a feature quantity for identifying whether the electromechanical transient simulation result is credible according to the power calculation;
[0059] The simulation method discrimination module is used to input the feature quantity into a classifier pre-generated offline. If the output result of the classifier is credible, the electromechanical transient simulation continues to be used and the electromagnetic transient simulation is terminated. Otherwise, the electromagnetic transient simulation is switched to perform subsequent simulation processes and the electromechanical transient simulation is terminated. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the feature quantity.
[0060] A computer program product comprises a computer program / instruction, wherein the computer program / instruction implements the above method when executed by a processor.
[0061] A computer device comprises a memory, a processor and a computer program stored in the memory and executable on the processor, wherein the processor executes the computer program to implement the method.
[0062] A computer-readable storage medium having a computer program / instruction stored thereon, which implements the method when executed by a processor.
[0063] Compared with the prior art, the present invention has the following beneficial effects: the present invention can adaptively analyze the calculation scenario, and identify whether the simulation process can be continued by the electromechanical transient simulation through the offline generated classifier, thereby realizing the adaptive switching of electromagnetic transient simulation and electromechanical transient simulation, taking into account the accuracy and speed of power grid simulation, and having engineering practicality. The classifier of the present invention performs electromagnetic transient simulation and electromechanical transient simulation in a very short time from the occurrence of the fault to the clearing of the fault for each calculation example of the expected fault set, and has high recognition accuracy and precision. BRIEF DESCRIPTION OF THE DRAWINGS
[0064] Figure 1 It is a flow chart of an embodiment of a method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario provided by the present invention;
[0065] Figure 2 It is a flow chart of an embodiment of the offline generation method of a classifier provided by the present invention. DETAILED DESCRIPTION
[0066] The technical solutions in the embodiments of the present invention will be described clearly and completely below in conjunction with the accompanying drawings in the embodiments of the present invention.
[0067] Embodiment 1
[0068] The embodiment of the present invention provides a method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario, such as Figure 1 As shown, including:
[0069] S101. Constructing an electromagnetic transient simulation system consistent with the electromechanical transient simulation system according to calculation data of the electromechanical transient simulation system.
[0070] This step specifically includes:
[0071] Obtain the calculation data of the electromechanical transient simulation system, including the power flow file, stability file and fault file of the electromechanical transient simulation system;
[0072] According to the parameters of nodes, lines, transformers in the power flow file and the parameters of dynamic components in the stability file, the parameters of the corresponding equipment in the electromagnetic transient simulation system are determined, and the topology of the electromagnetic transient simulation system is determined according to the network parameters in the stability file, thereby establishing the electromagnetic transient simulation system;
[0073] According to the fault information in the fault file, add corresponding timing switches in the electromagnetic transient simulation system;
[0074] The power flow calculation of the electromagnetic transient simulation system is performed based on the power flow file, and the relevant equipment in the electromagnetic transient simulation is initialized according to the power flow calculation results.
[0075] S102, respectively using an electromagnetic transient simulation system and an electromechanical transient simulation system to synchronously perform electromagnetic transient simulation and electromechanical transient simulation taking into account the electromagnetic transient process of the simulation object until a preset time period has passed after the fault is restored.
[0076] In the specific implementation, the electromagnetic transient simulation step is set to Δt, the electromechanical transient simulation step is set to ΔT, and the electromagnetic transient simulation system and the electromechanical transient simulation system are turned on at the same time; wherein, ΔT = k·Δt, k is a positive integer, and the default value is 200; starting from the first ΔT after the simulation starts, to a ΔT before the preset time period after the fault is restored, every time a ΔT is reached, according to the simulation results of the electromagnetic transient simulation system and the electromechanical transient simulation system in the previous ΔT time period, the computer electrical transient simulation reaches the compensation injection current of each simulation object when it is consistent with the electromagnetic transient simulation result, and the calculated compensation injection current is injected into the electromechanical transient simulation system, so as to perform the electromechanical transient simulation considering the electromagnetic transient process of the simulation object in the next ΔT time period. The interaction of electromechanical transient simulation and electromagnetic transient simulation results is realized through real-time network communication functions such as SOCKET.
[0077] The calculation steps of the compensation injection current specifically include:
[0078] According to the simulation results of the electromagnetic transient simulation system in the previous ΔT time period, calculate the active power transmitted to the grid through the grid-connected port of each simulation object in the electromagnetic transient simulation and reactive power The calculation of active power and reactive power can be specifically performed by obtaining the grid connection point voltage and injection current of each simulation object from the simulation results of the electromagnetic transient simulation system in the previous ΔT time period; performing fast Fourier transform (FFT) on the single-phase instantaneous values of the grid connection point voltage and injection current of each simulation object to obtain the ABC three-phase fundamental frequency phasor; extracting the positive sequence voltage and positive sequence current from the ABC three-phase fundamental frequency phasor based on the symmetrical component method, and calculating the active power and reactive power based on the positive sequence voltage and positive sequence current;
[0079] According to the simulation results of the electromechanical transient simulation system in the previous ΔT time period, the active power transmitted to the grid by each simulation object in the computer electrical transient simulation through the grid-connected port is and reactive power Active Power and reactive power It is directly calculated through the grid connection point voltage and injection current of the first simulation object;
[0080] Calculate the active power difference and reactive power difference of each simulation object in the electromagnetic transient simulation and the electromechanical transient simulation, and combine the grid connection point voltage of each simulation object in the electromechanical transient simulation to obtain the compensation injection current of each simulation object according to the following formula;
[0081]
[0082] In the formula, IR i,j , II i,j They represent the real and imaginary parts of the compensation injection current of the simulation object i at the jth ΔT moment, respectively. represents the active power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, It represents the reactive power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT moment in the electromechanical transient simulation. They represent the real part and imaginary part of the grid connection point voltage of simulation object i at the jth ΔT moment in the electromechanical transient simulation, and n is the total number of simulation objects.
[0083] Simulate to the specified time after the fault is restored, that is, after T aDeadline after time, T a is a preset time period. At each ΔT, the compensation injection current is injected into the electromechanical transient simulation system to perform electromechanical transient simulation considering the electromagnetic transient process of the simulation object.
[0084] S103, according to the simulation results of the electromagnetic transient simulation and the electromechanical transient simulation performed simultaneously, the power transmitted by the simulation object to the power grid through the grid-connected port is calculated, and a characteristic quantity for identifying whether the electromechanical transient simulation result is credible is obtained according to the power calculation.
[0085] The calculation steps of the feature quantity specifically include:
[0086] Active power and reactive power are obtained from the power. The active power and reactive power can be calculated according to the steps described in step S102. The active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment are calculated according to the following formulas:
[0087]
[0088] Where ΔP i,j , ΔQ i,j They represent the active power difference ratio and reactive power difference ratio of the simulation object i at the jth ΔT moment, respectively. They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT in the electromechanical transient simulation, n is the total number of simulation objects, and m is the total number of ΔTs from the start to the end of the simulation;
[0089] According to the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment, the active power intermediate characteristic quantity and reactive power intermediate characteristic quantity of each simulation object are calculated respectively according to the following formula;
[0090]
[0091] In the formula, D(ΔP i )、D(ΔQ i ) represent the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of the simulation object i, respectively. They represent the average active power difference ratio and average reactive power difference ratio of simulation object i respectively;
[0092] According to the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of each simulation object, the active power characteristic quantity and the reactive power characteristic quantity are calculated according to the following formula, which are used as characteristic quantities to identify whether the electromechanical transient simulation results are credible:
[0093]
[0094] In the formula, They represent the active power characteristic quantity and the reactive power characteristic quantity respectively.
[0095] S104, input the characteristic quantity into the classifier generated in advance offline. If the output result of the classifier is credible, continue to use the electromechanical transient simulation and terminate the electromagnetic transient simulation. Otherwise, switch to the electromagnetic transient simulation, perform the subsequent simulation process, and terminate the electromechanical transient simulation. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the characteristic quantity.
[0096] Among them, when the simulation is switched, if a DC commutation failure or a new energy grid disconnection event occurs in the previous electromagnetic transient simulation, the electromagnetic transient simulation is always performed.
[0097] Specifically, Figure 2 As shown, the classifier is generated in the following way:
[0098] Acquire calculation data of several electromechanical transient simulation systems with expected faults; and construct a corresponding electromagnetic transient simulation system consistent with the electromechanical transient simulation system according to the calculation data of each electromechanical transient simulation system; the specific steps are the same as step S101;
[0099] For each electromechanical transient simulation system, two electromechanical transient simulation processes Process1 and Process2 are started, and an electromagnetic transient simulation process Process3 is started in the corresponding electromagnetic transient simulation system until the fault is restored and a preset time period has passed, wherein Process1 performs a full-process electromechanical transient simulation, and Process2 performs an electromechanical transient simulation considering the electromagnetic transient simulation process Process3 of the simulation object. The specific simulation process is the same as the electromechanical transient simulation considering the electromagnetic transient process of the simulation object in step S102;
[0100] A feature quantity for identifying whether the electromechanical transient simulation result is credible is calculated based on the simulation result; the feature quantity calculation step is specifically the same as that in step S103;
[0101] After the fault is recovered and a preset time period has passed, Process3 is terminated, and Process2 is modified to no longer consider the electromagnetic transient simulation process Process3 of the simulation object, and the electromechanical transient simulation is continued until the simulation termination time;
[0102] According to the simulation results of Process2, the Extended Equal Area Criterion (EEAC) algorithm is used to calculate the first transient power angle stability margin η under fault conditions. hybridAccording to the simulation results of Process1, the EEAC algorithm is used to calculate the second transient power angle stability margin η under fault conditions. origin ;
[0103] According to η hybrid and η origin The margin difference of the simulation object is calculated as Δη=η hybrid -η origin |, if Δη>Δη th , it means that the electromechanical transient simulation results are unreliable, otherwise it means that the electromechanical transient simulation results are reliable, where Δη th is the margin difference threshold;
[0104] Through machine learning methods, such as SVM, etc., a nonlinear mapping relationship between the feature quantity and whether the electromechanical transient simulation result is credible is established as a classifier for identifying whether the electromechanical transient simulation result is credible.
[0105] Embodiment 2
[0106] An embodiment of the present invention provides a power system transient time domain simulation switching system for adaptive analysis and calculation scenarios, including:
[0107] A mirror simulation system building module is used to build an electromagnetic transient simulation system consistent with the electromechanical transient simulation system based on the calculation data of the electromechanical transient simulation system;
[0108] The electromechanical and electromagnetic synchronous simulation module is used to respectively use the electromagnetic transient simulation system and the electromechanical transient simulation system to synchronously perform electromagnetic transient simulation and electromechanical transient simulation considering the electromagnetic transient process of the simulation object until a preset time period has passed after the fault is restored;
[0109] A feature quantity information extraction module is used to calculate the power transmitted by the simulation object to the power grid through the grid-connected port according to the simulation results of the synchronous electromagnetic transient simulation and the electromechanical transient simulation, and obtain a feature quantity for identifying whether the electromechanical transient simulation result is credible according to the power calculation;
[0110] The simulation method discrimination module is used to input the characteristic quantity into the classifier pre-generated offline. If the output result of the classifier is credible, the electromechanical transient simulation continues to be used, and the electromagnetic transient simulation is terminated at the same time. Otherwise, it switches to the electromagnetic transient simulation, performs the subsequent simulation process, and terminates the electromechanical transient simulation at the same time. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the characteristic quantity. When the simulation is switched, if a DC commutation failure or a new energy grid disconnection event occurs in the previous electromagnetic transient simulation, the electromagnetic transient simulation is always performed.
[0111] The mirror simulation system building module specifically includes:
[0112] A simulation data acquisition unit, used to acquire calculation data of the electromechanical transient simulation system, including a power flow file, a stability file and a fault file of the electromechanical transient simulation system;
[0113] A simulation parameter determination unit is used to determine the parameters of the corresponding equipment in the electromagnetic transient simulation system according to the node, line, transformer parameters in the power flow file and the dynamic component parameters in the stability file, and to determine the topology of the electromagnetic transient simulation system according to the network parameters in the stability file, so as to establish the electromagnetic transient simulation system;
[0114] A fault information confirmation unit, used to add a corresponding timing switch in the electromagnetic transient simulation system according to the fault information in the fault file;
[0115] The initialization unit is used to perform power flow calculation of the electromagnetic transient simulation system based on the power flow file, and initialize the relevant equipment in the electromagnetic transient simulation according to the power flow calculation results.
[0116] The electromechanical and electromagnetic synchronous simulation module specifically includes:
[0117] A simulation setting unit, used to set the electromagnetic transient simulation step length to Δt, the electromechanical transient simulation step length to ΔT, and start the electromagnetic transient simulation system and the electromechanical transient simulation system at the same time; wherein ΔT=k·Δt, k is a positive integer;
[0118] A step timing unit is used to start from the first electromechanical transient simulation step after the simulation starts and end at an electromechanical transient simulation step before a preset time period after the fault is restored. Whenever an electromechanical transient simulation step ΔT is completed, the compensation injection current calculation unit is executed;
[0119] A compensation injection current calculation unit is used to calculate the compensation injection current of each simulation object when the computer electrical transient simulation is to be consistent with the electromagnetic transient simulation result according to the simulation results of the electromagnetic transient simulation system and the electromechanical transient simulation system in the previous ΔT time period;
[0120] The compensation current injection unit is used to inject the calculated compensation injection current into the electromechanical transient simulation system, so as to perform electromechanical transient simulation considering the electromagnetic transient process of the simulation object in the next ΔT time period.
[0121] Wherein, the compensation injection current calculation unit specifically includes:
[0122] The electromagnetic transient simulation power calculation subunit is used to calculate the active power and reactive power transmitted to the power grid through the grid-connected port of each simulation object in the electromagnetic transient simulation according to the simulation results of the electromagnetic transient simulation system in the previous ΔT time period; specifically, the grid-connected point voltage and injected current of each simulation object are obtained from the simulation results of the electromagnetic transient simulation system in the previous ΔT time period; the single-phase instantaneous values of the grid-connected point voltage and injected current of each simulation object are fast Fourier transformed to obtain the ABC three-phase fundamental frequency phasor; the positive sequence voltage and positive sequence current are extracted from the ABC three-phase fundamental frequency phasor based on the symmetrical component method, and the active power and reactive power are calculated according to the positive sequence voltage and positive sequence current;
[0123] The electromechanical transient simulation power calculation subunit is used to calculate the active power and reactive power transmitted to the power grid through the grid-connected port of each simulation object in the computer electrical transient simulation according to the simulation results of the electromechanical transient simulation system in the previous ΔT time period;
[0124] The compensation injection current calculation subunit is used to calculate the active power difference and reactive power difference of each simulation object in the electromagnetic transient simulation and the electromechanical transient simulation, and combine the grid connection point voltage of each simulation object in the electromechanical transient simulation to calculate the compensation injection current of each simulation object according to the following formula:
[0125]
[0126] In the formula, IR i,j ,II i,j They represent the real and imaginary parts of the compensation injection current of the simulation object i at the jth ΔT moment, respectively. represents the active power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, It represents the reactive power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT moment in the electromechanical transient simulation. They represent the real part and imaginary part of the grid connection point voltage of simulation object i at the jth ΔT moment in the electromechanical transient simulation, and n is the total number of simulation objects.
[0127] The feature quantity information extraction module specifically includes:
[0128] The power difference ratio calculation unit is used to calculate the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment according to the active power and reactive power in the power according to the following formula:
[0129]
[0130] Where ΔP i,j , ΔQ i,j They represent the active power difference ratio and reactive power difference ratio of the simulation object i at the jth ΔT moment, respectively. They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT in the electromechanical transient simulation, n is the total number of simulation objects, and m is the total number of ΔTs from the start to the end of the simulation;
[0131] The intermediate characteristic quantity calculation unit is used to calculate the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of each simulation object according to the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment according to the following formula;
[0132]
[0133] In the formula, D(ΔP i )、D(ΔQ i ) represent the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of the simulation object i, respectively. They represent the average active power difference ratio and average reactive power difference ratio of simulation object i respectively;
[0134] The characteristic quantity calculation unit is used to calculate the active power characteristic quantity and the reactive power characteristic quantity according to the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of each simulation object according to the following formula, as the characteristic quantity for identifying whether the electromechanical transient simulation result is credible:
[0135]
[0136] In the formula, They represent the active power characteristic quantity and the reactive power characteristic quantity respectively.
[0137] The system also includes a classifier offline generation module, which specifically includes:
[0138] A mirror simulation construction unit is used to obtain calculation data of a number of electromechanical transient simulation systems with anticipated faults, and to construct a corresponding electromagnetic transient simulation system consistent with the electromechanical transient simulation system according to the calculation data of each electromechanical transient simulation system;
[0139] The synchronous simulation unit is used to start two electromechanical transient simulation processes Process1 and Process2 for each electromechanical transient simulation system, and start an electromagnetic transient simulation process Process3 in the corresponding electromagnetic transient simulation system, until the fault is restored and a preset time period has passed, wherein Process1 performs the whole process electromechanical transient simulation, and Process2 performs the electromechanical transient simulation considering the electromagnetic transient simulation process Process3 of the simulation object;
[0140] A characteristic quantity calculation unit, used for calculating the characteristic quantity for identifying whether the electromechanical transient simulation result is credible according to the simulation result;
[0141] The simulation modification unit is used to terminate Process3 after a preset time period after recovery from the fault, and modify Process2 to no longer consider the electromagnetic transient simulation process Process3 of the simulation object, and continue the electromechanical transient simulation until the simulation termination time;
[0142] The stability margin calculation unit is used to calculate the first transient power angle stability margin η under fault conditions according to the Process2 simulation results. hybrid , according to the simulation results of Process1, calculate the second transient power angle stability margin η under fault origin ;
[0143] The simulation result credible result acquisition unit is bred according to η hybrid and η origin The margin difference of the simulation object is calculated as Δη=η hybrid -η origin |, if Δη>Δη th , it means that the electromechanical transient simulation results are unreliable, otherwise it means that the electromechanical transient simulation results are reliable, where Δη th is the margin difference threshold;
[0144] The classifier construction unit is used to establish a nonlinear mapping relationship between the characteristic quantity and whether the electromechanical transient simulation result is credible through a machine learning method, as a classifier for identifying whether the electromechanical transient simulation result is credible.
[0145] The system provided in the embodiment of the present invention can be used to execute the method provided in the first embodiment of the present invention, and has the corresponding functions and beneficial effects of the execution method.
[0146] It is worth noting that in the embodiment of the above system, the various units and modules included are only divided according to functional logic, but are not limited to the above division, as long as the corresponding functions can be achieved; in addition, the specific names of the functional units are only for the convenience of distinguishing each other, and are not used to limit the scope of protection of the present invention.
[0147] The embodiments described above are only illustrative, wherein the modules described as separate components may or may not be physically separated, and the components displayed as modules may or may not be physical modules, i.e., they may be located in one place, or they may be distributed on multiple network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the scheme of this embodiment. Those skilled in the art can clearly understand that each implementation method can be implemented by means of software plus a necessary general hardware platform, or of course, it can be implemented only by hardware, as long as the function or effect can be achieved.
[0148] Embodiment 3
[0149] The embodiment of the present invention provides a computer product, such as an app on a mobile phone, a tablet, an installation program on a computer, etc., which includes a computer program / instruction, and the computer program / instruction is executed by a processor to implement the method described in Example 1. The code of the computer executable program for performing the operation of the present invention can be written in one or more programming languages or a combination thereof, and the programming language includes an object-oriented programming language, such as Java, Smalltalk, C++, and also includes a conventional procedural programming language, such as "C" language or similar programming language. The program code can be executed completely on the user's computer, partially on the user's computer, as an independent software package, partially on the user's computer, partially on the remote computer, or completely on the remote computer or server. In the case of a remote computer, the remote computer can be connected to the user's computer through any type of network including a local area network (LAN) or a wide area network (WAN), or it can be connected to an external computer (for example, using an Internet service provider to connect through the Internet).
[0150] Embodiment 4
[0151] The embodiment of the present invention provides a computer device, which provides services for implementing the method of the above-mentioned embodiment 1. The device may include: a memory storing a computer executable program; a processor coupled to the memory; the processor calls the computer executable program stored in the memory to execute the steps in the method described in the embodiment 1.
[0152] The memory may include computer system readable media in the form of volatile memory, such as random access memory (RAM) and / or cache memory. The device may further include other removable / non-removable, volatile / non-volatile computer system storage media. By way of example only, the memory may be used to read and write non-removable, non-volatile magnetic media (commonly referred to as "hard drives"). A program / utility having a set (at least one) of program modules may be stored in, for example, the memory, such program modules including but not limited to an operating system, one or more application programs, other program modules, and program data, each of which or some combination may include an implementation of a network environment. The computer executable program of the program module typically performs the functions and / or methods in the embodiments described in the present invention.
[0153] The processor executes various functional applications and data processing by running the program stored in the memory, such as implementing the method provided in the first embodiment of the present invention.
[0154] The code of the computer executable program can be written in one or more programming languages or a combination thereof, including object-oriented programming languages such as Java, Smalltalk, C++, and conventional procedural programming languages such as "C" language or similar programming languages.
[0155] Embodiment 5
[0156] An embodiment of the present invention provides a storage medium including a computer executable program. When the computer executable program is executed by a computer processor, it is used to execute the method of the first embodiment.
[0157] The storage medium of the embodiment of the present invention can adopt any combination of one or more computer-readable media. Computer-readable media can be computer-readable signal media or computer-readable storage media. Computer-readable storage media can be, for example, but not limited to, electrical, magnetic, optical, electromagnetic, infrared, or semiconductor systems, devices or devices, or any combination of the above. More specific examples (non-exhaustive list) of computer-readable storage media include: electrical connections with one or more wires, portable computer disks, hard disks, random access memories (RAM), read-only memories (ROM), erasable programmable read-only memories (EPROM or flash memory), optical fibers, portable compact disk read-only memories (CD-ROM), optical storage devices, magnetic storage devices, or any suitable combination of the above. In this document, computer-readable storage media can be any tangible medium containing or storing a program, which can be used by an instruction execution system, device or device or used in combination with it.
[0158] Of course, the storage medium containing a computer executable program provided by an embodiment of the present invention, the computer executable program of which is not limited to the above method operations, can also execute related operations in the method provided by any embodiment of the present invention.
[0159] It should be understood that the above embodiments and descriptions only describe the principles, main features and advantages of the present invention. Without departing from the spirit and scope of the present invention, the present invention may be subject to various changes and improvements, and these changes and improvements all fall within the scope of protection of the present invention.
Claims
1. A method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario, characterized in that: include: According to the calculation data of the electromechanical transient simulation system, an electromagnetic transient simulation system consistent with the electromechanical transient simulation system is constructed; The electromagnetic transient simulation system and the electromechanical transient simulation system are respectively used to synchronously perform the electromagnetic transient simulation and the electromechanical transient simulation considering the electromagnetic transient process of the simulation object until the fault is restored and a preset time period has passed; According to the simulation results of electromagnetic transient simulation and electromechanical transient simulation carried out simultaneously, the power transmitted by the simulation object to the grid through the grid-connected port is calculated, and the characteristic quantity for identifying whether the electromechanical transient simulation result is credible is obtained according to the power calculation; The feature quantity is input into a classifier generated offline in advance. If the output result of the classifier is credible, the electromechanical transient simulation is continued and the electromagnetic transient simulation is terminated. Otherwise, the electromagnetic transient simulation is switched to perform subsequent simulation processes and the electromechanical transient simulation is terminated. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the feature quantity.
2. The method for switching power system transient time domain simulation in the adaptive analysis and calculation scenario according to claim 1 is characterized in that: During simulation switching, if a DC commutation failure or a new energy grid disconnection event occurs in the previous electromagnetic transient simulation, the electromagnetic transient simulation will always be performed.
3. The method for switching power system transient time domain simulation in the adaptive analysis and calculation scenario according to claim 1 is characterized in that: The step of constructing an electromagnetic transient simulation system consistent with the electromechanical transient simulation system based on the calculation data of the electromechanical transient simulation system specifically includes: Obtain the calculation data of the electromechanical transient simulation system, including the power flow file, stability file and fault file of the electromechanical transient simulation system; According to the parameters of nodes, lines, transformers in the power flow file and the parameters of dynamic components in the stability file, the parameters of the corresponding equipment in the electromagnetic transient simulation system are determined, and the topology of the electromagnetic transient simulation system is determined according to the network parameters in the stability file, thereby establishing the electromagnetic transient simulation system; According to the fault information in the fault file, add corresponding timing switches in the electromagnetic transient simulation system; The power flow calculation of the electromagnetic transient simulation system is performed based on the power flow file, and the relevant equipment in the electromagnetic transient simulation is initialized according to the power flow calculation results.
4. The method for switching power system transient time domain simulation in the adaptive analysis and calculation scenario according to claim 1 is characterized in that: The electromagnetic transient simulation system and the electromechanical transient simulation system are respectively used to synchronously perform electromagnetic transient simulation and electromechanical transient simulation considering the electromagnetic transient process of the simulation object, specifically including: Set the electromagnetic transient simulation step size to Δt, the electromechanical transient simulation step size to ΔT, and start the electromagnetic transient simulation system and the electromechanical transient simulation system at the same time; where ΔT = k·Δt, k is a positive integer; Starting from the first electromechanical transient simulation step after the simulation starts, to the electromechanical transient simulation step before the preset time period after the fault is restored, every time an electromechanical transient simulation step ΔT is completed, according to the simulation results of the electromagnetic transient simulation system and the electromechanical transient simulation system in the previous ΔT time period, the computer electrical transient simulation will reach the compensation injection current of each simulation object when it is consistent with the electromagnetic transient simulation results, and the calculated compensation injection current will be injected into the electromechanical transient simulation system, so as to perform electromechanical transient simulation considering the electromagnetic transient process of the simulation object in the next ΔT time period.
5. The method for switching power system transient time domain simulation in the adaptive analysis and calculation scenario according to claim 4 is characterized in that: The compensation injection current of each simulation object when the computer electrical transient simulation is to be consistent with the electromagnetic transient simulation result according to the simulation results of the electromagnetic transient simulation system and the electromechanical transient simulation system in the previous ΔT time period specifically includes: According to the simulation results of the electromagnetic transient simulation system in the previous ΔT time period, the active power and reactive power transmitted to the grid through the grid-connected port of each simulation object in the electromagnetic transient simulation are calculated; According to the simulation results of the electromechanical transient simulation system in the previous ΔT time period, the active power and reactive power transmitted to the grid through the grid-connected port by each simulation object in the computer electrical transient simulation; The active power difference and reactive power difference of each simulation object in the electromagnetic transient simulation and the electromechanical transient simulation are calculated, and the compensation injection current of each simulation object is calculated in combination with the grid connection point voltage of each simulation object in the electromechanical transient simulation.
6. The method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario according to claim 5 is characterized in that: The method of calculating the active power and reactive power transmitted to the power grid through the grid-connected port of each simulation object in the electromagnetic transient simulation according to the simulation result of the electromagnetic transient simulation system in the previous ΔT time period specifically includes: From the simulation results of the electromagnetic transient simulation system in the previous ΔT time period, the grid connection point voltage and injection current of each simulation object are obtained; Perform fast Fourier transformation on the single-phase instantaneous values of the grid connection point voltage and injection current of each simulation object to obtain the ABC three-phase fundamental frequency phasor; The positive-sequence voltage and positive-sequence current are extracted from the fundamental-frequency phasors of the three-phase ABC based on the symmetrical component method, and the active power and reactive power are calculated based on the positive-sequence voltage and positive-sequence current.
7. The method for switching power system transient time domain simulation in the adaptive analysis and calculation scenario according to claim 5 is characterized in that: The calculation formula of the compensation injection current is: In the formula, IR i,j ,II i,j They represent the real and imaginary parts of the compensation injection current of the simulation object i at the jth ΔT moment, respectively. represents the active power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, It represents the reactive power difference of the simulation object i in the electromagnetic transient simulation and the electromechanical transient simulation at the jth ΔT moment, They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT moment in the electromechanical transient simulation. They represent the real part and imaginary part of the grid connection point voltage of simulation object i at the jth ΔT moment in the electromechanical transient simulation, and n is the total number of simulation objects.
8. The method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario according to claim 1, characterized in that: The characteristic quantity for identifying whether the electromechanical transient simulation result is credible according to the power calculation specifically includes: According to the active power and reactive power in the power, the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment are calculated according to the following formula: Where ΔP i,j , ΔQ i,j They represent the active power difference ratio and reactive power difference ratio of the simulation object i at the jth ΔT moment, respectively. They represent the active power and reactive power transmitted to the grid through the grid-connected port by the simulation object i at the jth ΔT moment in the electromagnetic transient simulation. They represent the active power and reactive power transmitted to the grid through the grid-connected port of the simulation object i at the jth ΔT in the electromechanical transient simulation, n is the total number of simulation objects, and m is the total number of ΔTs from the start to the end of the simulation; According to the active power difference ratio and reactive power difference ratio of each simulation object at each ΔT moment, the active power intermediate characteristic quantity and reactive power intermediate characteristic quantity of each simulation object are calculated respectively according to the following formula; In the formula, D(ΔP i )、D(ΔQ i ) represent the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of the simulation object i, respectively. They represent the average active power difference ratio and average reactive power difference ratio of simulation object i respectively; According to the intermediate characteristic quantity of active power and the intermediate characteristic quantity of reactive power of each simulation object, the active power characteristic quantity and the reactive power characteristic quantity are calculated according to the following formula, which are used as characteristic quantities to identify whether the electromechanical transient simulation results are credible: In the formula, They represent the active power characteristic quantity and the reactive power characteristic quantity respectively.
9. The method for switching transient time-domain simulation of a power system in an adaptive analysis and calculation scenario according to claim 1, characterized in that: The classifier is generated in the following way: Obtain calculation data of several electromechanical transient simulation systems with anticipated faults; According to the calculation data of each electromechanical transient simulation system, a corresponding electromagnetic transient simulation system consistent with the electromechanical transient simulation system is constructed; For each electromechanical transient simulation system, two electromechanical transient simulation processes Process1 and Process2 are started, and an electromagnetic transient simulation process Process3 is started in the corresponding electromagnetic transient simulation system until the fault is restored and a preset time period has passed, wherein Process1 performs the whole process electromechanical transient simulation, and Process2 performs the electromechanical transient simulation considering the electromagnetic transient simulation process Process3 of the simulation object; According to the simulation results, characteristic quantities for identifying whether the electromechanical transient simulation results are credible are calculated; After the fault is recovered and a preset time period has passed, Process3 is terminated, and Process2 is modified to no longer consider the electromagnetic transient simulation process Process3 of the simulation object, and the electromechanical transient simulation is continued until the simulation termination time; According to the simulation results of Process2, the first transient power angle stability margin η under fault is calculated hybrid , according to the simulation results of Process1, calculate the second transient power angle stability margin η under fault origin ; According to η hybrid and η origin The margin difference of the simulation object is calculated as Δη=η hybrid -η origin |, if Δη>Δη th , it means that the electromechanical transient simulation results are unreliable, otherwise it means that the electromechanical transient simulation results are reliable, where Δη th is the margin difference threshold; Through machine learning methods, a nonlinear mapping relationship between the characteristic quantity and whether the electromechanical transient simulation results are credible is established as a classifier for identifying whether the electromechanical transient simulation results are credible.
10. A power system transient time domain simulation switching system for adaptive analysis and calculation scenarios, characterized in that: include: A mirror simulation system building module is used to build an electromagnetic transient simulation system consistent with the electromechanical transient simulation system based on the calculation data of the electromechanical transient simulation system; The electromechanical and electromagnetic synchronous simulation module is used to respectively use the electromagnetic transient simulation system and the electromechanical transient simulation system to synchronously perform electromagnetic transient simulation and electromechanical transient simulation considering the electromagnetic transient process of the simulation object until a preset time period has passed after the fault is restored; A feature quantity information extraction module is used to calculate the power transmitted by the simulation object to the power grid through the grid-connected port according to the simulation results of the synchronous electromagnetic transient simulation and the electromechanical transient simulation, and obtain a feature quantity for identifying whether the electromechanical transient simulation result is credible according to the power calculation; The simulation method discrimination module is used to input the feature quantity into a classifier pre-generated offline. If the output result of the classifier is credible, the electromechanical transient simulation continues to be used and the electromagnetic transient simulation is terminated. Otherwise, the electromagnetic transient simulation is switched to perform subsequent simulation processes and the electromechanical transient simulation is terminated. The classifier is used to identify whether the electromechanical transient simulation result is credible based on the feature quantity.
11. A computer program product comprising a computer program / instructions, characterized in that When the computer program / instructions are executed by a processor, the method according to any one of claims 1 to 9 is implemented.
12. A computer device comprising a memory, a processor, and a computer program stored in the memory and executable on the processor, wherein: The processor executes the computer program to implement the method according to any one of claims 1 to 9.
13. A computer-readable storage medium having a computer program / instruction stored thereon, characterized in that: The computer program / instructions, when executed by a processor, implement the method of any one of claims 1-9.