Electromagnetic Steady-State Solving Method and System for Synchronous Generators Based on Fourier Spectrum Method

Abstract

The present invention discloses a method and system for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method, which relates to the technical field of power systems. By using the operation data of the synchronous generator to construct a Fourier spectrum differential matrix, a Fourier spectrum differential matrix is generated; based on the operation data, an electromagnetic transient simulation is carried out to generate an initial electromagnetic transient time-domain model; based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model, a steady-state period solution is solved to generate electromagnetic steady-state solution data. When calculating the three-phase symmetrical steady-state period solution, the Fourier spectrum method can obtain accurate discrete points, and the Fourier spectrum method only needs very few discrete points to obtain an accurate series expression. The present invention constructs a Fourier spectrum differential matrix by using the Fourier spectrum method, and uses the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to solve the steady-state period, so that the steady-state period solution can be calculated with fewer discrete points, the solution speed is fast, and the solution time used is short.

Description

Technical Field

[0001] The present invention relates to the technical field of power systems, and in particular, to a method and system for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method. Background Art

[0002] China's power system is increasingly characterized by a high proportion of new energy, high power electronics, and high AC-DC coupling. The switching process of power electronic converters has a short time scale and also brings harmonics to the power system. In order to more precisely depict the transient process of the power system, it is necessary to perform electromagnetic transient simulation on it.

[0003] Traditional electromagnetic transient simulation requires the use of a small time step. Different from the equilibrium solution of electromechanical transients being an equilibrium point, the equilibrium solution of electromagnetic transients is a periodic solution. The existing method for solving the electromagnetic steady state of a synchronous generator is to use a zero-start method for solving, but this solving method takes a long calculation time to obtain a periodic solution. Summary of the Invention

[0004] The present invention provides a method and system for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method, which solves the technical problem that the existing method for solving the electromagnetic steady state of a synchronous generator uses a zero-start method for solving and takes a long calculation time to obtain a periodic solution.

[0005] A method for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method provided by the present invention includes:

[0006] Obtain the operation data corresponding to the synchronous generator, and use the operation data to construct a Fourier spectrum differential matrix to generate a Fourier spectrum differential matrix;

[0007] Perform electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time domain model;

[0008] Perform steady-state period solution based on the Fourier spectrum differential matrix and the initial electromagnetic transient time domain model to generate electromagnetic steady state solution data.

[0009] Optionally, the step of using the operation data to construct a Fourier spectrum differential matrix to generate a Fourier spectrum differential matrix includes:

[0010] Use the operation data to construct a differential-algebraic equation set to generate the differential-algebraic equation set corresponding to the synchronous generator;

[0011] Sample the periodic function corresponding to the differential-algebraic equation set multiple times within a preset domain to generate a discrete vector;

[0012] Substitute the discrete vector into the differential-algebraic equation set to construct the algebraic equation of the discrete point Fourier spectral method, and generate the algebraic equation of the Fourier spectral method;

[0013] Construct a Fourier spectral differential matrix based on the algebraic equation of the Fourier spectral method and the discrete vector, and generate a Fourier spectral differential matrix.

[0014] Optionally, the step of constructing a Fourier spectral differential matrix based on the algebraic equation of the Fourier spectral method and the discrete vector to generate a Fourier spectral differential matrix includes:

[0015] Perform discrete Fourier transform and inverse discrete Fourier transform on the algebraic equation of the Fourier spectral method to generate a transformed equation;

[0016] Perform derivative discrete point sampling on the transformed equation to generate a derivative discrete vector;

[0017] Construct a differential matrix using the derivative discrete vector and the discrete vector to generate a Fourier spectral differential matrix.

[0018] Optionally, the step of performing electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time-domain model includes:

[0019] Use the operation data to construct a model and generate a synchronous generator model;

[0020] Perform Park transformation on each state parameter in the synchronous generator model to generate an initial electromagnetic transient time-domain model;

[0021] The stator current equation corresponding to the initial electromagnetic transient time-domain model is:

[0022] ;

[0023] Where, is the direct-axis stator current; is the direct-axis open-circuit transient time constant; is the direct-axis open-circuit subtransient time constant; is the direct-axis transient electromotive force; is the direct-axis subtransient electromotive force; is the direct-axis transient reactance; is the direct-axis subtransient reactance; is the quadrature-axis stator current; is the quadrature-axis open-circuit transient time constant; is the quadrature-axis open-circuit subtransient time constant; is the quadrature-axis transient electromotive force; is the quadrature-axis subtransient electromotive force; is the quadrature-axis transient reactance; is the quadrature-axis subtransient reactance; is the stator resistance; is the direct-axis stator voltage; is the quadrature-axis stator voltage; is the rotor motion angle relative to the synchronous rotating coordinate system; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the zero-sequence current of the stator winding; is the zero-sequence inductance of the stator winding; is the zero-sequence voltage of the stator winding;

[0024] The generator electromagnetic power equation corresponding to the initial electromagnetic transient time-domain model is:

[0025] ;

[0026] where is the generator electromagnetic power; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the quadrature-axis stator current; is the quadrature-axis subtransient electromotive force; is the direct-axis stator current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance;

[0027] The generator electromagnetic torque equation corresponding to the initial electromagnetic transient time-domain model is:

[0028] ;

[0029] where is the generator electromagnetic torque; is the quadrature-axis stator current; is the quadrature-axis subtransient electromotive force; is the direct-axis stator current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance.

[0030] Optionally, the step of performing steady-state period solution based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to generate electromagnetic steady-state solution data includes:

[0031] Calculating the initial values of each state variable in the synchronous generator by using the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator to generate a target initial value data set;

[0032] Using the target initial data set to update the synchronous generator differential equation group in the initial electromagnetic transient time-domain model to generate a target electromagnetic transient time-domain model;

[0033] The target electromagnetic transient time domain model is solved for discrete points within a period using a Fourier spectrum method to generate a target discrete transposed vector;

[0034] Using the Fourier series coefficients corresponding to the target discrete transposed vector to update the Fourier series expression corresponding to the Fourier spectrum differential matrix, to generate a target Fourier series expression;

[0035] The Shannon sampling law and the target Fourier series expression are used to perform electromagnetic transient solution and steady-state periodic solution on the synchronous generator to generate electromagnetic steady-state solution data.

[0036] Optionally, the step of performing initial value calculation on each state variable in the synchronous generator by using the power flow data corresponding to the operating data and the power angle relationship data corresponding to the synchronous generator to generate a target initial value data set includes:

[0037] The power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator are used to perform The initial voltage and current values of the axis are calculated and generated The initial value of shaft voltage and Initial value of shaft current;

[0038] The flow data, the The initial value of the shaft voltage and the Substituting the initial value of the shaft current into a preset initial value formula of the state variable to perform initialization calculation on each state variable in the initial electromagnetic transient time domain model to generate an initial value data set;

[0039] Adopt the Initial value of shaft voltage, The shaft current initial value and the initial initial value data set are used to construct a target initial value data set.

[0040] Optionally, the step of using the Fourier spectrum method to solve the target electromagnetic transient time domain model for discrete points within a period to generate a target discrete transposed vector includes:

[0041] Initializing the target electromagnetic transient time domain model using a preset iteration accuracy and a preset iteration initial value to generate a time domain model;

[0042] Using the state equation function corresponding to each state variable of the time domain model within a period, constructing a discrete point Fourier spectrum method algebraic equation corresponding to the time domain model;

[0043] Calculate the product among the algebraic equation of the discrete point Fourier spectrum method, the Fourier spectrum differential transpose inverse matrix corresponding to the time domain model, and the state variable coefficient inverse matrix, generate an initial discrete transpose vector, and count the number of iterations;

[0044] Use the discrete vector corresponding to the initial discrete transpose vector and the historical discrete vector to judge the iteration accuracy, and generate accuracy judgment data;

[0045] When the accuracy judgment data is less than or equal to the preset iteration accuracy, use the initial discrete transpose vector at the current moment as the target discrete transpose vector;

[0046] When the accuracy judgment data is greater than the preset iteration accuracy, use the time domain model corresponding to the initial discrete transpose vector as the new time domain model, update the number of iterations by 1, and jump to execute the step of calculating the product among the state variable coefficient matrix, the discrete transpose vector, and the Fourier spectrum differential transpose matrix corresponding to the time domain model to generate the algebraic equation of the discrete point Fourier spectrum method corresponding to the time domain model.

[0047] The present invention also provides a synchronous generator electromagnetic steady-state solving system based on the Fourier spectrum method, including

[0048] A Fourier spectrum differential matrix generation module, configured to obtain the operation data corresponding to the synchronous generator, construct a Fourier spectrum differential matrix by using the operation data, and generate a Fourier spectrum differential matrix;

[0049] An initial electromagnetic transient time domain model generation module, configured to perform electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time domain model;

[0050] An electromagnetic steady-state solving data generation module, configured to perform steady-state period solving based on the Fourier spectrum differential matrix and the initial electromagnetic transient time domain model to generate electromagnetic steady-state solving data.

[0051] The present invention also provides an electronic device, including a memory and a processor. When the computer program stored in the memory is executed by the processor, the processor is caused to execute the steps of implementing the synchronous generator electromagnetic steady-state solving method based on the Fourier spectrum method as described in any one of the above.

[0052] The present invention also provides a computer-readable storage medium, on which a computer program is stored. When the computer program is executed, the synchronous generator electromagnetic steady-state solving method based on the Fourier spectrum method as described in any one of the above is implemented.

[0053] It can be seen from the above technical solutions that the present invention has the following advantages:

[0054] The present invention constructs a Fourier spectrum differential matrix by using the Fourier spectrum method. When calculating the three-phase symmetrical steady-state periodic solution, the Fourier spectrum method can obtain accurate discrete points, and only requires very few discrete points to obtain an accurate series expression. Therefore, by using the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model for steady-state periodic solution, the steady-state periodic solution can be calculated with fewer discrete points, providing a warm start for electromagnetic transient simulation, without the need to perform electromagnetic transient simulation until steady state with a large number of discrete points as in the traditional zero-start method. The present invention has a fast solution speed and short solution time used, solving the technical problem that the existing synchronous generator electromagnetic steady-state solution method uses the zero-start method for solution and requires a long calculation time to obtain the periodic solution. BRIEF DESCRIPTION OF THE DRAWINGS

[0055] In order to more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the following will briefly introduce the drawings required for use in the description of the embodiments or the prior art. Obviously, the drawings in the following description are only some embodiments of the present invention. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative labor.

[0056] Figure 1 It is a flowchart of the steps of a synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method provided in Embodiment 1 of the present invention;

[0057] Figure 2 It is a flowchart of the steps of a synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method provided in Embodiment 2 of the present invention;

[0058] Figure 3 It is a schematic diagram of the relationship between the dq rotating coordinate system and the abc three-phase stationary coordinate system provided in Embodiment 2 of the present invention;

[0059] Figure 4 It is a flowchart of a synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method provided in Embodiment 2 of the present invention;

[0060] Figure 5 It is a schematic diagram of the steady-state simulation result of the port voltage without harmonics of the Fourier spectrum method with 3 sampling points provided in Embodiment 2 of the present invention;

[0061] Figure 6 It is a schematic diagram of the steady-state simulation result of the port voltage without harmonics of the Fourier spectrum method with 10 sampling points provided in Embodiment 2 of the present invention;

[0062] Figure 7 It is a schematic diagram of the steady-state simulation result of the port voltage with (5th) harmonics of the Fourier spectrum method with 11 sampling points provided in Embodiment 2 of the present invention;

[0063] Figure 8 Schematic diagram of the steady-state simulation result of the port voltage containing (5th) harmonic of the Fourier spectrum method with 30 sampling points provided by the second embodiment of the present invention;

[0064] Figure 9 Block diagram of a synchronous generator electromagnetic steady-state solution system based on the Fourier spectrum method provided by the third embodiment of the present invention;

[0065] Figure 10 Block diagram of an electronic device provided by the fourth embodiment of the present invention. Detailed implementation manners

[0066] The embodiments of the present invention provide a synchronous generator electromagnetic steady-state solution method and system based on the Fourier spectrum method, which are used to solve the technical problem that the existing synchronous generator electromagnetic steady-state solution method uses a zero-start method for solution and requires a long calculation time to obtain a periodic solution.

[0067] In order to make the objectives, features, and advantages of the present invention more obvious and understandable, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the accompanying drawings in the embodiments of the present invention. Obviously, the embodiments described below are only a part of the embodiments of the present invention, rather than all the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present invention.

[0068] Embodiment 1

[0069] Please refer to Figure 1 , Figure 1 Flowchart of the steps of a synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method provided by the first embodiment of the present invention.

[0070] A synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method provided by the first example of the present invention includes:

[0071] Step 101: Obtain the operating data corresponding to the synchronous generator, and use the operating data to construct a Fourier spectrum differential matrix to generate a Fourier spectrum differential matrix.

[0072] In an embodiment of the present invention, operating data corresponding to a synchronous generator is acquired, and a differential-algebraic equation set is constructed using the operating data to generate a differential-algebraic equation set corresponding to the synchronous generator. The periodic function corresponding to the differential-algebraic equation set is sampled multiple times within a preset domain to generate a discrete vector. The discrete vector is substituted into the differential-algebraic equation set to construct an algebraic equation of the discrete-point Fourier spectrum method, generating an algebraic equation of the Fourier spectrum method. Based on the algebraic equation of the Fourier spectrum method and the discrete vector, a Fourier spectrum differential matrix is constructed, generating a Fourier spectrum differential matrix.

[0073] Step 102: Perform an electromagnetic transient simulation based on the operating data to generate an initial electromagnetic transient time-domain model.

[0074] In an embodiment of the present invention, a synchronous generator model is generated by constructing a model using the operating data. The state variables in the synchronous generator model are subjected to Park transformation to generate an initial electromagnetic transient time-domain model.

[0075] Step 103: Perform a steady-state period solution based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to generate electromagnetic steady-state solution data.

[0076] In an embodiment of the present invention, the initial values of the state variables in the synchronous generator are calculated using the power flow data corresponding to the operating data and the power angle relationship data corresponding to the synchronous generator to generate a target initial value data set. The synchronous generator differential equation set in the initial electromagnetic transient time-domain model is updated using the target initial data set to generate a target electromagnetic transient time-domain model. The target electromagnetic transient time-domain model is solved for discrete points within a period using the Fourier spectrum method to generate a target discrete transposed vector. The Fourier series expression corresponding to the Fourier spectrum differential matrix is updated using the Fourier series coefficients corresponding to the target discrete transposed vector to generate a target Fourier series expression. The electromagnetic transient of the synchronous generator is solved for the steady-state period using the Shannon sampling theorem and the target Fourier series expression to generate electromagnetic steady-state solution data.

[0077] In an embodiment of the present invention, by obtaining the operation data corresponding to a synchronous generator, a Fourier spectrum differential matrix is constructed using the operation data to generate a Fourier spectrum differential matrix; an electromagnetic transient simulation is performed based on the operation data to generate an initial electromagnetic transient time-domain model; and a steady-state period solution is obtained based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to generate electromagnetic steady-state solution data. When calculating the three-phase symmetric steady-state period solution, the Fourier spectrum method can obtain accurate discrete points, and the Fourier spectrum method only requires very few discrete points to obtain an accurate series expression. By constructing a Fourier spectrum differential matrix using the Fourier spectrum method and obtaining a steady-state period solution using the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model, the present invention can calculate the steady-state period solution with fewer discrete points, providing a warm start for electromagnetic transient simulation. There is no need to perform electromagnetic transient simulation until steady state using a large number of discrete points as in the traditional zero-start method. This solves the technical problem that the existing electromagnetic steady-state solution method for synchronous generators uses a zero-start method for solution and requires a long calculation time to obtain a period solution.

[0078] Embodiment 2

[0079] Please refer to Figure 2 , Figure 2 which is a flowchart of the steps of a method for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method provided by Embodiment 2 of the present invention.

[0080] Another method for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method provided by the second embodiment of the present invention includes:

[0081] Step 201: Obtain the operation data corresponding to the synchronous generator, and construct a Fourier spectrum differential matrix using the operation data to generate a Fourier spectrum differential matrix.

[0082] Further, step 201 may include the following sub-steps S11-S14:

[0083] S11: Construct a differential-algebraic equation system using the operation data to generate a differential-algebraic equation system corresponding to the synchronous generator.

[0084] S12: Sample the periodic function corresponding to the differential-algebraic equation system multiple times within a preset domain to generate a discrete vector.

[0085] S13: Substitute the discrete vector into the differential-algebraic equation system to construct an algebraic equation of the discrete-point Fourier spectrum method to generate an algebraic equation of the Fourier spectrum method.

[0086] S14: Construct a Fourier spectrum differential matrix based on the algebraic equation of the Fourier spectrum method and the discrete vector to generate a Fourier spectrum differential matrix.

[0087] Further, step S14 may include the following sub-steps S141-S143:

[0088] S141. Perform discrete Fourier transform and inverse discrete Fourier transform on the algebraic equation of the Fourier spectrum method to generate a transformed equation.

[0089] S142. Perform derivative discrete point sampling on the transformed equation to generate a derivative discrete vector.

[0090] S143. Use the derivative discrete vector and the discrete vector to construct a differential matrix to generate a Fourier spectrum differential matrix.

[0091] In the embodiment of the present invention, for a general component model, the general form of its differential-algebraic equation set, that is, the differential-algebraic equation set corresponding to a synchronous generator, is:

[0092] ;

[0093] Wherein, is the state variable; is the component port voltage; is the current injected by the component into the power grid; is the state variable coefficient matrix; is about , state equation function; is to calculate algebraic equation;

[0094] For the periodic function within the preset domain perform times of sampling, is the sampling period, and the discrete vector is obtained.

[0095] ;

[0096] Wherein, is the discrete vector; is the periodic function when the sampling period is 0; is the periodic function when the sampling period is ; is the periodic function when the sampling period is ; represents row column real matrix.

[0097] Define the Fourier spectrum differential matrix as wherein, represents row a real number matrix of columns, then the discrete vector of the derivative function corresponding to the derivative function of and the discrete vector have the following relationship:

[0098] ;

[0099] Substituting into the differential-algebraic equation system, the Fourier spectral method algebraic equation of discrete points can be obtained:

[0100] ;

[0101] where is the state variable coefficient matrix; is the transpose of the discrete vector i.e., the discrete transpose vector; is the Fourier spectral differential matrix the transpose of which is the Fourier spectral differential transpose matrix; is the voltage discrete matrix the transpose of which is the voltage discrete transpose matrix; is the current discrete matrix the transpose of which is the current discrete transpose matrix; T represents transpose; is the differential equation the Fourier spectral method algebraic equation of discrete points, where ; is the algebraic equation the Fourier spectral method algebraic equation of discrete points.

[0102] Next, generate the Fourier spectral differential matrix. Assume is a function with a period of and has the following Fourier series expression:

[0103] ;

[0104] ;

[0105] where is the th Fourier coefficient; is the total number of harmonics; is the Fourier coefficient vector; is the Fourier basis function vector; is the reference frequency; is the imaginary unit; is the Fourier basis function; represents conjugate; H represents conjugate transpose; is the conjugate of the th Fourier coefficient; T represents transpose.

[0106] The Fourier series expression of is denoted as , and we get:

[0107] ;

[0108] Among them, is the derivative of the Fourier series; is the th Fourier coefficient; is the total harmonic order; is the vector of Fourier coefficients; is the vector of Fourier basis functions; is the reference frequency; is the imaginary unit; is the Fourier basis function; represents conjugation; H represents conjugate transpose; is the conjugate of the th Fourier coefficient; represents transpose; is the derivative operator; is the vector of Fourier coefficients conjugate.

[0109] Sampling gives the discrete Fourier transform and inverse transform, that is, the transformation equation is;

[0110] ;

[0111] Among them, is the discrete vector; is the number of sampling times; is the vector of Fourier coefficients; H represents conjugate transpose; is the correction equation for the DC component; is the discrete matrix of Fourier basis functions, .

[0112] Then the derivative discrete point sampling is:

[0113] ;

[0114] Among them, is the discrete vector of the derivative function; is the vector of Fourier coefficients; represents conjugation; is the derivative operator; is the vector of Fourier coefficients conjugate; is the discrete matrix of Fourier basis functions, .

[0115] Since , , the derivative discrete vector can be obtained:

[0116] ;

[0117] Among them, is the derivative function discrete vector; is the discrete vector; is the number of sampling times; represents conjugate; is the derivative operator; is the correction equation for the DC component; T represents transpose; H represents conjugate transpose; is the discrete matrix of Fourier basis functions, .

[0118] After the above conversion is performed using the derivative discrete vector and the discrete vector, the Fourier spectral differential matrix obtained is:

[0119] ;

[0120] Among them, is the Fourier spectral differential matrix; is the number of sampling times; represents conjugate; is the derivative operator; is the correction equation for the DC component; T represents transpose; H represents conjugate transpose; is the discrete matrix of Fourier basis functions, .

[0121] Step 202: Perform electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time-domain model.

[0122] Furthermore, step 202 may include the following sub-steps S21 - S22:

[0123] S21: Use the operation data to construct a model to generate a synchronous generator model.

[0124] S22: Perform Park transformation on each state parameter in the synchronous generator model to generate an initial electromagnetic transient time-domain model.

[0125] In the embodiment of the present invention, the operation data is used to construct a model to generate a synchronous generator model. Since the angle between the synchronous generator rotating coordinate system and the three-phase stationary coordinate system (taking voltage as an example) is as Figure 3As shown in the figure, the conversion of each parameter in the rotating coordinate system and the corresponding parameters in the three-phase stationary coordinate system conforms to the Park transformation. When analyzing components, the 0-axis component is generally ignored. Therefore, through the Park transformation of each state parameter in the synchronous generator model, an initial electromagnetic transient time-domain model is constructed, and this initial electromagnetic transient time-domain model is used to solve the algebraic equation of the discrete-point Fourier spectrum method. Among them, the transformation matrix of the Park transformation and its inverse transformation coefficient matrix are as follows:

[0126] ;

[0127] ;

[0128] Among them, is the transformation matrix of the Park transformation; is the inverse transformation coefficient matrix of the Park transformation; is the rotor motion angle relative to the synchronous rotating coordinate system; is the rotor angular velocity relative to the synchronous rotating coordinate system.

[0129] Under the Park transformation, the variable in the coordinate system and the variable conversion equation in the

[0130] ;

[0131] Among them, is the general variable in the coordinate system, ; is the general variable in the coordinate system ;

[0132] The generator rotor motion equation is:

[0133] ;

[0134] Among them, is the rotor inertia time constant; T m is the mechanical torque; is the electromagnetic torque; is the rotor motion angle relative to the synchronous rotating coordinate system; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the rotor damping.

[0135] The synchronous generator mathematical model uses the motor parameters, and the generator convention is adopted on the stator side. Its electromagnetic transient rotor voltage equation (without considering saturation) is:

[0136] ;

[0137] wherein, is the direct-axis stator current; is the direct-axis open-circuit transient time constant; is the direct-axis open-circuit subtransient time constant; is the direct-axis transient electromotive force; is the direct-axis subtransient electromotive force; is the direct-axis transient reactance; is the direct-axis subtransient reactance; is the quadrature-axis stator current; is the quadrature-axis open-circuit transient time constant; is the quadrature-axis open-circuit subtransient time constant; is the quadrature-axis transient electromotive force; is the quadrature-axis subtransient electromotive force; is the quadrature-axis transient reactance; is the quadrature-axis subtransient reactance; is the excitation electromotive force;

[0138] Different from the electromechanical transient model, in the electromagnetic transient model, the stator current equation needs to consider the instantaneous change rate of the magnetic flux. Therefore, the stator current equation corresponding to the initial electromagnetic transient time-domain model:

[0139] ;

[0140] wherein, is the direct-axis stator current; is the direct-axis open-circuit transient time constant; is the direct-axis open-circuit subtransient time constant; is the direct-axis transient electromotive force; is the direct-axis subtransient electromotive force; is the direct-axis transient reactance; is the direct-axis subtransient reactance; is the quadrature-axis stator current; is the quadrature-axis open-circuit transient time constant; is the quadrature-axis open-circuit subtransient time constant; is the quadrature-axis transient electromotive force; is the quadrature-axis subtransient electromotive force; is the quadrature-axis transient reactance; is the quadrature-axis subtransient reactance; is the stator resistance; is the direct-axis stator voltage; is the quadrature-axis stator voltage; is the rotor motion angle relative to the synchronous rotating coordinate system; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the zero-sequence current of the stator winding; is the zero-sequence inductance of the stator winding; is the zero-sequence voltage of the stator winding.

[0141] The generator electromagnetic power equation and the generator electromagnetic torque equation corresponding to the initial electromagnetic transient time-domain model are as follows:

[0142] ;

[0143] ;

[0144] Wherein, is the generator electromagnetic power; is the generator electromagnetic torque; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the stator quadrature-axis current; is the quadrature-axis subtransient electromotive force; is the stator direct-axis current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance.

[0145] Step 203: Use the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator to calculate the initial values of each state variable in the synchronous generator, and generate a target initial value dataset.

[0146] Further, step 203 may include the following sub-steps S31-S33:

[0147] S31: Use the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator to perform calculation of the initial voltage value and the initial current value of the axis, and generate the initial voltage value of the axis and the initial current value of the

[0148] S32: Substitute the power flow data, the initial voltage value of the axis and the initial current value of the

[0149] axis into the preset state variable initial value formula to perform initialization calculation on each state variable in the initial electromagnetic transient time-domain model, and generate an initial initial value dataset. S33: Use the initial voltage value of the axis, the initial current value of the

[0150] axis and the initial initial value dataset to construct a target initial value dataset. In the embodiment of the present invention, according to the power flow data corresponding to the operation data and the power angle relationship of the synchronous generator, the initial voltage value of the axis and the initial current value of the

[0151] ;

[0152] ;

[0153] ;

[0154] ;

[0155] Among them, is the power angle at steady state; is the rotor angle at the initial moment; represents the steady-state voltage; is the steady-state voltage of the y-axis under the common rotating coordinate system; is the steady-state voltage of the x-axis under the common rotating coordinate system; represents the steady-state current; is the steady-state current of the y-axis under the common rotating coordinate system; is the steady-state current of the x-axis under the common rotating coordinate system; is the stator resistance; is the quadrature-axis reactance; is the power angle of the synchronous generator; is the steady-state voltage of the direct axis at the initial moment under the is the steady-state voltage of the quadrature axis at the initial moment under the is the steady-state current of the direct axis at the initial moment under the is the steady-state current of the quadrature axis at the initial moment under the is the steady-state voltage of the x-axis at the initial moment under the common rotating coordinate system; is the steady-state voltage of the y-axis at the initial moment under the common rotating coordinate system; is the steady-state current of the x-axis at the initial moment under the common rotating coordinate system; is the steady-state current of the y-axis at the initial moment under the common rotating coordinate system.

[0156] According to the voltage and current at the initial moment, the initial values of each state variable in the synchronous generator can be initialized in the following way, that is, substituting the power flow data, the initial value of the axis voltage and the initial value of the axis current into the following preset initial value formula of the state variable to initialize and calculate each state variable in the initial electromagnetic transient time-domain model, and generate an initial value data set. The preset initial value formula of the state variable is:

[0157] ;

[0158] Among them, is the rotor angular velocity relative to the synchronous rotating coordinate system at the initial moment; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance; is the quadrature-axis transient reactance; is the quadrature-axis subtransient reactance; is the initial value of the steady-state direct-axis subtransient electromotive force at the initial moment; is the initial value of the steady-state direct-axis subtransient electromotive force at the initial moment; is the initial value of the steady-state direct-axis transient electromotive force at the initial moment; is the initial value of the steady-state quadrature-axis transient electromotive force at the initial moment; is the initial value of the electromagnetic power; is the initial value of the mechanical power; is the initial value of the excitation electromotive force; is the direct-axis steady-state current in the is the quadrature-axis steady-state current in the is the direct-axis steady-state voltage in the is the quadrature-axis steady-state voltage in the

[0159] Step 204: Update the differential equation group of the synchronous generator in the initial electromagnetic transient time-domain model by using the target initial data set to generate a target electromagnetic transient time-domain model.

[0160] In the embodiment of the present invention, the initial values of the synchronous generator are initialized according to the power flow calculation results to obtain a target initial data set, and then the differential equation group of the synchronous generator in the initial electromagnetic transient time-domain model is updated by using the target initial data set, so that the differential terms of the differential equation group of the synchronous generator in the initial electromagnetic transient time-domain model are 0, and a target electromagnetic transient time-domain model is obtained.

[0161] Step 205: Solve the discrete points within a period of the target electromagnetic transient time-domain model by using the Fourier spectrum method to generate a target discrete transposed vector.

[0162] Further, step 205 may include the following sub-steps S41-S46:

[0163] S41: Initialize the target electromagnetic transient time-domain model by using a preset iteration accuracy and a preset iteration initial value to generate a time-domain model.

[0164] S42: Construct an algebraic equation of the discrete-point Fourier spectrum method corresponding to the time-domain model by using the state equation functions of the time-domain model corresponding to the state variables within a period.

[0165] S43. Calculate the product among the algebraic equation of the discrete point Fourier spectrum method, the Fourier spectrum differential transpose inverse matrix corresponding to the time-domain model, and the state variable coefficient inverse matrix, generate the initial discrete transpose vector, and count the number of iterations.

[0166] S44. Use the discrete vector corresponding to the initial discrete transpose vector and the historical discrete vector to judge the iteration accuracy, and generate accuracy judgment data.

[0167] S45. When the accuracy judgment data is less than or equal to the preset iteration accuracy, use the initial discrete transpose vector at the current moment as the target discrete transpose vector.

[0168] S46. When the accuracy judgment data is greater than the preset iteration accuracy, use the time-domain model corresponding to the initial discrete transpose vector as the new time-domain model, update the iteration count by 1, and jump to execute the step of calculating the product among the state variable coefficient matrix, the discrete transpose vector, and the Fourier spectrum differential transpose matrix corresponding to the time-domain model to generate the algebraic equation of the discrete point Fourier spectrum method corresponding to the time-domain model.

[0169] In the embodiment of the present invention, the Fourier spectrum method equation is written and the discrete points within the calculation period are solved. The flow chart of solving the synchronous generator electromagnetic transient model by the Fourier spectrum method is as Figure 4 shown. After calculating the Fourier spectrum differential matrix and the initial values of each variable at time 0 , a preset iteration accuracy and a preset iteration initial value are given, and is set to start iteration. The fixed-point iteration format is used for repeated iteration until accuracy convergence is achieved. That is, the target electromagnetic transient time-domain model is initialized with the preset iteration accuracy and the preset iteration initial value to generate the time-domain model. The state equation functions corresponding to each state variable within the period of the time-domain model are used to construct the algebraic equation of the discrete point Fourier spectrum method corresponding to the time-domain model. The algebraic equation of the discrete point Fourier spectrum method is:

[0170] ;

[0171] where is the state equation function with respect to , ; is the state variable; is the component port voltage; is the transpose of the discrete vector , i.e., the discrete transpose vector; is the transpose of the voltage discrete matrix , i.e., the voltage discrete transpose matrix; is the sampling period; denote row real number matrix of columns; denote transpose.

[0172] Then calculate the product among the algebraic equation of the discrete point Fourier spectrum method, the Fourier spectrum differential transpose inverse matrix corresponding to the time domain model, and the state variable coefficient inverse matrix, generate the initial discrete transpose vector and count the number of iterations. Among them, the initial discrete transpose vector is:

[0173] ;

[0174] Among them, is the algebraic equation of the discrete point Fourier spectrum method, where ; is the state variable coefficient inverse matrix; is the Fourier spectrum differential transpose inverse matrix corresponding to the time domain model.

[0175] Substitute the discrete vector corresponding to the initial discrete transpose vector and the historical discrete vector into the following formula to judge the iteration accuracy, and generate accuracy judgment data.

[0176] ;

[0177] Among them, is the discrete vector; is the discrete vector corresponding historical discrete vector; is the preset iteration accuracy.

[0178] When the accuracy judgment data is greater than the preset iteration accuracy, use the time domain model corresponding to the initial discrete transpose vector as the new time domain model, update the iteration number by 1, and jump to execute the step of calculating the product among the state variable coefficient matrix, the discrete transpose vector, and the Fourier spectrum differential transpose matrix corresponding to the time domain model to generate the algebraic equation of the discrete point Fourier spectrum method corresponding to the time domain model. When the accuracy judgment data is less than or equal to the preset iteration accuracy, use the initial discrete transpose vector at the current moment as the target discrete transpose vector.

[0179] Step 206: Update the Fourier series expression corresponding to the Fourier spectrum differential matrix with the Fourier series coefficients corresponding to the target discrete transpose vector to generate the target Fourier series expression.

[0180] In an embodiment of the present invention, the Fourier series coefficients corresponding to the discrete points calculated and solved through the target discrete transpose vector are used to obtain the target Fourier series expression of the interval solution, that is, the Fourier series expression corresponding to the Fourier spectral differential matrix is updated by using the Fourier series coefficients corresponding to the target discrete transpose vector to generate the target Fourier series expression. Thus, the variable values at any time on the entire interval can be obtained. The target Fourier series expression is as follows:

[0181] ;

[0182] Wherein, is the th Fourier coefficient; is the total number of harmonics; is the Fourier coefficient vector; is the Fourier basis function vector; is the reference frequency; is the imaginary unit; is the Fourier basis function; denotes conjugate; H denotes conjugate transpose; is the th conjugate of the Fourier coefficient; T denotes transpose.

[0183] Step 207: Use the Shannon sampling theorem and the target Fourier series expression to perform electromagnetic transient and steady-state periodic solution for the synchronous generator, and generate electromagnetic steady-state solution data.

[0184] In an embodiment of the present invention, when calculating the three-phase symmetrical steady-state periodic solution, the Fourier spectral method can obtain accurate discrete points, and the Fourier spectral method only needs very few discrete points to obtain an accurate series expression. After obtaining the target Fourier series expression, combined with the Shannon sampling theorem, perform electromagnetic transient and steady-state periodic solution for the synchronous generator. When the number of sampling points and the maximum number of harmonics meet the requirements, accurate results can be obtained, and electromagnetic steady-state solution data can be obtained.

[0185] In an embodiment of the present invention, as Figures 5 to 8 shown, the electromagnetic transient of the synchronous generator calculated by the Fourier spectral method without harmonics and the three-phase symmetrical steady-state results are as Figure 5 and Figure 6 shown, and the three-phase symmetrical steady-state results with (5th) harmonics are as Figure 7 and Figure 8As shown, the convergence accuracy is 10-14. The calculation results for different numbers of sampling points are given in the figure. It can be seen that when calculating the three-phase symmetrical steady-state periodic solution, the Fourier spectral method can obtain accurate discrete points, and the Fourier spectral method only requires fewer discrete points to obtain an accurate series expression. According to the Shannon sampling theorem, accurate results can be obtained when the number of sampling points and the maximum harmonic order satisfy the relationship. Therefore, in a synchronous generator, the Fourier spectral method can be used to calculate the steady-state periodic solution through fewer discrete points, without the need to use a large number of discrete points for electromagnetic transient simulation until steady state as in the traditional zero-start method.

[0186] Embodiment 3

[0187] Please refer to Figure 9 , Figure 9 which is the structural block diagram of a synchronous generator electromagnetic steady-state solving system based on the Fourier spectral method provided by Embodiment 3 of the present invention.

[0188] A synchronous generator electromagnetic steady-state solving system based on the Fourier spectral method provided by Embodiment 3 of the present invention includes:

[0189] A Fourier spectral differential matrix generation module 901, configured to obtain the operating data corresponding to the synchronous generator, construct a Fourier spectral differential matrix using the operating data, and generate a Fourier spectral differential matrix.

[0190] An initial electromagnetic transient time-domain model generation module 902, configured to perform electromagnetic transient simulation based on the operating data and generate an initial electromagnetic transient time-domain model.

[0191] An electromagnetic steady-state solving data generation module 903, configured to perform steady-state periodic solution based on the Fourier spectral differential matrix and the initial electromagnetic transient time-domain model, and generate electromagnetic steady-state solving data.

[0192] Optionally, the Fourier spectral differential matrix generation module 901 includes:

[0193] A differential-algebraic equation system generation module, configured to construct a differential-algebraic equation system using the operating data and generate a differential-algebraic equation system corresponding to the synchronous generator.

[0194] A discrete vector generation module, configured to sample the periodic function corresponding to the differential-algebraic equation system multiple times within a preset domain and generate a discrete vector.

[0195] A Fourier spectral method algebraic equation generation module, configured to substitute the discrete vector into the differential-algebraic equation system to construct a discrete-point Fourier spectral method algebraic equation and generate a Fourier spectral method algebraic equation.

[0196] The Fourier spectrum differential matrix generation sub-module is used to construct the Fourier spectrum differential matrix based on the algebraic equation of the Fourier spectrum method and the discrete vector, and generate the Fourier spectrum differential matrix.

[0197] Optionally, the Fourier spectrum differential matrix generation sub-module can perform the following steps:

[0198] Perform discrete Fourier transform and inverse discrete Fourier transform on the algebraic equation of the Fourier spectrum method to generate a transformed equation;

[0199] Perform derivative discrete point sampling on the transformed equation to generate a derivative discrete vector;

[0200] Use the derivative discrete vector and the discrete vector to construct a differential matrix to generate the Fourier spectrum differential matrix.

[0201] Optionally, the initial electromagnetic transient time-domain model generation module 902 can perform the following steps:

[0202] Use the operating data for model construction to generate a synchronous generator model;

[0203] Perform Park transformation on each state parameter in the synchronous generator model to generate the initial electromagnetic transient time-domain model;

[0204] The stator current equation corresponding to the initial electromagnetic transient time-domain model is:

[0205] ;

[0206] Where, is the direct-axis stator current; is the direct-axis open-circuit transient time constant; is the direct-axis open-circuit sub-transient time constant; is the direct-axis transient electromotive force; is the direct-axis sub-transient electromotive force; is the direct-axis transient reactance; is the direct-axis sub-transient reactance; is the quadrature-axis stator current; is the quadrature-axis open-circuit transient time constant; is the quadrature-axis open-circuit sub-transient time constant; is the quadrature-axis transient electromotive force; is the quadrature-axis sub-transient electromotive force; is the quadrature-axis transient reactance; is the quadrature-axis sub-transient reactance; is the stator resistance; is the direct-axis stator voltage; is the quadrature-axis stator voltage; is the rotor motion angle relative to the synchronous rotating coordinate system; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the zero-sequence current of the stator winding; is the zero-sequence inductance of the stator winding; is the zero-sequence voltage of the stator winding;

[0207] The generator electromagnetic power equation corresponding to the initial electromagnetic transient time-domain model is:

[0208] ;

[0209] where is the generator electromagnetic power; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the stator quadrature-axis current; is the quadrature-axis subtransient electromotive force; is the stator direct-axis current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance;

[0210] The generator electromagnetic torque equation corresponding to the initial electromagnetic transient time-domain model is:

[0211] ;

[0212] where is the generator electromagnetic torque; is the stator quadrature-axis current; is the quadrature-axis subtransient electromotive force; is the stator direct-axis current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance.

[0213] Optionally, the electromagnetic steady-state solution data generation module 903 includes:

[0214] A target initial value data set generation module, which is used to calculate the initial values of the state variables in the synchronous generator by using the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator, and generate a target initial value data set.

[0215] A target electromagnetic transient time-domain model generation module, which is used to update the synchronous generator differential equation set in the initial electromagnetic transient time-domain model by using the target initial data set, and generate a target electromagnetic transient time-domain model.

[0216] A target discrete transposed vector generation module, which is used to solve the discrete points within the period of the target electromagnetic transient time-domain model by using the Fourier spectrum method, and generate a target discrete transposed vector.

[0217] A target Fourier series expression generation module, which is used to update the Fourier series expression corresponding to the Fourier spectrum differential matrix with the Fourier series coefficients corresponding to the target discrete transposed vector, and generate a target Fourier series expression.

[0218] An electromagnetic steady-state solution data generation sub-module, which is used to perform electromagnetic transient solution steady-state periodic solution on a synchronous generator by using the Shannon sampling law and the target Fourier series expression, and generate electromagnetic steady-state solution data.

[0219] Optionally, the target initial value data set generation module can perform the following steps:

[0220] Use the power flow data corresponding to the operating data and the power angle relationship data corresponding to the synchronous generator to calculate the initial values of the voltage and current of the shaft, and generate the initial value of the shaft voltage and the initial value of the shaft current;

[0221] Substitute the power flow data, the initial value of the shaft voltage and the initial value of the shaft current into the preset initial value formula of the state variable to perform initialization calculation on each state variable in the initial electromagnetic transient time domain model, and generate an initial initial value data set;

[0222] Use the initial value of the shaft voltage, the initial value of the shaft current and the initial initial value data set to construct a target initial value data set.

[0223] Optionally, the target discrete transposed vector generation module can perform the following steps:

[0224] Initialize the target electromagnetic transient time domain model with a preset iteration accuracy and a preset iteration initial value to generate a time domain model;

[0225] Use the state equation functions corresponding to each state variable in the time domain model within a period to construct an algebraic equation of the discrete point Fourier spectrum method corresponding to the time domain model;

[0226] Calculate the product of the algebraic equation of the discrete point Fourier spectrum method, the Fourier spectrum differential transposed inverse matrix corresponding to the time domain model, and the state variable coefficient inverse matrix, generate an initial discrete transposed vector and count the number of iterations;

[0227] Use the discrete vector corresponding to the initial discrete transposed vector and the historical discrete vector to perform iteration accuracy judgment to generate accuracy judgment data;

[0228] When the accuracy judgment data is less than or equal to the preset iteration accuracy, use the initial discrete transposed vector at the current moment as the target discrete transposed vector;

[0229] When the accuracy judgment data is greater than the preset iteration accuracy, use the time-domain model corresponding to the initial discrete transpose vector as the new time-domain model, increment the iteration count by 1, and jump to execute the step of calculating the product between the state variable coefficient matrix, discrete transpose vector, and Fourier spectrum differential transpose matrix corresponding to the time-domain model to generate the algebraic equation of the discrete point Fourier spectrum method corresponding to the time-domain model.

[0230] Embodiment 4

[0231] Please refer to Figure 10 , Figure 10 which is a structural block diagram of an electronic device provided in Embodiment 4 of the present invention.

[0232] An electronic device according to an embodiment of the present invention includes a memory 1001 and a processor 1002, and a computer program is stored in the memory 1001; when the computer program is executed by the processor 1002, the processor 1002 is caused to execute the grounding grid topology determination method according to any of the above embodiments.

[0233] The memory 1001 may be an electronic memory such as a flash memory, EEPROM (electrically erasable programmable read-only memory), EPROM, hard disk, or ROM. The memory 1001 has a storage space 1003 for program code 1013 for executing any method step in the above method. For example, the storage space 1003 for program code may include respective program codes 1013 for implementing various steps in the above method. These program codes may be read from or written to one or more computer program products. These computer program products include program code carriers such as hard disks, compact discs (CDs), memory cards, or floppy disks. The program code may be compressed in an appropriate form. When these codes are run by a computing processing device, the computing processing device is caused to execute the respective steps in the synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method described above.

[0234] The embodiment of the present invention also provides a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, it implements the synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method according to any of the above embodiments.

[0235] Those skilled in the art can clearly understand that for the convenience and brevity of description, the specific working processes of the above-described systems, devices, and units can refer to the corresponding processes in the foregoing method embodiments, and will not be described herein again.

[0236] In several embodiments provided by the present application, it should be understood that the disclosed systems, devices, and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of units is only a logical function division. In actual implementation, there may be other division methods. For example, multiple units or components can be combined or integrated into another system, or some features can be ignored or not executed. Another point is that the displayed or discussed couplings or direct couplings or communication connections between each other can be through some interfaces. The indirect couplings or communication connections of devices or units can be in electrical, mechanical, or other forms.

[0237] The units described as separate components may or may not be physically separated. The components shown as units may or may not be physical units, that is, they can be located in one place, or they can be distributed to multiple network units. Some or all of the units can be selected according to actual needs to achieve the purpose of the solution of this embodiment.

[0238] In addition, in each embodiment of the present invention, the functional units can be integrated in a processing unit, or each unit can exist physically alone, or two or more units can be integrated in one unit. The above-mentioned integrated units can be implemented in the form of hardware or in the form of software functional units.

[0239] If the integrated unit is implemented in the form of a software functional unit and sold or used as an independent product, it can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, in essence, or the part that contributes to the prior art, or all or part of this technical solution, can be embodied in the form of a software product. This computer software product is stored in a storage medium and includes several instructions to enable a computer device (which can be a personal computer, a server, or a network device, etc.) to execute all or part of the steps of the methods in each embodiment of the present invention. The foregoing storage medium includes: USB flash drives, mobile hard disks, read-only memories (ROMs), random access memories (RAMs), magnetic disks, or optical discs, and other media that can store program codes.

[0240] The above embodiments are only used to illustrate the technical solutions of the present invention, rather than to limit them; although the present invention has been described in detail with reference to the foregoing embodiments, those of ordinary skill in the art should understand that they can still modify the technical solutions recorded in the foregoing embodiments, or perform equivalent replacements for some of the technical features; and these modifications or replacements do not cause the essence of the corresponding technical solutions to deviate from the spirit and scope of the technical solutions of each embodiment of the present invention.

Claims

1. A method for solving the electromagnetic steady state of a synchronous generator based on the Fourier spectrum method, characterized in that Including: Obtain the operation data corresponding to the synchronous generator, and use the operation data to construct a Fourier spectrum differential matrix to generate a Fourier spectrum differential matrix; Perform electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time-domain model; Perform steady-state period solution based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to generate electromagnetic steady-state solution data; The step of using the operation data to construct a Fourier spectrum differential matrix to generate a Fourier spectrum differential matrix includes: Use the operation data to construct a differential-algebraic equation set to generate the differential-algebraic equation set corresponding to the synchronous generator; Sample the periodic function corresponding to the differential-algebraic equation set multiple times within a preset domain to generate a discrete vector; Substitute the discrete vector into the differential-algebraic equation set to construct an algebraic equation of the discrete-point Fourier spectrum method to generate an algebraic equation of the Fourier spectrum method; Construct a Fourier spectrum differential matrix based on the algebraic equation of the Fourier spectrum method and the discrete vector to generate a Fourier spectrum differential matrix; The step of performing steady-state period solution based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to generate electromagnetic steady-state solution data includes: Calculate the initial values of each state variable in the synchronous generator using the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator to generate a target initial value data set; Update the synchronous generator differential equation set in the initial electromagnetic transient time-domain model using the target initial value data set to generate a target electromagnetic transient time-domain model; Perform discrete-point solution within a period of the target electromagnetic transient time-domain model using the Fourier spectrum method to generate a target discrete transpose vector; Update the Fourier series expression corresponding to the Fourier spectrum differential matrix using the Fourier series coefficients corresponding to the target discrete transpose vector to generate a target Fourier series expression; Perform electromagnetic transient steady-state period solution for the synchronous generator using the Shannon sampling theorem and the target Fourier series expression to generate electromagnetic steady-state solution data.

2. The electromagnetic steady-state solution method of a synchronous generator based on the Fourier spectrum method according to claim 1, wherein The step of constructing a Fourier spectrum differential matrix based on the algebraic equation of the Fourier spectrum method and the discrete vector to generate a Fourier spectrum differential matrix includes: Perform discrete Fourier transform and inverse discrete Fourier transform on the algebraic equation of the Fourier spectrum method to generate a transformed equation; Perform derivative discrete-point sampling on the transformed equation to generate a derivative discrete vector; Construct a differential matrix using the derivative discrete vector and the discrete vector to generate a Fourier spectrum differential matrix.

3. The electromagnetic steady-state solution method for a synchronous generator based on the Fourier spectrum method according to claim 1, wherein The step of performing electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time-domain model includes: Use the operation data to construct a model to generate a synchronous generator model; Perform Park transformation on each state parameter in the synchronous generator model to generate an initial electromagnetic transient time-domain model; The stator current equation corresponding to the initial electromagnetic transient time-domain model is: ； Wherein, is the direct-axis stator current; is the direct-axis open-circuit transient time constant; is the direct-axis open-circuit subtransient time constant; is the direct-axis transient electromotive force; is the direct-axis subtransient electromotive force; is the direct-axis transient reactance; is the direct-axis subtransient reactance; is the quadrature-axis stator current; is the quadrature-axis open-circuit transient time constant; is the quadrature-axis open-circuit subtransient time constant; is the quadrature-axis transient electromotive force; is the quadrature-axis subtransient electromotive force; is the quadrature-axis transient reactance; is the quadrature-axis subtransient reactance; is the stator resistance; is the direct-axis stator voltage; is the quadrature-axis stator voltage; is the rotor motion angle relative to the synchronous rotating coordinate system; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the zero-sequence current of the stator winding; is the zero-sequence inductance of the stator winding; is the zero-sequence voltage of the stator winding; The generator electromagnetic power equation corresponding to the initial electromagnetic transient time-domain model is: ； Among them, is the electromagnetic power of the generator; is the rotor angular velocity relative to the synchronous rotating coordinate system; is the stator quadrature-axis current; is the quadrature-axis subtransient electromotive force; is the stator direct-axis current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance; The generator electromagnetic torque equation corresponding to the initial electromagnetic transient time-domain model is: ； Among them, is the electromagnetic torque of the generator; is the stator quadrature-axis current; is the quadrature-axis subtransient electromotive force; is the stator direct-axis current; is the direct-axis subtransient electromotive force; is the quadrature-axis subtransient reactance; is the direct-axis subtransient reactance.

4. The electromagnetic steady-state solution method of a synchronous generator based on the Fourier spectrum method according to claim 1, characterized in that, The step of calculating the initial values of each state variable in the synchronous generator by using the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator to generate a target initial value data set includes: Using the power flow data corresponding to the said operating data and the power angle relationship data corresponding to the said synchronous generator to perform calculation of the initial voltage value and the initial current value of the axis to generate the initial value of the axis voltage and the Substitute the tidal current data, the initial value of the shaft voltage, and the initial value of the shaft current into the preset initial value formula of the state variables to perform initialization calculations on each state variable in the initial electromagnetic transient time-domain model, and generate an initial initial value data set; Using the said initial shaft voltage, the said initial shaft current, and the said initial dataset, construct a target initial dataset.

5. The electromagnetic steady-state solution method of a synchronous generator based on the Fourier spectrum method according to claim 1, wherein The step of using the Fourier spectrum method to solve the discrete points within a period of the target electromagnetic transient time-domain model to generate a target discrete transposed vector includes: Initializing the target electromagnetic transient time-domain model with a preset iteration accuracy and a preset initial value to generate a time-domain model; Using the state equation functions corresponding to each state variable within a period of the time-domain model to construct a discrete-point Fourier spectrum method algebraic equation corresponding to the time-domain model; Calculating the product among the discrete-point Fourier spectrum method algebraic equation, the Fourier spectrum differential transposed inverse matrix corresponding to the time-domain model, and the state variable coefficient inverse matrix to generate an initial discrete transposed vector and counting the number of iterations; Judging the iteration accuracy by using the discrete vector corresponding to the initial discrete transposed vector and the historical discrete vector to generate accuracy judgment data; When the accuracy judgment data is less than or equal to the preset iteration accuracy, using the initial discrete transposed vector at the current moment as the target discrete transposed vector; When the accuracy judgment data is greater than the preset iteration accuracy, using the time-domain model corresponding to the initial discrete transposed vector as a new time-domain model, incrementing the updated iteration count by 1, and jumping to execute the step of calculating the product among the state variable coefficient matrix, the discrete transposed vector, and the Fourier spectrum differential transposed matrix corresponding to the time-domain model to generate the discrete-point Fourier spectrum method algebraic equation corresponding to the time-domain model.

6. A synchronous generator electromagnetic steady-state solution system based on the Fourier spectrum method, characterized in that, Includes: A Fourier spectrum differential matrix generation module, configured to obtain operation data corresponding to a synchronous generator, and construct a Fourier spectrum differential matrix by using the operation data to generate a Fourier spectrum differential matrix; An initial electromagnetic transient time-domain model generation module, configured to perform electromagnetic transient simulation based on the operation data to generate an initial electromagnetic transient time-domain model; An electromagnetic steady-state solution data generation module, configured to perform steady-state period solution based on the Fourier spectrum differential matrix and the initial electromagnetic transient time-domain model to generate electromagnetic steady-state solution data; The Fourier spectrum differential matrix generation module includes: A differential-algebraic equation system generation module, configured to construct a differential-algebraic equation system corresponding to the synchronous generator by using the operation data to generate the differential-algebraic equation system corresponding to the synchronous generator; A discrete vector generation module, configured to sample a periodic function corresponding to the differential-algebraic equation system multiple times within a preset domain to generate a discrete vector; A Fourier spectrum method algebraic equation generation module, configured to substitute the discrete vector into the differential-algebraic equation system to construct a discrete-point Fourier spectrum method algebraic equation to generate a Fourier spectrum method algebraic equation; A Fourier spectrum differential matrix generation sub-module, configured to construct a Fourier spectrum differential matrix based on the Fourier spectrum method algebraic equation and the discrete vector to generate a Fourier spectrum differential matrix; The electromagnetic steady-state solution data generation module includes: A target initial value data set generation module, configured to calculate initial values of each state variable in the synchronous generator by using the power flow data corresponding to the operation data and the power angle relationship data corresponding to the synchronous generator, and generate a target initial value data set; A target electromagnetic transient time-domain model generation module, configured to update the synchronous generator differential equation set in the initial electromagnetic transient time-domain model by using the target initial value data set, and generate a target electromagnetic transient time-domain model; A target discrete transposed vector generation module, configured to solve discrete points within a period of the target electromagnetic transient time-domain model by using the Fourier spectrum method, and generate a target discrete transposed vector; A target Fourier series expression generation module, configured to update the Fourier series expression corresponding to the Fourier spectrum differential matrix by using the Fourier series coefficients corresponding to the target discrete transposed vector, and generate a target Fourier series expression; An electromagnetic steady-state solution data generation sub-module, configured to perform electromagnetic transient solution steady-state period solution on the synchronous generator by using the Shannon sampling theorem and the target Fourier series expression, and generate electromagnetic steady-state solution data.

7. An electronic device, characterized in that, It includes a memory and a processor. When the computer program stored in the memory is executed by the processor, the processor is caused to execute the steps of the synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method according to any one of claims 1 to 5.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that, When the computer program is executed, it implements the synchronous generator electromagnetic steady-state solution method based on the Fourier spectrum method according to any one of claims 1 to 5.

Images