Synchronous generator parameter solving method, computer storage medium and control equipment

Abstract

The invention relates to a synchronous generator parameter solving method, a computer storage medium and control equipment. The method comprises the following steps: acquiring reactance design parameters, a potential equation and a flux linkage equation of a synchronous generator; calculating a power grid voltage and a power grid current according to the reactance design parameters and a preset node current vector; respectively performing dq decomposition on the power grid voltage and the power grid current to obtain stator voltage component data and stator current component data; performing Park transformation on the potential equation and the flux linkage equation to obtain a Park equation; determining the electromotive force of the synchronous generator according to the reactance design parameters, the stator voltage component data, the stator current component data and a Park equation; and calculating electromotive force of the synchronous generator according to the stator voltage component data so as to determine output parameters of the synchronous generator according to the electromotive force. According to the invention, the angular velocity and related parameters of the synchronous generator can be efficiently solved, the solving efficiency is high, and the scheme realization difficulty is low.

Description

Technical Field

[0001] The present invention relates to the technical field of synchronous generators, and in particular to a synchronous generator parameter solving method, a computer storage medium and a control device. Background Art

[0002] Almost all AC generators in modern hydroelectric power stations, hydroelectric power stations and nuclear power stations are synchronous generators, which are the most important dynamic equipment in the power system. Although the original mathematical equation of the synchronous generator can describe its basic structure and operating characteristics, the parameters required by the original mathematical equation cannot be obtained directly. It requires practical testing and other work to obtain them. Not only is the acquisition efficiency low, but the acquisition method is also difficult. At present, there is an urgent need for a solution that can use simulation technology to efficiently solve the parameters of synchronous generators. Summary of the invention

[0003] The technical problem to be solved by the present invention is to provide a synchronous generator parameter solving method, a computer storage medium and a control device suitable for large-scale power grid node simulation calculation.

[0004] The technical solution adopted by the present invention to solve the technical problem is: constructing a synchronous generator parameter solution method for large-scale power grid node simulation calculation, the synchronous generator parameter solution method comprising:

[0005] Obtain the reactance design parameters of the synchronous generator, the original potential equation of the synchronous generator and the original flux equation of the synchronous generator;

[0006] Calculating the grid voltage and grid current according to the reactance design parameters and the preset node current vector;

[0007] Performing dq decomposition on the grid voltage and the grid current respectively to obtain stator voltage component data and stator current component data;

[0008] Performing Park transformation on the original potential equation of the synchronous generator and the original flux equation of the synchronous generator to obtain the Park equation;

[0009] The angular velocity of the synchronous generator is determined according to the reactance design parameters, the stator voltage component data, the stator current component data and the Park equation, and the electromotive force of the synchronous generator is calculated according to the stator voltage component data to determine the output parameters of the synchronous generator according to the electromotive force.

[0010] Preferably, the reactance design parameters include armature winding resistance, d-axis component of stator circuit self-inductance reactance, q-axis component of stator circuit self-inductance reactance, d-axis transient reactance, q-axis transient reactance, d-axis transient open-circuit time constant and q-axis transient open-circuit time constant;

[0011] The step of calculating the grid voltage and grid current according to the reactance design parameters and the preset node current vector comprises:

[0012] Determine a node admittance matrix according to the reactance design parameters, and calculate a grid voltage and a grid current according to the node current vector, the node admittance matrix and the node voltage matrix equation;

[0013] The node voltage matrix equation is expressed as: YU=I, where Y represents the node admittance matrix, U represents the node voltage vector, and I represents the node injection current vector.

[0014] Preferably, the step of determining the node admittance matrix according to the reactance design parameters comprises:

[0015] Obtain the equivalent circuit of the synchronous generator and the power system;

[0016] Based on the equivalent circuit of the synchronous generator and the power system, the reactance design parameters are equivalently transformed into the external power system network impedance, and the power system network impedance is integrated to obtain a node admittance matrix.

[0017] Preferably, the Parker equation includes the Parker potential equation and the Parker flux linkage equation;

[0018] The Parker potential equation is expressed as:

[0019]

[0020] The Parker flux equation is expressed as:

[0021]

[0022] Among them, u d 、u q and u0 represent the d-axis component, q-axis component and 0-axis component of the stator voltage respectively, u f represents the excitation voltage, R a Represents the armature winding resistance, R f Represents the excitation resistance, R D and R Q represents the d-axis resistance and q-axis resistance of the rotor equivalent damping winding, i d 、i q and i0 represent the d-axis component, q-axis component and 0-axis component of the instantaneous stator current respectively, i f represents the excitation current, i D and i Q represents the d-axis component and q-axis component of the rotor equivalent damping winding current, p represents the differential operator, ψ d , qand ψ0 represent the d-axis component, q-axis component and 0-axis component of the stator circuit flux, respectively. f , ψ D and ψ Q represents the d-axis component and q-axis component of the rotor circuit flux, ω represents the angular velocity, x d 、x q and x0 represent the d-axis component, q-axis component and 0-axis component of the stator circuit self-inductance reactance respectively, and x f represents the self-inductance reactance of the excitation winding, x D and x Q represents the d-axis component and q-axis component of the rotor circuit self-inductance reactance, respectively, af 、x aD 、x a Q 、x fa 、x fD 、x Da 、x Df and x Qa They represent the mutual inductance reactance between different windings.

[0023] Preferably, the stator voltage component data includes a d-axis component and a q-axis component of the stator voltage, and the stator current component data includes a d-axis component and a q-axis component of the stator current;

[0024] The step of determining the angular velocity of the synchronous generator according to the reactance design parameter, the stator voltage component data, the stator current component data and the Parker equation comprises:

[0025] The reactance design parameters, the stator voltage component data, and the stator current component data are substituted into the Parker equation to calculate the angular velocity.

[0026] Preferably, the output parameter of the synchronous generator includes the terminal current of the synchronous generator;

[0027] The expression of the terminal current is:

[0028] I t represents the terminal current, P t represents the preset end active power of the synchronous generator, Q t represents the preset end reactive power of the synchronous generator, E t represents the electromotive force of the synchronous generator.

[0029] Preferably, the output parameter of the synchronous generator also includes the power factor angle of the synchronous generator; the expression of the power factor angle is:

[0030] φ=cos ―1 [Pt / (E t I t )];

[0031] v represents the power factor angle.

[0032] Preferably, the output parameter of the synchronous generator also includes the power angle of the synchronous generator;

[0033] The expression of the power angle is:

[0034] δ i =tan ―1 [(x q I t cosφ―R a I t sinφ) / (E t +R a I t cosφ+x q I t sinφ)];

[0035] δ i represents the power angle, x q Represents the q-axis component of the stator circuit self-inductance, R a Represents the armature winding resistance.

[0036] The present invention also constructs a computer storage medium storing a computer program, which implements the steps of the method for solving the parameters of a synchronous generator described above when executed by a processor.

[0037] The present invention also constructs a control device, including a processor and a memory storing a computer program, wherein the processor implements the steps of the synchronous generator parameter solution method described above when executing the computer program.

[0038] The implementation of the present invention has the following beneficial effects: it provides a method for solving the parameters of a synchronous generator suitable for large-scale power grid node simulation calculations, and solves the angular velocity and related parameters of the synchronous generator through simulation technology, which not only has high solving efficiency but also has the advantage of low difficulty in implementing the solution. BRIEF DESCRIPTION OF THE DRAWINGS

[0039] The present invention will be further described below with reference to the accompanying drawings and embodiments, in which:

[0040] Figure 1 is a program flow chart of a method for solving parameters of a synchronous generator in some embodiments of the present invention;

[0041] Figure 2 is a schematic diagram of the structure of the original mathematical model of the synchronous generator in some embodiments of the present invention;

[0042] Figure 3 is an equivalent circuit of a synchronous generator and a power system in some embodiments of the present invention;

[0043] Figure 4 is a circuit structure block diagram of a control device in some embodiments of the present invention. DETAILED DESCRIPTION

[0044] In order to have a clearer understanding of the technical features, purposes and effects of the present invention, specific embodiments of the present invention are now described in detail with reference to the accompanying drawings.

[0045] It should be noted that the flowcharts shown in the accompanying drawings are only exemplary and do not necessarily include all the contents and operations / steps, nor do they have to be executed in the order described. For example, some operations / steps can be decomposed, while some operations / steps can be combined or partially combined, so the actual execution order may change according to actual conditions.

[0046] The block diagrams shown in the accompanying drawings are merely functional entities and do not necessarily correspond to physically independent entities. That is, these functional entities may be implemented in software form, or in one or more hardware modules or integrated circuits, or in different networks and / or processor devices and / or microcontroller devices.

[0047] Figure 1 It is a program flow chart of a method for solving synchronous generator parameters in some embodiments of the present invention. The method is suitable for large-scale grid node simulation calculation and can help staff solve the angular velocity, terminal current, power factor angle, power angle and other parameters of the synchronous generator. The method may include step S10, step S20, step S30, step S40, step S50 and step S60.

[0048] Step S10 includes: obtaining the reactance design parameters of the synchronous generator, the original potential equation of the synchronous generator, and the original flux equation of the synchronous generator.

[0049] In some embodiments, the reactance design parameters may include armature winding resistance, d-axis component of stator circuit self-inductance reactance, q-axis component of stator circuit self-inductance reactance, d-axis transient reactance, q-axis transient reactance, d-axis transient open-circuit time constant, and q-axis transient open-circuit time constant.

[0050] The armature winding resistance, the d-axis component of the stator circuit self-inductance reactance, and the q-axis component of the stator circuit self-inductance reactance are usually known parameters, so they can be determined by looking up the actual motor parameter table of the synchronous generator whose parameters are to be determined. Of course, they can also be measured using existing methods. Taking the d-axis component of the stator circuit self-inductance reactance and the q-axis component of the stator circuit self-inductance reactance as examples, they can be tested on the synchronous generator whose parameters are to be determined by LCR measurement method, high-frequency injection method, pulse voltage method, etc.

[0051] The d-axis transient reactance and q-axis transient reactance can be calculated by the following equations: x ′ d represents the d-axis transient reactance, x d represents the d-axis component of the stator circuit self-inductance, x ad represents the armature reaction reactance of the d-axis component, x f represents the self-inductance reactance of the excitation winding, x ′ q represents the q-axis transient reactance, x q represents the q-axis component of the stator circuit self-inductance, x ag represents the armature reaction reactance of the q-axis component, x g represents the armature reaction related reactance. ad 、x ag and x f Usually, the parameters are known and can be determined by looking up the actual motor parameter table of the synchronous generator whose parameters are to be determined. g It can be obtained by testing the synchronous generator with the parameters to be determined by LCR measurement method, etc. ad 、x ag and x f It can also be obtained by testing the synchronous generator whose parameters are to be determined. It should be noted that in the present invention, the d-axis represents the direct axis in the dq coordinate system, and the q-axis represents the quadrature axis in the dq coordinate system.

[0052] The d-axis transient open-circuit time constant and the q-axis transient open-circuit time constant can be calculated by the following equations: T ′ d0 represents the d-axis transient open-circuit time constant, T ′ q0 Represents the q-axis transient open-circuit time constant, R f Represents the field winding resistance, R g Represents the armature reaction related resistance, where R f The R can be determined by checking the actual motor parameter table of the synchronous generator whose parameters are to be determined, or by testing the synchronous generator whose parameters are to be determined by LCR measurement method or other methods. f and Rg .

[0053] In this step, the original potential equation of the synchronous generator and the original flux equation of the synchronous generator refer to the potential equation and the flux equation of the synchronous generator in the abc coordinate system, respectively.

[0054] In some embodiments, according to Figure 2 The original mathematical model of the synchronous generator shown in the figure shows that the potential equation in the abc coordinate system can be expressed as:

[0055]

[0056] Among them, u a 、u b and u c represents the voltage of phase a, phase b and phase c respectively, u f represents the excitation voltage, and They represent the time derivatives of the stator flux of phase a, phase b and phase c respectively, represents the time derivative of the field winding flux, and represents the derivative of the d-axis flux with respect to time and the derivative of the q-axis flux with respect to time of the rotor equivalent damping loop in the dq coordinate system, r represents the resistance of each phase, r f Represents the field winding resistance, r D and r Q represents the d-axis resistance and q-axis resistance of the rotor equivalent damping circuit in the dq coordinate system, i a 、i b and i c represents the current of stator phase a, phase b and phase c respectively, i f represents the excitation current, i D and i Q They represent the d-axis current and q-axis current of the rotor equivalent damping loop in the dq coordinate system respectively.

[0057] The magnetic flux equation in the abc coordinate system can be expressed as:

[0058]

[0059] Among them, ψ a , b and ψ c represents the flux of phase a, phase b and phase c respectively, ψ f represents the excitation flux, ψ D and ψ Q They represent the d-axis magnetic flux and the q-axis magnetic flux in the dq coordinate system respectively. qwIt represents the self-inductance or mutual inductance coefficient between the q winding and the w winding. q and w can be any letter from a, b, c, f, D and Q. L aa For example, L aa represents the self-inductance coefficient of phase a winding, and L af It represents the mutual inductance coefficient between the a-phase winding and the excitation winding.

[0060] Step S20 includes: calculating the grid voltage and grid current according to the reactance design parameters and the preset node current vector.

[0061] In some embodiments, the step of calculating the grid voltage and grid current based on the reactance design parameters and the preset node current vector may include: determining the node admittance matrix based on the reactance design parameters, and calculating the grid voltage and grid current based on the node current vector, the node admittance matrix and the node voltage matrix equation; wherein the node voltage matrix equation is expressed as: YU=I, Y represents the node admittance matrix, U represents the node voltage vector, and I represents the node injection current vector.

[0062] In this embodiment, the node voltage vector includes the voltage of each reference node in the power system, such as the grid voltage on the power system side. The node injection current vector includes the current of each reference node in the power system, such as the grid current on the power system side. Therefore, by substituting the node current vector and the node admittance matrix into the node voltage matrix equation, the node voltage vector can be calculated, and the grid voltage and grid current can be extracted. The node admittance matrix includes the admittance of each reference node in the power system. It should be noted that the node injection current vector is usually a parameter designed in advance by the staff based on the load that the synchronous generator with the parameters to be determined actually needs to drive, and is a known parameter.

[0063] In some embodiments, the step of determining the node admittance matrix based on the reactance design parameters may include: obtaining an equivalent circuit of the synchronous generator and the power system; based on the equivalent circuit of the synchronous generator and the power system, transforming the reactance design parameters into the external power system network impedance, and integrating the power system network impedance to obtain the node admittance matrix.

[0064] See also Figure 3 , Figure (a) corresponds to the equivalent circuit of the synchronous generator, X S represents the equivalent reactance of the synchronous generator, R a represents the equivalent resistance of the synchronous generator, and Figure (b) corresponds to the equivalent circuit of the power system. Represents the power system network impedance. This embodiment is based on Figure 3The circuit converts the reactance design parameters into the power system network impedance, and then converts the power system network impedance into the admittance form to obtain the node admittance matrix. It should be noted that the power system network impedance includes the impedance of each reference node in the power system, such as the grid impedance on the power system side.

[0065] Step S30 includes: performing dq decomposition on the grid voltage and the grid current respectively to obtain stator voltage component data and stator current component data. In the step, the grid voltage is decomposed by dq to obtain stator voltage component data related to the dq coordinate system, and the grid current is decomposed by dq to obtain stator current component data related to the dq coordinate system. Among them, the stator voltage component data may include the d-axis component, q-axis component and 0-axis component of the stator voltage, and the stator current component data may include the d-axis component, q-axis component and 0-axis component of the instantaneous stator current.

[0066] It should be noted that dq decomposition is also called dq transformation, which refers to converting the variables in the three-phase system into two orthogonal axis (d-axis and q-axis) components. The specific process can refer to the existing methods and will not be repeated here.

[0067] Step S40 includes: performing Park transformation on the original potential equation of the synchronous generator and the original flux equation of the synchronous generator to obtain the Park equation.

[0068] Specifically, the Parker equation may include the Parker potential equation and the Parker flux equation. The Parker potential equation obtained by performing the Parker transformation on the potential equation in the abc coordinate system can be expressed as:

[0069]

[0070] Among them, u d 、u q and u0 represent the d-axis component, q-axis component and 0-axis component of the stator voltage respectively, u f represents the excitation voltage, R a Represents the armature winding resistance, R f Represents the excitation resistance, R D and R Q represents the d-axis resistance and q-axis resistance of the rotor equivalent damping winding, i d 、i q and i0 represent the d-axis component, q-axis component and 0-axis component of the instantaneous stator current respectively, i f represents the excitation current, i D and i Q represents the d-axis component and q-axis component of the rotor equivalent damping winding current, p represents the differential operator, ψ d , qand ψ0 represent the d-axis component, q-axis component and 0-axis component of the stator circuit flux, respectively. f , ψ D and ψ Q They represent the d-axis component and q-axis component of the rotor circuit flux respectively, and ω represents the angular velocity.

[0071] The Parker magnetic flux equation obtained by performing Parker transformation on the magnetic flux equation in the abc coordinate system can be expressed as:

[0072]

[0073] Among them, x d 、x q and x0 represent the d-axis component, q-axis component and 0-axis component of the stator circuit self-inductance reactance respectively, and x f represents the self-inductance reactance of the excitation winding, x D and x Q represents the d-axis component and q-axis component of the rotor circuit self-inductance reactance, respectively, af 、x aD 、x a Q 、x fa 、x fD 、x Da 、x Df and x Qa They represent the mutual inductance between different windings. Specifically, x af represents the mutual inductive reactance from the armature winding to the field winding, x aD represents the mutual inductance reactance of the armature winding to the rotor circuit d-axis component, x aQ represents the mutual inductance reactance of the q-axis component of the armature winding to the rotor circuit, x fa represents the mutual inductance from the field winding to the armature winding, x fD represents the mutual inductive reactance of the d-axis component from the field winding to the rotor circuit, x Da represents the mutual inductive reactance from the d-axis component of the rotor circuit to the armature winding, x Df represents the mutual inductive reactance from the d-axis component of the rotor circuit to the excitation winding, x Qa It represents the mutual inductive reactance from the q-axis component of the rotor circuit to the armature winding.

[0074] Step S50 includes: determining the angular velocity of the synchronous generator according to the reactance design parameters, the stator voltage component data, the stator current component data and the Park equation.

[0075] In some embodiments, the step of determining the angular velocity of the synchronous generator based on the reactance design parameters, stator voltage component data, stator current component data and the Park equation may include: substituting the reactance design parameters, stator voltage component data, and stator current component data into the Park equation to calculate the excitation voltage and excitation current, and then calculating the angular velocity of the synchronous generator.

[0076] In this embodiment, as long as the reactance design parameters, stator voltage component data, and stator current component data are substituted into the Parker equation, the angular velocity ω can be calculated by solving the matrix.

[0077] Step S60 includes: calculating the electromotive force of the synchronous generator according to the stator voltage component data, so as to determine the output parameters of the synchronous generator according to the electromotive force. In this step, the expression of the electromotive force can be: u d and u q The stator voltage component data includes the d-axis component and the q-axis component of the stator voltage, E t Represents the electromotive force of the synchronous generator.

[0078] In some embodiments, the output parameter of the synchronous generator may include the terminal current of the synchronous generator. Accordingly, the expression of the terminal current may be: I t Indicates terminal current, P t Indicates the preset end active power of the synchronous generator, Q t Indicates the preset end reactive power of the synchronous generator.

[0079] In some embodiments, the output parameter of the synchronous generator may also include the power factor angle of the synchronous generator. Accordingly, the expression of the power factor angle may be: φ = cos ―1 [P t / (E t I t )], φ represents the power factor angle.

[0080] In some embodiments, the output parameter of the synchronous generator may also include the power angle of the synchronous generator. Accordingly, the expression of the power angle may be: δ i =tan ―1 [(x q I t cosφ―R a I t sinφ) / (E t +R a I t cosφ+x q I t sinφ)],δ i represents the power angle, xq Represents the q-axis component of the stator circuit self-inductance, R a Represents the armature winding resistance.

[0081] Since the output of the synchronous generator changes dynamically, the stator voltage component data and the stator current component data will change with the angular velocity. If it is necessary to determine the real-time changes of parameters such as angular velocity, terminal current, power factor angle, and power angle, the node current vector can be updated based on the latest solved angular velocity, and return to step S20, so that the latest angular velocity and related output parameters can be further solved. Furthermore, in some embodiments, after step S60, step S70 can be included: updating the node current vector according to the angular velocity, and returning to step S20. In this way, data iteration can be achieved, thereby solving the dynamic changes of the angular velocity and related output parameters of the synchronous generator.

[0082] It can be understood that the present invention obtains the reactance design parameters, the original potential equation and the original flux equation of the synchronous generator, and then calculates the grid voltage and the grid current according to the reactance design parameters and the preset node current vector, and then performs dq decomposition on the grid voltage and the grid current respectively to obtain the stator voltage component data and the stator current component data, and then performs Park transformation on the original potential equation and the original flux equation of the synchronous generator to obtain the Park equation, and then determines the electromotive force of the synchronous generator according to the reactance design parameters, the stator voltage component data, the stator current component data and the Park equation, and finally determines the output parameters of the synchronous generator according to the electromotive force, and realizes the solution of the angular velocity, terminal current, power factor angle, power angle and other parameters of the synchronous generator through simulation technology, which not only has high solution efficiency, but also has the advantage of low difficulty in implementing the solution.

[0083] The present invention further provides a computer storage medium storing a computer program, which, when executed by a processor, implements the steps of the method for solving the parameters of a synchronous generator provided in an embodiment of the present invention.

[0084] like Figure 4 As shown, the present invention further provides a control device, including a processor and a memory storing a computer program, and the processor implements the steps of the synchronous generator parameter solution method provided in an embodiment of the present invention when executing the computer program.

[0085] In this specification, each embodiment is described in a progressive manner, and each embodiment focuses on the differences from other embodiments. The same or similar parts between the embodiments can be referred to each other. For the device disclosed in the embodiment, since it corresponds to the method disclosed in the embodiment, the description is relatively simple, and the relevant parts can be referred to the method part.

[0086] Professionals may further appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, computer software, or a combination of the two. In order to clearly illustrate the interchangeability of hardware and software, the composition and steps of each example have been generally described in the above description according to function. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians may use different methods to implement the described functions for each specific application, but such implementation should not be considered to be beyond the scope of the present invention.

[0087] The steps of the method or algorithm described in conjunction with the embodiments disclosed herein may be implemented directly using hardware, a software module executed by a processor, or a combination of the two. The software module may be placed in a random access memory (RAM), a memory, a read-only memory (ROM), an electrically programmable ROM, an electrically erasable programmable ROM, a register, a hard disk, a removable disk, a CD-ROM, or any other form of storage medium known in the art.

[0088] It can be understood that the above embodiments only express the preferred implementation modes of the present invention, and the descriptions thereof are relatively specific and detailed, but they cannot be understood as limiting the patent scope of the present invention. It should be pointed out that, for ordinary technicians in this field, the above technical features can be freely combined without departing from the concept of the present invention, and several deformations and improvements can be made, which all belong to the protection scope of the present invention. Therefore, all equivalent changes and modifications made to the scope of the claims of the present invention should belong to the coverage of the claims of the present invention.

Claims

1. A method for solving synchronous generator parameters for large-scale power grid node simulation calculation, characterized in that: The synchronous generator parameter solving method comprises: Obtain the reactance design parameters of the synchronous generator, the original potential equation of the synchronous generator and the original flux equation of the synchronous generator; Calculating the grid voltage and grid current according to the reactance design parameters and the preset node current vector; Performing dq decomposition on the grid voltage and the grid current respectively to obtain stator voltage component data and stator current component data; Performing Park transformation on the original potential equation of the synchronous generator and the original flux equation of the synchronous generator to obtain the Park equation; The angular velocity of the synchronous generator is determined according to the reactance design parameters, the stator voltage component data, the stator current component data and the Park equation, and the electromotive force of the synchronous generator is calculated according to the stator voltage component data to determine the output parameters of the synchronous generator according to the electromotive force.

2. The method for solving the parameters of a synchronous generator according to claim 1, characterized in that: The reactance design parameters include armature winding resistance, d-axis component of stator circuit self-inductance reactance, q-axis component of stator circuit self-inductance reactance, d-axis transient reactance, q-axis transient reactance, d-axis transient open-circuit time constant and q-axis transient open-circuit time constant; The step of calculating the grid voltage and grid current according to the reactance design parameters and the preset node current vector comprises: Determine a node admittance matrix according to the reactance design parameters, and calculate a grid voltage and a grid current according to the node current vector, the node admittance matrix and the node voltage matrix equation; The node voltage matrix equation is expressed as: YU=I, where Y represents the node admittance matrix, U represents the node voltage vector, and I represents the node injection current vector.

3. The method for solving the parameters of a synchronous generator according to claim 2, characterized in that: The step of determining the node admittance matrix according to the reactance design parameters comprises: Obtain the equivalent circuit of the synchronous generator and the power system; Based on the equivalent circuit of the synchronous generator and the power system, the reactance design parameters are equivalently transformed into the external power system network impedance, and the power system network impedance is integrated to obtain a node admittance matrix.

4. The method for solving the parameters of a synchronous generator according to claim 2, characterized in that: The Parker equations include the Parker potential equation and the Parker flux linkage equation; The Parker potential equation is expressed as: The Parker flux equation is expressed as: Among them, u d 、u q and u0 represent the d-axis component, q-axis component and 0-axis component of the stator voltage respectively, u f represents the excitation voltage, R a Represents the armature winding resistance, R f Represents the excitation resistance, R D and R Q represents the d-axis resistance and q-axis resistance of the rotor equivalent damping winding, i d 、i q and i0 represent the d-axis component, q-axis component and 0-axis component of the instantaneous stator current respectively, i f represents the excitation current, i D and i Q represents the d-axis component and q-axis component of the rotor equivalent damping winding current, p represents the differential operator, ψ d , q and ψ0 represent the d-axis component, q-axis component and 0-axis component of the stator circuit flux, respectively. f , ψ D and ψ Q represents the d-axis component and q-axis component of the rotor circuit flux, ω represents the angular velocity, x d 、x q and x0 represent the d-axis component, q-axis component and 0-axis component of the stator circuit self-inductance reactance respectively, and x f represents the self-inductance reactance of the excitation winding, x D and x Q represents the d-axis component and q-axis component of the rotor circuit self-inductance reactance, respectively, af 、x aD 、x a Q 、x fa 、x fD 、x Da 、x Df and x Qa They represent the mutual inductance reactance between different windings.

5. The method for solving the parameters of a synchronous generator according to claim 4, characterized in that: The stator voltage component data includes a d-axis component and a q-axis component of the stator voltage, and the stator current component data includes a d-axis component and a q-axis component of the stator current; The step of determining the angular velocity of the synchronous generator according to the reactance design parameter, the stator voltage component data, the stator current component data and the Parker equation comprises: The reactance design parameters, the stator voltage component data, and the stator current component data are substituted into the Parker equation to calculate the angular velocity.

6. The method for solving the parameters of a synchronous generator according to any one of claims 2 to 5, characterized in that: The output parameter of the synchronous generator includes the terminal current of the synchronous generator; The expression of the terminal current is: I t represents the terminal current, P t represents the preset end active power of the synchronous generator, Q t represents the preset end reactive power of the synchronous generator, E t represents the electromotive force of the synchronous generator.

7. The method for solving the parameters of a synchronous generator according to claim 6, characterized in that: The output parameters of the synchronous generator also include the power factor angle of the synchronous generator; The expression of the power factor angle is: φ=what ―1 [P t / (AND t THE t )]; φ represents the power factor angle.

8. The method for solving the parameters of a synchronous generator according to claim 7, characterized in that: The output parameters of the synchronous generator also include the power angle of the synchronous generator; The expression of the power angle is: δ i =tan ―1 [(x q I t cosφ―R a I t sinφ) / (E t +R a I t cosφ+x q I t sinφ)]; δ i represents the power angle, x q Represents the q-axis component of the stator circuit self-inductance, R a Represents the armature winding resistance.

9. A computer storage medium storing a computer program, characterized in that: The computer program implements the steps of the method for solving the parameters of a synchronous generator according to any one of claims 1 to 8 when executed by a processor.

10. A control device, characterized in that: It comprises a processor and a memory storing a computer program, and the processor implements the steps of the method for solving the parameters of a synchronous generator as described in any one of claims 1 to 8 when executing the computer program.

Images