Emergency generator tripping and load shedding control method and system for alternating-current and direct-current power transmission system

Abstract

The invention discloses an emergency generator tripping and load shedding control method and system for an AC / DC power transmission system, and mainly relates to the technical field of high-voltage DC power transmission. Comprising the following steps: establishing a high-voltage direct-current power transmission system model based on a dynamic phasor model, and carrying out modeling optimization on transient stability emergency control of a power system; converting the transient stability constraint into an integral form by using a constraint aggregation method; under a solving framework of a direct sequential method, converting the transient stability emergency control model of the power system into a nonlinear programming problem and a time domain simulation problem; a nonlinear programming problem is converted into a standard form of an optimization algorithm, and a Kriging function proxy model optimization algorithm framework based on dynamic updating is designed to solve the problem. The method has the beneficial effects that the calculation efficiency can be improved on the premise that the calculation precision is not obviously influenced.

Description

Technical Field

[0001] The present invention relates to the technical field of direct current (DC) power transmission in power grids, and in particular to a method and system for controlling emergency generator and load shedding in an AC / DC power transmission system. Background Art

[0002] As the scale of power grids continues to expand, high-voltage direct current (HVDC) transmission systems have been widely used as an economical and reliable cross-regional power transmission method. However, if AC / DC systems fail to commutate due to a grid fault, which in turn triggers a large-scale cascading failure, it will have a huge impact on the safe operation of the grid. After the AC / DC grid suffers an asymmetric fault, how to calculate reasonable generator and load shedding control strategies is crucial to ensuring the safe and stable operation of the grid.

[0003] Traditional emergency generator and load shedding models often use a quasi-steady-state HVDC system DC model. However, these quasi-steady-state models are inaccurate when simulating asymmetric faults in the AC system and inverter-side commutation failures. Related research shows that asymmetric faults are more likely to cause commutation failures. Three consecutive commutation failures will trigger DC blocking protection, leading to large-scale power flow transfer in the AC / DC hybrid system and severe cascading failures. In existing academic research and engineering applications, researchers have used electromagnetic transient models to improve the accuracy of HVDC system calculations. However, these models are computationally inefficient and cannot meet the computational time requirements of emergency control algorithms.

[0004] Therefore, there is an urgent need for a method and system for generating an emergency generator and load shedding control strategy for an AC / DC transmission system to solve the above problems. Summary of the Invention

[0005] The object of the present invention is to provide an AC / DC power transmission system emergency generator and load shedding control method and system, which can improve the calculation efficiency without significantly affecting the calculation accuracy.

[0006] To achieve the above-mentioned purpose, the present invention is implemented through the following technical solutions:

[0007] In one aspect, a method for controlling emergency generator and load shedding in an AC / DC power transmission system is provided, preferably comprising the following steps:

[0008] S1: Establish a HVDC transmission system based on a dynamic phasor model and optimize the transient stability emergency control of the power system. Establish a dynamic optimization model based on differential-algebraic equations describing the AC / DC power system and transient stability as constraints, with the minimum amount of generator and load shedding as the objective function.

[0009] S2: Using the path constraint aggregation method, the multi-dimensional transient stability constraints are transformed into a one-dimensional integral form;

[0010] S3: According to step S2, the modeling optimization model of the power system transient stability emergency control in step S1 is converted into a nonlinear programming problem and an AC / DC system time domain simulation problem;

[0011] S4: Convert the nonlinear programming problem in step S3 into a standard form of an optimization algorithm based on a radial basis function Kriging function surrogate model.

[0012] S5: The time domain simulation problem of the AC / DC system in step S3 is used to establish a composite sequence network simulation model considering the momentum phasor model.

[0013] S6: Use Latin hypercube sampling method to sample and establish the Kriging function proxy model of the AC / DC system in the original step S4.

[0014] S7: Execute the dynamic proxy model optimization strategy and use the Memetic algorithm to find the optimal solution of the nonlinear model in step S3. If the convergence criterion is not met, repeat step S6 until convergence. Preferably, the step S1 is specifically as follows:

[0015]

[0016] Among them, Ф(u) is the amount of cutting machine and load, <x> k (t)、 <y> k (t) is the state variable and algebraic variable in the dynamic phasor model of the HVDC transmission system; the differential equation f D and the algebraic equation g D It is a dynamic phasor model that describes the HVDC system; <x (0) > k represents the initial value of the dynamic phasor model, which can be obtained through power flow calculation; t∈[0, T], T is the simulation planting time; x is the state variable vector of the generator, HVDC transmission system and related controllers, x0 is the initial value of x; y is the algebraic variable vector including the bus node voltage; u is the control variable vector including the generator cutting u G , load shedding u L ; f is a differential equation, including the generator rotor equation and the HVDC system dynamic equation; G is a set of inequality constraints, including power angle constraints and shutdown angle constraints; G and Indicates the lower and upper bounds of the constraint.

[0017] Preferably, the step S2 is specifically as follows:

[0018]

[0019] Among them, the variable θ(t|u) is defined as:

[0020]

[0021] Preferably, the step S3 is specifically as follows:

[0022] Nonlinear programming problem:

[0023] min u Ф(u)

[0024]

[0025] AC / DC system time domain simulation problems:

[0026]

[0027] The differential equation in formula (5) is converted into a difference equation by numerical differentiation.

[0028] Preferably, the step S4 is specifically as follows:

[0029] min u Ф(u)+σh(u)

[0030]

[0031] Among them, σ is the penalty factor, which is usually a large number.

[0032] Preferably, it also includes:

[0033] Design of hybrid simulation interface for DC system with electromechanical transient and improved switching function;

[0034] t k Time and t k The converter AC bus voltage at time +Δt is When calculating the DC system, linear interpolation is used to obtain τ k Voltage at the moment:

[0035]

[0036] The current injected into the AC system by the converter can be expressed as the average injected current within a large simulation step Δt:

[0037]

[0038] Where superscript "s" is "+", "-", or "0", representing positive-sequence, negative-sequence, and zero-sequence voltages, respectively. m = Δt / Δτ represents the number of interpolation points in the small-step system.

[0039] Modeling and simulation of the AC / DC network model in asymmetric operation state are performed. Based on the improved switching function, the injection current of each sequence of the DC converter at the interface bus is calculated, and the voltage of each node of the AC network is solved.

[0040] Preferably, the voltage of each node in the AC network is obtained by:

[0041] Under asymmetric faults, the positive sequence equivalent network of the AC / DC system is:

[0042]

[0043] Under asymmetric faults, the negative sequence equivalent network of the AC / DC system is:

[0044]

[0045] Under asymmetric faults, the zero-sequence equivalent network of the AC / DC system is:

[0046]

[0047] The fault boundary condition equation is:

[0048]

[0049] The subscript "c" can be "r" and "i", which represent the rectifier node and the inverter node respectively; + , Y - , Y 0 They represent the positive and negative zero-sequence node admittance matrices of the AC and DC systems respectively; B represents the boundary conditions of different faults; and They represent the positive and negative zero-sequence currents equivalently injected by the DC converter respectively; Indicates the positive and negative zero-sequence voltages at the fault point; and Represent the positive and negative zero-sequence currents at the fault point respectively.

[0050] Preferably, the step S4 includes:

[0051] Based on the Kriging function surrogate model, Latin hypercube sampling method is adopted to establish the Kriging function surrogate model of the original AC / DC system.

[0052] The basic form of the Kriging function is:

[0053]

[0054] Among them, N s is the sample size, u i represents the i-th sampling point, ω i is the weight coefficient, ω is the weight vector, z(x) is the random error of x, which satisfies the mean of 0 and the variance of σ 2 Normal distribution.

[0055] Preferably, the step S7 includes: executing a dynamic proxy model optimization strategy, using a Memetic algorithm to optimize and solve the model described in formula (5) to obtain the load shedding amount, until convergence. It is characterized in that:

[0056]

[0057] When k>1, it is necessary to calculate the relative error of the true objective function value corresponding to the possible optimal solution calculated for the kth time and the k-1th time, and judge whether the relative error meets the given convergence criterion. If it satisfies formula (14), the loop is stopped, then This is the optimal solution to the original optimization problem; if it is not satisfied, execute step S6.

[0058] On the other hand, a control system based on an AC / DC power transmission system emergency generator and load shedding control method is provided, preferably comprising:

[0059] HVDC system modeling module, used to establish a HVDC system based on a dynamic phasor model and to establish an emergency generator and load shedding model for AC and DC power systems;

[0060] Process constraint conversion module, used to convert transient stability constraints into integral form;

[0061] A control model generation module is used to convert the emergency generator and load shedding model of the AC / DC power system into a nonlinear programming model and an AC / DC system time-domain simulation model;

[0062] The problem-solving module is used to solve nonlinear programming problems and has a built-in dynamic optimization algorithm based on the Kriging function surrogate model.

[0063] Compared with the prior art, the beneficial effects of the present invention are:

[0064] 1. It can generate emergency control strategies when asymmetric faults occur in the AC system, avoiding the defect that the quasi-steady-state model cannot accurately describe the dynamic characteristics of the HVDC transmission system under asymmetric operating conditions, and improving the reliability of the emergency control strategies under asymmetric faults.

[0065] 2. By reducing the difficulty of solving the optimization problem and avoiding the calculation of complex gradient information, this algorithm improves the effectiveness and computational efficiency of solving transient stability emergency control problems. BRIEF DESCRIPTION OF THE DRAWINGS

[0066] Figure 1 is a flow chart of the method of the present invention;

[0067] Figure 2 It is a schematic diagram of the interactive relationship between AC and DC simulation data of the present invention;

[0068] Figure 3 This is a schematic diagram of a network with equivalent values ​​of each sequence under an asymmetric short circuit fault of the present invention;

[0069] Figure 4 is a 22-node system topology diagram including a high-voltage direct current transmission system of the present invention;

[0070] Figure 5 This is a schematic diagram of the time domain simulation verification of the emergency control strategy of the 22-node system of the present invention;

[0071] Figure 6 1. It is a schematic diagram of time domain simulation verification of the emergency control strategy of the present invention in different DC models;

[0072] Figure 7 It is a schematic diagram of the system structure of the present invention. DETAILED DESCRIPTION

[0073] Below in conjunction with specific embodiment, further set forth the present invention.Should be understood that these embodiments are only used to illustrate the present invention and are not used in limiting the scope of the present invention.In addition, should be understood that after reading the content taught by the present invention, those skilled in the art can make various changes or modifications to the present invention, and these equivalent forms fall within the scope limited by the application equally.

[0074] In the present invention, terms such as "upper", "lower", "left", "right", "front", "back", "vertical", "horizontal", "side", "bottom", etc. indicating directions or positional relationships are based on the directions or positional relationships shown in the accompanying drawings. They are merely relational terms determined for the convenience of describing the structural relationships of the various parts or elements of the present invention, and do not specifically refer to any part or element in the present invention, and should not be understood as limiting the present invention.

[0075] In the present invention, terms such as "fixed connection," "connected," and "connection" should be interpreted broadly to mean a fixed connection, an integral connection, or a detachable connection; a direct connection or an indirect connection through an intermediary. Relevant researchers or technicians in this field may determine the specific meanings of these terms in the present invention based on specific circumstances, and they should not be construed as limitations of the present invention.

[0076] Example:

[0077] like Figure 1 As shown, this embodiment provides an AC / DC power transmission system emergency generator and load shedding control method, comprising the following steps:

[0078] S1: Considering the HVDC transmission system based on the dynamic phasor model, the power system transient stability emergency control can be modeled as a dynamic optimization problem of formula (1):

[0079]

[0080]

[0081] Among them, Ф(u) is the amount of cutting machine and load, <x> k (t)、 <y> k (t) is the state variable and algebraic variable in the dynamic phasor model of the HVDC transmission system; the differential equation f D and the algebraic equation g D It is a dynamic phasor model that describes the HVDC system; <x (0) > k represents the initial value of the dynamic phasor model, which can be obtained through power flow calculation; t∈[0, T], T is the simulation planting time; x is the state variable vector of the generator, HVDC transmission system and related controllers, x0 is the initial value of x; y is the algebraic variable vector including the bus node voltage; u is the control variable vector including the generator cutting u G , load shedding u L ; f is a differential equation, including the generator rotor equation and the HVDC system dynamic equation; G is a set of inequality constraints, including power angle constraints and shutdown angle constraints; G and Indicates the lower and upper bounds of the constraints;

[0082] S2: Using the path constraint aggregation method, the transient stability constraint can be transformed into the integral form of the formula:

[0083]

[0084] Among them, the variable θ(t|u) is defined as:

[0085]

[0086] S3: The original optimization problem (1) is transformed into a nonlinear programming problem:

[0087] min u Ф(u)

[0088]

[0089] AC / DC system time domain simulation problems:

[0090]

[0091] S4: The above optimization problem (4) is transformed into the standard form of the optimization algorithm based on the radial basis function surrogate model. The transient stability constraint in the optimization is added as a penalty function to the objective function, thereby further transforming the emergency control model into a nonlinear optimization problem:

[0092] min u Ф(u)+σh(u)

[0093]

[0094] Among them, σ is the penalty factor, which is usually a large number. Compared with the original optimization problem, the optimization scale of formula (5) is significantly reduced, and the calculation of the transient stability constraint h(u) includes the time domain simulation process of the AC / DC system.

[0095] S5: Design of hybrid simulation interface of electromechanical transient and DC system with improved switching function. Electromechanical transient of alternating iteration method requires multiple iterations in one time step calculation, such as Figure 2 As shown, the backward dashed arrow indicates that the HVDC transmission system calculates and provides the injection current to the AC system based on the voltage information provided by the AC side, and then re-performs the electromechanical transient simulation within this large step size until the current time step simulation meets the convergence criterion;

[0096] t k Time and t k The converter AC bus voltage at time +Δt is When calculating the DC system, linear interpolation is used to obtain τ k Voltage at the moment:

[0097]

[0098] The current injected into the AC system by the converter can be expressed as the average injected current within a large simulation step Δt:

[0099]

[0100] Where superscript "s" is "+", "-", or "0", representing positive-sequence, negative-sequence, and zero-sequence voltages, respectively. m = Δt / Δτ represents the number of interpolation points in the small-step system.

[0101] Modeling and simulation of the AC / DC network model in asymmetric operation state, according to the improved switching function, the sequence injection current of the DC converter at the interface bus is obtained; by solving the network equations and fault boundary conditions simultaneously, the voltage of each node of the AC network can be obtained, such as Figure 3 (a), (b) and (c) show specifically:

[0102] Under asymmetric faults, the positive sequence equivalent network of the AC / DC system is:

[0103]

[0104] Under asymmetric faults, the negative sequence equivalent network of the AC / DC system is:

[0105]

[0106] Considering the wiring method of the converter transformer of the HVDC transmission system, the zero-sequence component of the DC system cannot flow into the AC system. From the AC port, the DC system is equivalent to a zero-sequence impedance. Therefore, under an asymmetric fault, the zero-sequence equivalent network of the AC and DC systems is:

[0107]

[0108] The fault boundary condition equation is:

[0109]

[0110] The subscript "c" can be "r" and "i", which represent the rectifier node and the inverter node respectively; + , Y - , Y 0 They represent the positive and negative zero-sequence node admittance matrices of the AC and DC systems respectively; B represents the boundary conditions of different faults; and They represent the positive and negative zero-sequence currents equivalently injected by the DC converter respectively; Indicates the positive and negative zero-sequence voltages at the fault point; and Represent the positive and negative zero-sequence currents of the fault point respectively;

[0111] S6: Use Latin hypercube sampling method to sample and establish the Kriging function proxy model of the original AC / DC system.

[0112] The basic form of the Kriging function is:

[0113]

[0114] Among them, N s is the sample size, u i represents the i-th sampling point, w i is the weight coefficient, is the weight vector, z(x) is the random error of x, which satisfies the mean of 0 and the variance of σ 2 Normal distribution.

[0115] S7: Execute the dynamic agent model optimization strategy and use the Memetic algorithm to optimize the model described in formula (5) to obtain the load shedding amount until convergence. The convergence criterion is:

[0116]

[0117] When k>1, it is necessary to calculate the relative error of the true objective function value corresponding to the possible optimal solution calculated for the kth time and the k-1th time, and judge whether the relative error meets the given convergence criterion. If it satisfies formula (14), the loop is stopped, then This is the optimal solution to the original optimization problem; if it is not satisfied, execute step S6.

[0118] like Figure 4 The figure shows a schematic diagram of a 22-node grid structure in this embodiment. The above transmission system model was built using the MATLAB software platform. Six scenarios, including single-phase ground fault, two-phase interphase fault, and two-phase ground fault, were applied to lines 19-21 and 16-20, respectively. The implementation results are shown based on case data:

[0119] Operating environment:

[0120] Intel Core i3-10105 CPU 3.70GHz, 16GB RAM, Microsoft Windows10X64MATLAB 2020b

[0121] Implementation results: The simulation results of the generator power angle curve with / without emergency control measures under different fault scenarios are given, such as Figure 5 (a), (b), (c), (d), (e), and (f) show the simulation results of the generator power angle curve with and without emergency control measures under different fault scenarios. The solid line represents the relative power angle of the synchronous generator after the emergency control measures are applied.

[0122] The dashed line represents the relative power angle of the synchronous generators without emergency control measures. In scenarios 3, 5, and 6, the system can maintain transient stability through relay protection action (disconnecting the faulty line). However, in scenarios 1, 2, and 4, the system power angle will lose synchronization. Without emergency control measures, relying solely on the first line of defense (disconnecting the faulty line) cannot guarantee system transient stability. Therefore, the simulation results show that a more accurate model is beneficial for improving the reliability of the transient stability emergency control strategy.

[0123] The emergency control strategies generated based on the quasi-steady-state model are applied to the systems with different DC models, such as Figure 6 The simulation results show that in the two AC asymmetric fault scenarios mentioned above, the emergency control strategy generated based on the quasi-steady-state model can prevent transient instability of the system power angle in the quasi-steady-state model, but fails in the AC / DC system based on the dynamic phasor model.

[0124] According to the results of this case, it can be seen that the proposed electromechanical transient-improved switching function model hybrid simulation framework can simulate and calculate the system under different faults, which makes up for the defect that the quasi-steady-state model of the high-voltage direct current transmission system cannot accurately simulate asymmetric faults. The algorithm proposed in this chapter can improve the reliability of the emergency control strategy under asymmetric faults.

[0125] like Figure 7 As shown, this embodiment also provides an AC / DC power transmission system emergency generator and load shedding control system, including:

[0126] HVDC system modeling module, used to establish a HVDC system based on a dynamic phasor model and to establish an emergency generator and load shedding model for AC and DC power systems;

[0127] Process constraint conversion module, used to convert transient stability constraints into integral form;

[0128] A control model generation module is used to convert the emergency generator and load shedding model of the AC / DC power system into a nonlinear programming model and an AC / DC system time-domain simulation model;

[0129] The problem-solving module is used to solve nonlinear programming problems and has a built-in optimization algorithm based on the Kriging function surrogate model.

[0130] The above is a specific description of the preferred implementation of the present invention, but the present invention is not limited to the described embodiments. Those skilled in the art can make various equivalent modifications or substitutions to the transaction features between nodes without violating the spirit of the present invention. These equivalent modifications or substitutions are all included in the scope defined by the claims of this application.< / y> < / x> < / y> < / x>

Claims

1. A method for controlling emergency generator and load shedding in an AC / DC power transmission system, characterized in that: The following steps are involved: S1: Establish a dynamic optimization model for a HVDC transmission system based on a dynamic phasor model, using differential-algebraic equations describing the dynamic behavior of the AC / DC power system and transient power angle stability as constraints, with the minimum amount of generator and load shedding as the objective function. S2: Using the constraint aggregation method, the multi-dimensional transient stability constraints are transformed into a one-dimensional integral form; S3: According to step S2, the power system transient stability emergency control model in step S1 is converted into a nonlinear programming problem and an AC / DC system time domain simulation problem; S4: Convert the nonlinear programming problem in step S3 into a standard form of an optimization algorithm based on a Kriging function surrogate model. S5: Establish a composite sequence network simulation model considering the momentum phasor model to solve the AC / DC system time domain simulation problem in step S3. S6: Use Latin hypercube sampling method to sample and establish the Kriging function proxy model of the AC / DC system in the original step S4. S7: Execute the dynamic agent model optimization strategy and use the Memetic algorithm to solve the optimal solution of the nonlinear model in step S3. If the convergence criterion is not met, repeat step S6 until convergence.

2. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: The step S1 is specifically as follows: minute u F(u) Among them, Ф(u) is the amount of cutting machine and load, <x> k (t)、 <y> k (t) is the state variable and algebraic variable in the dynamic phasor model of the HVDC transmission system; the differential equation f D and the algebraic equation g D It is a dynamic phasor model that describes the HVDC system; <x (0) > k represents the initial value of the dynamic phasor model, which can be obtained through power flow calculation; t∈[0, T], T is the simulation planting time; x is the state variable vector of the generator, HVDC transmission system and related controllers, x0 is the initial value of x; y is the algebraic variable vector including the bus node voltage; u is the control variable vector including the generator cutting u G , load shedding u L ; f is a differential equation, including the generator rotor equation and the HVDC system dynamic equation; G is a set of inequality constraints, including power angle constraints and shutdown angle constraints; G and Indicates the lower and upper bounds of the constraint.< / y> < / x> 3. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: The step S2 is specifically as follows: Among them, the variable θ(t|u) is defined as:

4. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: The step S3 is specifically as follows: Nonlinear programming problem: minute u F(u) AC / DC system time domain simulation problems: The differential equation in formula (5) is converted into a difference equation by numerical differentiation.

5. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: The step S4 is specifically as follows: min u F(u)+σh(u) Among them, σ is the penalty factor, which is usually a large number.

6. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: Also includes: The design of the hybrid simulation interface of AC / DC system with electromechanical transient-improved switching function is as follows: t k Time and t k The converter AC bus voltage at time +Δt is When calculating the DC system, linear interpolation is used to obtain τ k Voltage at the moment: The current injected into the AC system by the converter can be expressed as the average injected current within a large simulation step Δt: Where superscript "s" is "+", "-", or "0", representing positive-sequence, negative-sequence, and zero-sequence voltages, respectively. m = Δt / Δτ represents the number of interpolation points in the small-step system. The AC / DC network model in asymmetric operating state is modeled and simulated. The injected current of the DC converter at the interface bus is calculated based on the improved switching function, and is used to solve the voltage of each node in the AC network.

7. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: The step S6 is used to solve the voltage of each node of the AC network, specifically: Under asymmetric faults, the positive sequence equivalent network of the AC / DC system is: Under asymmetric faults, the negative sequence equivalent network of the AC / DC system is: Under asymmetric faults, the zero-sequence equivalent network of the AC / DC system is: The fault boundary condition equation is: The subscript "c" can be "r" and "i", which represent the rectifier node and the inverter node respectively; + , Y - , Y 0 They represent the positive and negative zero-sequence node admittance matrices of the AC and DC systems respectively; B represents the boundary conditions of different faults; and They represent the positive and negative zero-sequence currents equivalently injected by the DC converter respectively; Indicates the positive and negative zero-sequence voltages at the fault point; and Represent the positive and negative zero-sequence currents at the fault point respectively.

8. The method for controlling emergency generator and load shedding in an AC / DC power transmission system according to claim 1, characterized in that: The step S6 is specifically as follows: For the Kriging function-based proxy model in step S4, Latin hypercube sampling method is adopted to sample and establish the Kriging function proxy model of the original AC / DC system. The basic form of the Kriging function is: Among them, N s is the sample size, u i represents the i-th sampling point, ω i is the weight coefficient, ω is the weight vector, z(x) is the random error of x, which satisfies the mean of 0 and the variance of σ 2 Normal distribution.

9. According to the method for controlling emergency generator and load shedding in an AC / DC power transmission system of claim 1, in step S7, a dynamic agent model optimization strategy is executed, and the Memetic algorithm is used to optimize and solve the model described in equation (5) to obtain the generator and load shedding amount until convergence. It is characterized in that: When k>1, it is necessary to calculate the relative error of the true objective function value corresponding to the possible optimal solution calculated for the kth time and the k-1th time, and judge whether the relative error meets the given convergence criterion. If it satisfies formula (14), the loop is stopped, then This is the optimal solution to the original optimization problem; if it is not satisfied, execute step S6.

10. A control system based on the AC / DC power transmission system emergency generator and load shedding control method according to claim 1, characterized in that: include: HVDC system modeling module, used to establish a HVDC system based on a dynamic phasor model and to establish an emergency generator and load shedding model for AC and DC power systems; Process constraint conversion module, used to convert transient stability constraints into integral form; A control model generation module is used to convert the power system transient stability emergency control model into a nonlinear programming model and an AC / DC system time-domain simulation model; The problem-solving module is used to solve nonlinear programming problems and has a built-in optimization algorithm based on the Kriging function surrogate model.

Images