CN115935879A - Modeling method and device for electromechanical transient six-order mathematical model of distributed phase modulator - Google Patents
Modeling method and device for electromechanical transient six-order mathematical model of distributed phase modulator Download PDFInfo
- Publication number
- CN115935879A CN115935879A CN202210291359.3A CN202210291359A CN115935879A CN 115935879 A CN115935879 A CN 115935879A CN 202210291359 A CN202210291359 A CN 202210291359A CN 115935879 A CN115935879 A CN 115935879A
- Authority
- CN
- China
- Prior art keywords
- axis
- equation
- winding
- transient
- obtaining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E60/00—Enabling technologies; Technologies with a potential or indirect contribution to GHG emissions mitigation
Landscapes
- Control Of Eletrric Generators (AREA)
Abstract
The invention discloses a modeling method and a device for an electromechanical transient six-order chemical model of a distributed phase modulator, which comprise the following steps: obtaining a q-axis rotor winding voltage equation, a q-axis exciting winding current and a damping winding current according to a park equation, a d-axis transient electromotive force and a sub-transient electromotive force of the double-axis exciting generator; obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and a q-axis rotor winding voltage equation; substituting the q-axis transient reactance and the sub-transient reactance into a differential equation of a q-axis rotor winding to obtain a first equation of a q-axis state variable; obtaining a q-axis damping winding voltage equation and a damping winding resistance; obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation; a second equation for the q-axis state variable is obtained based on the fundamental assumptions in the conventional generator model and the differential equation for the q-axis damping winding. The problem of the demand of the current simulation analysis on an electromechanical transient practical model is solved.
Description
Technical Field
The application relates to the technical field of mathematical modeling of motors and electric appliances, in particular to a modeling method and device of an electromechanical transient six-order mathematical model for a distributed phase modulator.
Background
The national clean energy war provides a wide market prospect for the vigorous development of new energy, however, the voltage support capability and the frequency stability of the power grid are obviously reduced because the voltage regulation and frequency modulation capability and the high and low voltage ride through capability of the new energy unit are far inferior to those of the conventional unit.
At present, extra-high voltage direct current has become a main mode of new energy power transmission. In an extra-high direct current transmission end system, with the gradual increase of new energy ratio, a conventional generator set faces a severe test in the aspect of maintaining the stable operation of the system. Firstly, the conventional generator is limited by the self-static stability limit, and the rotational inertia of the whole system is insufficient after the new energy proportion is greatly improved, so that the system is difficult to keep transient stability when serious fault impact occurs; secondly, the conventional generator is limited by the minimum exciting current, the phase advancing capability is limited, and the capability of inhibiting the transient state and the steady-state overvoltage of the system during the direct current fault and recovery period is insufficient; thirdly, the dynamic stability problem caused by insufficient system damping is difficult to solve through PSS optimization of a small number of conventional units.
Meanwhile, in an extra-high voltage direct current receiving end system, the number of local conventional units is gradually reduced due to the vigorous development of large-scale direct current feed-in and distributed energy, the dynamic reactive power reserve and the rotational inertia of the system are reduced, the transient stability and voltage stability risk of the system are aggravated, and higher requirements are provided for the voltage supporting capacity and the capacity of maintaining the transient stability of a single unit. For example, the voltage regulating capability and the short-circuit current supplying capability of the generator should be further improved to suppress transient overvoltage or transient undervoltage; the frequency modulation capability is more excellent, and the frequency modulation capability is used for compensating the defects caused by the reduction of system inertia and the reduction of frequency modulation units; the power transmission device has better stability, can absorb a large amount of fault impact power, and prevents the power from being transferred to a weak alternating current channel to cause the loss of stability of an alternating current power grid.
Therefore, a double-shaft excitation generator with ultra-strong stability has become a hot point of research in recent years. The double-shaft excitation generator is characterized in that the d shaft and the q shaft are respectively provided with an excitation winding, and the direction of excitation magnetomotive force can be changed at will by adjusting the magnitude and the direction of the excitation current of the d shaft and the q shaft, so that the power angle of the generator can be controlled within a reasonable range all the time in the transient process, and cannot be out of step due to the limitation of a static stability limit. Meanwhile, through the coordination control of the d axis and the q axis, a synthetic magnetic field rotating relative to the rotor can be obtained to realize asynchronous operation, and the active impact of a power grid is absorbed and converted into the kinetic energy of the rotor in the transient process, so that the impact resistance of the system is greatly improved; in addition, the double-shaft excitation generator can also realize active and reactive independent control, the d-shaft excitation winding can reversely flow through to provide phase advance reactive power more than twice of rated power, and meanwhile, the motor outputs required active power through the control of the q-shaft excitation winding to maintain synchronous operation. The transient dynamic stability, voltage stability and frequency stability of the system can be greatly improved. However, at present, a practical electromechanical transient model of the dual-axis excitation distributed phase modulator for simulation analysis of a large power grid is not established.
Disclosure of Invention
In order to solve the problems and research the dynamic characteristics of a double-shaft excitation distributed phase modulator in a large power grid and the influence of the double-shaft excitation distributed phase modulator on a power system, the application provides a modeling method of an electromechanical transient state sextuple order model for the distributed phase modulator, which comprises the following steps:
acquiring q-axis transient reactance, q-axis secondary dynamic reactance, d-axis transient electromotive force and d-axis secondary transient electromotive force of the double-shaft excitation generator;
obtaining a q-axis rotor winding voltage equation, a q-axis excitation winding current and a q-axis damping winding current according to a park equation, a d-axis transient electromotive force and a d-axis sub-transient electromotive force of the double-axis excitation generator;
obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and the q-axis rotor winding voltage equation;
substituting the q-axis transient state reactance and the q-axis secondary state reactance into a differential equation of the q-axis rotor winding to obtain a first equation of a q-axis state variable of the double-axis excitation generator;
acquiring a q-axis damping winding voltage equation and q-axis damping winding resistance; obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation;
and acquiring a second equation of the q-axis state variable of the double-shaft excitation generator according to basic assumptions in a traditional generator model and a differential equation of the q-axis damping winding.
Further, the q-axis transient reactance and the q-axis secondary dynamic reactance of the dual-axis excitation generator are obtained according to a park equation and a q-axis operational reactance of the dual-axis excitation generator, and specifically include:
substituting the park equation of the double-shaft excitation generator after conversion into q-axis operational reactance, and obtaining q-axis transient reactance X 'of the double-shaft excitation generator according to the initial value theorem' q And q-axis sub-dynamic reactance X ″) q 。
Further, obtaining a q-axis rotor winding voltage equation, a q-axis exciting winding current and a q-axis damping winding current according to a park equation, a d-axis transient electromotive force and a d-axis sub-transient electromotive force of the double-shaft exciting generator, comprising:
determining d-axis transient electromotive force E 'according to induced electromotive force in a traditional generator model' d And d-axis sub-transient electromotive force E ″) d ;
Obtaining a q-axis rotor winding voltage equation, a q-axis excitation winding current and a q-axis damping winding current according to the park equation;
using d-axis transient electromotive force E' d And d-axis sub-transient electromotive force E ″ d Replacing q-axis excitation winding magnetic linkage and q-axis damping winding magnetic linkage in the q-axis excitation winding current and the q-axis damping winding current to obtain a pass E' d And E ″) d The q-axis field winding current and q-axis damping winding current are shown.
Further, obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and the q-axis rotor winding voltage equation comprises:
substituting basic assumptions in a traditional generator model and q-axis excitation winding current into the q-axis rotor winding voltage equation to obtain a q-axis rotor winding differential equation.
Further, obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation, comprising:
obtaining a q-axis damping winding voltage equation according to the park equation;
according to q-axis no-load open circuit transient time constant T' q0 And q-axis no-load open-circuit sub-transient time constant T q0 Obtaining q-axis damping winding resistance;
and substituting the q-axis damping winding resistance into a q-axis damping winding voltage equation to obtain a differential equation of the q-axis damping winding.
Further, obtaining a second equation of a q-axis state variable of the dual-axis excitation generator according to basic assumptions in a traditional generator model and a differential equation of the q-axis damping winding comprises:
substituting basic assumptions in a traditional generator model into a differential equation of the q-axis damping winding to obtain a second equation for obtaining a q-axis state variable of the double-shaft excitation generator.
Further, the method also comprises the following steps:
and the first equation and the second equation are combined to generate a q-axis state variable equation of the double-shaft excitation generator.
The invention also provides a modeling device of the electromechanical transient six-order mathematical model for the distributed phase modulator, which comprises the following components:
the parameter acquisition unit is used for acquiring q-axis transient reactance, q-axis secondary dynamic reactance, d-axis transient electromotive force and d-axis secondary transient electromotive force of the double-shaft excitation generator;
the q-axis rotor winding voltage equation acquisition unit is used for acquiring a q-axis rotor winding voltage equation, q-axis excitation winding current and q-axis damping winding current according to a park equation, d-axis transient electromotive force and d-axis secondary transient electromotive force of the double-axis excitation generator;
the first differential equation obtaining unit is used for obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and the voltage equation of the q-axis rotor winding;
the first equation obtaining unit is used for substituting the q-axis transient state reactance and the q-axis secondary state reactance into a differential equation of the q-axis rotor winding to obtain a first equation of a q-axis state variable of the double-shaft excitation generator;
the second differential equation obtaining unit is used for obtaining a q-axis damping winding voltage equation and q-axis damping winding resistance; obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation;
and the second equation acquisition unit is used for acquiring a second equation of the q-axis state variable of the double-shaft excitation generator according to basic assumption in a traditional generator model and a differential equation of the q-axis damping winding.
Further, the rotor winding voltage equation obtaining unit includes:
an electromotive force determining subunit for determining d-axis transient electromotive force E 'according to induced electromotive force in the conventional generator model' d And d-axis sub-transient electromotive force E ″ d ;
The voltage equation obtaining subunit is used for obtaining a rotor winding voltage equation, q-axis excitation winding current and q-axis damping winding current according to the park equation;
a winding current obtaining subunit for using d-axis transient electromotive force E' d And d-axis sub-transient electromotive force E ″ d Replacing q-axis excitation winding magnetic linkage and q-axis damping winding magnetic linkage in the q-axis excitation winding current and the q-axis damping winding current to obtain a pass E' d And E ″) d The q-axis field winding current and q-axis damping winding current are shown.
Further, the method also comprises the following steps:
and the variable equation generating unit is used for generating a q-axis state variable equation of the double-shaft excitation generator by combining the first equation and the second equation.
The invention also provides an electronic device, which comprises a memory, a processor and a computer program stored on the memory and capable of running on the processor, wherein the processor implements the steps of any one of the above methods when executing the computer program.
The invention also provides a readable storage medium having stored thereon a computer program which, when executed by a processor, carries out the steps of the method of any one of the above.
The modeling method and the device for the electromechanical transient six-order chemical model of the distributed phase modulator provided by the invention solve the problem of the requirement of the current simulation analysis on the electromechanical transient practical model.
Drawings
Fig. 1 is a schematic flowchart of a modeling method for an electromechanical transient six-order chemical model of a distributed phase modulator according to an embodiment of the present invention;
fig. 2 is a schematic diagram of a reference direction of an electrical quantity of a double-shaft excitation generator according to an embodiment of the invention;
FIG. 3 is an equivalent circuit diagram of a q-axis open-circuit transient time constant according to an embodiment of the present invention;
FIG. 4 is an equivalent circuit diagram of a q-axis open-circuit transient time constant in accordance with an embodiment of the present invention;
fig. 5 is a schematic structural diagram of a modeling apparatus for an electromechanical transient six-order mathematical model of a distributed phase modulator according to an embodiment of the present invention.
Detailed Description
In the following description, numerous specific details are set forth in order to provide a thorough understanding of the present application. This application is capable of implementation in many different ways than those herein set forth and of similar import by those skilled in the art without departing from the spirit of this application and is therefore not limited to the specific implementations disclosed below.
Fig. 1 is a schematic flow chart of a modeling method for an electromechanical transient sextant chemical model of a distributed phase modulator according to an embodiment of the present invention, and the method according to the embodiment of the present invention is described in detail below with reference to fig. 1.
First, a brief description is given of the original park model of the generator, which involves 18 original parameters, and the acquisition of these original parameters is very difficult, so that a properly simplified practical model is generally used in the power system analysis. According to the difference of the simplified mode, 11 or 12 practical parameters are contained in the practical model. The two-shaft excitation generator differs from the conventional generator only in that one q-axis excitation winding is substituted for one damping winding, and the reference directions of the electrical quantities are as shown in fig. 2.
Wherein, f: a d-axis excitation winding;
g: a q-axis excitation winding;
d is a D-axis damping winding;
q: a q-axis damping winding;
d. q, 0: a coordinate system formed by transforming three phases of the stators a, b and c through park;
ψ f : d-axis excitation winding flux linkage;
ψ g : q-axis excitation winding flux linkage;
ψ D : d-axis damping winding flux linkage;
ψ Q : q-axis damping winding flux linkage;
ψ d : stator winding d-axis flux linkage;
ψ q : stator winding q-axis flux linkage
ψ a、b、c : stator a, b, c phase winding flux linkage;
u fd d-axis excitation winding voltage
u fq : q-axis field winding voltage
u a、b、c : stator a, b, c phase voltage;
i a、b、c : stator a, b, c phase currents;
i D : d-axis damping winding current;
i Q q-axis damping winding current;
in fig. 2, the respective electrical quantity reference directions satisfy the following regulations:
(1) The positive current of the phase modulator stator winding generates a flux linkage which is opposite to the positive direction of the flux linkage (namely, the positive current generates a negative flux linkage); the positive current of each winding of the rotor generates the magnetic flux linkage with the same positive direction as the magnetic flux linkage (the positive current generates the positive magnetic flux linkage).
(2) The stator windings follow generator convention; the rotor windings follow the motor convention.
(3) The q axis leads the d axis by 90 degrees.
(4) The positive direction of the flux linkage of the stator winding of the generator and the current flowing out of the generator (flowing from the tail end XYZ to the head end ABC) accord with the right-hand spiral rule.
Step S101, obtaining q-axis transient reactance, q-axis secondary dynamic reactance, d-axis transient electromotive force and d-axis secondary transient electromotive force of the double-shaft excitation generator.
Computing q-axis reactanceSubstituting into a park equation of the double-shaft excitation generator to obtain q-axis transient reactance X 'of the double-shaft excitation generator according to the initial value theorem' q And q-axis sub-dynamic reactance X ″) q 。
The park equation model for a two-shaft excitation generator is as follows:
wherein, X d : equivalent d inductive reactance of the stator winding;
X q : equivalent q inductive reactance of the stator winding;
X f : d-axis excitation winding inductive reactance;
X g : q-axis excitation winding inductive reactance;
X D : d-axis damping winding inductive reactance;
X Q : q-axis damping winding inductive reactance;
X af : the stator winding and the d-axis excitation winding are mutually inductively resisted;
X ag : the stator winding and the q-axis excitation winding are mutually inductively resisted;
X aD : the stator winding and the d-axis damping winding are mutually inductively resisted;
X aQ : the stator winding and the q-axis damping winding are mutually inductively resistant;
X fD : the d-axis excitation winding and the d-axis damping winding are mutually inductively coupled;
X gQ : the q-axis excitation winding and the q-axis damping winding are mutually inductively resisted;
U d : stator winding d-axis voltage;
i d : stator winding d-axis current;
r a : a stator winding resistance;
U q : a stator winding q-axis voltage;
i q : stator winding q-axis current;
r q : a stator winding q-axis resistor;
i f : d-axis field winding current;
r f : d-axis excitation winding resistance;
i g : q-axis field winding current;
r g : q-axis field winding resistance;
U D : d-axis damping winding voltage;
r D : d-axis damping winding resistance;
U Q : q-axis damping winding voltage;
r Q : q-axis damping winding resistance.
ω: the angular velocity.
In the expression, D, f and D represent D-axis windings, Q, g and Q represent Q-axis windings, the D-axis state equation of the double-shaft excitation generator is completely the same as that of the traditional generator, and the Q-axis state equation is derived.
Due to omega psi q And ω ψ d Much greater than p psi d And p psi q Therefore, in the practical generator model derivation process, the transformer electromotive force of the stator winding is ignored, i.e. p psi is considered d =pψ q =0。
The Q-axis operational reactance may be defined as follows:
from the latter two equations of equation (1)
Wherein p is a differential operator
When q-axis excitation winding voltage U fq If =0, it can be obtained according to the formula (2)
Combined (4) and (5) to obtain
According to the formula (1) to obtain
ψ q =-X q i q +X ag i g +X aQ i Q (7)
Substituting (6) and (7) into (3) to obtain
Obtained by the theorem of initial values
When the Q winding is not considered, there is r Q If = ∞ is substituted into formula (8) to obtain
Obtained by the theorem of initial values
And S102, obtaining a rotor winding voltage equation, a q-axis exciting winding current and a q-axis damping winding current according to the park equation, the d-axis transient electromotive force and the d-axis transient electromotive force.
D-axis transient electromotive force E 'is determined according to induced electromotive force in a traditional generator model' d And d-axis sub-transient electromotive force E ″ d (ii) a Obtaining a rotor winding voltage equation, q-axis excitation winding current and q-axis damping winding current according to the park equation; use ofd-axis transient electromotive force E' d And d-axis sub-transient electromotive force E ″ d Replacing q-axis excitation winding magnetic linkage and q-axis damping winding magnetic linkage in the q-axis excitation winding current and the q-axis damping winding current to obtain a pass E' d And E ″) d The q-axis field winding current and q-axis damping winding current are shown.
To convert the rotor variables to the stator side for analysis and metrology, the following practical variables were introduced.
T′ q0 And T ″) q0 The definition and equivalent circuit of (2) are shown in fig. 3 and 4.
In FIG. 3, X g1 Leakage reactance of q-axis excitation winding, r g For q-axis excitation winding resistance
In FIG. 4, X Q1 Is leakage reactance of q-axis damping winding, r Q Is a q-axis damping winding resistance.
Wherein:
T′ q0 : q-axis no-load open circuit transient time constant;
T″ q0 : q-axis no-load open-circuit sub-transient time constant.
With reference to the definition of induced electromotive force in a conventional generator model in PSD-BPA, the following practical variables are defined:
wherein:
E fq : q-axis field winding electromotive force;
E d : d-axis induced electromotive force;
E′ d : d-axis transient induced electromotive force;
E″ d : d-axis sub-transient induced electromotive force.
The equation for the rotor winding voltage obtained from equation (1) is as follows
pψ g =u fq -r g i (14)
T′ q0 pE′ d =E fq +X ag i g (15)
From the last two equations of equation (1) can be derived
Namely, it is
To eliminate the flux linkage variable in the equation, the flux linkage ψ may be set g And psi Q From E' d And E ″) d And (4) showing. According to formula (13) to obtain
Substituting (18) into (17) to obtain
And step S103, obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and the q-axis rotor winding voltage equation.
Substituting basic assumptions in a traditional generator model into the q-axis rotor winding voltage equation to obtain a q-axis rotor winding differential equation.
Reference to the basic assumptions in the conventional Generator model in PSD-BPA
Substituting the primary formulae in (19) and (20) into (15) can obtain E' d I.e. the differential equation of the q-axis rotor winding.
And step S104, substituting the q-axis transient state reactance and the q-axis secondary state reactance into a differential equation of the q-axis rotor winding to obtain a first equation of a q-axis state variable of the double-shaft excitation generator.
Substituting (9) and (11) into (21) to obtain
The above equation is the first equation of the q-axis state variable of the two-axis excitation generator.
Step S105, obtaining a q-axis damping winding voltage equation and q-axis damping winding resistance; and obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation.
Obtaining a q-axis damping winding voltage equation according to the park equation; according to q-axis no-load open-circuit transient time constant T' q0 And q-axis no-load open-circuit sub-transient time constant T q0 Obtaining q-axis damping winding resistance; and substituting the q-axis damping winding resistance into a q-axis damping winding voltage equation to obtain a differential equation of the q-axis damping winding. The q-axis damping winding voltage equation can also be obtained according to the equation (1) as follows
pψ Q =-r Q i Q (23)
Is obtainable from the formula (12)
Substituting equation (24) into equation (23) and multiplying by equation on both sidesCan obtain the product
Substitution of the second formula in (20) into (25) yields the result for E ″ " d A differential equation of (2).
Substituting the second formula in (19) into (26) to obtain
And step S106, acquiring a second equation of the q-axis state variable of the double-axis excitation generator according to basic assumptions in a traditional generator model and a differential equation of the q-axis damping winding.
The above equation is a second equation of the q-axis state variable of the two-axis excitation generator.
The first equation and the second equation are combined to generate a q-axis state variable equation of the double-shaft excitation generator.
Equations (22) and (28) are the q-axis state variable equation of the double-shaft excitation generator; the d-axis state equation and the rotor motion equation of the double-shaft excitation generator are the same as those of the traditional generator and are not deduced.
The specific application examples are as follows:
in order to make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions in the embodiments of the present invention will be clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are some, but not all, embodiments of the present invention. Thus, the following detailed description of the embodiments of the present invention, presented in the figures, is not intended to limit the scope of the invention, as claimed, but is merely representative of selected embodiments of the invention. All other embodiments, which can be obtained by a person skilled in the art without inventive step based on the embodiments of the present invention, are within the scope of protection of the present invention.
For simplicity of calculation, simplified assumptions were made in PSD-BPA and PSASP for the mutual inductance between the d and q axis windings, respectively.
The basic assumption A: PSD-BPA hypothesis
Meanwhile, each electrical quantity of PSD-BPA is defined as follows:
basic assumption B: PSASP hypothesis
Meanwhile, the PSASP electric quantities are defined as follows:
example 1:
a PSD-BPA double-shaft excitation phase modulator/generator six-order practical model based on hypothesis A is as follows:
example 2:
the six-order utility model of the PSASP two-axis excitation phase modulator/generator based on hypothesis B is as follows:
in the above formula X 1 Is the stator leakage reactance.
Based on the same inventive concept, the present invention also provides a modeling apparatus 500 for an electromechanical transient six-order mathematical model of a distributed phase modulator, as shown in fig. 5, including:
a parameter obtaining unit 510, configured to obtain a q-axis transient reactance, a q-axis secondary dynamic reactance, a d-axis transient electromotive force, and a d-axis secondary transient electromotive force of the dual-axis excitation generator;
a q-axis rotor winding voltage equation obtaining unit 520, configured to obtain a q-axis rotor winding voltage equation, a q-axis excitation winding current, and a q-axis damping winding current according to a park equation, a d-axis transient electromotive force, and a d-axis sub-transient electromotive force of the dual-axis excitation generator;
a first differential equation obtaining unit 530, configured to obtain a differential equation of the q-axis rotor winding according to a traditional generator model and the q-axis rotor winding voltage equation;
a first equation obtaining unit 540, configured to substitute the q-axis transient reactance and the q-axis sub-dynamic reactance into a differential equation of the q-axis rotor winding to obtain a first equation of a q-axis state variable of the dual-axis excitation generator;
a second differential equation obtaining unit 550, configured to obtain a q-axis damping winding voltage equation and a q-axis damping winding resistance; obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation;
the second equation obtaining unit 560 is configured to obtain a second equation of the q-axis state variable of the dual-axis excitation generator according to a basic assumption in a conventional generator model and a differential equation of the q-axis damping winding.
Further, the rotor winding voltage equation obtaining unit includes:
an electromotive force determining subunit, configured to determine a d-axis transient electromotive force E 'according to the induced electromotive force in the conventional generator model' d And d-axis sub-transient electromotive force E ″) d ;
A voltage equation obtaining subunit, configured to obtain a q-axis rotor winding voltage equation, a q-axis field winding current, and a q-axis damping winding current according to the park equation
A winding current obtaining subunit for obtaining the d-axis transient electromotive force E' d And d-axis sub-transient electromotive force E ″ d Replacing q-axis excitation winding magnetic linkage and q-axis damping winding magnetic linkage in the q-axis excitation winding current and the q-axis damping winding current to obtain a pass E' d And E ″) d The q-axis field winding current and q-axis damping winding current are shown.
Further, the method also comprises the following steps:
and the variable equation generating unit is used for generating a q-axis state variable equation of the double-shaft excitation generator by combining the first equation and the second equation.
Compared with the prior art, the modeling method of the six-order electromechanical transient practical model of the generator/phase modulator/motor based on double-shaft excitation considers the influence of q-axis excitation voltage on the model, can be directly adopted by large power grid simulation software such as PSD-BPA (phase-sensitive Power-Process-array) and PSASP (Power System analysis software analysis) and the like, and solves the problem of the requirement of the current simulation analysis on the electromechanical transient practical model of the double-shaft excitation distributed phase modulator.
As will be appreciated by one skilled in the art, embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flow diagrams and/or block diagrams, and combinations of flows and/or blocks in the flow diagrams and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer implemented process such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Finally, it should be noted that: although the present invention has been described in detail with reference to the foregoing embodiments, it should be understood by those skilled in the art that various changes may be made and equivalents may be substituted for elements thereof without departing from the spirit and scope of the invention.
Claims (12)
1. A method of modeling an electromechanical transient sextant mathematical model for a distributed phase modulator, comprising:
acquiring q-axis transient reactance, q-axis secondary dynamic reactance, d-axis transient electromotive force and d-axis secondary transient electromotive force of the double-shaft excitation generator;
obtaining a q-axis rotor winding voltage equation, a q-axis excitation winding current and a q-axis damping winding current according to a park equation, a d-axis transient electromotive force and a d-axis sub-transient electromotive force of the double-axis excitation generator;
obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and the q-axis rotor winding voltage equation;
substituting the q-axis transient state reactance and the q-axis secondary state reactance into a differential equation of the q-axis rotor winding to obtain a first equation of a q-axis state variable of the double-shaft excitation generator;
acquiring a q-axis damping winding voltage equation and q-axis damping winding resistance; obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation;
and acquiring a second equation of the q-axis state variable of the double-shaft excitation generator according to basic assumptions in a traditional generator model and a differential equation of the q-axis damping winding.
2. The method according to claim 1, wherein the q-axis transient reactance and the q-axis sub-dynamic reactance of the dual-axis excitation generator are obtained according to a park equation and a q-axis operational reactance of the dual-axis excitation generator, and specifically comprises:
substituting the park equation of the double-shaft excitation generator after conversion into q-axis operational reactance, and obtaining q-axis transient reactance X 'of the double-shaft excitation generator according to the initial value theorem' q And q-axis sub-dynamic reactance X ″) q 。
3. The method of claim 1, wherein obtaining the q-axis rotor winding voltage equation and the q-axis field winding current and the q-axis damping winding current from a park equation, a d-axis transient electromotive force and a d-axis sub-transient electromotive force of the dual-axis field generator comprises:
d-axis transient electromotive force E 'is determined according to induced electromotive force in a traditional generator model' d And d-axis sub-transient electromotive force E ″) d ;
Obtaining a q-axis rotor winding voltage equation, a q-axis excitation winding current and a q-axis damping winding current according to the park equation;
using d-axis transient electromotive force E' d And d-axis sub-transient electromotive force E ″ d Replacing q-axis excitation winding magnetic linkage and q-axis damping winding magnetic linkage in the q-axis excitation winding current and the q-axis damping winding current to obtain a pass E' d And E ″) d The q-axis field winding current and q-axis damping winding current are shown.
4. The method of claim 1, wherein obtaining a differential equation for the q-axis rotor winding from a conventional generator model and the q-axis rotor winding voltage equation comprises:
substituting basic assumptions in a traditional generator model and q-axis excitation winding current into the q-axis rotor winding voltage equation to obtain a q-axis rotor winding differential equation.
5. The method of claim 1, wherein obtaining a differential equation for the q-axis damping winding based on the q-axis damping winding resistance and the q-axis damping winding voltage equation comprises:
obtaining a q-axis damping winding voltage equation according to the park equation;
according to q-axis no-load open-circuit transient time constant T' q0 And q-axis no-load open-circuit sub-transient time constant T q0 Obtaining q-axis damping winding resistance;
and substituting the q-axis damping winding resistance into a q-axis damping winding voltage equation to obtain a differential equation of the q-axis damping winding.
6. The method according to claim 1, wherein obtaining a second equation of a q-axis state variable of a two-axis excited generator from a fundamental assumption in a conventional generator model and a differential equation of the q-axis damping winding comprises:
substituting basic assumptions in a traditional generator model into a differential equation of the q-axis damping winding to obtain a second equation for obtaining a q-axis state variable of the double-shaft excitation generator.
7. The method of claim 1, further comprising:
the first equation and the second equation are combined to generate a q-axis state variable equation of the double-shaft excitation generator.
8. A modeling apparatus for an electromechanical transient sextant mathematical model of a distributed phase modulator, comprising:
the parameter acquisition unit is used for acquiring q-axis transient reactance, q-axis secondary dynamic reactance, d-axis transient electromotive force and d-axis secondary transient electromotive force of the double-shaft excitation generator;
the q-axis rotor winding voltage equation acquisition unit is used for acquiring a q-axis rotor winding voltage equation, q-axis excitation winding current and q-axis damping winding current according to a park equation, d-axis transient electromotive force and d-axis sub-transient electromotive force of the double-axis excitation generator;
the first differential equation obtaining unit is used for obtaining a differential equation of the q-axis rotor winding according to a traditional generator model and the voltage equation of the q-axis rotor winding;
the first equation obtaining unit is used for substituting the q-axis transient state reactance and the q-axis secondary state reactance into a differential equation of the q-axis rotor winding to obtain a first equation of a q-axis state variable of the double-shaft excitation generator;
the second differential equation obtaining unit is used for obtaining a q-axis damping winding voltage equation and a q-axis damping winding resistance; obtaining a differential equation of the q-axis damping winding according to the q-axis damping winding resistance and the q-axis damping winding voltage equation;
and the second equation acquisition unit is used for acquiring a second equation of the q-axis state variable of the double-shaft excitation generator according to basic assumption in a traditional generator model and a differential equation of the q-axis damping winding.
9. The apparatus of claim 8, wherein the rotor winding voltage equation obtaining unit comprises:
an electromotive force determining subunit, configured to determine a d-axis transient electromotive force E 'according to the induced electromotive force in the conventional generator model' d And d-axis sub-transient electromotive force E ″ d ;
The voltage equation obtaining subunit is used for obtaining a rotor winding voltage equation, q-axis excitation winding current and q-axis damping winding current according to the park equation;
a winding current obtaining subunit for using d-axis transient electromotive force E' d And d-axis sub-transient electromotive force E ″ d Replacing q-axis excitation winding magnetic linkage and q-axis damping winding magnetic linkage in the q-axis excitation winding current and the q-axis damping winding current to obtain a pass E' d And E ″) d The q-axis field winding current and q-axis damping winding current are shown.
10. The apparatus of claim 8, further comprising:
and the variable equation generating unit is used for generating a q-axis state variable equation of the double-shaft excitation generator by combining the first equation and the second equation.
11. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, characterized in that the processor implements the steps of the method according to any one of claims 1 to 7 when executing the computer program.
12. A readable storage medium, on which a computer program is stored, which, when being executed by a processor, carries out the steps of the method according to any one of claims 1 to 7.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210291359.3A CN115935879A (en) | 2022-03-23 | 2022-03-23 | Modeling method and device for electromechanical transient six-order mathematical model of distributed phase modulator |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202210291359.3A CN115935879A (en) | 2022-03-23 | 2022-03-23 | Modeling method and device for electromechanical transient six-order mathematical model of distributed phase modulator |
Publications (1)
Publication Number | Publication Date |
---|---|
CN115935879A true CN115935879A (en) | 2023-04-07 |
Family
ID=86654745
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202210291359.3A Pending CN115935879A (en) | 2022-03-23 | 2022-03-23 | Modeling method and device for electromechanical transient six-order mathematical model of distributed phase modulator |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN115935879A (en) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117349576A (en) * | 2023-12-05 | 2024-01-05 | 华能澜沧江水电股份有限公司 | Method and system for calculating overvoltage of generator power-on braking rotor |
-
2022
- 2022-03-23 CN CN202210291359.3A patent/CN115935879A/en active Pending
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN117349576A (en) * | 2023-12-05 | 2024-01-05 | 华能澜沧江水电股份有限公司 | Method and system for calculating overvoltage of generator power-on braking rotor |
CN117349576B (en) * | 2023-12-05 | 2024-03-01 | 华能澜沧江水电股份有限公司 | Method and system for calculating overvoltage of generator power-on braking rotor |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Mo et al. | Evaluation of virtual synchronous machines with dynamic or quasi-stationary machine models | |
CN105281350B (en) | Micro-capacitance sensor frequency control method and system | |
Prasad et al. | Mathematical modeling and position-sensorless algorithm for stator-side field-oriented control of rotor-tied DFIG in rotor flux reference frame | |
Long et al. | Asymmetrical low-voltage ride through of brushless doubly fed induction generators for the wind power generation | |
CN104836503B (en) | Double-fed aerogenerator stator short-circuit current calculation method under unbalanced fault | |
CN107947186B (en) | Starting control method for accessing high-voltage direct-current transmission inversion side to synchronous phase modulator | |
CN104579080A (en) | Torque pulsation inhibition method for permanent magnet synchronous motor | |
CN106998074B (en) | Control method and system for modular multi-level flexible direct current converter station | |
CN107124126A (en) | The no phase-locked loop current control method and device of a kind of double fed induction generators | |
Maharjan et al. | Real-time simulation for active and reactive power control of doubly fed induction generator | |
CN108573094A (en) | The method for building up and system of the VBR electromagnetic transient simulation models of synchronous generator | |
CN115935879A (en) | Modeling method and device for electromechanical transient six-order mathematical model of distributed phase modulator | |
Liu et al. | Recent advances of control technologies for brushless doubly-fed generators | |
Banerjee et al. | Control architecture for a switched doubly fed machine propulsion drive | |
CN102739147A (en) | Back electromotive force harmonic compensation control method of permanent magnetic synchronous motor | |
Pimple et al. | New direct torque control of DFIG under balanced and unbalanced grid voltage | |
Zhang et al. | Stator-current-based MRAS observer for the sensorless control of the brushless doubly-fed induction machine | |
Chen et al. | Modeling and current control of a double salient permanent magnet generator (DSPMG) | |
Tu et al. | Short-circuit sequence network model of DFIG with different slips | |
dos Santos Serra et al. | Implementation of a virtual synchronous machine to improve the dynamic response of inverters | |
Prasad et al. | Position-sensorless algorithm for stator-flux-oriented control of single VSC-based DFIG | |
Yang et al. | Direct power control of brushless doubly-fed reluctance generator based on sliding mode control | |
Ofualagba | Analysis of the dynamic characteristics of an isolated self-excited induction generator driven by a wind-turbine | |
Onah et al. | Dynamic Modelling and Simulation of Salient Pole Synchronous Motor Using Embedded Matlab | |
Hasanova et al. | STARTING MODE OF SYNCHRONOUS MACHINES WITH MASSIVE ROTORS |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination |