CN113241779B - Stability analysis method and device for direct-drive wind power plant grid-connected system - Google Patents

Stability analysis method and device for direct-drive wind power plant grid-connected system Download PDF

Info

Publication number
CN113241779B
CN113241779B CN202110580826.XA CN202110580826A CN113241779B CN 113241779 B CN113241779 B CN 113241779B CN 202110580826 A CN202110580826 A CN 202110580826A CN 113241779 B CN113241779 B CN 113241779B
Authority
CN
China
Prior art keywords
grid
direct
voltage
equivalent
network side
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.)
Active
Application number
CN202110580826.XA
Other languages
Chinese (zh)
Other versions
CN113241779A (en
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.)
Guangdong Power Grid Co Ltd
Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd
Original Assignee
Guangdong Power Grid Co Ltd
Electric Power Dispatch Control Center of Guangdong Power Grid 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 Guangdong Power Grid Co Ltd, Electric Power Dispatch Control Center of Guangdong Power Grid Co Ltd filed Critical Guangdong Power Grid Co Ltd
Priority to CN202110580826.XA priority Critical patent/CN113241779B/en
Publication of CN113241779A publication Critical patent/CN113241779A/en
Application granted granted Critical
Publication of CN113241779B publication Critical patent/CN113241779B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

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/24—Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241—The oscillation concerning frequency
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00—Administration; Management
    • G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063—Operations research, analysis or management
    • G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Energy or water supply
    • 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
    • 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
    • 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
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Health & Medical Sciences (AREA)
  • Physics & Mathematics (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Power Engineering (AREA)
  • Development Economics (AREA)
  • Marketing (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Tourism & Hospitality (AREA)
  • Educational Administration (AREA)
  • General Business, Economics & Management (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Game Theory and Decision Science (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Control Of Eletrric Generators (AREA)

Abstract

The invention discloses a stability analysis method and device for a direct-drive wind power plant grid-connected system, and relates to the technical field of stability analysis of power systems. The method comprises the following steps: acquiring a state equation and an algebraic equation of a grid-connected system, and performing simultaneous combination by taking a grid-side voltage as a carrier; representing a phase-locked loop and a network side control link in the form of equivalent impedance and an equivalent voltage source to obtain an equivalent circuit and equivalent impedance of a direct-drive wind power plant; combining the equivalent impedance to obtain a stability margin expression of the grid-connected system; obtaining the normalized sensitivity of different parameters in the grid-connected system to the stability margin of the system according to the stability margin expression; and quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin, and screening out dominant influence factors. According to the method, the stability margin of the direct-drive wind power plant grid-connected system can be deduced by constructing an analytical expression of the equivalent impedance of the direct-drive fan and adopting a Nyquist stability criterion, so that the influence degree of different parameters on the system stability is quantized.

Description

Stability analysis method and device for direct-drive wind power plant grid-connected system
Technical Field
The invention relates to the technical field of power system stability analysis, in particular to a stability analysis method of a direct-drive wind power plant grid-connected system.
Background
With the rapid development of renewable energy sources such as wind power and the like, the wind power has a wide development prospect, but the power electronization trend of a power system is more obvious due to the improvement of the capacity of a fan, and new challenges are brought to the stable operation of the power system. The direct-drive fan is one of the mainstream fan types, and the direct-drive fan is easy to cause the problem of oscillation instability when being connected to a power grid, such as the phenomenon of subsynchronous oscillation in the Hami area of Xinjiang 2015.
Some scholars have achieved certain results around how to solve the oscillation instability. In the existing research methods, such as characteristic value analysis, impedance analysis, time domain simulation, etc., the system stability can be judged only by calculating characteristic values or impedances or after simulation, and the stability degree of the direct-drive wind farm grid-connected system by different parameters cannot be intuitively reflected. In addition, for the dominant influence factors of the direct-drive wind turbine, the existing research has not reached a consensus. In order to quantify the influence degree of different parameters on the stability of the grid-connected system of the direct-driven wind farm and identify main influence factors, further discussion is needed.
Disclosure of Invention
The invention aims to provide a stability analysis method of a direct-drive wind power plant grid-connected system, which is characterized in that the stability margin of the direct-drive wind power plant grid-connected system is deduced by constructing an analytical expression of equivalent impedance of a direct-drive fan and adopting Nyquist stability criterion so as to quantify the influence degree of different parameters on the system stability.
In order to achieve the above object, an embodiment of the present invention provides a method for analyzing stability of a direct-drive wind farm grid-connected system, including:
acquiring a state equation and an algebraic equation of a grid-connected system, and performing simultaneous combination by taking the network side voltage as a carrier to obtain an equivalent network side voltage equation considering a phase-locked loop and a network side control link;
representing the phase-locked loop and the network side control link in the form of equivalent impedance and an equivalent voltage source to obtain an RLC equivalent circuit and corresponding equivalent impedance of a direct-drive wind power plant;
combining the equivalent impedance and the Nyquist stability criterion to obtain a stability margin expression of the grid-connected system;
obtaining the normalized sensitivity of different parameters in the grid-connected system to the stability margin of the system according to the stability margin expression;
and quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin, and screening out dominant influence factors.
Preferably, the state equation and the algebraic equation of the grid-connected system are obtained, and the grid-side voltage is taken as a carrier to perform simultaneous combination, so as to obtain an equivalent grid-side voltage equation considering the phase-locked loop and the grid-side control link, wherein the adopted expression comprises:
Figure BDA0003085968620000021
Figure BDA0003085968620000022
Figure BDA0003085968620000023
Figure BDA0003085968620000024
Figure BDA0003085968620000025
wherein u is g s* Is the net side equivalent voltage; u. u g Is the net side voltage; i.e. i g A grid side current; i.e. i m Is the machine side current; u. of m Is the machine side voltage; s is a laplace operator; superscript s represents the grid dq coordinate system; subscripts d and q denote d and q-axis components, respectively; subscript 0 represents the initial value of the corresponding variable; theta pll Is a phase lock angle; u shape g Is the net side voltage amplitude; l is a radical of an alcohol g Inductance of the incoming line reactor; u shape dc Is the direct current side voltage of the converter; q gpu Is a per unit value of reactive power; c dc Is a DC bus capacitor; k is ppll Is the proportionality coefficient of the phase-locked loop; k is ipll Is the integral coefficient of the phase-locked loop; k is pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k is ig The integral coefficient of the current inner loop of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; k idc Is the integral coefficient of the outer ring of the direct current voltage; u shape gb Is a reference value of the network side voltage; I.C. A gb Is the reference value of the current at the network side; u shape dcb Is a reference value of the direct current voltage; s gb Is the reference value of the network side power.
Preferably, the phase-locked loop and the network side control link are characterized in the form of an equivalent impedance and an equivalent voltage source, and an RLC equivalent circuit and a corresponding equivalent impedance of the direct-drive wind farm are obtained, wherein the expression of the equivalent impedance is as follows:
Figure BDA0003085968620000031
the expression of the equivalent voltage source is as follows:
Figure BDA0003085968620000032
wherein R is l Is a line resistance; l is l Is a line inductance; l is a radical of an alcohol T Is a transformer inductor; l is p Is an equivalent inductance of the direct-drive fan; n is the number of fans; k is ig The integral coefficient of a current inner loop of the grid-side converter is obtained; u shape dc Is the direct current side voltage of the converter; u shape gb At network side voltageA reference value; I.C. A gb Is a reference value of the current on the network side; c dc Is a DC bus capacitor; c p Is the equivalent capacitance of the direct-drive fan; k pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; k idc The integral coefficient of the outer ring of the direct current voltage is obtained; omega is angular frequency; omega s Is the rotation angular velocity of the synchronous magnetic field; u shape dcb Is a reference value of the direct current voltage; s gb Is a reference value of the network side power; u shape dc0 Is a reference value of the direct current voltage; j is an imaginary unit.
Preferably, the equivalent impedance and the nyquist stability criterion are combined to obtain a stability margin expression of the grid-connected system, where the stability margin expression of the grid-connected system is:
Figure BDA0003085968620000033
Figure BDA0003085968620000034
wherein, K g Representing an amplitude margin; γ represents a phase margin; r is p Representing an equivalent impedance; r is s Representing the system equivalent resistance; r is l Is a line resistance; l is a radical of an alcohol l Is a line inductance; u shape gb Is a reference value of the network side voltage; k pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; i is gb Is a reference value of the current on the network side; s. the gb Is a reference value of the network side power; l is a radical of an alcohol p Is an equivalent inductance of the direct-drive fan; l is s Is the equivalent inductance of the system; and n is the number of the fans.
Preferably, the normalized sensitivity of different parameters in the grid-connected system to the system stability margin is obtained according to the stability margin expression, where the expression of the normalized sensitivity of different parameters to the system stability margin is:
Figure BDA0003085968620000041
wherein a represents the parameters of the grid-connected system, m is the number of the parameters, and H is the amplitude margin or the phase margin.
The embodiment of the invention provides a stability analysis device of a direct-drive wind power plant grid-connected system, which comprises:
the simultaneous combination module is used for acquiring a state equation and an algebraic equation of a grid-connected system, and performing simultaneous combination by taking the network side voltage as a carrier to obtain an equivalent network side voltage equation considering a phase-locked loop and a network side control link;
the characterization module is used for characterizing the phase-locked loop and the network side control link in the form of equivalent impedance and an equivalent voltage source to obtain an RLC equivalent circuit and corresponding equivalent impedance of the direct-drive wind power plant;
the stability margin obtaining module is used for combining the equivalent impedance and the Nyquist stability criterion to obtain a stability margin expression of the grid-connected system;
the sensitivity acquisition module is used for acquiring the normalized sensitivity of different parameters in the grid-connected system to the stability margin of the system according to the stability margin expression;
and the influence degree screening module is used for quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin and screening out leading influence factors.
Preferably, the expression adopted by the simultaneous combination module includes:
Figure BDA0003085968620000042
Figure BDA0003085968620000043
Figure BDA0003085968620000051
Figure BDA0003085968620000052
Figure BDA0003085968620000053
wherein u is g s* Is the net side equivalent voltage; u. u g Is the net side voltage; i.e. i g A grid side current; i.e. i m Is the machine side current; u. of m Is the machine side voltage; s is a Laplace operator; superscript s represents the grid dq coordinate system; subscripts d and q denote d and q-axis components, respectively; subscript 0 represents the initial value of the corresponding variable; theta.theta. pll Is a phase lock angle; u shape g Is the net side voltage amplitude; l is a radical of an alcohol g The inductance of the incoming line reactor; u shape dc Is the direct current side voltage of the converter; q gpu Is a per unit value of reactive power; c dc Is a direct current bus capacitor; k is ppll Is the proportionality coefficient of the phase-locked loop; k ipll Is the integral coefficient of the phase-locked loop; k pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k ig The integral coefficient of a current inner loop of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; k is idc Is the integral coefficient of the outer ring of the direct current voltage; u shape gb Is a reference value of the network side voltage; i is gb Is a reference value of the current on the network side; u shape dcb Is a reference value of the direct current voltage; s gb Is the reference value of the network side power.
Preferably, in the characterization module, the expression of the equivalent impedance is:
Figure BDA0003085968620000054
the expression of the equivalent voltage source is as follows:
Figure BDA0003085968620000055
wherein R is l Is a line resistance; l is l Is a line inductance; l is T Is a transformer inductor; l is p The equivalent inductance of the direct-drive fan is obtained; n is the number of fans; k ig The integral coefficient of the current inner loop of the grid-side converter is obtained; u shape dc Is the converter DC side voltage; u shape gb Is a reference value of the network side voltage; i is gb Is a reference value of the current on the network side; c dc Is a DC bus capacitor; c p Is the equivalent capacitance of the direct-drive fan; k pg The proportionality coefficient of the current inner ring of the grid-side converter is obtained; k idc The integral coefficient of the outer ring of the direct current voltage is obtained; omega is angular frequency; omega s Is the rotation angular velocity of the synchronous magnetic field; u shape dcb Is a reference value of the direct current voltage; s gb The reference value is the network side power; u shape dc0 Is a reference value of the direct current voltage; j is an imaginary unit.
Preferably, in the stability margin obtaining module, the stability margin expression of the grid-connected system is as follows:
Figure BDA0003085968620000061
Figure BDA0003085968620000062
wherein, K g Representing an amplitude margin; γ represents a phase margin; r is p Representing an equivalent impedance; r is s Representing the equivalent resistance of the system; r is l Is a line resistance; l is a radical of an alcohol l Is a line inductance; u shape gb Is a reference value of the network side voltage; k pg The proportionality coefficient of the current inner ring of the grid-side converter is obtained; k is pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; i is gb Is a reference value of the current on the network side; s gb Is a reference value of the network side power; l is p The equivalent inductance of the direct-drive fan is obtained; l is s Is the equivalent inductance of the system; and n is the number of the fans.
Preferably, in the sensitivity acquisition module, an expression of normalized sensitivity of the different parameters to the system stability margin is as follows:
Figure BDA0003085968620000063
wherein a represents the parameters of the grid-connected system, m is the number of the parameters, and H is the amplitude margin or the phase margin.
Compared with the prior art, the invention has the following beneficial effects:
according to the method and the device for analyzing the stability of the direct-drive wind power plant grid-connected system, the equivalent grid-side voltage considering the control of the phase-locked loop and the grid-side converter is redefined from the aspect of circuit properties through a simultaneous grid-side circuit combining equation and a control equation, the grid-side converter control is represented in the form of equivalent impedance and an equivalent voltage source, and an analytical expression of the equivalent impedance is deduced; calculating the amplitude margin and the phase margin of the direct-drive wind power plant grid-connected system according to the Nyquist stability criterion and the equivalent impedance; and measuring the influence degree of different parameters on the grid-connected stability margin by utilizing the normalized sensitivity so as to adjust the parameters and realize the stable operation of the grid-connected system.
Drawings
In order to more clearly illustrate the technical solution of the present invention, the drawings required to be used in the embodiments will be briefly described below, and obviously, the drawings in the following description are only some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a schematic flow chart of a stability analysis method of a direct-drive wind farm grid-connected system according to an embodiment of the present invention;
FIG. 2 is a simplified model diagram of a direct-drive wind farm grid-connected system according to an embodiment of the present invention;
fig. 3 is a structural diagram of a direct-drive fan according to an embodiment of the present invention;
FIG. 4 is a diagram of an equivalent circuit of a direct-drive wind farm RLC according to an embodiment of the present invention;
fig. 5 is an equivalent circuit diagram of a direct-drive wind farm grid-connected system according to an embodiment of the present invention;
fig. 6 is a schematic structural diagram of a stability analysis device of a direct-drive wind farm grid-connected system according to an embodiment of the present invention;
fig. 7 is a schematic structural diagram of a computer terminal device according to an embodiment of the present invention.
Detailed Description
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 only a part of the embodiments of the present invention, and not all of the embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
It should be understood that the step numbers used herein are for convenience of description only and are not intended as limitations on the order in which the steps are performed.
It is to be understood that the terminology used in the description of the invention herein is for the purpose of describing particular embodiments only and is not intended to be limiting of the invention. As used in the specification of the present invention and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise.
The terms "comprises" and "comprising" indicate the presence of the described features, integers, steps, operations, elements, and/or components, but do not preclude the presence or addition of one or more other features, integers, steps, operations, elements, components, and/or groups thereof.
The term "and/or" refers to and includes any and all possible combinations of one or more of the associated listed items.
Referring to fig. 1, fig. 1 is a schematic flow chart of a stability analysis method of a direct-drive wind farm grid-connected system according to an embodiment of the present invention. The method for analyzing the stability of the direct-drive wind power plant grid-connected system provided by the embodiment comprises the following steps:
and S110, acquiring a state equation and an algebraic equation of the grid-connected system, and performing simultaneous combination by taking the grid-side voltage as a carrier to obtain an equivalent grid-side voltage equation considering the phase-locked loop and a grid-side control link.
In a specific embodiment, a simplified model diagram of a direct-drive wind farm grid-connected system is shown in fig. 2, wherein, ignoring the current collection lines in the direct-drive wind farm, all the wind turbines are connected in parallel to a common collection bus, i.e. the whole wind farm is replaced by the same type of low-capacity wind turbine, and the wind farm is interconnected with the grid through a transformer and line impedance.
Because the grid-side converter is controlled by constant direct-current voltage, the direct-current voltage dynamic can be ignored, so that a direct-current part of the grid-side converter is equivalently represented by a direct-current voltage source, and the structure diagram of the direct-drive fan is shown in fig. 3.
In one embodiment, the equivalent impedance estimation process for a direct drive wind turbine is as follows.
Considering the phase-locked loop dynamics, the phase-locked angle theta in the s-domain pll With stator q-axis voltage u s gq The relationship between them is as follows.
Figure BDA0003085968620000081
Wherein, K ppll And K ipll Proportional and integral coefficients, U, of the phase-locked loop, respectively g And (4) the grid side voltage amplitude value is obtained, s is a Laplace operator, and superscript s represents a grid dq coordinate system.
And the grid dq coordinate system variable and the phase-locked loop dq coordinate system variable have the following relation:
Figure BDA0003085968620000082
wherein the superscript c represents the phase-locked loop dq coordinate system; subscripts d and q denote d, q-axis components, respectively; x represents a voltage or a current.
Regarding the direct drive wind turbine grid-side converter (GSC) and its control, under dq rotation coordinate system, considering the PLL, the grid-side voltage of the direct drive wind turbine grid-side converter (GSC) is as follows:
Figure BDA0003085968620000091
wherein u is g 、i g And e g Respectively, the network side voltage, the network side current and the converter alternating current side voltage; l is a radical of an alcohol g Inductance, omega, of wire-in reactors s Is the synchronous magnetic field rotation angular velocity.
The equation of the control output voltage of the grid-side converter is shown as follows:
Figure BDA0003085968620000092
wherein u is gdpu And u gqpu D and q components of the voltage per unit value of the network side are respectively; i.e. i gdpu And i gqpu D and q components of the per unit value of the network side current are respectively; e.g. of a cylinder gdpu And e gqpu D and q components of the per unit voltage value of the AC side of the converter are respectively; u shape dcref And Q gref Direct-current voltage and reactive power reference values are respectively; u shape dcpu And Q gpu The DC voltage and the reactive power per unit value are respectively; k is pg And K ig The proportional coefficient and the integral coefficient of a current inner ring of the grid-side converter are respectively; k is pdc And K idc Proportional and integral coefficients of the direct-current voltage outer ring are respectively; k pq And K iq Respectively, the proportional and integral coefficients of the reactive power outer loop.
The reactive power expression is as follows:
Figure BDA0003085968620000093
regarding the direct current link, considering the direct current dynamics between the direct-drive wind turbine grid-side converter (GSC) and the machine-side converter (MSC), the direct current voltage dynamic equation is as follows:
U dc ＝(P m -P g )/(sC dc U dc ) (6)
wherein, U dc Is the direct current side voltage of the converter; c dc Is a DC bus capacitor; pm and P g Output and input power of a direct-drive fan machine side converter and a grid side converter respectively are as follows:
Figure BDA0003085968620000094
Figure BDA0003085968620000095
wherein i m And u m Respectively the machine side current and the voltage.
In the s domain, the following formulae (7) to (8) are substituted for formula (6) to give the following formulae.
Figure BDA0003085968620000101
In general, the proportional-integral coefficients of the dq-axis outer rings are equal, i.e., K pdc ＝K pq ，K idc ＝K iq . According to the formulas (2) and (4), the control output voltage of the grid-side converter can be transferred to a grid dq coordinate system, and then the formulas (3), (5) and (9) are replaced, and linearization is carried out, so that the following formula is obtained.
Figure BDA0003085968620000102
Wherein, U gb And I gb Reference values of the network side voltage and the current are respectively; u shape dcb And S gb The reference values are direct-current voltage and network side power reference values; the subscript 0 indicates the initial value of the corresponding variable.
Figure BDA0003085968620000103
Figure BDA0003085968620000104
Figure BDA0003085968620000105
Figure BDA0003085968620000106
As can be seen from the equation (10), the equation has no net side current component on the left side of the equation equal sign, and there is a coupling effect between the net side and the machine side current and voltage, so the left side of the equation can be characterized as the controlled voltage source U s* g In the form of (1). Therefore, redefining the equivalent grid-side voltage after taking into account the phase-locked loop and the grid-side converter control can be expressed as follows:
Figure BDA0003085968620000111
and S120, representing a phase-locked loop and a network side control link in the form of equivalent impedance and an equivalent voltage source, and obtaining an RLC equivalent circuit and corresponding equivalent impedance of the direct-drive wind power plant.
In the embodiment of the present invention, the expression between the grid-side voltage and the current is obtained from equations (11) and (15) as follows:
Figure BDA0003085968620000112
based on Thevenin's theorem, the formula (16) can be sorted and combined, and represented in the form of equivalent impedance, capacitance and inductance, and simultaneously, dq axis components are combined and converted into a stationary coordinate system, and because the difference omega between the stationary coordinate system and the dq rotating coordinate system is omega s I.e. with (s-j ω) s ) Instead of s, the result is 1/(s-j ω) s ) 2 Term at divided fundamental frequency f s ＝ω s The frequency bands except the frequency band around/2 pi are very small, the related terms can be ignored, and in addition, in a phase-locked loop coordinate system, u is c gd0 ≈1，u c gq0 0, the following expression is obtained:
Figure BDA0003085968620000113
wherein,
Figure BDA0003085968620000114
I s g ＝i s gd +ji s gq 。
Figure BDA0003085968620000115
Figure BDA0003085968620000116
R GSC and C GSC Namely the equivalent impedance and the equivalent capacitance introduced by the grid-side converter. Therefore, considering the number n of wind turbines, the RLC equivalent circuit diagram of the direct-drive wind farm in the stationary coordinate system is shown in fig. 4. In the figure, R l And L l Respectively a line resistor and a line inductor; l is T Is a transformer inductor.
Obviously, at any angular frequency ω, the line and transformer inductances are much greater than L g And therefore L can be ignored g The equivalent resistance and the equivalent voltage source of the direct-drive wind power plant are respectively as follows:
Figure BDA0003085968620000121
Figure BDA0003085968620000122
and S130, combining the equivalent impedance and the Nyquist stability criterion to obtain a stability margin expression of the grid-connected system.
In the embodiment of the invention, the stability margin expression of the grid-connected system is obtained by combining the analytical expression of the equivalent impedance and the Nyquist stability criterion, and the influence of different parameters on the stability margin of the system can be intuitively reflected.
The direct-drive wind power plant grid connection is divided into two subsystems: the system comprises a direct-drive wind power plant subsystem and a power grid subsystem. The grid subsystem may typically use an ideal voltage source U s And impedance Z s (Z s ＝R s +jX s ) Series representation, and equivalent voltage source is used for the direct-drive wind power plant subsystem
Figure BDA0003085968620000123
And impedance Z p (Z p ＝R p +jX p ) In series, as shown in FIG. 5. Thus, in the s domain, the current flowing from the direct-drive wind farm to the grid is as follows.
Figure BDA0003085968620000124
In the general case of the above-mentioned,
Figure BDA0003085968620000125
and Z p (s) is stable, so the stability of the direct drive wind farm sourced current depends on the stability of the closed loop transfer function G(s). According to the control theory, by judging the open-loop transfer function Z s (s)/Z p If(s) is stable, it can be determined whether G(s) is stable.
Whereas the stability of the system can be quantified by the amplitude margin and the phase margin. For the open loop transfer function Z s (s)/Z p (s) the amplitude margin is defined as follows.
Figure BDA0003085968620000126
Wherein, ω is g The following expression is satisfied:
Figure BDA0003085968620000127
namely X s (jω g )/R s ＝X p (jω g )/R p And satisfy R s And R p Different sign or X s (jω g ) And X p (jω g ) And (4) carrying out opposite sign.
The phase margin is defined as follows:
Figure BDA0003085968620000131
wherein, ω is c The following expression is satisfied:
Figure BDA0003085968620000132
substituting the impedance expressions of the direct-drive wind power plant and the power grid to obtain approximate values of the amplitude margin and the phase margin, wherein the approximate values are as follows:
Figure BDA0003085968620000133
Figure BDA0003085968620000134
and according to the expression of the amplitude margin and the phase margin, the action condition of each parameter in the direct-drive wind power plant on the stability margin can be visually reflected.
And S140, obtaining the normalized sensitivity of different parameters in the grid-connected system to the stability margin of the system according to the stability margin expression.
S150, quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin, and screening out dominant influence factors.
In the embodiment of the invention, a normalization concept is introduced, and the normalized sensitivity of different parameters to the stability margin of the interconnected system is solved as follows:
Figure BDA0003085968620000135
wherein a represents parameters of a direct-drive wind power plant grid-connected system, m is the number of the parameters, and H is amplitude margin or phase margin.
The specific steps for screening the dominant influence factors are as follows:
(1) Forming a multi-element complex function K containing different factors based on a system equivalent circuit g (K pg ,K pdc ,R l ,n)，γ(K pg ,K pdc ,L l ,L T N) analytical expressions;
(2) Calculating the normalized sensitivity of the amplitude margin and the phase margin to each parameter one by one according to an equation (29);
(3) And comparing the normalized sensitivity of different parameters to find the parameter which influences the system stability to the maximum.
Although the larger the amplitude margin and the phase margin, the better the stability of the interconnected system, the too large margin may cause the response speed of the system to become slow, and therefore, the stability analysis method provided by the invention can be used for adjusting parameters to realize the stable operation of the system.
Referring to fig. 6, fig. 6 is a schematic structural diagram of a stability analysis device of a direct-drive wind farm grid-connected system according to an embodiment of the present invention. The stability analysis device for the direct-drive wind power plant grid-connected system provided by the embodiment comprises:
the simultaneous combination module 210 is configured to obtain a state equation and an algebraic equation of a grid-connected system, perform simultaneous combination by using a network-side voltage as a carrier, and obtain an equivalent network-side voltage equation considering a phase-locked loop and a network-side control link;
the characterization module 220 is used for characterizing the phase-locked loop and the network side control link in the form of equivalent impedance and an equivalent voltage source to obtain an RLC equivalent circuit and corresponding equivalent impedance of the direct-drive wind power plant;
the stability margin obtaining module 230 is configured to obtain a stability margin expression of the grid-connected system by combining the equivalent impedance and the nyquist stability criterion;
the sensitivity obtaining module 240 is configured to obtain normalized sensitivities of different parameters in the grid-connected system to the system stability margin according to the stability margin expression;
and the influence degree screening module 250 is used for quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin and screening out dominant influence factors.
For specific limitations of the stability analysis device of the direct-drive wind power plant grid-connected system, reference may be made to the above limitations of the stability analysis method of the direct-drive wind power plant grid-connected system, and details are not repeated here. All or part of each module in the stability analysis device of the direct-drive wind power plant grid-connected system can be realized through software, hardware and combination thereof. The modules can be embedded in a hardware form or independent from a processor in the computer device, and can also be stored in a memory in the computer device in a software form, so that the processor can call and execute operations corresponding to the modules.
Referring to fig. 7, an embodiment of the invention provides a computer terminal device, which includes one or more processors and a memory. The storage is coupled to the processor and is used for storing one or more programs, and when the one or more programs are executed by the one or more processors, the one or more processors enable the one or more processors to implement the stability analysis method of the direct-drive wind farm grid-connected system in any one of the embodiments.
The processor is used for controlling the overall operation of the computer terminal device so as to complete all or part of the steps of the stability analysis method of the direct-drive wind power plant grid-connected system. The memory is used to store various types of data to support the operation at the computer terminal device, which data may include, for example, instructions for any application or method operating on the computer terminal device, as well as application-related data. The Memory may be implemented by any type or combination of volatile and non-volatile Memory devices, such as Static Random Access Memory (SRAM), electrically Erasable Programmable Read-Only Memory (EEPROM), erasable Programmable Read-Only Memory (EPROM), programmable Read-Only Memory (PROM), read-Only Memory (ROM), magnetic Memory, flash Memory, magnetic disk, or optical disk.
In an exemplary embodiment, the computer terminal Device may be implemented by one or more Application Specific 1 integrated circuits (AS 1C), digital Signal Processors (DSP), digital Signal Processing Devices (DSPD), programmable Logic Devices (PLD), field Programmable Gate Arrays (FPGA), controllers, microcontrollers, microprocessors or other electronic components, and is configured to perform the above-mentioned stabilization analysis method for the wind farm grid-connected system, and achieve the technical effects consistent with the above-mentioned method.
In another exemplary embodiment, a computer readable storage medium including program instructions is further provided, which when executed by a processor, implement the steps of the stability analysis method of the direct drive wind farm grid-connected system in any one of the above embodiments. For example, the computer readable storage medium may be the above memory including program instructions, and the program instructions may be executed by a processor of a computer terminal device to complete the above method for stability analysis of a direct-drive wind farm grid-connected system, and achieve the technical effects consistent with the above method.
While the foregoing is directed to the preferred embodiment of the present invention, it will be understood by those skilled in the art that various changes and modifications may be made without departing from the spirit and scope of the invention.

Claims (2)

1. A stability analysis method for a direct-drive wind power plant grid-connected system is characterized by comprising the following steps:
acquiring a state equation and an algebraic equation of a grid-connected system, and performing simultaneous combination by taking grid-side voltage as a carrier to obtain an equivalent grid-side voltage equation considering a phase-locked loop and a grid-side control link;
representing the phase-locked loop and the network side control link in the form of equivalent impedance and an equivalent voltage source to obtain an RLC equivalent circuit and corresponding equivalent impedance of a direct-drive wind power plant;
combining the equivalent impedance and the Nyquist stability criterion to obtain a stability margin expression of the grid-connected system;
obtaining the normalized sensitivity of different parameters in the grid-connected system to the stability margin of the grid-connected system according to the stability margin expression;
quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin, and screening out dominant influence factors;
the method comprises the following steps of obtaining a state equation and an algebraic equation of a grid-connected system, carrying out simultaneous combination by taking a network side voltage as a carrier, and obtaining an equivalent network side voltage equation considering a phase-locked loop and a network side control link, wherein the adopted expression comprises the following steps:
Figure FDA0003808107670000011
Figure FDA0003808107670000012
Figure FDA0003808107670000013
Figure FDA0003808107670000014
Figure FDA0003808107670000015
u g s* is the net side equivalent voltage; u. of g Is the net side voltage; i.e. i g A grid side current; i.e. i m Is the machine side current; u. of m Is the machine side voltage; s is a laplace operator; the superscript s represents the grid dq coordinate system; subscripts d and q denote d and q-axis components, respectively; subscript 0 represents the initial value of the corresponding variable; theta.theta. pll Is a phase lock angle; u shape g Is the net side voltage amplitude; l is a radical of an alcohol g The inductance of the incoming line reactor; u shape dc Is the converter DC side voltage; q gpu Is a reactive power per unit value; c dc Is a DC bus capacitor; k ppll Is the proportionality coefficient of the phase-locked loop; k ipll Is the integral coefficient of the phase-locked loop; k pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k ig The integral coefficient of the current inner loop of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; k is idc The integral coefficient of the outer ring of the direct current voltage is obtained; u shape gb Is a reference value of the network side voltage; I.C. A gb Is a reference value of the current on the network side; u shape dcb Is a reference value of the direct current voltage; s. the gb The reference value is the network side power;
the expression of the equivalent impedance is as follows:
Figure FDA0003808107670000021
the expression of the equivalent voltage source is as follows:
Figure FDA0003808107670000022
wherein R is l Is a line resistance; l is l Is a line inductance; l is T Is a transformer inductor; l is p The equivalent inductance of the direct-drive fan is obtained; n is the number of fans; k ig The integral coefficient of the current inner loop of the grid-side converter is obtained; u shape dc Is the converter DC side voltage; u shape gb Is a reference value of the network side voltage; i is gb Is a reference value of the current on the network side; c dc Is a direct current bus capacitor; c p Is equivalent to a direct-drive fanA capacitor; k pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; k idc Is the integral coefficient of the outer ring of the direct current voltage; omega is angular frequency; omega s Is the rotation angular velocity of the synchronous magnetic field; u shape dcb Is a reference value of the direct current voltage; s gb The reference value is the network side power; u shape dc0 Is the initial value of the direct current voltage; j is an imaginary unit;
the stability margin expression of the grid-connected system is as follows:
Figure FDA0003808107670000023
Figure FDA0003808107670000024
wherein, K g Representing an amplitude margin; γ represents a phase margin; r p Representing an equivalent impedance; r s Representing the system equivalent resistance; r is l Is a line resistance; l is l Is a line inductance; u shape gb Is a reference value of the network side voltage; k pg The proportionality coefficient of the current inner ring of the grid-side converter is obtained; k pdc Is the proportionality coefficient of the direct current voltage outer ring; i is gb Is a reference value of the current on the network side; s gb Is a reference value of the network side power; l is a radical of an alcohol p The equivalent inductance of the direct-drive fan is obtained; l is a radical of an alcohol s Is the equivalent inductance of the system; n is the number of fans;
the expression of the normalized sensitivity is as follows:
Figure FDA0003808107670000031
wherein a represents the parameters of the grid-connected system, m is the number of the parameters, and H is the amplitude margin or the phase margin.
2. The utility model provides a stability analysis device of direct-drive wind-powered electricity generation field grid-connected system which characterized in that includes:
the simultaneous combination module is used for acquiring a state equation and an algebraic equation of a grid-connected system, and performing simultaneous combination by taking the network side voltage as a carrier to obtain an equivalent network side voltage equation considering a phase-locked loop and a network side control link;
the characterization module is used for characterizing the phase-locked loop and the network side control link in the form of equivalent impedance and an equivalent voltage source to obtain an RLC equivalent circuit and corresponding equivalent impedance of the direct-drive wind power plant;
the stability margin acquisition module is used for combining the equivalent impedance and the Nyquist stability criterion to obtain a stability margin expression of the grid-connected system;
the sensitivity acquisition module is used for acquiring the normalized sensitivity of different parameters in the grid-connected system to the stability margin of the system according to the stability margin expression;
the influence degree screening module is used for quantifying the influence degree of different parameters on the normalized sensitivity of the system stability margin and screening out dominant influence factors;
wherein, the expression that the said simultaneous combination module adopts includes:
Figure FDA0003808107670000032
Figure FDA0003808107670000033
Figure FDA0003808107670000041
Figure FDA0003808107670000042
Figure FDA0003808107670000043
wherein u is g s* Is the net side equivalent voltage; u. u g Is the net side voltage; i all right angle g A grid side current; i.e. i m Is the machine side current; u. of m Is the machine side voltage; s is a laplace operator; the superscript s represents the grid dq coordinate system; subscripts d and q denote d and q-axis components, respectively; subscript 0 represents the initial value of the corresponding variable; theta pll Is a phase lock angle; u shape g Is the net side voltage amplitude; l is g The inductance of the incoming line reactor; u shape dc Is the converter DC side voltage; q gpu Is a per unit value of reactive power; c dc Is a DC bus capacitor; k is ppll Is the proportionality coefficient of the phase-locked loop; k ipll Is the integral coefficient of the phase-locked loop; k pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k ig The integral coefficient of the current inner loop of the grid-side converter is obtained; k pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; k idc The integral coefficient of the outer ring of the direct current voltage is obtained; u shape gb Is a reference value of the network side voltage; i is gb Is the reference value of the current at the network side; u shape dcb Is a reference value of the direct current voltage; s gb The reference value is the network side power;
in the characterization module, the expression of the equivalent impedance is as follows:
Figure FDA0003808107670000044
the expression of the equivalent voltage source is as follows:
Figure FDA0003808107670000045
wherein R is l Is a line resistance; l is l Is a line inductance; l is T Is a transformer inductor; l is p Is an equivalent inductance of the direct-drive fan; n is the number of fans; k ig The integral coefficient of the current inner loop of the grid-side converter is obtained; u shape dc Is the direct current side voltage of the converter; u shape gb Is a reference value of the network side voltage; i is gb Is a reference value of the current on the network side; c dc Is a DC busA capacitor; c p Is the equivalent capacitance of the direct-drive fan; k is pg The proportionality coefficient of the current inner ring of the grid-side converter is obtained; k idc The integral coefficient of the outer ring of the direct current voltage is obtained; omega is angular frequency; omega s Is the rotation angular velocity of the synchronous magnetic field; u shape dcb Is a reference value of the direct current voltage; s. the gb The reference value is the network side power; u shape dc0 Is the initial value of the direct current voltage; j is an imaginary unit;
in the stability margin obtaining module, a stability margin expression of the grid-connected system is as follows:
Figure FDA0003808107670000051
Figure FDA0003808107670000052
wherein, K g Representing an amplitude margin; γ represents a phase margin; r p Representing an equivalent impedance; r s Representing the equivalent resistance of the system; r l Is a line resistance; l is a radical of an alcohol l Is a line inductance; u shape gb Is a reference value of the network side voltage; k is pg The proportionality coefficient of a current inner ring of the grid-side converter is obtained; k is pdc The proportionality coefficient of the direct-current voltage outer ring is obtained; I.C. A gb Is a reference value of the current on the network side; s gb The reference value is the network side power; l is a radical of an alcohol p The equivalent inductance of the direct-drive fan is obtained; l is s Is the equivalent inductance of the system; n is the number of fans;
in the sensitivity obtaining module, the expression of the normalized sensitivity of the different parameters to the system stability margin is as follows:
Figure FDA0003808107670000053
wherein a represents the parameters of the grid-connected system, m is the number of the parameters, and H is the amplitude margin or the phase margin.
CN202110580826.XA 2021-05-26 2021-05-26 Stability analysis method and device for direct-drive wind power plant grid-connected system Active CN113241779B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110580826.XA CN113241779B (en) 2021-05-26 2021-05-26 Stability analysis method and device for direct-drive wind power plant grid-connected system

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110580826.XA CN113241779B (en) 2021-05-26 2021-05-26 Stability analysis method and device for direct-drive wind power plant grid-connected system

Publications (2)

Publication Number Publication Date
CN113241779A CN113241779A (en) 2021-08-10
CN113241779B true CN113241779B (en) 2022-10-28

Family

ID=77139015

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110580826.XA Active CN113241779B (en) 2021-05-26 2021-05-26 Stability analysis method and device for direct-drive wind power plant grid-connected system

Country Status (1)

Country Link
CN (1) CN113241779B (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113848720B (en) * 2021-10-09 2023-11-03 国核电力规划设计研究院有限公司 System stability analysis method based on full-power direct-drive fan control inertia

Family Cites Families (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
EP2124311A1 (en) * 2008-05-23 2009-11-25 ABB Research LTD Time delay compensation in power system control
WO2012118768A1 (en) * 2011-03-01 2012-09-07 Williams Michael Lamar Method and apparatus for controlling electrical loads to provide positive damping of power grid oscillation
CN112054549B (en) * 2020-09-08 2022-11-04 华北电力大学 Method for improving stability of direct-current power transmission system
CN112260264B (en) * 2020-09-16 2023-04-14 北京理工大学 Method and device for judging stability of land double-fed wind field through traditional direct current sending system

Also Published As

Publication number Publication date
CN113241779A (en) 2021-08-10

Similar Documents

Publication Publication Date Title
Li et al. Wind in weak grids: Low-frequency oscillations, subsynchronous oscillations, and torsional interactions
CN109449958B (en) A kind of double-fed blower grid-connected system method for analyzing stability
Nian et al. Stability and power quality enhancement strategy for DFIG system connected to harmonic grid with parallel compensation
CN109103903A (en) A kind of judgment method causing sub-synchronous oscillation for straight drive blower
CN109586337B (en) Frequency domain modeling-based VSC grid-connected system subsynchronous oscillation risk assessment method
Zhang et al. Impedance modeling and SSR analysis of DFIG using complex vector theory
CN109921421B (en) Method for establishing harmonic current output model of doubly-fed wind turbine generator
CN108880300B (en) Double-fed fan rectifier impedance calculation method based on double closed-loop control
CN112615393B (en) Vector fitting-based direct-drive wind turbine generator controller parameter identification method and device
CN109066725B (en) Direct-drive fan equivalent modeling method for subsynchronous oscillation simulation
CN107069800A (en) A kind of method for building up of double-fed fan motor short circuit current flow model
CN112217203B (en) Direct-drive wind power plant subsynchronous oscillation equivalence method capable of keeping model structure
CN110086207A (en) A kind of grid-connected converter Control method, apparatus of energy storage and computer storage medium
CN113241779B (en) Stability analysis method and device for direct-drive wind power plant grid-connected system
CN111697618A (en) Impedance modeling and stability analysis method for wide frequency domain full model of large-scale wind power plant
CN103490423B (en) Wind-farm reactive instruction distribution method
Lin et al. Vector-variable based modeling method for large-scale wind farm considering collector system dynamics
CN105552951B (en) A kind of DFIG system control methods based on repetition sliding formwork
CN109088436B (en) VSC multi-scale transient modeling method
CN111769574B (en) Direct-drive wind field sub/super frequency oscillation suppression method and device based on energy dissipation
CN110311409B (en) Improved double-ring DFIG low voltage ride through control strategy under unbalanced voltage
Gao et al. Improved extended kalman filter based dynamic equivalent method of DFIG wind farm cluster
CN113482852A (en) Control method, terminal and storage medium for permanent magnet direct-drive wind power generation converter
Liu et al. Research on power flow algorithm for power system including wind farm
Sun et al. General impedance model of DFIG for wide-range-frequency oscillation studies

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
GR01 Patent grant
GR01 Patent grant