PSD-based power system transient stability automatic batch calculation and analysis method

Abstract

The application discloses a kind of based on PSD's power system transient stability automatic batch calculation analysis method, belong to the technical field of power system electromechanical transient stability analysis.The method is realized based on PSD-BPA platform, and core includes line data intelligent analysis and check, standardized fault scene automatic generation, transient batch calculation, adaptive convergence adjustment, intelligent stability determination and post-processing, report export six big steps.The application can automatically generate the fault set in line with the "power system safety and stability guide", support the targeted investigation of user-defined fault subset, improve the stability evaluation accuracy through the dynamic stability determination algorithm based on time-domain energy attenuation characteristics, have adaptive convergence adjustment capability, realize the full-process automation closed loop from data verification to report generation, greatly improve the efficiency and accuracy of power grid transient stability analysis, provide efficient technical support for power grid planning and operation safety analysis.

Description

Technical Field

[0001] This invention belongs to the field of electromechanical transient stability analysis technology of power systems, specifically relating to an automatic batch calculation and analysis method for transient stability of power systems based on PSD. Background Technology

[0002] Transient stability calculation is a core component of power grid planning and dispatching. Currently, relevant units in my country's power system operation, planning, design, and testing widely use the PSD power system analysis software developed by the China Electric Power Research Institute for transient stability calculations. However, in the analysis of annual power grid operation modes and safety and stability control, various maintenance methods and fault types need to be considered. This software lacks an automatic function for judging transient simulation results, leading to a large amount of repetitive manual operation during the analysis process and extremely low work efficiency.

[0003] Existing transient-stable batch computing technologies for the PSD-BPA software platform have the following shortcomings: First, batch calculations often adopt a "full scan" mode, which makes it difficult for users to flexibly customize fault subsets based on specific concerns such as new energy gathering areas, heavy load sections, and power plant transmission lines. This results in a large amount of computing resources being wasted on non-critical faults, making it impossible to conduct in-depth targeted investigations of specific weak links in the power grid, and leading to long analysis cycles and poor flexibility.

[0004] Second, traditional stability assessment methods rely on static criteria such as whether the maximum power angle difference at a single moment exceeds a threshold, which poses a serious risk of misjudgment and omission: for systems experiencing violent but controllable power swings under strong disturbances (i.e., large oscillations at the front end of the waveform, but the energy decay at the back end tends to be stable), it is very easy to misjudge as instability; for slow dynamic instability processes that are initially stable but gradually diverge in the later stage, the maximum power angle criterion in the short time window is very easy to miss, which cannot meet the requirements of high-precision power grid stability assessment. Summary of the Invention

[0005] Based on the above, the present invention aims to propose an automatic batch calculation and analysis method for transient stability of power systems based on PSD, so as to solve the problems of wasted computing resources and insufficient stability judgment accuracy in the prior art.

[0006] To achieve the above-mentioned technical objectives, the technical solution adopted by this invention is as follows: This invention proposes an automatic batch calculation and analysis method for transient stability of power systems based on PSD, comprising the following steps: S1. Obtain the line node information to be calculated, parse the PSD power flow data file, verify the existence and topological connection relationship of the target line node, distinguish the line type and complete the standard naming, and generate standardized line basic information. S2. Based on the fault parameters specified in the "Guidelines for the Safety and Stability of Power Systems", and combined with the voltage level and topology of the line, typical fault types are automatically generated; at the same time, user-defined fault constraint conditions are supported to construct custom fault subsets. S3. Based on the selected operating mode, the set of maintenance lines, and the set of faults, construct a transient simulation task queue corresponding to the Cartesian product of the set of maintenance lines and the set of faults, and execute the transient simulation tasks in batches; modify the power flow data file .dat and the transient stability file .swi in turn for the task queue to construct the power grid model under the corresponding maintenance mode and the fault cards that need to be verified; when an abnormal pop-up window of the PSD program is detected during the simulation, automatically restart and re-execute the current calculation case; S4. Monitor the calculation log of the PSD program in real time. When a calculation non-convergence anomaly is detected, trigger the convergence adjustment mechanism, automatically locate and fine-tune the active and reactive power values ​​of the specified partition in the power flow data file, write the modification back to the power flow file and call the PSD kernel to calculate again. Repeat the process until the calculation converges or the preset maximum number of adjustments is reached. Mark and log the final non-converged case. S5. Read the simulation result file and perform quantitative judgment on the system's power angle stability, voltage stability, frequency stability and dynamic stability respectively. For cases judged as unstable or critically stable in power angle, voltage and frequency, automatically redraw the corresponding characteristic curves and export them for storage. For cases judged as weakly damped or unstable in dynamic stability, prompt the user to pay special attention or prompt the user to enter the second level of Prony analysis. S6. Based on simulation and judgment results, generate a standardized analysis report containing multi-dimensional stability assessment data, and generate special detailed records for working conditions that do not meet stability requirements.

[0007] Furthermore, in step S1, the methods for obtaining line node information include batch import by the user through an Excel spreadsheet or manual entry. Based on the PSD power flow data file parsed from the line nodes, the system quickly retrieves the target line and verifies the existence of the target node and the line topology connection relationship. If the line does not exist, a prompt is issued to the user. If the line exists, the system automatically distinguishes the line type. The line type includes single-circuit lines and double-circuit lines erected on the same pole.

[0008] Furthermore, in step S2, the fault parameters specified in the "Guidelines for the Safety and Stability of Power Systems" include fault type, fault clearing time, and parameters of reclosing logic. User-defined fault constraints include fault location area, voltage level, fault nature, whether it involves important sections, and whether it involves power plant transmission lines.

[0009] Furthermore, in step S3, based on the selected operating mode, line nodes and generated fault types, users can specify the lines and fault types to be verified, or perform batch calculations on all line maintenance conditions and fault types to achieve flexible scheduling of targeted fault verification and full batch fault calculation. The verification lines include lines in areas with high penetration of new energy sources, heavy-load interconnection lines, and power plant transmission lines.

[0010] Furthermore, in step S4, the convergence adjustment mechanism alleviates the problems of insufficient local voltage support or excessive reactive circulating current by changing the reactive power distribution of the system, thereby achieving automatic repair of power flow ills.

[0011] Furthermore, in step S5, the determination of the stability of the work angle specifically includes: Extract the power angle curve time series data of the specified generator group, first filter out small and normal fluctuations by using the band-limited threshold zero-crossing counting method, and capture significant oscillations exceeding the preset amplitude; The power angle stability is initially verified by calculating the number of fluctuations within a set amplitude range. If the number of fluctuations exceeds the set threshold, dynamic stability judgment based on the time-domain energy decay characteristics is further performed. The dynamic stability assessment includes: identifying all local peaks and troughs in the power angle time series, constructing upper and lower envelopes respectively, and obtaining the envelope slope through least squares linear fitting; if the envelope slope is less than a threshold... β Determine if the system's power angle is stable; if the envelope slope is greater than the threshold. β Determine if the system is instable at the power angle; if it is in [- β，β The power angle curve is then exported between [ ] for manual verification.

[0012] Furthermore, in step S5, the determination of voltage stability specifically includes: Based on the per-unit value, monitor the recovery process of each bus voltage after fault clearance, locate the continuous intervals where the voltage is below 0.75pu, and calculate the duration of each over-limit interval through time series difference; if the over-limit duration of any bus voltage exceeds the preset maximum allowable duration, the voltage under this condition is determined to be unstable, otherwise it is determined to be stable.

[0013] Furthermore, in step S5, the determination of frequency stability specifically includes: The frequency-time curve is sampled according to the preset calculated cycle. If the maximum frequency obtained by sampling is greater than 51Hz, it is judged as transient high cycle and frequency instability. If the sampling frequency is within the range of 49Hz-51Hz, it is judged as frequency stability. If there is a continuous interval with sampling frequency below 49Hz and positioning frequency below 47Hz and the duration is calculated, if the duration exceeds the preset maximum allowable duration, it is judged as transient low cycle and frequency instability.

[0014] Furthermore, in step S5, the determination of dynamic stability specifically includes: Extract the system damping ratio from the simulation output file .OUT. If the damping ratio is ≥0.015, the system is considered dynamically stable. If the damping ratio is between 0.01 and 0.015, the system is considered weakly damped and marked as a key area of ​​concern. If the damping ratio is <0.01, a prompt is triggered and the user is guided to perform a second-level Prony analysis.

[0015] Furthermore, in step S6, the standardized analysis report includes multi-dimensional stability assessment data, which includes the comprehensive system stability assessment results of power angle stability, voltage stability, frequency stability, and dynamic stability.

[0016] By adopting the above technical solution, the present invention has the following beneficial effects compared with the prior art: 1. This invention realizes the standardized automatic generation of fault sets, strictly follows the "Guidelines for the Safety and Stability of Power Systems" to automatically generate fault cards, completely replaces the manual operation mode of writing each item one by one, greatly reduces labor costs, completely avoids omissions and errors in manual writing, and ensures the standardization and accuracy of fault sets.

[0017] 2. This invention supports the construction of custom fault subsets, enabling flexible scheduling of full scanning and targeted troubleshooting. Users can customize fault sets for key links such as new energy aggregation areas, heavy-load sections, and power plant transmission lines, significantly reducing invalid calculations for non-critical faults, saving computing resources, and shortening the analysis cycle for specific scenarios.

[0018] 3. This invention proposes a dynamic stability assessment method based on time-domain energy decay characteristics. It abandons the traditional static criterion of a single threshold and determines the system damping characteristics by the slope of the envelope of the power angle curve. It can accurately identify controllable power swing and slow dynamic instability processes under strong disturbances, greatly reduce the risk of misjudgment and omission, and significantly improve the accuracy of stability assessment.

[0019] 4. This invention has an adaptive convergence adjustment mechanism, which can automatically identify power flow non-convergence anomalies and complete parameter fine-tuning and recalculation, repair power flow ill-conditions under complex maintenance methods, ensure continuous and stable execution of batch calculations, and, combined with the automatic restart function of abnormal pop-up windows, realize unattended, fully automated calculation closed loop.

[0020] 5. The entire process of this invention can be achieved through a graphical interface for human-computer interaction, from importing and verifying line data, configuring faults, batch calculations to intelligent stability assessment and report generation, without the need for manual intervention. This significantly improves the efficiency of power grid transient stability analysis and provides accurate and efficient technical support for power grid planning, operation mode analysis and safety and stability control. Attached Figure Description

[0021] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the accompanying drawings used in the description of the embodiments will be briefly introduced below.

[0022] Figure 1 This is a simplified flowchart of the present invention.

[0023] Figure 2 Fault settings interface; Figure 3 This is the interface for transient batch calculations. Detailed Implementation

[0024] The present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be particularly noted that the following embodiments are for illustrative purposes only and do not limit the scope of the invention. Similarly, the following embodiments are only some, not all, embodiments of the present invention, and all other embodiments obtained by those skilled in the art without creative effort are within the scope of protection of the present invention.

[0025] This embodiment provides an automatic batch calculation and analysis method for transient stability of power systems based on PSD, implemented using the PSD-BPA power system analysis software platform. The overall process is as follows: Figure 1 As shown, the specific implementation steps are as follows: Step 1: Intelligent Analysis and Verification of Line Data Users can input the line node information to be calculated through a graphical interface, either by batch importing from an Excel spreadsheet or by manual entry. This includes basic data such as the names of the stations at both ends of the line, the number of circuits, and the voltage level. The program automatically parses the PSD power flow data file (.dat) based on the line node information, quickly retrieves the target line, and verifies the existence of the target node and the line topology connection relationship. If the target line node does not exist, a prompt message is immediately displayed to the user. If the line exists, it automatically distinguishes the line type, such as single-circuit line or double-circuit line erected on the same pole. Subsequently, it automatically identifies the line information and names the line in a standardized manner, providing standard line basic information for subsequent fault generation.

[0026] Step 2: Automatic generation of standardized fault scenarios The program incorporates all the specifications of the "Guidelines for the Safety and Stability of Power Systems." Based on parameters such as fault types, fault clearing times, and reclosing logic specified in the guidelines, and combined with information such as line voltage and topology obtained in step one, it automatically generates typical fault types as defined in the "Guidelines for the Safety and Stability of Power Systems" and completes the automatic generation of fault card codes. Simultaneously, the program provides a custom configuration entry point, allowing users to anticipate accidents based on other factors (such as fault location area, voltage level, fault nature, whether it involves important sections, and whether it involves power plant transmission lines), constructing a custom fault subset. The fault setting interface is as follows... Figure 2 As shown, this meets users' needs for targeted investigation of specific weak links in the power grid.

[0027] Step 3: Transient Batch Calculation Based on the selected operating mode, line nodes, and generated fault types, users can select the lines to be verified (such as high-penetration areas of new energy sources, heavy-load tie lines, and power plant transmission lines) and fault types (such as calculating only three-phase short circuits and single-phase failure to operate) in the graphical interface. Alternatively, users can choose to perform batch calculations for all line maintenance conditions and fault types. The transient batch calculation interface is shown below. Figure 3 As shown, the program generates a transient simulation task queue corresponding to the Cartesian product of the selected set of maintenance lines and fault sets, and performs batch execution of transient simulation tasks. For each case in the task queue, it automatically modifies the power flow data file (.dat) and the transient stability file (.swi) sequentially, completing the automatic construction of the power grid model and the fault cards to be verified under the corresponding maintenance mode. This achieves a semantic-level conversion from "graphical selection" to "machine-executable code," and then calls the PSD kernel to execute the transient simulation. During the simulation execution, the program monitors the PSD program's running status in real time. If an unknown PSD error pop-up window is encountered, the program automatically restarts the PSD program after a preset waiting time and re-executes the current case, ensuring the continuity of batch calculations.

[0028] Step 4: Adaptive Convergence Adjustment During the transient stability calculation of the simulated line maintenance mode, the program monitors the calculation log output of the PSD program in real time. When the PSD calculation program detects an abnormal state of "calculation non-convergence," the convergence adjustment mechanism is automatically triggered. Based on the user-predefined line partition, the program automatically locates the active and reactive power values ​​of the specified partition in the power flow data file and makes fine adjustments according to a preset step size to change the reactive power distribution of the system and alleviate the problems of insufficient local voltage support or excessive reactive power circulation current.

[0029] After parameter adjustments are completed, the program writes the modified data back to the power flow file and re-calls the PSD kernel for calculation. This adjustment process can be repeated until the calculation successfully converges or the maximum number of adjustments preset by the user is reached. If convergence fails after reaching the maximum number of adjustments, the program marks the calculation as "Power Flow Non-Convergence" and records a complete log on the program's main interface for subsequent manual review.

[0030] Step 5: Intelligent Stability Assessment and Post-processing After the simulation calculation is completed, the program automatically reads the simulation result file and performs quantitative intelligent judgment on the power grid operating conditions from four dimensions: power angle stability, voltage stability, frequency stability, and dynamic stability, as detailed below: 5.1 Determination of Stability of Work Angle After the calculation is completed, the simulation result file (.swx) is read. For the specified generator group, the power angle curve time series data is extracted. A band-limited threshold zero-crossing counting method is used, setting a non-zero amplitude band (upper limit - lower limit > 0) threshold to effectively filter small and normal fluctuations, and capture significant oscillations exceeding the preset amplitude. By calculating the number of fluctuations within the set amplitude range, the power angle stability is initially verified. If the number of fluctuations exceeds the set threshold, the power angle attenuation characteristics of the generator group under this operating condition are evaluated.

[0031] This invention abandons the traditional method of static threshold criterion for maximum power angle difference at a single moment, and proposes a dynamic stability judgment method based on time-domain energy decay characteristics. This method divides the simulation waveform into "transient impact period" and "late evolution period", and focuses on monitoring the envelope slope of the power angle swing amplitude during the late evolution period.

[0032] The dynamic stability assessment method includes: identifying all local peaks and troughs in the power angle time series, constructing upper and lower envelopes respectively, and obtaining the envelope slope through least squares linear fitting; if the envelope slope is less than 0, the system is considered stable in the power angle; if the envelope slope is greater than 0, the system is considered unstable in the power angle; if the envelope slope approaches 0, the system is considered to be in a critically stable state. The specific execution process is as follows: The program segments and identifies the waveform, first identifying the first swing peak to determine the end point of the transient impact. Then, it tracks the amplitude trends of subsequent peaks and calculates the ratio of the amplitudes of adjacent peaks and the slope of their envelopes.

[0033] Identifying extreme points: First, smooth the power angle time series data δ(t), then identify all local maxima (peaks) and local minima (troughs) from δ(t). Let the amplitude and time of the i-th peak be (δ_peak_i, t_peak_i), and the amplitude and time of the i-th trough be (δ_valley_i, t_valley_i).

[0034] Constructing the envelope: Upper envelope: Connects or fits all peak points (t_peak_i, δ_peak_i).

[0035] Lower envelope: Connect or fit all valley points (t_valley_i, δ_valley_i).

[0036] Calculate the slope: Perform a linear fit on the constructed envelope (e.g., using the least squares method), and the slope of the fitted line is the slope of the envelope.

[0037] Using linear regression statistics, assuming a linear fit on the first N peaks, the upper envelope can be expressed as δ_peak(t) ≈ k_upper * t + b_upper. The slope k_upper of the envelope can be calculated using the following least squares formula:

[0038] get

[0039] Similarly, the slope k_lower of the lower envelope can be calculated using the trough points.

[0040] k_upper<- β This indicates that the oscillation amplitude is decaying, the system has positive damping, and is stable.

[0041] k_upper> β This indicates that the oscillation amplitude has increased, the system has negative damping, and is unstable.

[0042] k_upper:[- β，β Export the power angle curve for manual verification.

[0043] For the above example of curve evaluation, redraw the waveform of the work angle curve and export it to the specified folder for manual review.

[0044] 5.2 Determination of Voltage Stability Using per-unit (pu) values ​​as a benchmark, the recovery process of each bus voltage after fault clearance is monitored to identify continuous intervals where the voltage is below 0.75 pu. For each continuous over-limit segment, the duration Δt = t_end - t_start is calculated using time series differential calculation, and the maximum allowable duration T_max = 1 second is set. If Δt > T_max exists for any bus, the voltage under that operating condition is determined to be unstable; otherwise, the voltage is determined to be stable.

[0045] 5.3 Determination of Frequency Stability Monitor the frequency-time curve under each working condition, sample the frequency-time values of the set calculation cycle, and make judgments according to the following rules: If the maximum sampled frequency f_max > 51 Hz, it is determined that the system is in a transient high-frequency state and the frequency is unstable; If all sampled frequencies fi satisfy 49 Hz < fi < 51 Hz, it is determined that the system frequency is stable; If there is a sampled frequency fi < 49 Hz, the program automatically locates the continuous interval where the frequency is lower than 47 Hz, and calculates the duration Δt = t_end - t_star through time series difference. Set the maximum allowable duration T_max = 0.5 s. If there is Δt > T_max, it is determined that the system is in a transient low-frequency state and the frequency is unstable; if not, it is determined that the frequency is stable. For the working conditions determined to have unstable frequencies, the program automatically redraws the frequency curve waveform diagram and exports it to the specified folder for manual review.

[0046] 5.4 Judgment of dynamic stability Extract the system damping ratio ζ from the simulation output file.OUT, and make judgments according to the following rules: If ζ ≥ 0.015, it is determined that the working condition is dynamically stable; if 0.01 ≤ ζ < 0.015, it is determined as weakly damped and marked as a working condition that needs to be focused on; if ζ < 0.01, a prompt is triggered to guide the user to manually enter the second-level Prony analysis to conduct in-depth dynamic stability assessment.

[0047] Step 6. Report export Based on the simulation calculation results and stability judgment results, the program automatically generates a multi-dimensional stability analysis report. The report adopts a standardized format, presenting the key index values of the system comprehensive stability assessment of power angle stability state, voltage stability, frequency stability, and dynamic stability. For any working condition with unqualified stability in any category, a special detailed record is automatically generated, providing complete data support for power grid safety and stability analysis and control measure formulation.

[0048] Furthermore, the functional units in the embodiments of the present invention can be integrated into one processing unit, or each unit can exist physically separately, or two or more units can be integrated into one unit. The integrated unit can be implemented in hardware or as a software functional unit. If the integrated unit is implemented as a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of the technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to cause a computer device (which may be a personal computer, server, or network device, etc.) or processor to execute all or part of the steps of the methods of the various embodiments of the present invention. The aforementioned storage medium includes various media capable of storing program code, such as USB flash drives, portable hard drives, read-only memory (ROM), random access memory (RAM), magnetic disks, or optical disks.

[0049] The above description is only a part of the embodiments of the present invention and does not limit the scope of protection of the present invention. Any equivalent device or equivalent process transformation made based on the content of the present invention specification and drawings, or direct or indirect application in other related technical fields, are similarly included within the patent protection scope of the present invention.

Claims

1. A method for automatic batch calculation and analysis of transient stability in power systems based on PSD, characterized in that, include: S1. Obtain the line node information to be calculated, parse the PSD power flow data file, verify the existence and topological connection relationship of the target line node, distinguish the line type and complete the standard naming, and generate standardized line basic information. S2. Based on the fault parameters specified in the "Guidelines for the Safety and Stability of Power Systems", and combined with the voltage level and topology of the line, typical fault types are automatically generated; at the same time, user-defined fault constraint conditions are supported to construct custom fault subsets. S3. Based on the selected operating mode, the set of maintenance lines and the set of faults, construct a transient simulation task queue corresponding to the Cartesian product of the set of maintenance lines and the set of faults, and perform batch execution of transient simulation tasks; modify the power flow data file .dat and the transient stability file .swi in turn for the task queue to construct the power grid model and the fault cards that need to be verified under the corresponding maintenance mode. S4. Monitor the calculation log of the PSD program in real time. When a calculation non-convergence anomaly is detected, trigger the convergence adjustment mechanism, automatically locate and fine-tune the active and reactive power values ​​of the specified partition in the power flow data file, write the modification back to the power flow file and call the PSD kernel to calculate again. Repeat the process until the calculation converges or the preset maximum number of adjustments is reached. Mark and log the final non-converged case. S5. Read the simulation result file and perform quantitative judgment on the system's power angle stability, voltage stability, frequency stability and dynamic stability respectively. For cases judged as unstable or critically stable in power angle, voltage and frequency, automatically redraw the corresponding characteristic curves and export them for storage. For cases judged as weakly damped or unstable in dynamic stability, prompt the user to pay special attention or prompt the user to enter the second level of Prony analysis. S6. Based on simulation and judgment results, generate a standardized analysis report containing multi-dimensional stability assessment data, and generate special detailed records for working conditions that do not meet stability requirements.

2. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S1, the methods for obtaining line node information include batch import by the user through an Excel spreadsheet or manual entry. Based on the PSD power flow data file parsed from the line nodes, the system quickly retrieves the target line and verifies the existence of the target node and the line topology connection relationship. If the line does not exist, a prompt is issued to the user. If the line exists, the system automatically distinguishes the line type. The line type includes single-circuit lines and double-circuit lines erected on the same pole.

3. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S2, the fault parameters specified in the "Guidelines for the Safety and Stability of Power Systems" include fault type, fault clearing time, and parameters of reclosing logic; User-defined fault constraints include fault location area, voltage level, fault nature, whether it involves important sections, and whether it involves power plant transmission lines.

4. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S3, based on the selected operating mode, line nodes and generated fault types, users can specify the lines and fault types to be verified, or perform batch calculations on all line maintenance conditions and fault types to achieve flexible scheduling of targeted fault verification and full batch fault calculation. The verification lines include lines in areas with high penetration of new energy sources, heavy-load interconnection lines, and power plant transmission lines.

5. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S4, the convergence adjustment mechanism alleviates the problems of insufficient local voltage support or excessive reactive circulating current by changing the reactive power distribution of the system.

6. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S5, the determination of the stability of the work angle specifically includes: Extract the power angle curve time series data of the specified generator group, first filter out small and normal fluctuations by using the band-limited threshold zero-crossing counting method, and capture significant oscillations exceeding the preset amplitude; The power angle stability is initially verified by calculating the number of fluctuations within a set amplitude range. If the number of fluctuations exceeds the set threshold, dynamic stability judgment based on the time-domain energy decay characteristics is further performed. The dynamic stability assessment includes: identifying all local peaks and troughs in the power angle time series, constructing upper and lower envelopes respectively, and obtaining the envelope slope through least squares linear fitting; if the envelope slope is less than a threshold... β Determine if the system's power angle is stable; if the envelope slope is greater than the threshold. β Determine if the system is instable at the power angle; if it is in [- β，β The power angle curve is then exported between [ ] for manual verification.

7. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S5, the determination of voltage stability specifically includes: Based on the per-unit value, monitor the recovery process of each bus voltage after fault clearance, locate the continuous intervals where the voltage is below 0.75pu, and calculate the duration of each over-limit interval through time series difference; if the over-limit duration of any bus voltage exceeds the preset maximum allowable duration, the voltage under this condition is determined to be unstable, otherwise it is determined to be stable.

8. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S5, the determination of frequency stability specifically includes: The frequency-time curve is sampled according to the preset calculated cycle. If the maximum frequency obtained by sampling is greater than 51Hz, it is judged as transient high cycle and frequency instability. If the sampling frequency is within the range of 49Hz-51Hz, it is judged as frequency stability. If there is a continuous interval with sampling frequency below 49Hz and positioning frequency below 47Hz and the duration is calculated, if the duration exceeds the preset maximum allowable duration, it is judged as transient low cycle and frequency instability.

9. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S5, the determination of dynamic stability specifically includes: Extract the system damping ratio from the simulation output file .OUT. If the damping ratio is ≥0.015, the system is considered dynamically stable. If the damping ratio is between 0.01 and 0.015, the system is considered weakly damped and marked as a key area of ​​concern. If the damping ratio is <0.01, a prompt is triggered and the user is guided to perform a second-level Prony analysis.

10. The automatic batch calculation and analysis method for transient stability of power systems based on PSD according to claim 1, characterized in that, In step S6, the standardized analysis report includes multi-dimensional stability assessment data, which includes the comprehensive system stability assessment results of power angle steady state, voltage stability, frequency stability, and dynamic stability.

Images