CN117039996A - Transient stability analysis and adjustment method and device for full-power wind turbine generator - Google Patents

Transient stability analysis and adjustment method and device for full-power wind turbine generator Download PDF

Info

Publication number
CN117039996A
CN117039996A CN202310987654.7A CN202310987654A CN117039996A CN 117039996 A CN117039996 A CN 117039996A CN 202310987654 A CN202310987654 A CN 202310987654A CN 117039996 A CN117039996 A CN 117039996A
Authority
CN
China
Prior art keywords
transient stability
full
wind turbine
power
pll
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
Application number
CN202310987654.7A
Other languages
Chinese (zh)
Inventor
孙大卫
刘辉
宋嘉炜
宋鹏
白恺
吴林林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Corp of China SGCC
State Grid Jibei Electric Power Co Ltd
Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Jibei Electric Power Co Ltd
Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by State Grid Corp of China SGCC, State Grid Jibei Electric Power Co Ltd, Electric Power Research Institute of State Grid Jibei Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202310987654.7A priority Critical patent/CN117039996A/en
Publication of CN117039996A publication Critical patent/CN117039996A/en
Pending legal-status Critical Current

Links

Classifications

    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The application provides a transient stability analysis and adjustment method and device for a full-power wind turbine, wherein the method comprises the following steps: acquiring a transient stability analysis model under a limit condition based on an electromagnetic power angle theory of the full-power wind turbine generator; determining electromagnetic power angle characteristic curves before and after the occurrence of the fault working condition; calculating limit cutting time when the acceleration area and the deceleration area in the electromagnetic power angle characteristic curve are equal; evaluating the transient stability margin value of the full-power wind turbine generator according to the limit cutting time; and in response to the transient stability margin value not meeting the requirement, increasing the transient stability margin value of the full-power wind turbine generator by increasing the limit cutting time. The method obtains an expression of transient stability criterion-fault limit excision time of the full-power wind turbine generator under the voltage deep drop working condition. And then according to the variables involved in the expression, the objective of improving the transient stability margin value and further improving the transient stability of the system is achieved by adjusting the physical quantity corresponding to the variables.

Description

Transient stability analysis and adjustment method and device for full-power wind turbine generator
Technical Field
The application relates to the technical field of wind power generation stability, in particular to a transient stability analysis and adjustment method and device for a full-power wind turbine generator.
Background
At present, the transient stability analysis basic methods of the traditional power system are divided into three main types: numerical solution, equal area law and direct method, respectively. It should be noted that in the use of the above-mentioned various analysis methods, the power system is assumed to have different degrees, and when the transient stability mechanism of the wind power system is analyzed by using the traditional power system method, the assumption under different methods needs to be considered for the wind power grid-connected system, so that no related study exists in the aspect of transient stability analysis and adjustment of the full-power wind turbine, and therefore, a scheme capable of implementing analysis and adjustment of the transient stability of the full-power wind turbine is needed.
Disclosure of Invention
In view of the above, the present application provides a method and an apparatus for transient stability analysis adjustment of a full power wind turbine generator to solve at least one of the above-mentioned problems.
In order to achieve the above purpose, the present application adopts the following scheme:
according to a first aspect of the present application, an embodiment of the present application provides a transient stability analysis adjustment method for a full-power wind turbine, where the method includes: acquiring a transient stability analysis model under a limit condition based on an electromagnetic power angle theory of the full-power wind turbine generator; determining electromagnetic power angle characteristic curves before and after the occurrence of the fault working condition; calculating an acceleration area and a deceleration area based on the electromagnetic power angle characteristic curve; calculating a limit cutting angle and a corresponding limit cutting time when the acceleration area and the deceleration area are equal; evaluating the transient stability margin value of the full-power wind turbine according to the limit cutting time; and in response to the transient stability margin value not meeting the requirement, increasing the transient stability margin value of the full-power wind turbine by increasing the limit cutting time.
Preferably, the transient stability analysis model in the method of the embodiment of the present application is:
in the above formula: t (T) pll Is a mechanical-like inertia that is analogous to the mechanical inertia in the synchronous generator rotor equation of motion; t (T) m Is a mechanical-like torque that is analogous to the mechanical torque in the synchronous generator rotor equation of motion; t (T) e Is an electromagnetic-like torque that is analogous to the electromagnetic torque in the synchronous generator rotor equation of motion; d (D) pll Is an electrical-like damping that is analogous to the amount of electrical damping in the synchronous generator rotor motion equation; omega pll The phase-locked angular velocity is output by the phase-locked loop; omega g The angular speed is the corresponding angular speed of the power grid at the power frequency of 50 Hz; omega pll -ω g ＝Δω pll ，δ pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, in the above method according to the embodiment of the present application, calculating the limit cut angle when the acceleration area and the deceleration area are equal includes: and carrying out iterative solution by using the following formula to obtain an expression of the limit cutting angle:
in the above, delta cr For extreme cut angle, delta A Is the electromagnetic power angle delta at the point A J Is the electromagnetic power angle at the J point, T m In the form of a mechanical-like torque,for the pre-fault class electromagnetic torque +.>As the electromagnetic torque after failure, delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, in the above method according to the embodiment of the present application, the limit cutting time corresponding to the limit cutting angle is calculated by the following formula:
in the above, t cr For extreme excision time, delta cr For extreme cut angle, delta A Is the electromagnetic power angle at the point A, T pll Is of similar mechanical inertia, T m In the form of a mechanical-like torque,to be the electromagnetic torque of the pre-fault class omega g For the corresponding angular speed of the power frequency of 50Hz of the power grid, S I For accelerating area, delta pll ＝θ pll -Δ pcc ，θ pll For phase-locked loop controlled output phase, delta pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, in the method according to the embodiment of the present application, estimating the transient stability margin value of the full-power wind turbine according to the limit excision time includes: the larger the limit cutting time of the fault is, the larger the transient stability margin value of the corresponding full-power wind turbine generator is estimated.
Preferably, in the method according to the embodiment of the present application, in response to the transient stability margin value not meeting a requirement, the increasing the transient stability margin value of the full-power wind turbine by increasing the limit cutting time includes: by reducing the phase-locked loop bandwidth within a set range; and/or reducing the voltage sag depth; and/or to increase grid strength conditions; and/or the adjusting system is set to enable the full-power wind turbine generator to have a phase angle abrupt change value with a set negative value under the fault working condition; and/or reducing the electrical distance from the fan port to the fault point in a set range to improve the transient stability margin value of the full-power wind turbine generator.
According to a second aspect of the present application, an embodiment of the present application further provides a transient stability analysis adjustment device of a full-power wind turbine, where the device includes: the model acquisition unit is used for acquiring a transient stability analysis model under a limit condition based on the electromagnetic power angle theory of the full-power wind turbine generator; the fault electromagnetic power angle characteristic curve determining unit is used for determining electromagnetic power angle characteristic curves before and after the fault working condition occurs; an acceleration/deceleration area calculation unit for calculating an acceleration area and a deceleration area amount based on the electromagnetic work angle characteristic curve; a limit cut-off time calculation unit configured to calculate a limit cut-off angle and a corresponding limit cut-off time when the acceleration area and the deceleration area are equal in amount; the margin value evaluation unit is used for evaluating the transient stability margin value of the full-power wind turbine generator according to the limit cutting time; and the margin value adjusting unit is used for increasing the transient stability margin value of the full-power wind turbine generator by increasing the limit cutting time in response to the transient stability margin value not meeting the requirement.
Preferably, the transient stability analysis model in the device of the embodiment of the present application is:
in the above formula: t (T) pll Is a mechanical-like inertia that is analogous to the mechanical inertia in the synchronous generator rotor equation of motion; t (T) m Is a similar mechanical torque, which is analogous to the rotor movement of synchronous generatorMechanical torque in the journey; t (T) e Is an electromagnetic-like torque that is analogous to the electromagnetic torque in the synchronous generator rotor equation of motion; d (D) pll Is an electrical-like damping that is analogous to the amount of electrical damping in the synchronous generator rotor motion equation; omega pll The phase-locked angular velocity is output by the phase-locked loop; omega g The angular speed is the corresponding angular speed of the power grid at the power frequency of 50 Hz; omega pll -ω g ＝Δω pll ，δ pll ＝Δ pll -Δ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, the limit cut-off time calculation unit in the above-described apparatus of the embodiment of the present application, when the acceleration area and the deceleration area are equal in amount, includes:
and carrying out iterative solution by using the following formula to obtain an expression of the limit cutting angle:
in the above, delta cr For extreme cut angle, delta A Is the electromagnetic power angle delta at the point A J Is the electromagnetic power angle at the J point, T m In the form of a mechanical-like torque,for the pre-fault class electromagnetic torque +.>As the electromagnetic torque after failure, delta pll ＝Δ pll -Δ pcc ，δ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, in the device according to the embodiment of the present application, the limit cutting time calculating unit calculates the limit cutting time corresponding to the limit cutting angle by the following formula:
in the above, t cr For extreme excision time, delta cr For extreme cut angle, delta A Is the electromagnetic power angle at the point A, T pll Is of similar mechanical inertia, T m In the form of a mechanical-like torque,to be the electromagnetic torque of the pre-fault class omega g For the corresponding angular speed of the power frequency of 50Hz of the power grid, S I For accelerating area, delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, the evaluation unit for the margin value in the device according to the embodiment of the present application evaluates the transient stability margin value of the full-power wind turbine generator according to the limit cutting time, including: the larger the limit cutting time of the fault is, the larger the transient stability margin value of the corresponding full-power wind turbine generator is estimated.
Preferably, in the foregoing apparatus according to the embodiment of the present application, in response to the transient stability margin value not meeting a requirement, the margin value adjustment unit is specifically configured to: by reducing the phase-locked loop bandwidth within a set range; and/or reducing the voltage sag depth; and/or to increase grid strength conditions; and/or the adjusting system is set to enable the full-power wind turbine generator to have a phase angle abrupt change value with a set negative value under the fault working condition; and/or reducing the electrical distance from the fan port to the fault point in a set range to improve the transient stability margin value of the full-power wind turbine generator.
According to a third aspect of the present application, there is also provided an electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, the processor implementing the steps of the above method when executing the computer program.
According to a fourth aspect of the present application, embodiments of the present application also provide a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the above-described method.
According to a fifth aspect of the application, embodiments of the application also provide a computer program product comprising computer programs/instructions which, when executed by a processor, implement the steps of the above method.
According to the transient stability analysis and adjustment method and device for the full-power wind turbine, the electromagnetic power angle theory is utilized to provide a criterion for transient stability of the full-power wind turbine, namely, an expression of transient stability criterion-fault limit removal time of the full-power wind turbine under the voltage deep drop working condition is obtained by carrying out iterative solution and integral solution on a transient stability analysis model formula under the limit condition based on the electromagnetic power angle theory. And then according to the variables involved in the expression, the objective of improving the transient stability margin value and further improving the transient stability of the system is achieved by adjusting the physical quantity corresponding to the variables.
Drawings
In order to more clearly illustrate the embodiments of the application or the technical solutions in the prior art, the drawings that are required in the embodiments or the description of the prior art will be briefly described, it being obvious that the drawings in the following description are only some embodiments of the application, and that other drawings may be obtained according to these drawings without inventive effort for a person skilled in the art. In the drawings:
FIG. 1 is a schematic flow chart of a transient stability analysis and adjustment method for a full-power wind turbine generator provided by an embodiment of the application;
FIG. 2 is a schematic diagram of an electromagnetic angle characteristic curve after failure according to an embodiment of the present application;
FIG. 3 is a graph of phase-locked loop bandwidth versus limit cut time provided by an embodiment of the present application;
FIG. 4 is a graph of fault condition variables versus limit cut-off time provided by an embodiment of the present application;
FIG. 5 is a graph of reactance versus limit cut-off time for a fan port to a point of failure provided by an embodiment of the present application;
fig. 6 is a schematic structural diagram of a transient stability analysis and adjustment device of a full-power wind turbine generator provided by an embodiment of the present application;
fig. 7 is a schematic structural diagram of an electronic device according to an embodiment of the present application.
Detailed Description
For the purpose of making the objects, technical solutions and advantages of the embodiments of the present application more apparent, the embodiments of the present application will be described in further detail with reference to the accompanying drawings. The exemplary embodiments of the present application and their descriptions herein are for the purpose of explaining the present application, but are not to be construed as limiting the application.
Fig. 1 is a schematic flow chart of a transient stability analysis and adjustment method for a full-power wind turbine, which is provided by the embodiment of the application, and includes the following steps:
step S101: and obtaining a transient stability analysis model under the limit condition based on the electromagnetic power angle theory of the full-power wind turbine generator. In this embodiment, the limit condition indicated in this step refers to a scenario that no phase angle mutation exists under the working condition that no balance point exists after the full-power wind turbine generator fails.
Preferably, the expression of the transient stability analysis model under the limit condition is obtained based on the electromagnetic power angle theory of the full-power wind turbine generator set as follows (1):
wherein: t (T) pll Is a mechanical-like inertia that is analogous to the mechanical inertia in the synchronous generator rotor equation of motion; t (T) m Is a mechanical-like torque that is analogous to the mechanical torque in the synchronous generator rotor equation of motion; t (T) e Is an electromagnetic-like torque that is analogous to the electromagnetic torque in the synchronous generator rotor equation of motion; d (D) pll Is an electrical-like damping that is analogous to the amount of electrical damping in the synchronous generator rotor motion equation; omega pll The phase-locked angular velocity is output by the phase-locked loop; omega g Is electric powerThe corresponding angular speed of the network power frequency is 50 Hz; omega pll -ω g ＝Δω pll ，δ pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Further, T pll 、T m 、T e And D pll Can be represented by the following formulas (2) to (5), respectively:
T m ＝I d X+I q R (3)
T e ＝|U pcc |sinδ pll (4)
wherein the method comprises the steps ofPhase-locked loop PI parameter I of full-power wind turbine generator d 、I q Current vector I respectively output by wind turbine generator wt The d and q axis corresponding components after dq decomposition are X, R, i U, respectively, the sum of reactance values and resistance values of the box transformer, the main transformer and the corresponding transmission lines between the fan port and the grid connection point pcc And I is the modulus of the voltage vector of the grid-connected point.
The transient stability analysis model is developed from the electromagnetic power angle concept of the full-power converter, and reveals that the transient stability mechanism of the full-power wind turbine generator can be compared with the theory used by the traditional power system for analyzing the transient stability mechanism of the synchronous generator; after the model is built, the model change before and after the fault can be analyzed by using the equal area rule in the subsequent step, so that the transient stability mechanism of the full-power wind turbine generator can be revealed, and corresponding adjustment is provided to improve the transient stability.
Step S102: and determining electromagnetic power angle characteristic curves before and after the fault working condition occurs.
The transient stability analysis and adjustment method provided by the embodiment of the application is applicable to systems with different voltage drop degrees. Taking the case of deeper voltage drop, the electromagnetic power angle characteristic curve after the fault can be shown in fig. 2 under the working condition that no balance point exists after the fault and no phase angle mutation exists.
As can be seen from FIG. 2, the full-power wind turbine generator is firstly operated at the balance point A, at this time, the sudden fault causes the voltage depth to drop, the operation point is suddenly changed to the point B, and when the fault working condition has no factor inducing the phase angle suddenly changed, i.e. no phase angle suddenly changed, the electromagnetic torque-like characteristic curve after the fault is as shown in the graphThe curve shows that during a fault the wind power system will be along +.>Curve operation, namely that during fault period, the similar mechanical torque is always larger than the similar electromagnetic torque, so that the rotating speed is continuously increased, and the electromagnetic power angle delta is increased pll The wind power system cannot restore the balance of the similar mechanical torque and the similar electromagnetic torque by itself under the working condition, and faults need to be quickly removed to restore the balance of the similar mechanical torque and the similar electromagnetic torque. Under the working condition, if the fault is cut when the wind power system operates to the point C, at the moment, the operating point of the wind power system is suddenly changed to the point D, and the instantaneous mechanical torque is smaller than the electromagnetic torque, but the rotating speed is still larger than the rotating speed corresponding to the point A of the balance point at the moment, so that the wind power system still can be along the T before the fault e - The characteristic is slowly accelerated for a period of time, when S is present I ＝S II (wherein S I To accelerate the area S II For deceleration area), after the wind power system operates to the point E, the electromagnetic-like torque at the point E is larger than the mechanical-like torque, and the damping inertia of the system exists, so that the wind power system is finally restored to the point A of the balance point, and the system is stable at the moment, wherein the coordinate of the point A is A= (delta omega) pll ，δ pll )＝(0，δ A )。
Step S103: calculating the acceleration area and the deceleration area based on the electromagnetic angle characteristic curve, namely calculating S in the electromagnetic angle characteristic curve of the figure 2 I And S is II Is a part of the surface area of the substrate.
Step S104: and calculating a limit cutting angle and a corresponding limit cutting time when the acceleration area and the deceleration area are equal. The limit cut-off time is a very important parameter in transient stability analysis and control of the power system, and represents the longest fault duration that the power system can bear without losing stability after suffering a certain fault, and as can be seen from the description, when the acceleration area and the deceleration area are equal, the wind power system will be stable, so that the limit cut-off time is obtained, only the limit cut-off angle when the acceleration area and the deceleration area are equal needs to be calculated, and the limit cut-off time is obtained according to the limit cut-off angle.
In this embodiment, the transient stability analysis model, that is, the formula (1), may be subjected to integral transformation, and then the following is obtained:
due to self delta A →δ E In the process, T e The curve before and after the fault is divided into two curves, so the formula (6) is split, and the damping quantity D is given pll When=0, it is possible to obtain:
due to the left side of the equationThe above formula (7) can be simplified as:
obviously when delta E ＝δ J At delta, when A To delta J Between which there is a limit cutting angle delta cr The following is established:
in the above formulas (6) - (9), δ cr Is the limit cutting angle; delta A 、δ C 、δ E 、δ J Electromagnetic power angles at A, C, E, J points respectively; t (T) pll Is a similar mechanical inertia; d (D) pll Is an electrical-like damping; t (T) m Is a mechanical-like torque; t (T) e Is electromagnetic torque;is the electromagnetic torque before failure; />The electromagnetic torque is the post-fault type electromagnetic torque; Δω pll ＝ω pll -ω g ，ω pll The phase-locked angular velocity is output by the phase-locked loop; omega g The angular speed is the corresponding angular speed of the power grid at the power frequency of 50 Hz; delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage. The left side of the equal sign of the above formula (9) represents the acceleration area, and the right side represents the deceleration area.
Therefore, the limit cutting angle delta can be obtained by carrying out iterative solution on the formula (9) cr Further, the limit cutting angle delta is obtained cr 。
The limit cut angle delta obtained by the above formula (9) is then used cr The expression of (2) is substituted into the transient stability analysis model, and the corresponding time quantity, namely limit cutting time, can be obtained by integral solution, wherein the limit cutting time t cr The expression of (2) is as follows:
in the above, t cr For extreme excision time, delta cr For extreme cut angle, delta A Is the electromagnetic power angle at the point A, T pll Is of similar mechanical inertia, T m In the form of a mechanical-like torque,to be the electromagnetic torque of the pre-fault class omega g For the corresponding angular speed of the power frequency of 50Hz of the power grid, S I For accelerating area, delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Step S105: and evaluating the transient stability margin value of the full-power wind turbine generator according to the limit cutting time.
If the fault removal time is less than or equal to the limit removal time, the full-power wind turbine generator can be restored to be stable after the fault is removed, and the corresponding limit removal time is the transient stability criterion of the full-power wind turbine generator on the premise of theoretical calculation based on the equal area rule. The corresponding evaluation criteria were: the larger the limit cutting time of the fault is, the larger the transient stability margin value of the corresponding full-power wind turbine generator is estimated.
Step S106: and in response to the transient stability margin value not meeting the requirement, increasing the transient stability margin value of the full-power wind turbine by increasing the limit cutting time.
From the limit cutting time t cr The calculation flow of (2) is reversely deduced, and the salient elements mainly influencing the transient stability margin can be divided into the following five elements: phase locked loop bandwidth, fault conditions and electrical distance from the fan port to the fault point.
How to adjust the five elements to improve the transient stability margin value of the full-power wind turbine generator is described below.
(1) Phase-locked loop bandwidth
When the phase-locked loop is in the same fault condition, the action starting time of the phase-locked loop is different and is related to the value of the limit cutting time, the corresponding time scale is represented by the phase-locked loop bandwidth variable, and a relation graph of the phase-locked loop bandwidth and the limit cutting time is provided, wherein the relation graph is shown in figure 3. As can be seen from fig. 3, when the voltage drop depth is the same, the larger the bandwidth of the phase-locked loop is, the smaller the value of the limit cut-off time is, and the smaller the transient stability margin is, so that the bandwidth frequency of the phase-locked loop can be properly reduced to expand the transient stability margin when a fault occurs.
(2) Failure condition
The analysis of the fault condition in this embodiment may be equivalent to the following three elements: 1. the voltage drop depth caused by faults; 2. the strength of the power grid, namely the magnitude of the short circuit ratio value; 3. the magnitude of the phase angle abrupt change radian caused by the fault. When the bandwidth of the phase-locked loop is fixed to be 50Hz, the three variables are combined with the limit cutting time t cr The relationship of (2) is also shown for analysis in FIG. 4.
From the analysis of fig. 4, it can be known that when the power grid strength is certain and no phase angle abrupt change scene exists, the deeper the voltage drop, the smaller the corresponding fault limit cutting time, the smaller the transient stability margin, and the transient stability is not facilitated; the strong power grid is more beneficial than the weak power grid in expanding transient stability margin, and is more beneficial particularly when the voltage drop depth is smaller, and the transient stability margin is most obvious under a certain phase angle abrupt change. For the aspect of phase angle mutation, when a system fault causes a phase angle mutation quantity with a certain negative value, the expansion of the transient stability margin is beneficial, but when the phase angle mutation quantity is larger and exceeds a certain quantity, the corresponding transient stability margin is smaller than the transient stability margin in the situation of no phase angle mutation, and as the phase angle mutation quantity is increased, the transient stability margin is reduced.
(3) Electrical distance from fan port to failure point
The present application mainly discusses X (the sum Z of the impedance of the corresponding transmission lines) ∑ Reactance component of (c) value pair T m (mechanical torque-like) influence. When the faults are assumed to uniformly occur at grid connection points, the bandwidth of the phase-locked loop is set to be 50Hz, the voltage drops to a value of 0.2pu due to the occurrence of the faults, and when no phase angle abrupt change scene exists, a fan port is drawn until the faults occurThe characteristic curve of reactance at the point of origin versus limit cut-off time is shown in fig. 5. As can be seen from FIG. 5, the limit cut-off time of the fault decreases with the electrical distance from the fan port to the fault point, regardless of whether the grid condition is a strong grid or a weak grid, and the analysis of the mechanism is mainly due to the influence of the electrical distance from the fan port to the fault point on T m When the electrical distance from the fan port to the fault occurrence point becomes larger, T m The value becomes larger, correspondingly, the acceleration area of the working condition is increased and the deceleration area is reduced under the same condition, so that the transient stability margin is continuously reduced, and the conclusion is consistent with the conclusion of the figure 5.
As can be seen from the above, in the embodiment, step S106 may increase the transient stability margin value of the full-power wind turbine generator set by: by reducing the phase-locked loop bandwidth within a set range; and/or reducing the voltage sag depth; and/or to increase grid strength conditions; and/or the adjusting system is set to enable the full-power wind turbine generator to have a phase angle abrupt change value with a set negative value under the fault working condition; and/or reducing the electrical distance from the fan port to the fault point in a set range to improve the transient stability margin value of the full-power wind turbine generator. That is, the transient stability margin value of the full-power wind turbine generator can be improved by adjusting one or more of the five influencing elements, so that the transient stability of the system is improved.
According to the application, a criterion is provided for transient stability of the full-power wind turbine generator by utilizing an electromagnetic power angle theory, namely, an expression of transient stability criterion-fault limit removal time of the full-power wind turbine generator under a voltage deep drop working condition is obtained by carrying out iterative solution and integral solution on a transient stability analysis model formula under a limit condition based on the electromagnetic power angle theory. And then according to the variables involved in the expression, the objective of improving the transient stability margin value and further improving the transient stability of the system is achieved by adjusting the physical quantity corresponding to the variables.
Fig. 6 is a schematic structural diagram of a transient stability analysis and adjustment device for a full-power wind turbine according to an embodiment of the present application, where the device includes: the model acquisition unit 610, the failure curve determination unit 620, the acceleration/deceleration area calculation unit 630, the limit cut-off time calculation unit 640, the margin value evaluation unit 650, and the margin value adjustment unit 660 are connected in this order therebetween.
The model obtaining unit 610 is configured to obtain a transient stability analysis model under a limit condition based on an electromagnetic power angle theory of the full-power wind turbine.
The fault curve determining unit 620 is configured to determine electromagnetic angle characteristic curves before and after occurrence of the fault condition.
The acceleration/deceleration area calculating unit 630 is configured to calculate an acceleration area and a deceleration area amount based on the electromagnetic angle characteristic curve.
The limit cut-off time calculation unit 640 is configured to calculate a limit cut-off angle and a corresponding limit cut-off time when the acceleration area and the deceleration area amount are equal.
The margin value evaluation unit 650 is used for evaluating the transient stability margin value of the full-power wind turbine generator according to the limit cutting time.
The margin value adjusting unit 660 is configured to increase the transient stability margin value of the full-power wind turbine generator set by increasing the limit cut time in response to the transient stability margin value not meeting a requirement.
Preferably, the stability analysis model acquired by the model acquisition unit 610 is as follows
In the above formula: t (T) pll Is a mechanical-like inertia that is analogous to the mechanical inertia in the synchronous generator rotor equation of motion; t (T) m Is a mechanical-like torque that is analogous to the mechanical torque in the synchronous generator rotor equation of motion; t (T) e Is an electromagnetic-like torque that is analogous to the electromagnetic torque in the synchronous generator rotor equation of motion; d (D) pll Is an electrical-like damping that is analogous to the amount of electrical damping in the synchronous generator rotor motion equation; omega pll The phase-locked angular velocity is output by the phase-locked loop; omega g The angular speed is the corresponding angular speed of the power grid at the power frequency of 50 Hz; omega pll -ω g ＝Δω pll ，δ pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, the limit cut time calculating unit 640, when calculating the limit cut angle when the acceleration area and the deceleration area are equal in amount, includes:
and carrying out iterative solution by using the following formula to obtain an expression of the limit cutting angle:
in the above, delta cr For extreme cut angle, delta A Is the electromagnetic power angle delta at the point A J Is the electromagnetic power angle at the J point, T m In the form of a mechanical-like torque,for the pre-fault class electromagnetic torque +.>As the electromagnetic torque after failure, delta pll ＝θ pll -θ pcc ，δ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage. Preferably, the limit cut time calculating unit 640 calculates the limit cut time corresponding to the limit cut angle by the following equation:
in the above, t cr For extreme excision time, delta cr For extreme cut angle, delta A Is the electromagnetic power angle at the point A, T pll Is of similar mechanical inertia, T m In the form of a mechanical-like torque,to be the electromagnetic torque of the pre-fault class omega g For the corresponding angular speed of the power frequency of 50Hz of the power grid, S I For accelerating area, delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
Preferably, the above-mentioned margin value evaluation unit 650 evaluates the transient stability margin value of the full power wind turbine generator according to the limit cut time, including: the larger the limit cutting time of the fault is, the larger the transient stability margin value of the corresponding full-power wind turbine generator is estimated.
Preferably, in response to the transient stability margin value not meeting the requirement, the margin value adjusting unit 660 is specifically configured to: by reducing the phase-locked loop bandwidth within a set range; and/or reducing the voltage sag depth; and/or to increase grid strength conditions; and/or the adjusting system is set to enable the full-power wind turbine generator to have a phase angle abrupt change value with a set negative value under the fault working condition; and/or reducing the electrical distance from the fan port to the fault point in a set range to improve the transient stability margin value of the full-power wind turbine generator.
The detailed description of each unit may be referred to the corresponding description in the foregoing method embodiments, and will not be repeated here.
According to the transient stability analysis and adjustment device for the full-power wind turbine, provided by the application, the electromagnetic power angle theory is utilized to provide a criterion for transient stability of the full-power wind turbine, namely, the transient stability analysis model formula under the limit condition is subjected to iterative solution and integral solution based on the electromagnetic power angle theory to obtain an expression of the transient stability criterion-fault limit excision time of the full-power wind turbine under the voltage deep drop working condition. And then according to the variables involved in the expression, the objective of improving the transient stability margin value and further improving the transient stability of the system is achieved by adjusting the physical quantity corresponding to the variables.
Fig. 7 is a schematic diagram of an electronic device according to an embodiment of the present application. The electronic device shown in fig. 7 is a general-purpose data processing apparatus comprising a general-purpose computer hardware structure including at least a processor 801 and a memory 802. The processor 801 and the memory 802 are connected by a bus 803. The memory 802 is adapted to store one or more instructions or programs executable by the processor 801. The one or more instructions or programs are executed by the processor 801 to implement the steps in the source network load store interactive method described above.
The processor 801 may be a separate microprocessor or a collection of one or more microprocessors. Thus, the processor 801 performs the process of processing data and controlling other devices by executing the commands stored in the memory 802, thereby executing the method flow of the embodiment of the present application as described above. The bus 803 connects the above-described components together, while connecting the above-described components to a display controller 804 and a display device and an input/output (I/O) device 805. Input/output (I/O) devices 805 may be a mouse, keyboard, modem, network interface, touch input device, somatosensory input device, printer, and other devices known in the art. Typically, input/output (I/O) devices 805 are connected to the system through input/output (I/O) controllers 806.
The memory 802 may store software components such as an operating system, communication modules, interaction modules, and application programs, among others. Each of the modules and applications described above corresponds to a set of executable program instructions that perform one or more functions and methods described in the embodiments of the application.
The embodiment of the application also provides a computer readable storage medium, on which a computer program is stored, which is executed by a processor to realize the steps of the transient stability analysis adjustment method of the full-power wind turbine generator.
The embodiment of the application also provides a computer program product, which comprises a computer program/instruction, wherein the computer program/instruction realizes the steps of the transient stability analysis and adjustment method of the full-power wind turbine generator set when being executed by a processor.
In summary, the transient stability analysis adjustment method and device for the full-power wind turbine provided by the application provide a criterion for transient stability of the full-power wind turbine by utilizing the electromagnetic power angle theory, namely, the transient stability analysis model formula under the limit condition is subjected to iterative solution and integral solution based on the electromagnetic power angle theory to obtain an expression of the transient stability criterion-fault limit excision time of the full-power wind turbine under the voltage deep drop working condition. And then according to the variables involved in the expression, the objective of improving the transient stability margin value and further improving the transient stability of the system is achieved by adjusting the physical quantity corresponding to the variables.
Preferred embodiments of the present application are described above with reference to the accompanying drawings. The many features and advantages of the embodiments are apparent from the detailed specification, and thus, it is intended by the appended claims to cover all such features and advantages of the embodiments which fall within the true spirit and scope thereof. Further, since numerous modifications and changes will readily occur to those skilled in the art, it is not desired to limit the embodiments of the application to the exact construction and operation illustrated and described, and accordingly, all suitable modifications and equivalents may be resorted to, falling within the scope thereof.
It will be appreciated by those skilled in the art that 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 flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations 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.
The foregoing description of the embodiments has been provided for the purpose of illustrating the general principles of the application, and is not meant to limit the scope of the application, but to limit the application to the particular embodiments, and any modifications, equivalents, improvements, etc. that fall within the spirit and principles of the application are intended to be included within the scope of the application.

Claims (10)

1. The transient stability analysis and adjustment method for the full-power wind turbine generator is characterized by comprising the following steps of:
acquiring a transient stability analysis model under a limit condition based on an electromagnetic power angle theory of the full-power wind turbine generator;
determining electromagnetic power angle characteristic curves before and after the occurrence of the fault working condition;
calculating an acceleration area and a deceleration area based on the electromagnetic power angle characteristic curve;
calculating a limit cutting angle and a corresponding limit cutting time when the acceleration area and the deceleration area are equal;
evaluating the transient stability margin value of the full-power wind turbine according to the limit cutting time;
and in response to the transient stability margin value not meeting the requirement, increasing the transient stability margin value of the full-power wind turbine by increasing the limit cutting time.
2. The method for analyzing and adjusting the transient stability of the full-power wind turbine generator set according to claim 1, wherein the transient stability analysis model is as follows:
in the above formula: t (T) pll Is a mechanical-like inertia that is analogous to the mechanical inertia in the synchronous generator rotor equation of motion; t (T) m Is a mechanical-like torque that is analogous to the mechanical torque in the synchronous generator rotor equation of motion; t (T) e Is an electromagnetic-like torque that is analogous to the electromagnetic torque in the synchronous generator rotor equation of motion; d (D) pll Is an electrical-like damping that is analogous to the amount of electrical damping in the synchronous generator rotor motion equation; omega pll The phase-locked angular velocity is output by the phase-locked loop; omega g The angular speed is the corresponding angular speed of the power grid at the power frequency of 50 Hz; omega pll -ω g ＝Δω pll ，δ pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
3. The method for transient stability analysis adjustment of a full power wind turbine according to claim 2, wherein calculating a limit cut angle when the acceleration area and the deceleration area are equal comprises:
and carrying out iterative solution by using the following formula to obtain an expression of the limit cutting angle:
in the above, delta cr For extreme cut angle, delta A Is the electromagnetic power angle delta at the point A J Is the electromagnetic power angle at the J point, T m In the form of a mechanical-like torque,for the pre-fault class electromagnetic torque +.>As the electromagnetic torque after failure, delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
4. The transient stability analysis adjustment method of a full power wind turbine according to claim 3, wherein the limit cut time corresponding to the limit cut angle is calculated by the following formula:
in the above, t cr For extreme excision time, delta cr For extreme cut angle, delta A Is the electromagnetic power angle at the point A, T pll Is of similar mechanical inertia, T m In the form of a mechanical-like torque,for the pre-fault class electromagnetic torque omega g For the corresponding angular speed of the power frequency of 50Hz of the power grid, S I For accelerating area, delta pll ＝θ pll -θ pcc ，θ pll For phase locked loop controlled output phase, θ pcc The synchronous phase angle of the grid-connected point is the corresponding phase of the grid-connected point voltage.
5. The method for transient stability analysis adjustment of a full-power wind turbine according to claim 1, wherein the evaluating the transient stability margin value of the full-power wind turbine according to the limit cut-off time comprises: the larger the limit cutting time of the fault is, the larger the transient stability margin value of the corresponding full-power wind turbine generator is estimated.
6. The method of claim 5, wherein in response to the transient stability margin value not meeting a requirement, the increasing the transient stability margin value of the full power wind turbine by increasing the limit cut time comprises:
by reducing the phase-locked loop bandwidth within a set range; and/or reducing the voltage sag depth; and/or to increase grid strength conditions; and/or the adjusting system is set to enable the full-power wind turbine generator to have a phase angle abrupt change value with a set negative value under the fault working condition; and/or reducing the electrical distance from the fan port to the fault point in a set range to improve the transient stability margin value of the full-power wind turbine generator.
7. The transient stability analysis and adjustment device of the full-power wind turbine generator is characterized by comprising:
the model acquisition unit is used for acquiring a transient stability analysis model under a limit condition based on the electromagnetic power angle theory of the full-power wind turbine generator;
the fault curve determining unit is used for determining electromagnetic power angle characteristic curves before and after the fault working condition occurs;
an acceleration/deceleration area calculation unit for calculating an acceleration area and a deceleration area amount based on the electromagnetic work angle characteristic curve;
a limit cut-off time calculation unit configured to calculate a limit cut-off angle and a corresponding limit cut-off time when the acceleration area and the deceleration area are equal in amount;
the margin value evaluation unit is used for evaluating the transient stability margin value of the full-power wind turbine generator according to the limit cutting time;
and the margin value adjusting unit is used for increasing the transient stability margin value of the full-power wind turbine generator by increasing the limit cutting time in response to the transient stability margin value not meeting the requirement.
8. 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 of any of claims 1 to 6 when the computer program is executed by the processor.
9. A computer readable storage medium, on which a computer program is stored, characterized in that the computer program, when being executed by a processor, implements the steps of the method according to any one of claims 1 to 6.
10. A computer program product comprising computer programs/instructions which, when executed by a processor, implement the steps of the method of any of claims 1 to 6.
CN202310987654.7A 2023-08-07 2023-08-07 Transient stability analysis and adjustment method and device for full-power wind turbine generator Pending CN117039996A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310987654.7A CN117039996A (en) 2023-08-07 2023-08-07 Transient stability analysis and adjustment method and device for full-power wind turbine generator

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310987654.7A CN117039996A (en) 2023-08-07 2023-08-07 Transient stability analysis and adjustment method and device for full-power wind turbine generator

Publications (1)

Publication Number Publication Date
CN117039996A true CN117039996A (en) 2023-11-10

Family

ID=88642328

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310987654.7A Pending CN117039996A (en) 2023-08-07 2023-08-07 Transient stability analysis and adjustment method and device for full-power wind turbine generator

Country Status (1)

Country Link
CN (1) CN117039996A (en)

Similar Documents

Publication Publication Date Title
CN109193752B (en) Method for optimizing low-frequency oscillation control parameters of virtual inertia-containing doubly-fed wind turbine grid-connected system
US9300142B2 (en) Method for emulation of synchronous machine
CN108649847A (en) Motor PI attitude conirol methods based on frequency method and fuzzy control
CN115313524B (en) Photovoltaic power generation grid-connected control method and system based on network-building type converter
CN115313419B (en) Direct-drive wind power system control method and system based on network-building type converter
CN110932319B (en) Method and system for inhibiting subsynchronous oscillation of doubly-fed wind turbine generator
Errami et al. Design of a nonlinear backstepping control strategy of grid interconnected wind power system based PMSG
CN106786673B (en) The suppressing method and device of double-fed blower compensated transmission system subsynchronous resonance
CN107370181B (en) Grid connection control method and system
CN111478365B (en) Optimization method and system for control parameters of virtual synchronizer of direct-drive wind turbine generator
JP4350998B2 (en) Transient stability analysis method for induction generator
CN114841005A (en) Method for evaluating equivalent inertia of asynchronous motor in inertia response stage
CN105450119B (en) Accelerated method of the induction machine of V/f controls in weak magnetic area
CN112787325A (en) Quantitative evaluation method for transient synchronization stability of permanent magnet synchronous wind driven generator based on Lyapunov direct method
CN112350346A (en) Method for establishing transient stability simulation model of power system
CN116632866A (en) Hybrid energy storage self-adaptive inertia VSG control method for liquid flow super-capacity lithium battery
CN111654039A (en) Method and system for judging sub/super synchronous oscillation stability of double-fed wind power grid-connected system
Amutha et al. Effect of modeling of induction generator based wind generating systems on determining CCT
CN117039996A (en) Transient stability analysis and adjustment method and device for full-power wind turbine generator
CN112886611B (en) Subsynchronous oscillation suppression method for direct-drive fan grid-connected system
Papangelis et al. Design principles of wind turbine inertia emulators for autonomous systems
CN114243762A (en) Analysis and control method for fan grid connection
CN112260225B (en) Method and system for calculating three-phase short-circuit current of doubly-fed wind turbine generator based on combined model
Abulizi et al. Research of current control strategies for doubly-fed wind power generation system
CN108011368A (en) The control method of complicated electric power system maximum angle of oscillation is limited based on EEAC

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