Stability evaluation method, device and equipment of power system, storage medium and program product

Abstract

The invention relates to a stability evaluation method and device of a power system, equipment, a storage medium and a program product. The method comprises the following steps: constructing a power system model carrying slow-varying disturbance according to a plurality of related parameters of a synchronous generator in the power system, analyzing the power system model, determining a target boundary of the power system model, evaluating the stability of the power system according to the target boundary, and determining the stability of the power system. The method does not depend on a specific balance point, and only needs to perform dynamic characteristic analysis on the power system model to determine that the disturbance change rate of the power system model is within a certain range and determine the corresponding boundary, so that the stability of the power system is simply and quickly determined according to the boundary.

Description

Technical Field

[0001] The present application relates to the technical field of power systems, and in particular to a method, device, equipment, storage medium and program product for evaluating the stability of a power system. Background Art

[0002] In traditional power systems dominated by thermal power units, most disturbances that have a significant impact on system stability are instantaneous, such as load switching, instantaneous faults, and line disconnection. However, with the access of a high proportion of new energy to the power system, disturbances that have a significant impact on power system stability include not only instantaneous disturbances, but also fast-changing continuous disturbances and slow-changing continuous disturbances.

[0003] Compared with instantaneous disturbances, when there are slow-changing continuous disturbances, the balance point in the power system has enough response time to change with the continuous disturbance. Therefore, the balance point of the power system may deviate from the original balance point. Then, if the power system stability analysis is still based on the original balance point, the error of the power system stability analysis will increase, which will have a non-negligible impact on the system stability.

[0004] Currently, there is no method for evaluating the stability of power systems with slow-changing and persistent disturbances. Summary of the invention

[0005] Based on this, it is necessary to provide a method, device, equipment, storage medium and program product that can realize stability assessment of power systems with slow-changing disturbances in response to the above technical problems.

[0006] In a first aspect, the present application provides a method for evaluating the stability of a power system, comprising:

[0007] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0008] Analyze the power system model and determine the target boundary of the power system model;

[0009] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0010] In one embodiment, analyzing the power system model to determine the target boundary of the power system model includes:

[0011] Based on the linear matrix inequality relationship and the power system model, determining a first preset matrix of the power system within the operation domain;

[0012] Determine the continuity constants and related parameters within the operation domain based on the preset continuity conditions and power system model;

[0013] According to any two data points in the power system model, obtaining a second preset matrix corresponding to the power system model;

[0014] A target boundary of the power system model is determined according to the first preset matrix, the continuous constant, the relevant parameters in the operation domain and the second preset matrix.

[0015] In one embodiment, determining the target boundary of the power system model according to the first preset matrix, the continuous constant, the maximum value in the operating domain, and the second preset matrix includes:

[0016] Obtaining a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of a first preset matrix;

[0017] Obtaining a second matrix minimum eigenvalue of a second preset matrix;

[0018] The target boundary of the power system model is determined according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value in the operation domain and the minimum eigenvalue of the second matrix.

[0019] In one embodiment, the above-mentioned determination of the continuous constant and related parameters in the operation domain based on the preset continuous condition and the power system model includes:

[0020] If the power system model meets the preset continuity conditions, then the relevant parameters in the operation domain are determined according to the power system model;

[0021] The continuous constant is determined according to the relevant parameters in the operation domain and the relevant parameters of each state variable in the operation domain.

[0022] In one embodiment, the step of obtaining a second preset matrix corresponding to the power system model according to any two data points in the power system model includes:

[0023] Obtain the difference between the power system models corresponding to the two data points;

[0024] According to the difference between the power system models corresponding to the two data points, a second preset matrix corresponding to the power system model is determined.

[0025] In one embodiment, the method further comprises:

[0026] Performing coordinate transformation on the power system model to obtain the power system model after coordinate transformation;

[0027] Analyze the power system model and determine the target boundaries of the power system model, including:

[0028] Analyze the power system model after coordinate transformation and determine the target boundary of the power system model.

[0029] In one embodiment, the method further comprises:

[0030] Applying preset control to the power system model to obtain the power system model after the preset control is applied; the preset control includes integral proportional control and negative feedback control;

[0031] Analyze the power system model and determine the target boundaries of the power system model, including:

[0032] The power system model after applying preset control is analyzed to determine the target boundary of the power system model.

[0033] In a second aspect, the present application also provides a stability assessment device for a power system, comprising:

[0034] A construction module is used to construct a power system model carrying slow-varying disturbances according to multiple relevant parameters of the synchronous generator in the power system; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0035] An analysis module, used to analyze the power system model and determine the target boundary of the power system model;

[0036] The evaluation module is used to evaluate the stability of the power system according to the target boundary and determine the stability of the power system.

[0037] In a third aspect, the present application further provides a computer device, including a memory and a processor, wherein the memory stores a computer program, and when the processor executes the computer program, the following steps are implemented:

[0038] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0039] Analyze the power system model and determine the target boundary of the power system model;

[0040] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0041] In a fourth aspect, the present application further provides a computer-readable storage medium having a computer program stored thereon, wherein when the computer program is executed by a processor, the following steps are implemented:

[0042] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0043] Analyze the power system model and determine the target boundary of the power system model;

[0044] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0045] In a fifth aspect, the present application further provides a computer program product, including a computer program, which implements the following steps when executed by a processor:

[0046] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0047] Analyze the power system model and determine the target boundary of the power system model;

[0048] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0049] The above-mentioned power system stability assessment method, device, equipment, storage medium and program product, the above-mentioned method provides a stability assessment method for a power system with slow-changing disturbances. A major difficulty that slow-changing disturbances bring to stability analysis is that the balance point will continue to change with the disturbance, and the analysis of a given balance point is difficult to apply. The above-mentioned method does not rely on a specific balance point, and only needs to perform a dynamic characteristic analysis on the power system model to determine that the disturbance change rate of the power system model is within a certain range, and determine the corresponding boundary, so as to simply and quickly determine the stability of the power system based on the boundary. BRIEF DESCRIPTION OF THE DRAWINGS

[0050] In order to more clearly illustrate the technical solutions in the embodiments of the present application or related technologies, the drawings required for use in the embodiments of the present application or related technical descriptions will be briefly introduced below. Obviously, the drawings described below are only some embodiments of the present application. For ordinary technicians in this field, other related drawings can be obtained based on these drawings without paying creative work.

[0051] FIG1 is an application environment diagram of a method for evaluating the stability of a power system in one embodiment;

[0052] Figure 2 is a flow chart of a method for evaluating the stability of a power system in one embodiment;

[0053] Figure 3 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0054] Figure 4 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0055] Figure 5 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0056] Figure 6 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0057] Figure 7 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0058] Figure 8 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0059] Fig. 9 is a flow chart of a method for evaluating the stability of a power system in another embodiment;

[0060] FIG10 is a structural block diagram of a device for evaluating the stability of a power system in one embodiment;

[0061] FIG. 11 is a diagram showing the internal structure of a computer device in one embodiment. DETAILED DESCRIPTION

[0062] In order to make the purpose, technical solution and advantages of the present application more clearly understood, the present application is further described in detail below in conjunction with the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present application and are not used to limit the present application.

[0063] In traditional power systems dominated by thermal power units, most of the disturbances that have a significant impact on system stability are instantaneous, such as load switching, instantaneous faults, line removal, etc. The power system is easily deviated from the equilibrium state and enters the non-equilibrium state when affected by disturbances. The stability analysis of the power system mainly focuses on whether the power system can return to the original equilibrium state or reach a new steady state after being disturbed.

[0064] Generally speaking, the traditional power system can be regarded as an autonomous system. With the access of a high proportion of renewable energy to the power system, the power system will face continuous disturbances caused by the volatility of renewable energy, and the operating conditions and balance points will continue to change. In addition, continuous power fluctuations will bring fast-changing continuous disturbances and slow-changing continuous disturbances to the power system, which is essentially different from the impact of instantaneous disturbances and also brings new challenges to the stability analysis of the power system.

[0065] The output changes of new energy sources such as wind power can be decomposed into slow, fast and sudden (ramp) components. The safe and stable operation of the power system requires the real-time balance between power generation and demand. Therefore, under the background of high proportion of new energy access, this balance will be affected by events in three time frames: fast-changing continuous disturbances, slow-changing continuous disturbances and instantaneous sudden disturbances. Compared with fast-changing continuous disturbances, under slow-changing continuous disturbances, the system balance point has enough response time to change with the continuous disturbance, and may deviate from the original balance point, which increases the error of the stability analysis around the original balance point and has a non-negligible impact on the stability of the system.

[0066] However, currently, there is no method for evaluating the stability of a power system with slow-changing and continuous disturbances. This application aims to solve this problem.

[0067] After the above introduction to the background technology of a method for evaluating the stability of an electric power system provided by an embodiment of the present application, the implementation environment involved in the method for evaluating the stability of an electric power system provided by an embodiment of the present application will be briefly described below. The method for evaluating the stability of an electric power system provided by an embodiment of the present application can be applied to the implementation environment shown in Figure 1, which includes a server 104, and the server 104 can be implemented with an independent server 104 or a server cluster composed of multiple servers 104, and the data storage system 102 can store the data that the server 104 needs to process. The data storage system 102 can be integrated on the server 104, or it can be placed on a cloud or other network server. Among them, the server 104 can obtain multiple relevant parameters of the electric power system, and then, according to multiple relevant parameters of the synchronous generator in the electric power system, a model of the electric power system carrying slow-changing disturbances is constructed, and then, the electric power system model is analyzed to determine the target boundary of the electric power system model, and finally, the stability of the electric power system is evaluated according to the target boundary of the electric power system to determine the stability of the electric power system.

[0068] In other possible implementations, the power system stability assessment method provided in the embodiments of the present application may also be applied to terminals, which may be but are not limited to various personal computers, laptops, smart phones, tablet computers, Internet of Things devices, and portable wearable devices.

[0069] After introducing the application scenarios of the power system stability assessment method provided in the embodiments of the present application, the following focuses on introducing the power system stability assessment method described in the present application.

[0070] In one embodiment, Figure 2 As shown, a method for evaluating the stability of a power system is provided, and the method is applied to the server in FIG1 as an example for explanation, and includes the following steps:

[0071] S201. Construct a power system model carrying slow-varying disturbances according to multiple relevant parameters of synchronous generators in the power system.

[0072] Among them, the relevant parameters include the d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator.

[0073] In the embodiment of the present application, multiple relevant parameters of the synchronous generator in the power system can be obtained first, and a power system model carrying slow-varying disturbances can be constructed based on the multiple relevant parameters of the synchronous generator and the basic model of the internal structure of the generator. Optionally, the power system model carrying slow-varying disturbances can be expressed by the following formula (1):

[0074]

[0075] Among them, is the generator voltage phase angle, is the rotor angular frequency, is the potential amplitude in the generator q axis, M is the inertia constant, D is the damping coefficient, is the q-axis open-loop transient time constant, is the d-axis reactance, is the d-axis transient reactance, is the line reactance, Simplify the total circuit reactance for the system, and are the mechanical power and excitation voltage of the synchronous generator, respectively. , They are disturbances related to active and reactive output of new energy respectively.

[0076] S202: Analyze the power system model to determine the target boundary of the power system model.

[0077] Among them, the target boundary is used to characterize the stability of the power system model. The target boundary refers to the dynamic final consistency bound of the power system model.

[0078] In the embodiment of the present application, after the power system model carrying the slow-varying disturbance is constructed as described above, the dynamic characteristics of the power system model may be analyzed to determine the target boundary of the power system model.

[0079] S203. Evaluate the stability of the power system according to the target boundary to determine the stability of the power system.

[0080] In an embodiment of the present application, after the target boundary of the power system model is determined as described above, the stability of the power system can be evaluated according to the target boundary to determine the stability of the power system. Optionally, a boundary threshold can be predetermined. If the target boundary is greater than the boundary threshold, it means that the power system model is unstable. If the target boundary is not greater than the boundary threshold, it means that the power system model is stable. Alternatively, the higher the target boundary, the looser the requirements of the power system stability on the nature of the disturbance, which also means that under the same disturbance, the smaller the final consistent bound of the system, that is, the stronger the slow-changing stability of the system itself.

[0081] The stability assessment method of the power system provided in the embodiment of the present application constructs a power system model carrying slow-varying disturbances based on multiple relevant parameters of the synchronous generator in the power system, analyzes the power system model, determines the target boundary of the power system model, assesses the stability of the power system based on the target boundary, and determines the stability of the power system, wherein the relevant parameters include the d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator. The above method provides a stability assessment method for a power system with slow-varying disturbances. A major difficulty that slow-varying disturbances bring to stability analysis is that the balance point will change continuously with the disturbance, and the analysis of a given balance point is difficult to apply. The above method does not rely on a specific balance point, and only needs to perform a dynamic characteristic analysis on the power system model to determine that the disturbance change rate of the power system model is within a certain range, and determine the corresponding boundary, thereby simply and quickly determining the stability of the power system based on the boundary.

[0082] In one embodiment, Figure 2 Based on the embodiment shown, the process of determining the target boundary of the power system model can be described as follows: Figure 3 As shown, the above S202 “analyzing the power system model to determine the target boundary of the power system model” includes:

[0083] S301. Determine a first preset matrix of the power system in an operation domain based on a linear matrix inequality relationship and a power system model.

[0084] The linear matrix inequality relation refers to the linear matrix inequality (LMI).

[0085] The operation domain of the power system refers to , where r is the radius of the system operating domain. The operating domain of the power system is pre-set. For example, the voltage of the power system can vary by plus or minus 30% of the threshold voltage.

[0086] Among them, the first preset matrix is ​​a positive definite real symmetric matrix.

[0087] In the embodiment of the present application, a linear matrix inequality (LMI) and a power system model can be used to solve the problem and obtain a first preset matrix of the power system in the operation domain. , solve the following inequality (2) to obtain the first preset matrix P:

[0088]

[0089] in, , That is The Jacobian matrix of The preset balance point ,have , z is the translation distance of the coordinate translation of formula (1), , I is the identity matrix.

[0090] Alternatively, the power system model (1) can be expressed in a compact form as shown in the following formula (3):

[0091]

[0092] The state variable x can be expressed by the following formula (4):

[0093]

[0094] Among them, the function It can be expressed by the following formula (5):

[0095]

[0096] in, , They correspond to the three nonlinear equations in formula (1).

[0097] in, It can be expressed by the following formula (6):

[0098]

[0099] S302: Determine the continuity constant and related parameters within the operation domain based on the preset continuity condition and the power system model.

[0100] The preset continuity condition refers to the Lipschitz continuity condition.

[0101] The relevant parameters in the operation domain are The Jacobian matrix The maximum value of .

[0102] In an embodiment of the present application, when the power system model meets the preset continuity conditions, the relevant parameters in the operating domain can be determined according to the power system model, and then the continuous constants can be determined according to the relevant parameters in the operating domain and the relevant parameters of each state variable in the operating domain.

[0103] Optionally, a method for determining a continuous constant and a maximum value within the operating domain is provided below, such as Figure 4 As shown, the above S302 "determine the continuous constant and related parameters in the operation domain based on the preset continuous conditions and the power system model" includes:

[0104] S3021. If the power system model meets the preset continuity conditions, determine the relevant parameters within the operation domain according to the power system model.

[0105] In the embodiment of the present application, when the above-mentioned power system model meets the preset continuity condition, it is determined The Jacobian matrix The maximum value of The Jacobian matrix The maximum value of is taken as the relevant parameter in the operating domain.

[0106] S3022. Determine a continuous constant according to a maximum value in the operation domain and relevant parameters of each state variable in the operation domain.

[0107] The relevant parameters of each state variable in the operation domain refer to the maximum absolute value of each state variable in the operation domain. .

[0108] In the embodiment of the present application, after the maximum value in the operating domain is determined as above, the relevant parameters of each state variable in the operating domain can be further obtained, and the continuous constant can be determined based on the maximum value in the operating domain and the relevant parameters of each state variable in the operating domain.

[0109] Optionally, if the power system model satisfies the Lipschitz condition in the operating domain, its Lipschitz constant is denoted as L, then exist In the case of Satisfies the following formula (7):

[0110]

[0111] Then for the power system model shown in formula (1), its Lipschitz constant can be calculated by the following formula (8):

[0112] (8);

[0114] in, is the maximum absolute value of each state variable in the operating domain, is the integral control parameter, refers to the proportional control parameter, It refers to the negative feedback control parameter.

[0115] It should be noted that the Lipschitz constant is used to characterize the continuous characteristics of the power system dynamics.

[0116] S303: Obtain a second preset matrix corresponding to the power system model according to any two data points in the power system model.

[0117] In an embodiment of the present application, two data points can be arbitrarily selected within the operating domain of the power system model, and then the values ​​corresponding to the two data points in the power system model can be determined, and then the difference between the values ​​corresponding to the two data points in the power system model can be obtained, and based on the difference between the values ​​corresponding to the two data points in the power system model, the second preset matrix corresponding to the power system model can be determined.

[0118] Optionally, a method for determining a second preset matrix corresponding to a power system model is provided below, such as Figure 5 As shown, the above S303 "obtaining a second preset matrix corresponding to the power system model according to any two data points in the power system model" includes:

[0119] S3031. Obtain the difference between the power system models corresponding to the two data points.

[0120] In the embodiment of the present application, two data points x and x can be randomly selected in the operation domain of the power system model. , these two data points x and Input them into the power system model respectively to obtain the corresponding value of data point x in the power system model , and at the same time, determine the data points The corresponding value in the power system model , then determine the value corresponding to the data point x in the power system model With data points The corresponding value in the power system model The difference between .

[0121] S3032. Determine a second preset matrix corresponding to the power system model according to the difference between the power system models corresponding to the two data points.

[0122] The second preset matrix is ​​a positive definite real symmetric matrix.

[0123] In the embodiment of the present application, the value corresponding to the data point x in the power system model is obtained as follows: With data points The corresponding value in the power system model The difference between Afterwards, you can Taylor expansion is performed to obtain the following formula (9):

[0124]

[0125] in, is a positive definite real symmetric matrix.

[0126] S304: Determine a target boundary of the power system model according to the first preset matrix, the continuous constant, relevant parameters in the operation domain, and the second preset matrix.

[0127] In an embodiment of the present application, after the first preset matrix, the continuous constant, the relevant parameters within the operating domain and the second preset matrix are obtained as mentioned above, the maximum eigenvalue of the first matrix in the first preset matrix and the maximum eigenvalue of the second matrix in the first preset matrix can be obtained first, and then the minimum eigenvalue of the second matrix in the second preset matrix can be obtained. Then, according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value within the operating domain and the minimum eigenvalue of the second matrix, the target boundary of the power system model is determined.

[0128] Optionally, a method for determining the target boundary of the power system model is provided below, such as Figure 6 As shown, the above S304 "determine the target boundary of the power system model according to the first preset matrix, the continuous constant, the relevant parameters in the operation domain and the second preset matrix" includes:

[0129] S3041. Obtain a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of a first preset matrix.

[0130] In the embodiment of the present application, after the first preset matrix P is obtained as described above, the first maximum eigenvalue of the first preset matrix is ​​obtained. , and obtain the first matrix minimum eigenvalue of the first preset matrix .

[0131] S3042. Obtain a second minimum eigenvalue of a second preset matrix.

[0132] In the embodiment of the present application, after the second preset matrix S is obtained as described above, the second minimum eigenvalue of the second preset matrix in the operation domain D is obtained. .

[0133] S3043. Determine the target boundary of the power system model according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value within the operating domain, and the minimum eigenvalue of the second matrix.

[0134] In the embodiment of the present application, after obtaining the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value in the operating domain, and the minimum eigenvalue of the second matrix as described above, the first parameter can be determined according to the product of the minimum eigenvalue of the first matrix and the minimum eigenvalue of the second matrix, and the second parameter can be determined according to the maximum eigenvalue of the first matrix and the continuous constant, and the minimum eigenvalue of the second matrix can be determined as the third parameter, and the fourth parameter can be determined according to the maximum value in the operating domain, the maximum eigenvalue of the second matrix, and the continuous constant, and the fifth parameter can be determined according to the maximum eigenvalue of the first matrix and the continuous constant. Then, the target boundary of the power system model is determined according to the first parameter, the second parameter, the third parameter, the fourth parameter, and the fifth parameter.

[0135] Below, the determination method of the first parameter, the second parameter, the third parameter, the fourth parameter and the fifth parameter is provided as shown in formula (10):

[0136]

[0137] in, , , , and They refer to the first parameter, the second parameter, the third parameter, the fourth parameter and the fifth parameter respectively. is the maximum eigenvalue of the first matrix, refers to the minimum eigenvalue of the second matrix, L refers to a continuous constant, is the maximum eigenvalue of the first matrix, It refers to the maximum value within the operating domain.

[0138] Optionally, after the first parameter, the second parameter, the third parameter, the fourth parameter and the fifth parameter are determined, the initial state of the power system model is The following formula (11) is satisfied:

[0139]

[0140] For the disturbance change rate of the power system model The following formula (12) is satisfied:

[0141]

[0142] Then, according to the above formula (10), the power system model is is uniformly bounded, and its final uniform bound can be expressed by the following formula (13):

[0143]

[0144] Where b is the final consistent bound of the power system model, is the disturbance parameter of the power system model.

[0145] Optionally, the perturbation parameter in the above final consistent bound can be differentiated to obtain the target boundary of the power system model shown in the following formula (14):

[0146]

[0147] in, is the target boundary of the power system model.

[0148] From the above formula (14), we can see that the target boundary The target boundary is proportional to the requirement of the disturbance change rate (12). The higher it is, the looser the requirements of the system stability on the nature of the disturbance. It also means that under the same disturbance, the final consistent bound of the system is smaller, that is, the slow-changing stability of the system itself is stronger.

[0149] The method for determining the target boundary of the power system model provided in the embodiment of the present application does not rely on the specific balance point of the power system model, but abstracts the concept of the system balance point. It is not necessary to specifically obtain the balance point information of the power system model, but only requires that the balance point itself and its rate of change with disturbance are within a certain range. Ultimately, the distance between the actual trajectory of the power system model and the balance point, that is, the target boundary, can be evaluated, and the stability of the power system model can be evaluated based on the target boundary.

[0150] In one embodiment, Figure 2-Figure 6 Based on the embodiment shown, Figure 7 As shown, the above method also includes:

[0151] S204: performing coordinate transformation on the power system model to obtain the power system model after coordinate transformation.

[0152] The coordinate transformation may be a coordinate translation transformation.

[0153] In the embodiment of the present application, after the power system model is obtained as described above, the power system model can be subjected to coordinate translation transformation. Optionally, the power system model can be translated horizontally. , the equilibrium point of the power system model is translated to the origin to obtain the translated slow-varying system as shown in the following formula (15):

[0154]

[0155] Among them, the state variable Indicates the distance between the power system state and the balance point.

[0156] The above S202 “analyzing the power system model and determining the target boundary of the power system model” includes:

[0157] S202: Analyze the power system model after the coordinate transformation to determine the target boundary of the power system model.

[0158] In the application embodiment, after the power system model carrying slow-varying disturbances is constructed as described above, the power system model can be subjected to coordinate transformation to obtain the power system model after coordinate transformation, and then analysis can be performed based on the power system model after coordinate transformation to determine the target boundary of the power system model.

[0159] The coordinate transformation method for the power system model provided in the embodiment of the present application provides a data basis for subsequent analysis of the power system model after the coordinate transformation to obtain the target boundary of the power system model.

[0160] In one embodiment, Figure 2-Figure 6 Based on the embodiment shown, Figure 8 As shown, the above method also includes:

[0161] S205 , applying preset control to the power system model to obtain the power system model after applying the preset control.

[0162] Among them, the preset control includes integral proportional control and negative feedback control.

[0163] In the embodiment of the present application, after the power system model is obtained as described above, the mechanical power of the synchronous generator in formula (1) can be calculated as Adopt integral proportional control (PI control), and the excitation voltage of the synchronous generator Adopt negative feedback control to obtain the mechanical power of PI control , and the excitation voltage applied with negative feedback control As shown in the following formula (16):

[0164]

[0165] in, and To control the set value, is the initial amplitude of the potential in the generator q axis.

[0166] The above S202 “analyzing the power system model and determining the target boundary of the power system model” includes:

[0167] S202: Analyze the power system model after applying preset control to determine the target boundary of the power system model.

[0168] In the embodiment of the application, after the power system model carrying slow-varying disturbances is constructed as described above, the mechanical power in the power system model can be and excitation voltage Apply preset control respectively to obtain the mechanical power after applying preset control and excitation voltage , and then analyze the power system model after applying the preset control to determine the target boundary of the power system model.

[0169] Optionally, you can also first calculate the mechanical power in the power system model and excitation voltage Apply preset control respectively to obtain the mechanical power after applying preset control and excitation voltage , and then based on the mechanical power after applying the preset control and excitation voltage A coordinate transformation is performed to obtain a power system model with the coordinate transformation and preset control applied, and then an analysis is performed based on the power system model with the coordinate transformation and preset control applied to determine a target boundary of the power system model.

[0170] The method for applying control to the power system model provided in the embodiment of the present application enhances the stability of the system and provides a data basis for subsequent analysis of the power system model after coordinate transformation to obtain the target boundary of the power system model.

[0171] In one embodiment, Fig. 9 As shown, a stability assessment method for a power system is also provided, including:

[0172] S10, constructing a power system model carrying slow-varying disturbances according to multiple relevant parameters of synchronous generators in the power system;

[0173] S11, applying a preset control to the power system model to obtain the power system model after the preset control is applied;

[0174] S12, performing coordinate transformation on the power system model after applying preset control, to obtain the power system model after coordinate transformation;

[0175] S13, determining a first preset matrix of the power system within the operation domain based on the linear matrix inequality relationship and the power system model;

[0176] S14. If the power system model satisfies the preset continuity condition, then determine the relevant parameters within the operation domain according to the power system model;

[0177] S15, determining a continuous constant according to relevant parameters in the operation domain and relevant parameters of each state variable in the operation domain;

[0178] S16, obtaining the difference between the power system models corresponding to the two data points;

[0179] S17, determining a second preset matrix corresponding to the power system model according to the difference between the power system models corresponding to the two data points;

[0180] S18, obtaining a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of a first preset matrix;

[0181] S19, obtaining a second minimum eigenvalue of a second preset matrix;

[0182] S20, determining a target boundary of the power system model according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value and the second minimum eigenvalue in the operation domain;

[0183] S21. Evaluate the stability of the power system according to the target boundary and determine the stability of the power system.

[0184] The above method provides a stability assessment method for a power system with slow-varying disturbances. A major difficulty that slow-varying disturbances bring to stability analysis is that the balance point will change continuously with the disturbance, and the analysis of a given balance point is difficult to apply. The above method does not rely on a specific balance point. It only needs to perform a dynamic characteristic analysis on the power system model to determine that the disturbance change rate of the power system model is within a certain range and determine the corresponding boundaries, so as to simply and quickly determine the stability of the power system based on the boundaries.

[0185] It should be understood that, although the various steps in the flowcharts involved in the above-mentioned embodiments are displayed in sequence according to the indication of the arrows, these steps are not necessarily executed in sequence according to the order indicated by the arrows. Unless there is a clear explanation in this article, the execution of these steps does not have a strict order restriction, and these steps can be executed in other orders. Moreover, at least a part of the steps in the flowcharts involved in the above-mentioned embodiments can include multiple steps or multiple stages, and these steps or stages are not necessarily executed at the same time, but can be executed at different times, and the execution order of these steps or stages is not necessarily to be carried out in sequence, but can be executed in turn or alternately with other steps or at least a part of the steps or stages in other steps.

[0186] Based on the same inventive concept, the embodiment of the present application also provides a stability assessment device for an electric power system for implementing the stability assessment method for an electric power system involved above. The implementation scheme for solving the problem provided by the device is similar to the implementation scheme recorded in the above method, so the specific limitations in the embodiments of one or more stability assessment devices for electric power systems provided below can refer to the limitations of the stability assessment method for electric power systems above, and will not be repeated here.

[0187] In an exemplary embodiment, as shown in FIG10 , a stability assessment device for a power system is provided, comprising: a construction module 10, a first determination module 11 and a second determination module 12, wherein:

[0188] The construction module 10 is used to construct a power system model carrying slow-varying disturbances based on multiple relevant parameters of the synchronous generator in the power system; the relevant parameters include the d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator.

[0189] The first determination module 11 is used to analyze the power system model and determine the target boundary of the power system model.

[0190] The second determination module 12 is used to evaluate the stability of the power system according to the target boundary and determine the stability of the power system.

[0191] In an exemplary embodiment, the first determining module 11 includes: a first determining unit, a second determining unit, an acquiring unit and a third determining unit, wherein:

[0192] A first determining unit, specifically configured to determine a first preset matrix of the power system within an operation domain based on a linear matrix inequality relationship and a power system model;

[0193] A second determination unit, specifically configured to determine the continuous constant and related parameters within the operation domain based on a preset continuous condition and a power system model;

[0194] An acquisition unit, specifically configured to acquire a second preset matrix corresponding to the power system model according to any two data points in the power system model;

[0195] The third determination unit is specifically used to determine the target boundary of the power system model according to the first preset matrix, the continuous constant, the relevant parameters in the operation domain and the second preset matrix.

[0196] In an exemplary embodiment, the above-mentioned third determination unit is specifically used to obtain the first matrix maximum eigenvalue and the first matrix minimum eigenvalue of the first preset matrix; obtain the second matrix minimum eigenvalue of the second preset matrix; and determine the target boundary of the power system model based on the first matrix maximum eigenvalue, the first matrix minimum eigenvalue, the continuous constant, the maximum value within the operating domain, and the second matrix minimum eigenvalue.

[0197] In an exemplary embodiment, the above-mentioned second determination unit is specifically used to determine relevant parameters within the operating domain according to the power system model if the power system model meets the preset continuity conditions; and determine the continuous constants according to the relevant parameters within the operating domain and the relevant parameters of each state variable within the operating domain.

[0198] In an exemplary embodiment, the acquisition unit is further configured to acquire a difference between power system models corresponding to two data points; and determine a second preset matrix corresponding to the power system model according to the difference between the power system models corresponding to the two data points.

[0199] In an exemplary embodiment, the above-mentioned device further includes: a transformation module, which is used to perform coordinate transformation on the power system model to obtain the power system model after coordinate transformation;

[0200] The first determination module is used to analyze the power system model after the coordinate transformation and determine the target boundary of the power system model.

[0201] In an exemplary embodiment, the above-mentioned device further includes: a control module, which is used to apply a preset control to the power system model to obtain the power system model after the preset control is applied; the preset control includes integral proportional control and negative feedback control;

[0202] The first determination module is used to analyze the power system model after the preset control is applied, and determine the target boundary of the power system model.

[0203] Each module in the above-mentioned power system stability assessment device can be implemented in whole or in part by software, hardware and a combination thereof. Each of the above-mentioned modules can be embedded in or independent of a processor in a computer device in the form of hardware, or can be stored in a memory in a computer device in the form of software, so that the processor can call and execute the operations corresponding to each of the above modules.

[0204] In an exemplary embodiment, a computer device is provided, which may be a server, and its internal structure diagram may be shown in FIG11. The computer device includes a processor, a memory, an input / output interface (Input / Output, referred to as I / O) and a communication interface. The processor, the memory and the input / output interface are connected via a system bus, and the communication interface is connected to the system bus via the input / output interface. The processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The database of the computer device is used to store power system model data. The input / output interface of the computer device is used to exchange information between the processor and an external device. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, a stability assessment method for a power system is implemented.

[0205] Those skilled in the art will understand that the structure shown in FIG. 11 is merely a block diagram of a partial structure related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than shown in the figure, or combine certain components, or have a different arrangement of components.

[0206] In an exemplary embodiment, a computer device is provided, including a memory and a processor, wherein a computer program is stored in the memory, and when the processor executes the computer program, the following steps are implemented:

[0207] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0208] Analyze the power system model and determine the target boundary of the power system model;

[0209] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0210] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0211] Based on the linear matrix inequality relationship and the power system model, determining a first preset matrix of the power system within the operation domain;

[0212] Determine the continuity constants and related parameters within the operation domain based on the preset continuity conditions and power system model;

[0213] According to any two data points in the power system model, obtaining a second preset matrix corresponding to the power system model;

[0214] A target boundary of the power system model is determined according to the first preset matrix, the continuous constant, the relevant parameters in the operation domain and the second preset matrix.

[0215] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0216] Obtaining a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of a first preset matrix;

[0217] Obtaining a second matrix minimum eigenvalue of a second preset matrix;

[0218] The target boundary of the power system model is determined according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value in the operation domain and the minimum eigenvalue of the second matrix.

[0219] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0220] If the power system model meets the preset continuity conditions, then the relevant parameters in the operation domain are determined according to the power system model;

[0221] The continuous constant is determined according to the relevant parameters in the operation domain and the relevant parameters of each state variable in the operation domain.

[0222] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0223] Obtain the difference between the power system models corresponding to the two data points;

[0224] According to the difference between the power system models corresponding to the two data points, a second preset matrix corresponding to the power system model is determined.

[0225] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0226] Performing coordinate transformation on the power system model to obtain the power system model after coordinate transformation;

[0227] Analyze the power system model and determine the target boundaries of the power system model, including:

[0228] Analyze the power system model after coordinate transformation and determine the target boundary of the power system model.

[0229] In one embodiment, when the processor executes the computer program, the processor further implements the following steps:

[0230] Applying preset control to the power system model to obtain the power system model after the preset control is applied; the preset control includes integral proportional control and negative feedback control;

[0231] Analyze the power system model and determine the target boundaries of the power system model, including:

[0232] The power system model after applying preset control is analyzed to determine the target boundary of the power system model.

[0233] In one embodiment, a computer readable storage medium is provided, on which a computer program is stored, and when the computer program is executed by a processor, the following steps are implemented:

[0234] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0235] Analyze the power system model and determine the target boundary of the power system model;

[0236] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0237] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0238] Based on the linear matrix inequality relationship and the power system model, determining a first preset matrix of the power system within the operation domain;

[0239] Determine the continuity constants and related parameters within the operation domain based on the preset continuity conditions and power system model;

[0240] According to any two data points in the power system model, obtaining a second preset matrix corresponding to the power system model;

[0241] A target boundary of the power system model is determined according to the first preset matrix, the continuous constant, the relevant parameters in the operation domain and the second preset matrix.

[0242] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0243] Obtaining a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of a first preset matrix;

[0244] Obtaining a second matrix minimum eigenvalue of a second preset matrix;

[0245] The target boundary of the power system model is determined according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value in the operation domain and the minimum eigenvalue of the second matrix.

[0246] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0247] If the power system model meets the preset continuity conditions, then the relevant parameters in the operation domain are determined according to the power system model;

[0248] The continuous constant is determined according to the relevant parameters in the operation domain and the relevant parameters of each state variable in the operation domain.

[0249] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0250] Obtain the difference between the power system models corresponding to the two data points;

[0251] According to the difference between the power system models corresponding to the two data points, a second preset matrix corresponding to the power system model is determined.

[0252] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0253] Performing coordinate transformation on the power system model to obtain the power system model after coordinate transformation;

[0254] Analyze the power system model and determine the target boundaries of the power system model, including:

[0255] Analyze the power system model after coordinate transformation and determine the target boundary of the power system model.

[0256] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0257] Applying preset control to the power system model to obtain the power system model after the preset control is applied; the preset control includes integral proportional control and negative feedback control;

[0258] Analyze the power system model and determine the target boundaries of the power system model, including:

[0259] The power system model after applying preset control is analyzed to determine the target boundary of the power system model.

[0260] In one embodiment, a computer program product is provided, comprising a computer program, which, when executed by a processor, implements the following steps:

[0261] According to multiple relevant parameters of synchronous generators in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator;

[0262] Analyze the power system model and determine the target boundary of the power system model;

[0263] The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

[0264] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0265] Based on the linear matrix inequality relationship and the power system model, determining a first preset matrix of the power system within the operation domain;

[0266] Determine the continuity constants and related parameters within the operation domain based on the preset continuity conditions and power system model;

[0267] According to any two data points in the power system model, obtaining a second preset matrix corresponding to the power system model;

[0268] A target boundary of the power system model is determined according to the first preset matrix, the continuous constant, the relevant parameters in the operation domain and the second preset matrix.

[0269] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0270] Obtaining a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of a first preset matrix;

[0271] Obtaining a second matrix minimum eigenvalue of a second preset matrix;

[0272] The target boundary of the power system model is determined according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value in the operation domain and the minimum eigenvalue of the second matrix.

[0273] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0274] If the power system model meets the preset continuity conditions, then the relevant parameters in the operation domain are determined according to the power system model;

[0275] The continuous constant is determined according to the relevant parameters in the operation domain and the relevant parameters of each state variable in the operation domain.

[0276] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0277] Obtain the difference between the power system models corresponding to the two data points;

[0278] According to the difference between the power system models corresponding to the two data points, a second preset matrix corresponding to the power system model is determined.

[0279] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0280] Performing coordinate transformation on the power system model to obtain the power system model after coordinate transformation;

[0281] Analyze the power system model and determine the target boundaries of the power system model, including:

[0282] Analyze the power system model after coordinate transformation and determine the target boundary of the power system model.

[0283] In one embodiment, when the computer program is executed by a processor, the following steps are also implemented:

[0284] Applying preset control to the power system model to obtain the power system model after the preset control is applied; the preset control includes integral proportional control and negative feedback control;

[0285] Analyze the power system model and determine the target boundaries of the power system model, including:

[0286] The power system model after applying preset control is analyzed to determine the target boundary of the power system model.

[0287] A person of ordinary skill in the art can understand that all or part of the processes in the above-mentioned embodiment method can be completed by instructing the relevant hardware through a computer program, and the computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above-mentioned methods. Among them, any reference to the memory, database or other medium used in the embodiments provided in the present application can include at least one of non-volatile memory and volatile memory. Non-volatile memory can include read-only memory (ROM), magnetic tape, floppy disk, flash memory, optical memory, high-density embedded non-volatile memory, resistive random access memory (ReRAM), magnetic random access memory (MRAM), ferroelectric random access memory (FRAM), phase change memory (PCM), graphene memory, etc. Volatile memory can include random access memory (RAM) or external cache memory, etc. As an illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM). The database involved in each embodiment provided in this application may include at least one of a relational database and a non-relational database. Non-relational databases may include distributed databases based on blockchains, etc., but are not limited to this. The processor involved in each embodiment provided in this application may be a general-purpose processor, a central processing unit, a graphics processor, a digital signal processor, a programmable logic device, a data processing logic device based on quantum computing, an artificial intelligence (AI) processor, etc., but are not limited to this.

[0288] The technical features of the above embodiments may be combined arbitrarily. To make the description concise, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, they should be considered to be within the scope of this application.

[0289] The above-described embodiments only express several implementation methods of the present application, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the scope of the present application. It should be pointed out that, for a person of ordinary skill in the art, several variations and improvements can be made without departing from the concept of the present application, and these all belong to the protection scope of the present application. Therefore, the protection scope of the present application shall be subject to the attached claims.

Claims

1. A method for evaluating the stability of a power system, characterized in that: The method comprises: According to multiple relevant parameters of the synchronous generator in the power system, a power system model carrying slow-varying disturbances is constructed; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator; Analyzing the power system model to determine a target boundary of the power system model; The stability of the power system is evaluated according to the target boundary to determine the stability of the power system.

2. The method according to claim 1, characterized in that: The analyzing the power system model to determine the target boundary of the power system model includes: Based on the linear matrix inequality relationship and the power system model, determining a first preset matrix of the power system in the operation domain; Determining a continuous constant and related parameters within the operation domain based on a preset continuous condition and the power system model; According to any two data points in the power system model, obtaining a second preset matrix corresponding to the power system model; A target boundary of the power system model is determined according to the first preset matrix, the continuous constant, relevant parameters within the operating domain, and the second preset matrix.

3. The method according to claim 2, characterized in that The step of determining the target boundary of the power system model according to the first preset matrix, the continuous constant, the maximum value in the operation domain, and the second preset matrix includes: Obtaining a first matrix maximum eigenvalue and a first matrix minimum eigenvalue of the first preset matrix; Obtaining a second matrix minimum eigenvalue of the second preset matrix; The target boundary of the power system model is determined according to the maximum eigenvalue of the first matrix, the minimum eigenvalue of the first matrix, the continuous constant, the maximum value within the operating domain and the minimum eigenvalue of the second matrix.

4. The method according to claim 2, characterized in that: The determining of the continuous constant and the related parameters in the operation domain based on the preset continuous condition and the power system model includes: If the power system model satisfies the preset continuity condition, determining relevant parameters within the operation domain according to the power system model; The continuous constant is determined according to relevant parameters in the operation domain and relevant parameters of each state variable in the operation domain.

5. The method according to claim 2, characterized in that: The step of obtaining a second preset matrix corresponding to the power system model according to any two data points in the power system model includes: Obtaining a difference between the power system models corresponding to the two data points; According to the difference between the power system models corresponding to the two data points, a second preset matrix corresponding to the power system model is determined.

6. The method according to any one of claims 1 to 5, characterized in that: The method further comprises: Performing coordinate transformation on the power system model to obtain a power system model after coordinate transformation; The analyzing the power system model to determine the target boundary of the power system model includes: The power system model after the coordinate transformation is analyzed to determine the target boundary of the power system model.

7. The method according to any one of claims 1 to 5, characterized in that: The method further comprises: Applying a preset control to the power system model to obtain the power system model after the preset control is applied; the preset control includes integral proportional control and negative feedback control; The analyzing the power system model to determine the target boundary of the power system model includes: The power system model after applying the preset control is analyzed to determine the target boundary of the power system model.

8. A stability assessment device for a power system, characterized in that: The device comprises: A construction module is used to construct a power system model carrying slow-varying disturbances according to multiple relevant parameters of a synchronous generator in the power system; the relevant parameters include d-axis reactance, d-axis transient reactance, total circuit reactance, q-axis open-loop transient time constant and control coefficient in the synchronous generator; An analysis module, used to analyze the power system model and determine a target boundary of the power system model; An evaluation module is used to evaluate the stability of the power system according to the target boundary to determine the stability of the power system.

9. A computer device comprising a memory and a processor, wherein the memory stores a computer program, wherein: When the processor executes the computer program, the steps of the method according to any one of claims 1 to 7 are implemented.

10. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 7 are implemented.

Images