CN113991704A - Method and device for identifying inertia and primary frequency modulation parameters based on synchronous generator set - Google Patents

Method and device for identifying inertia and primary frequency modulation parameters based on synchronous generator set Download PDF

Info

Publication number
CN113991704A
CN113991704A CN202111283166.5A CN202111283166A CN113991704A CN 113991704 A CN113991704 A CN 113991704A CN 202111283166 A CN202111283166 A CN 202111283166A CN 113991704 A CN113991704 A CN 113991704A
Authority
CN
China
Prior art keywords
primary frequency
inertia
frequency
synchronous motor
frequency modulation
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.)
Granted
Application number
CN202111283166.5A
Other languages
Chinese (zh)
Other versions
CN113991704B (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.)
State Grid Corp of China SGCC
State Grid Tianjin Electric Power Co Ltd
Original Assignee
State Grid Corp of China SGCC
State Grid Tianjin 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 Tianjin Electric Power Co Ltd filed Critical State Grid Corp of China SGCC
Priority to CN202111283166.5A priority Critical patent/CN113991704B/en
Publication of CN113991704A publication Critical patent/CN113991704A/en
Application granted granted Critical
Publication of CN113991704B publication Critical patent/CN113991704B/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
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/10—Power transmission or distribution systems management focussing at grid-level, e.g. load flow analysis, node profile computation, meshed network optimisation, active network management or spinning reserve management
    • 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
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • 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

Landscapes

  • Engineering & Computer Science (AREA)
  • Power Engineering (AREA)
  • Control Of Ac Motors In General (AREA)

Abstract

The invention provides a method and a device for identifying inertia and primary frequency modulation parameters based on a synchronous generator set, and relates to the technical field of frequency stability control of an electric power system, wherein the method comprises the following steps: acquiring a system power basic value, a voltage value and a current value of a PCC (point of charge control) measuring point, acquiring the electromagnetic power variation of an equivalent unit at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing an inertia estimation equation of the synchronous motor, and acquiring an inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor; the method comprises the steps of obtaining steady state deviation of the frequency of the equivalent unit and the electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing a target function, and obtaining a primary frequency modulation parameter by the primary frequency control coefficient through an iteration method based on an inertia estimated value of a synchronous motor. The method and the device provided by the invention can solve the problem that the control design of the new energy power electronic interface is difficult when the key parameters of the synchronous unit are unknown in the prior art.

Description

Method and device for identifying inertia and primary frequency modulation parameters based on synchronous generator set
Technical Field
The invention relates to the technical field of power system frequency stability control, in particular to a method and a device for identifying parameters based on inertia and primary frequency modulation of a synchronous generator set.
Background
In order to improve the effective inertia and the frequency stability of the dual-high system, the control mode and the control parameters of a new energy and energy storage power electronic interface need to be reasonably designed. In most of the existing researches, the control mode and parameter (such as virtual inertia and damping coefficient in virtual synchronous control) design are based on the known conditions of the frequency dynamic related parameters of the synchronous unit in the researched area. However, in practical engineering, the real inertia and the primary frequency control parameter of the synchronous unit may not be directly obtained. For renewable energy control, the synchronous set behaves as a "black box". Therefore, in the double-high system, when the frequency dynamic related parameters of the synchronous generator set are unknown, how to reasonably design a new energy interface is a main motivation of the invention. The main idea is to respectively identify the inertia of the unit and the primary frequency control parameter according to different stages of the synchronous unit after being disturbed by power.
The basic idea of synchronous machine inertia estimation is based on a synchronous machine swing equation, a variable-order polynomial frequency fitting algorithm is provided in the prior art, the precision of a traditional estimation method based on frequency fitting is improved, the prior art also identifies a state space model in a data driving mode through online inertia estimation, and inertia is estimated according to unit step response of the state space model. Aiming at a system comprising a synchronous machine and a virtual synchronous machine, the equivalent inertia of the system is estimated by adopting a typical second-order system fitting frequency dynamic method. Notably, most existing inertia estimation approaches utilize a Phase Measurement Unit (PMU) to obtain the frequency signal, which increases the cost of parameter identification. In addition, the identification of the primary frequency modulation parameter is not considered in the research work, and the identification also has great influence on the control design of the new energy interface.
Disclosure of Invention
In view of the above, the present invention provides a method and an apparatus for identifying an inertia and a primary frequency modulation parameter of a synchronous generator set, so as to solve the problem in the prior art that a control design of a new energy power electronic interface is difficult when a key parameter of the synchronous generator set is unknown.
The invention provides a method for identifying inertia and primary frequency modulation parameters based on a synchronous generator set, which comprises the following steps:
acquiring a system power basic value, a voltage value and a current value of a PCC (point of charge control) measuring point, acquiring the electromagnetic power variation of an equivalent unit at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing an inertia estimation equation of a synchronous motor, and acquiring an inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor;
the method comprises the steps of obtaining steady state deviation of the frequency of the equivalent unit and the electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining a primary frequency modulation parameter by the primary frequency control coefficient through an iteration method based on the inertia estimated value of the synchronous motor.
Preferably, the synchronous motor inertia estimation equation is as follows:
Figure BDA0003331993680000021
Δ P (0+) -the equivalent unit electromagnetic power variation at the disturbance moment;
Sbase-a system power base value;
a1-the rate of change of frequency at the moment of disturbance.
Preferably, the step of obtaining the steady state deviation of the equivalent unit frequency and the electromagnetic power variation of the equivalent unit in the steady state to obtain the primary frequency control coefficient, constructing the objective function, and obtaining the primary frequency modulation parameter by using the iterative method for the primary frequency control coefficient based on the inertia estimated value of the synchronous motor includes:
the primary frequency control coefficient is obtained by adopting the following formula:
Figure BDA0003331993680000031
Δωsg(∞) -steady state deviation of the equivalent unit frequency;
Δ P (∞) -the variation of the electromagnetic power of the equivalent unit in a steady state;
the objective function is:
Figure BDA0003331993680000032
omega-experimental frequency curve
Figure BDA0003331993680000033
-estimating a model frequency curve
ωi-data points of the experimental frequency curve;
Figure BDA0003331993680000034
-estimating data points of the model frequency;
m-total number of frequency data
TGTime constant of speed regulator
FHPHigh pressure turbine coefficient
TRH-reheat time constant
TCHTime constant of the gas cell
Preferably, the step of obtaining the system power basic value and the disturbance power, constructing an inertia estimation equation of the synchronous motor, and obtaining the inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor includes:
the frequency change rate at the moment of disturbance is obtained by the following steps:
obtaining the frequency of the PCC measurement point after disturbance;
obtaining an estimated frequency change rate RoCoF based on the frequency after the PCC measurement point disturbance, and judging whether the frequency change rate RoCoF exceeds a first threshold value;
if so, acquiring the frequency change rate at the moment of disturbance by using a polynomial fitting method;
and if not, acquiring an estimated frequency change rate RoCoF based on the frequency after the PCC measurement point disturbance, and judging whether the frequency change rate RoCoF exceeds a threshold value.
Preferably, the step of constructing an objective function, and obtaining a primary frequency modulation parameter by using a primary frequency control coefficient using an iterative method based on the inertia estimated value of the synchronous motor, includes:
setting primary frequency modulation parameters, wherein the sequential frequency modulation parameters comprise TG、FHP、TRH、TCH；
Obtaining an estimated model frequency based on the primary frequency modulation parameter
Figure BDA0003331993680000041
And a primary frequency modulation parameter step length;
judging whether the step length of the primary frequency modulation parameter is smaller than a second threshold value;
if so, outputting the primary frequency modulation parameter;
if not, based on the primary frequency modulation parameter and the primary frequency modulationUpdating the primary frequency modulation parameter by the parameter step length, and re-executing the frequency of the estimation model obtained based on the primary frequency modulation parameter
Figure BDA0003331993680000042
And step of primary frequency modulation parameter step length.
On the other hand, the invention provides a device for identifying inertia and primary frequency modulation parameters based on a synchronous generator set, preferably, the device comprises:
an inertia estimation module: the system comprises a synchronous motor, a PCC measuring point, a synchronous motor inertia estimation equation and a control system, wherein the synchronous motor is used for acquiring a system power basic value, a voltage value and a current value of the PCC measuring point, acquiring the equivalent unit electromagnetic power variation at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing the synchronous motor inertia estimation equation and acquiring an inertia estimation value of the synchronous motor based on the synchronous motor inertia estimation equation;
a primary frequency modulation parameter acquisition module: the method is used for obtaining steady state deviation of the equivalent unit frequency and electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining a primary frequency modulation parameter by using the primary frequency control coefficient through an iteration method based on the inertia estimated value of the synchronous motor.
The embodiment of the invention has the following beneficial effects: the invention provides a method and a device for identifying inertia and primary frequency modulation parameters based on a synchronous generator set, wherein the method comprises the following steps: acquiring a system power basic value, a voltage value and a current value of a PCC (point of charge control) measuring point, acquiring the electromagnetic power variation of an equivalent unit at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing an inertia estimation equation of the synchronous motor, and acquiring an inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor; the method comprises the steps of obtaining steady state deviation of the frequency of the equivalent unit and the electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing a target function, and obtaining a primary frequency modulation parameter by the primary frequency control coefficient through an iteration method based on an inertia estimated value of a synchronous motor. The method and the device provided by the invention can solve the problem that the control design of the new energy power electronic interface is difficult when the key parameters of the synchronous unit are unknown in the prior art.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are 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 flowchart of a method for identifying inertia and primary frequency modulation parameters based on a synchronous generator set according to an embodiment of the present invention;
fig. 2 is a flow chart for acquiring a frequency change rate at a disturbance moment based on a synchronous generator set inertia and primary frequency modulation parameter identification method according to an embodiment of the present invention;
fig. 3 is a flow chart of primary frequency modulation parameter acquisition based on a synchronous generator set inertia and primary frequency modulation parameter identification method provided by the embodiment of the invention;
FIG. 4 is a block diagram of a system based on a synchronous generator set according to an embodiment of the present invention;
fig. 5(a) is a sample frequency curve calculated based on the method for identifying the inertia and the primary frequency modulation parameters of the synchronous generator set according to the embodiment of the present invention;
fig. 5(b) is a fitting curve of 0 to 4 seconds of example disturbance based on the identification method of the inertia and the primary frequency modulation parameters of the synchronous generator set according to the embodiment of the present invention;
fig. 5(c) is a primary frequency modulation parameter process based on the synchronous generator set inertia and primary frequency modulation parameter identification method provided by the embodiment of the present invention;
fig. 5(d) is a frequency dynamic identification result diagram based on the synchronous generator set inertia and primary frequency modulation parameter identification method provided by the embodiment of the invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. 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.
At present, the basic idea of synchronous machine inertia estimation is based on a synchronous machine oscillation equation, most of existing inertia estimation means are that a Phase Measurement Unit (PMU) is used for obtaining a frequency signal, which increases the cost of parameter identification, and identification of primary frequency modulation parameters is not considered in the prior art, which also has a great influence on the control design of a new energy interface.
To facilitate understanding of the embodiment, a method for identifying parameters based on inertia and primary frequency modulation of a synchronous generator set disclosed in the embodiment of the present invention is first described in detail.
The first embodiment is as follows:
with reference to fig. 1, the first embodiment provides a method for identifying inertia and primary frequency modulation parameters based on a synchronous generator set, including:
acquiring a system power basic value, a voltage value and a current value of a PCC (point of charge control) measuring point, acquiring the electromagnetic power variation of an equivalent unit at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing an inertia estimation equation of a synchronous motor, and acquiring an inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor;
the method comprises the steps of obtaining steady state deviation of the frequency of the equivalent unit and the electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining a primary frequency modulation parameter by the primary frequency control coefficient through an iteration method based on the inertia estimated value of the synchronous motor.
Preferably, the synchronous motor inertia estimation equation is as follows:
Figure BDA0003331993680000071
Δ P (0+) -the equivalent unit electromagnetic power variation at the disturbance moment;
Sbase-a system power base value;
a1-the rate of change of frequency at the moment of disturbance.
Further, in the embodiment provided by the present invention, when unbalanced power is applied (t ═ 0+), it is considered that the mechanical power remains unchanged and the electromagnetic power disturbance Δ P is known at this time;
preferably, the step of obtaining the steady state deviation of the equivalent unit frequency and the electromagnetic power variation of the equivalent unit in the steady state to obtain the primary frequency control coefficient, constructing the objective function, and obtaining the primary frequency modulation parameter by using the iterative method for the primary frequency control coefficient based on the inertia estimated value of the synchronous motor includes:
the primary frequency control coefficient is obtained by adopting the following formula:
Figure BDA0003331993680000072
Δωsg(∞) -steady state deviation of the equivalent unit frequency;
Δ P (∞) -the variation of the electromagnetic power of the equivalent unit in a steady state;
further, the steady state deviation delta omega of the equivalent unit frequencysgThe (∞) and the steady-state equivalent unit electromagnetic power variation delta P (∞) pass through the current of the PCC measuring point in the steady stateObtaining a value and a voltage value;
in the embodiment provided by the invention, the steady state deviation of the equivalent unit frequency and the electromagnetic power variation of the equivalent unit in the steady state can be obtained through the voltage and the current of the PCC measuring point;
the objective function is:
Figure BDA0003331993680000081
omega-experimental frequency curve
Figure BDA0003331993680000082
-estimating a model frequency curve
ωi-data points of the experimental frequency curve;
Figure BDA0003331993680000083
-estimating data points of the model frequency;
m-total number of frequency data
TGTime constant of speed regulator
FHPHigh pressure turbine coefficient
TRH-reheat time constant
TCHTime constant of the gas cell
Further, in the embodiment provided by the present invention, the expression of the transfer function g(s) is:
Figure BDA0003331993680000084
the objective function can be obtained based on the transfer function and the estimated value of the inertia of the synchronous motor;
preferably, the step of obtaining the system power basic value and the disturbance power, constructing an inertia estimation equation of the synchronous motor, and obtaining the inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor includes:
with reference to fig. 2, the following steps are adopted to obtain the frequency change rate at the moment of disturbance:
obtaining the frequency of the PCC measurement point after disturbance;
obtaining an estimated frequency change rate RoCoF based on the frequency after the PCC measurement point disturbance, and judging whether the frequency change rate RoCoF exceeds a first threshold value;
if so, acquiring the frequency change rate at the moment of disturbance by using a polynomial fitting method;
and if not, acquiring an estimated frequency change rate RoCoF based on the frequency after the PCC measurement point disturbance, and judging whether the frequency change rate RoCoF exceeds a threshold value.
In an embodiment of the present invention, the frequency change rate rocef is d Δ ω at a time t ═ 0+sg/dt；
For a conventional synchronous machine, the inertia estimation is mainly based on the roll equation:
Figure BDA0003331993680000091
where H is the synchronous unit inertia, ωsgIs the per unit value, P, of the rotor frequency of the synchronous machinem、PeNamed values, S, of prime mover input mechanical power and unit electromagnetic power, respectivelybaseIs a system power base value;
wherein, the unit frequency can be expressed as:
ωsg＝antn+an-1tn-1+...+a1t+a0；
when the rotor frequency of the synchronous motor does not change suddenly, the RoCoF is approximate to a1However, since the synchronous motor suddenly changes, the frequency change rate a at the moment of disturbance occurs1Is not RoCoF;
with reference to fig. 3, preferably, the step of constructing an objective function, and obtaining a primary frequency modulation parameter by using a primary frequency control coefficient using an iterative method based on the inertia estimated value of the synchronous motor includes:
is provided with aA secondary frequency modulation parameter, the sequential frequency modulation parameter comprises TG、FHP、TRH、TCH；
Obtaining an estimated model frequency based on the primary frequency modulation parameter
Figure BDA0003331993680000092
And a primary frequency modulation parameter step length;
judging whether the step length of the primary frequency modulation parameter is smaller than a second threshold or exceeds a preset range;
if so, outputting the primary frequency modulation parameter;
if not, updating the primary frequency modulation parameter based on the primary frequency modulation parameter and the primary frequency modulation parameter step length, and re-executing the frequency of the estimation model obtained based on the primary frequency modulation parameter
Figure BDA0003331993680000101
And step of primary frequency modulation parameter step length.
Example two:
the embodiment of the invention provides a synchronous generator set inertia and primary frequency modulation parameter identification device, which comprises:
an inertia estimation module: the system comprises a synchronous motor, a PCC measuring point, a synchronous motor inertia estimation equation and a control system, wherein the synchronous motor is used for acquiring a system power basic value, a voltage value and a current value of the PCC measuring point, acquiring the equivalent unit electromagnetic power variation at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing the synchronous motor inertia estimation equation and acquiring an inertia estimation value of the synchronous motor based on the synchronous motor inertia estimation equation;
a primary frequency modulation parameter acquisition module: the method is used for obtaining steady state deviation of the equivalent unit frequency and electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining a primary frequency modulation parameter by using the primary frequency control coefficient through an iteration method based on the inertia estimated value of the synchronous motor.
Example three:
the third embodiment of the invention is a specific embodiment provided by the invention:
in order to verify the accuracy of the inertia and the identification method of the synchronous unit, an IEEE 5 machine 14 node system simulation model is built, power disturbance is applied to a #13 node, and the parameters of the synchronous unit in the area are shown in the table 1:
TABLE 1 IEEE 5 machine 14 node synchronizer parameters
Capacity (MW) Inertia moment 1/R TG TCH TRH FHP
Unit 1 300 6.5 13.33 0.2 0.37 10.5 0.28
Unit 2 200 6.5 20 0.22 0.36 12 0.39
Unit 3 200 6.175 20 0.27 0.41 6 0.23
Unit 4 100 6.175 16.67 0.3 0.48 14 0.32
Unit 5 100 6.175 16.67 0.3 0.48 14 0.32
1) Inertia estimation method verification
The simulation working condition is as follows: five machines are operated with load, and a rated disturbance power of 300MW is applied at a node # 13. The node frequency is obtained by phase locking the voltage data of the #13 node, polynomial fitting is carried out on a 4s frequency curve after disturbance occurs, and the least square method is used for fitting the 4s frequency result after power disturbance to show that: the coefficient a1 of the first order of the fitted curve is 0.0165; with reference to fig. 5(a) and 5(b), according to the load power curve, the equivalent unit electromagnetic power variation Δ P (0+) at the moment of disturbance is 180MW, the system power base value is 900MW, and the synchronous motor inertia estimation result is:
Figure BDA0003331993680000111
compared with the true value 6.3556, the relative error of the estimated value is only 4.6%.
Continuously using the inertia estimated value, carrying out other frequency modulation parameter estimation of the synchronous unit under the same simulation working condition, and selecting an initial value as TG＝0.2、TCH＝0.3、TRH＝0.3、FHPThe parameter identification process and results are shown in fig. 7. The final governor estimate is: t isG＝0.5、TCH＝0.5、TRH＝9.0417、FHP is 0.3365. As can be seen from fig. 5(c) and 5(d), even if the initial value deviates greatly from the true value, the convergence is very fast according to the proposed estimation method. Equivalent unit parameters of the identification area are shown in table 2.
TABLE 2 equivalent unit identification parameters
Capacity (MW) Inertia moment 1/R TG TCH TRH FHP
Equivalent unit 900 6.06 11.75 0.5 0.5 9.0417 0.3365
The invention has the following advantages:
1) the identification of the inertia of the equivalent synchronous unit and the primary frequency control parameter can be carried out through the dynamic process of frequency multi-stage recovery after the analysis system is disturbed when the number of synchronous units in the region is too large or the unit parameters are difficult to obtain, so that a frequency analysis model of the regional equivalent synchronous unit is established.
2) The voltage and power of the PCC point can be directly measured in situ, and then the frequency can be obtained through a phase-locked loop (PLL) technology. By adopting the estimation method without PMU, the cost of parameter identification can be greatly reduced.
Unless specifically stated otherwise, the relative steps, numerical expressions, and values of the components and steps set forth in these embodiments do not limit the scope of the present invention.
The flowchart and block diagrams in the figures illustrate the architecture, functionality, and operation of possible implementations of systems, methods and computer program products according to various embodiments of the present invention. In this regard, each block in the flowchart or block diagrams may represent a module, segment, or portion of code, which comprises one or more executable instructions for implementing the specified logical function(s). It should also be noted that, in some alternative implementations, the functions noted in the block may occur out of the order noted in the figures. For example, two blocks shown in succession may, in fact, be executed substantially concurrently, or the blocks may sometimes be executed in the reverse order, depending upon the functionality involved. It will also be noted that each block of the block diagrams and/or flowchart illustration, and combinations of blocks in the block diagrams and/or flowchart illustration, can be implemented by special purpose hardware-based systems which perform the specified functions or acts, or combinations of special purpose hardware and computer instructions.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In addition, in the description of the embodiments of the present invention, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (6)

1. A method for identifying parameters based on inertia and primary frequency modulation of a synchronous generator set is characterized by comprising the following steps:
acquiring a system power basic value, a voltage value and a current value of a PCC (point of charge control) measuring point, acquiring the electromagnetic power variation of an equivalent unit at the moment of disturbance based on the voltage value and the current value of the PCC measuring point, constructing an inertia estimation equation of a synchronous motor, and acquiring an inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor;
the method comprises the steps of obtaining steady state deviation of the frequency of the equivalent unit and the electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining a primary frequency modulation parameter by the primary frequency control coefficient through an iteration method based on the inertia estimated value of the synchronous motor.
2. The method of claim 1,
the synchronous motor inertia estimation equation is as follows:
Figure FDA0003331993670000011
Hest-a synchronous machine inertia estimate;
Δ P (0+) -the equivalent unit electromagnetic power variation at the disturbance moment;
Sbase-a system power base value;
a1-the rate of change of frequency at the moment of disturbance.
3. The method according to claim 1, wherein the step of obtaining steady state deviation of the equivalent unit frequency and the variation of the equivalent unit electromagnetic power in the steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining the primary frequency control coefficient by using an iterative method to obtain a primary frequency modulation parameter based on the inertia estimated value of the synchronous motor comprises the steps of:
the primary frequency control coefficient is obtained by adopting the following formula:
Figure FDA0003331993670000021
Δωsg(∞) -steady state deviation of the equivalent unit frequency;
Δ P (∞) -the variation of the electromagnetic power of the equivalent unit in a steady state;
the objective function is:
Figure FDA0003331993670000022
omega-experimental frequency curve
Figure FDA0003331993670000023
-estimating a model frequency curve
ωi-data points of the experimental frequency curve;
Figure FDA0003331993670000024
-estimating data points of the model frequency;
m-total number of frequency data
TGTime constant of speed regulator
FHPHigh pressure turbine coefficient
TRH-reheat time constant
TCH-the time constant of the gas cell.
4. The method of claim 1, wherein the steps of obtaining a system power basic value and a disturbance power, constructing an inertia estimation equation of the synchronous motor, and obtaining an inertia estimation value of the synchronous motor based on the inertia estimation equation of the synchronous motor comprise:
the frequency change rate at the moment of disturbance is obtained by the following steps:
obtaining the frequency of the PCC measurement point after disturbance;
obtaining an estimated frequency change rate RoCoF based on the frequency after the PCC measurement point disturbance, and judging whether the frequency change rate RoCoF exceeds a first threshold value;
if so, acquiring the frequency change rate at the moment of disturbance by using a polynomial fitting method;
and if not, acquiring an estimated frequency change rate RoCoF based on the frequency after the PCC measurement point disturbance, and judging whether the frequency change rate RoCoF exceeds a threshold value.
5. The method of claim 1, wherein the step of constructing an objective function, based on the estimated value of inertia of the synchronous machine, and using an iterative method for the primary frequency control coefficients to obtain the primary frequency modulation parameters comprises:
setting primary frequency modulation parameters, wherein the sequential frequency modulation parameters comprise TG、FHP、TRH、TCH；
Obtaining an estimated model frequency based on the primary frequency modulation parameterRate of change
Figure FDA0003331993670000031
And a primary frequency modulation parameter step length;
judging whether the step length of the primary frequency modulation parameter is smaller than a second threshold value;
if so, outputting the primary frequency modulation parameter;
if not, updating the primary frequency modulation parameter based on the primary frequency modulation parameter and the primary frequency modulation parameter step length, and re-executing the frequency of the estimation model obtained based on the primary frequency modulation parameter
Figure FDA0003331993670000032
And step of primary frequency modulation parameter step length.
6. The method for identifying the inertia and primary frequency modulation parameters based on the synchronous generator set is characterized by comprising the following steps of:
an inertia estimation module: the synchronous motor inertia estimation method is used for acquiring a system power basic value and disturbance power, constructing a synchronous motor inertia estimation equation and acquiring an inertia estimation value of a synchronous motor based on the synchronous motor inertia estimation equation;
a primary frequency modulation parameter acquisition module: the method is used for obtaining steady state deviation of the equivalent unit frequency and electromagnetic power variation of the equivalent unit in a steady state to obtain a primary frequency control coefficient, constructing an objective function, and obtaining a primary frequency modulation parameter by using the primary frequency control coefficient through an iteration method based on the inertia estimated value of the synchronous motor.
CN202111283166.5A 2021-11-01 2021-11-01 Synchronous generator set inertia and primary frequency modulation parameter identification method and device Active CN113991704B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111283166.5A CN113991704B (en) 2021-11-01 2021-11-01 Synchronous generator set inertia and primary frequency modulation parameter identification method and device

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111283166.5A CN113991704B (en) 2021-11-01 2021-11-01 Synchronous generator set inertia and primary frequency modulation parameter identification method and device

Publications (2)

Publication Number Publication Date
CN113991704A true CN113991704A (en) 2022-01-28
CN113991704B CN113991704B (en) 2024-07-02

Family

ID=79745355

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111283166.5A Active CN113991704B (en) 2021-11-01 2021-11-01 Synchronous generator set inertia and primary frequency modulation parameter identification method and device

Country Status (1)

Country Link
CN (1) CN113991704B (en)

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2339144C1 (en) * 2007-07-19 2008-11-20 Юрий Владимирович Шаров Method for improving dynamic stability and damping oscillations of electric-power systems and device for its realisation
CN110323763A (en) * 2019-07-18 2019-10-11 河海大学 A kind of inverter frequency modulation method of comprehensive rotatory inertia simulation and primary frequency modulation control
CN111444596A (en) * 2020-03-12 2020-07-24 国家电网公司东北分部 Method and device for detecting frequency safety after power disturbance of power system
CN111756039A (en) * 2020-07-09 2020-10-09 西安交通大学 New energy power system inertia estimation method based on probability statistics
CN112072674A (en) * 2020-07-21 2020-12-11 中国电力科学研究院有限公司 Method and system for determining inertia of power system
CN112290590A (en) * 2020-10-20 2021-01-29 国网四川省电力公司电力科学研究院 PMSG fan control method based on super-capacitor virtual inertia and short-term primary frequency modulation
CN112736988A (en) * 2020-12-24 2021-04-30 西安交通大学 New energy power system inertia estimation method considering influence of primary frequency modulation and droop control
CN113469332A (en) * 2021-06-28 2021-10-01 上海电机学院 Virtual synchronous generator inertia damping self-adaptive control method based on fuzzy nerves

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
RU2339144C1 (en) * 2007-07-19 2008-11-20 Юрий Владимирович Шаров Method for improving dynamic stability and damping oscillations of electric-power systems and device for its realisation
CN110323763A (en) * 2019-07-18 2019-10-11 河海大学 A kind of inverter frequency modulation method of comprehensive rotatory inertia simulation and primary frequency modulation control
CN111444596A (en) * 2020-03-12 2020-07-24 国家电网公司东北分部 Method and device for detecting frequency safety after power disturbance of power system
CN111756039A (en) * 2020-07-09 2020-10-09 西安交通大学 New energy power system inertia estimation method based on probability statistics
CN112072674A (en) * 2020-07-21 2020-12-11 中国电力科学研究院有限公司 Method and system for determining inertia of power system
CN112290590A (en) * 2020-10-20 2021-01-29 国网四川省电力公司电力科学研究院 PMSG fan control method based on super-capacitor virtual inertia and short-term primary frequency modulation
CN112736988A (en) * 2020-12-24 2021-04-30 西安交通大学 New energy power system inertia estimation method considering influence of primary frequency modulation and droop control
CN113469332A (en) * 2021-06-28 2021-10-01 上海电机学院 Virtual synchronous generator inertia damping self-adaptive control method based on fuzzy nerves

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
刘方蕾;毕天姝;闫家铭;王凡;胥国毅;: "基于PMU同步测量的分区惯量估计方法", 华北电力大学学报(自然科学版), no. 03, 30 May 2020 (2020-05-30) *

Also Published As

Publication number Publication date
CN113991704B (en) 2024-07-02

Similar Documents

Publication Publication Date Title
CN103473478B (en) Power Network Transient Stability appraisal procedure based on energy function
CN101119095A (en) Large destabilization real-time simulation system based on nonlinear robust power system stabilizer
CN109861296B (en) Inertia and damping identification method for island virtual synchronous generator system
CN106972484A (en) Fired power generating unit excitation system dynamic property online evaluation method and apparatus under many scenes
CN112511056A (en) Robust generator dynamic state estimation method based on phasor measurement
Davarikia et al. Real-time integrity indices in power grid: a synchronization coefficient based clustering approach
CN105224728B (en) A kind of Power Network Transient Stability energy function analysis method and system containing detailed generator model
Karrari et al. Identification of Heffron-Phillips model parameters for synchronous generators using online measurements
CN108574290B (en) Method, device, terminal and readable storage medium for positioning oscillation source of forced oscillation
Rigatos et al. Fault diagnosis of gas-turbine power units with the Derivative-free nonlinear Kalman Filter
CN113991704B (en) Synchronous generator set inertia and primary frequency modulation parameter identification method and device
Verij Kazemi et al. Adaptive frequency control support of a DFIG based on second‐order derivative controller using data‐driven method
CN108536116A (en) A kind of test method and system of governing system
Barrios-Gomez et al. RoCoF calculation using low-cost hardware in the loop: Multi-area Nordic power system
CN111293693A (en) Doubly-fed wind turbine converter control parameter identification method based on extended Kalman filtering
CN114069705B (en) Method and system for judging large disturbance stability of grid-connected voltage source type converter
CN112910006B (en) Universal electromagnetic transient modeling method for direct-drive wind turbine generator
Zerdali Strong tracking extended Kalman filter based speed and load torque estimations of induction motor
Arastou et al. Dynamic state and parameter estimation of the improved Heffron-Phillips model using a fast UKF-based algorithm and a novel rotor angle measurement approach
CN108512470B (en) Synchronous generator switched energization control method based on switch control
Wang et al. SDAE-based probabilistic stability analysis of wind integrated power systems
JP4757782B2 (en) Method and program for calculating constant of load model of subordinate system
Mao et al. State-of-charge estimation of lithium-ion polymer battery based on sliding mode observer
CN117895529A (en) Frequency response analysis method based on wind turbine generator equivalent inertia characterization model
Zuccaro et al. Identification of linear data-driven models for large-scale power systems

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