CN116191478B - Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors - Google Patents
Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors Download PDFInfo
- Publication number
- CN116191478B CN116191478B CN202310439751.2A CN202310439751A CN116191478B CN 116191478 B CN116191478 B CN 116191478B CN 202310439751 A CN202310439751 A CN 202310439751A CN 116191478 B CN116191478 B CN 116191478B
- Authority
- CN
- China
- Prior art keywords
- power system
- equivalent model
- frequency response
- dynamic equivalent
- test
- 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
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
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F3/00—Input arrangements for transferring data to be processed into a form capable of being handled by the computer; Output arrangements for transferring data from processing unit to output unit, e.g. interface arrangements
- G06F3/01—Input arrangements or combined input and output arrangements for interaction between user and computer
- G06F3/02—Input arrangements using manually operated switches, e.g. using keyboards or dials
- G06F3/023—Arrangements for converting discrete items of information into a coded form, e.g. arrangements for interpreting keyboard generated codes as alphanumeric codes, operand codes or instruction codes
- G06F3/027—Arrangements for converting discrete items of information into a coded form, e.g. arrangements for interpreting keyboard generated codes as alphanumeric codes, operand codes or instruction codes for insertion of the decimal point
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06N—COMPUTING ARRANGEMENTS BASED ON SPECIFIC COMPUTATIONAL MODELS
- G06N3/00—Computing arrangements based on biological models
- G06N3/12—Computing arrangements based on biological models using genetic models
- G06N3/126—Evolutionary algorithms, e.g. genetic algorithms or genetic programming
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/04—Constraint-based CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/06—Multi-objective optimisation, e.g. Pareto optimisation using simulated annealing [SA], ant colony algorithms or genetic algorithms [GA]
-
- 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]
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Engineering & Computer Science (AREA)
- Biophysics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Evolutionary Biology (AREA)
- General Physics & Mathematics (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Biomedical Technology (AREA)
- Data Mining & Analysis (AREA)
- Physiology (AREA)
- Genetics & Genomics (AREA)
- Artificial Intelligence (AREA)
- Human Computer Interaction (AREA)
- Computational Linguistics (AREA)
- Power Engineering (AREA)
- Evolutionary Computation (AREA)
- General Health & Medical Sciences (AREA)
- Molecular Biology (AREA)
- Computing Systems (AREA)
- Mathematical Physics (AREA)
- Software Systems (AREA)
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention provides a method for evaluating equivalent inertia of a plurality of asynchronous motors and modeling frequency response, belonging to the technical field of frequency stability analysis of an electric power system; the problem that when a plurality of asynchronous motors are connected into a power system, inertia evaluation and frequency response modeling of the power system are difficult is solved; the method comprises the following steps: step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set; step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form; step three: fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm; step four: testing whether the parameters of the established dynamic equivalent model meet the test indexes or not on the test set; the invention is applied to stability analysis of an electric power system comprising a plurality of asynchronous motors.
Description
Technical Field
The invention provides a method for evaluating equivalent inertia of a plurality of asynchronous motors and modeling frequency response, and belongs to the technical field of frequency stability analysis of power systems.
Background
In recent years, more and more new energy power generation systems are connected into a power system through a power electronic converter, and the ratio of the photovoltaic power station to the wind turbine generator in the new energy is more than 80%. The photovoltaic power station does not have rotational inertia, and the power electronic converter makes the inertia of the wind turbine generator difficult to effectively utilize. Therefore, a large number of conventional generator sets are replaced by new energy sources, which causes the moment of inertia of the power system to gradually decrease, causing concerns about frequency stability.
Under the above circumstances, the asynchronous motor directly connected to the power grid has the capability of participating in the frequency response, and can improve primary frequency modulation performance, so that the asynchronous motor has an increasingly important role in the frequency response.
However, when evaluating the inertia of an asynchronous motor and establishing a mechanism model of the frequency response of an electric power system, detailed parameters of the asynchronous motor are often required, for example, an equivalent inertia evaluation method of the asynchronous motor in an inertia response phase disclosed in the patent application No. 2022105164271, but it is very difficult to implement the equivalent inertia evaluation method for a large system having hundreds of asynchronous motors.
Disclosure of Invention
The invention provides a method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors, which aims to solve the problem that the inertia evaluation and the frequency response modeling of a power system are difficult when the plurality of asynchronous motors are connected into the power system.
In order to solve the technical problems, the invention adopts the following technical scheme: a method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors comprises the following steps:
step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set;
step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form;
step three: fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm;
step four: and testing whether the parameters of the established dynamic equivalent model meet the test indexes or not on the test set.
The change rule of the frequency deviation of the dynamic equivalent model established in the second step is as follows:
in the above formula:Din order to be a damping coefficient,His an inertial time constant of the power system;
ΔP d for load disturbance, the expression of load disturbance is:
in the above formula:P d,mag the amplitude of the load disturbance is a step function;
ΔP g for the output variation of the synchronous generator, the expression of the output variation of the synchronous generator is as follows:
in the above formula:T R for the reheat time constant, the temperature of the refrigerant is,K m as a function of the mechanical power gain factor,F H for the high-pressure turbine coefficient,Ris a governor coefficient.
In the third step, parameters of the dynamic equivalent model are fitted from the test set through a genetic algorithm, wherein the parameters of the dynamic equivalent model comprise: speed regulator coefficientRInertial time constant of power systemHDamping coefficientDCoefficient of mechanical power gainK m High pressure turbine coefficientF H And reheat time constantT R 。
The expression of the optimization target of the genetic algorithm is as follows:
in the above formula:Jin order to adapt the function of the degree of adaptation,mas the number of sets of data to be used for fitting,n i for the number of samples of the power system frequency in each set of data,data for true frequency offset, +.>Is the frequency offset of the dynamic equivalent model.
Constraint conditions to be met by the optimization target of the genetic algorithm are as follows:
the expression of the test index in the fourth step is as follows:
in the above formula:J test for the test index, used to characterize the average fitting error,lfor the number of test sets, the standard for test passing is that the test index is less than 1%.
Compared with the prior art, the invention has the following beneficial effects: the method for evaluating the equivalent inertia of the plurality of asynchronous motors and modeling the frequency response aims to establish a dynamic equivalent model, adopts a genetic algorithm to fit parameters of the dynamic equivalent model from historical data, is convenient and quick, does not need detailed asynchronous motor parameters, and can be suitable for evaluating the inertia and modeling the frequency response of a large-scale power system comprising the plurality of asynchronous motors.
Drawings
The invention is further described below with reference to the accompanying drawings:
FIG. 1 is a flow chart of the method of the present invention;
FIG. 2 is a schematic diagram of an IEEE9 node system with three asynchronous motors in accordance with an embodiment of the present invention;
FIG. 3 is a schematic illustration of a Simulink modeling architecture of the IEEE9 node system with three asynchronous motors of FIG. 2;
FIG. 4 is a schematic diagram of a transfer function of a dynamic equivalent model of the frequency response of the power system according to the present invention;
FIG. 5 is a schematic diagram of the power system frequency response under three different sets of perturbations used for fitting in accordance with the present invention;
FIG. 6 is a graph showing the frequency response of the power system under four different sets of perturbations used in the test according to the present invention.
Detailed Description
The invention aims to provide a method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors driven by data, which comprises the following implementation steps of:
step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set;
step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form;
step three: fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm;
step four: and testing the established dynamic equivalent model parameters on a test set, wherein the standard of test passing is that the test index is less than 1%.
The method of the invention is further described below with reference to specific examples.
Embodiment one: the schematic diagram of the power system to be tested is shown in fig. 2, the Simulink model of the power system to be tested is shown in fig. 3, and the Simulink model is a frequency response historical data source of the power system to be tested.
Historical frequency response data of a power system comprising a plurality of asynchronous motors is obtained, and the historical data is divided into three fitting sets and four testing sets.
The dynamic equivalent model of the frequency response of the power system is shown in fig. 4, a dynamic equivalent model of the power system to be analyzed is established, and the dynamic equivalent model adopts a low-order power system model in a transfer function form. In this model, the law of variation of the frequency deviation is:
in the method, in the process of the invention,Dthe damping coefficient is generally 0 to 2,Hthe inertia time constant of the power system is generally 3-9 s.
ΔP d For load disturbance, the expression of load disturbance is:
in the method, in the process of the invention,P d,mag the amplitude of the load disturbance is a step function.
ΔP g For the output variation of the synchronous generator, the expression of the output variation of the synchronous generator is
In the method, in the process of the invention,T R for the reheat time constant, it is usually between 6 to 14s,K m the mechanical power gain factor is mainly influenced by the power factor and the rotation reserve, and is generally between 0.4 and 1,F H the high-pressure turbine coefficient is generally between 0.15 and 0.4,Rthe speed regulator coefficient is generally 0.04-1.
Fitting parameters of a dynamic equivalent model from the test set through a genetic algorithm, wherein the parameters of the dynamic equivalent model comprise: speed regulator coefficientRInertial time constant of power systemHDamping coefficientDCoefficient of mechanical power gainK m High pressure turbine coefficientF H And reheat time constantT R 。
The optimization targets of the genetic algorithm are:
in the method, in the process of the invention,Jin order to adapt the function of the degree of adaptation,mas the number of sets of data to be used for fitting,n i for the number of samples of the power system frequency in each set of data,data for true frequency offset, +.>Is the frequency offset of the dynamic equivalent model. The constraint conditions to be satisfied by the above optimization problem are:
i.e. all parameters are guaranteed to be within reasonable limits.
The genetic algorithm fits on three fitting sets, and the frequency response of the obtained equivalent dynamic model and the equivalent dynamic model pair of the power system to be analyzed are as shown in fig. 5.
And testing the established dynamic equivalent model parameters on a test set, wherein the test indexes are as follows:
in the method, in the process of the invention,J test for the test index, used to characterize the average fitting error,lfor the number of test sets, the standard for test passing is that the test index is less than 1%.
The equivalent dynamic model and the test set result of the power system to be analyzed are shown in fig. 6.
According to the invention, the inertia of the power system can be evaluated and a frequency response model can be established under the condition that the detailed parameters of the asynchronous motor are unknown, so that the frequency stability of the power system can be effectively analyzed.
The specific structure of the invention needs to be described that the connection relation between the component modules adopted by the invention is definite and realizable, and besides the specific description in the embodiment, the specific connection relation can bring about corresponding technical effects, and on the premise of not depending on execution of corresponding software programs, the technical problems of the invention are solved, the types of the components, the modules and the specific components, the connection modes of the components and the expected technical effects brought by the technical characteristics are clear, complete and realizable, and the conventional use method and the expected technical effects brought by the technical characteristics are all disclosed in patents, journal papers, technical manuals, technical dictionaries and textbooks which can be acquired by a person in the field before the application date, or the prior art such as conventional technology, common knowledge in the field, and the like, so that the provided technical scheme is clear, complete and the corresponding entity products can be reproduced or obtained according to the technical means.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and not for limiting the same; although the invention has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical scheme described in the foregoing embodiments can be modified or some or all of the technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit of the invention.
Claims (2)
1. A method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors is characterized by comprising the following steps: the method comprises the following steps:
step one: acquiring frequency response data of a power system to be analyzed through a power grid SCADA system, and dividing historical data into a fitting set and a testing set;
step two: establishing a dynamic equivalent model of the power system to be analyzed, wherein the dynamic equivalent model adopts a low-order power system model in a transfer function form;
step three: fitting parameters of the dynamic equivalent model from the fitting set through a genetic algorithm, wherein the parameters of the dynamic equivalent model comprise: speed regulator coefficientRInertial time constant of power systemHDamping coefficientDCoefficient of mechanical power gainK m High pressure turbine coefficientF H And reheat time constantT R ;
The expression of the optimization target of the genetic algorithm is as follows:
in the above formula:Jin order to adapt the function of the degree of adaptation,mas the number of sets of data to be used for fitting,n i for the number of samples of the power system frequency in each set of data,data for true frequency offset, +.>Frequency offset for dynamic equivalent model;
constraint conditions to be met by the optimization target of the genetic algorithm are as follows:
step four: testing whether the parameters of the established dynamic equivalent model meet the test indexes or not on the test set;
the expression of the test index in the fourth step is as follows:
in the above formula:J test for the test index, used to characterize the average fitting error,lfor the number of test sets, the standard for test passing is that the test index is less than 1%.
2. A method for evaluating equivalent inertia and modeling frequency response of a plurality of asynchronous motors according to claim 1, wherein: the change rule of the frequency deviation of the dynamic equivalent model established in the second step is as follows:
in the above formula:Din order to be a damping coefficient,His an inertial time constant of the power system;
in the above formula:P d,mag is the amplitude of the load disturbance;
ΔP g for the output variation of the synchronous generator, the expression of the output variation of the synchronous generator is as follows:
in the above formula:T R for the reheat time constant, the temperature of the refrigerant is,K m as a function of the mechanical power gain factor,F H for the high-pressure turbine coefficient,Ris a governor coefficient.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310439751.2A CN116191478B (en) | 2023-04-23 | 2023-04-23 | Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202310439751.2A CN116191478B (en) | 2023-04-23 | 2023-04-23 | Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors |
Publications (2)
Publication Number | Publication Date |
---|---|
CN116191478A CN116191478A (en) | 2023-05-30 |
CN116191478B true CN116191478B (en) | 2023-07-11 |
Family
ID=86452381
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202310439751.2A Active CN116191478B (en) | 2023-04-23 | 2023-04-23 | Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN116191478B (en) |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113972654A (en) * | 2021-11-01 | 2022-01-25 | 国网江苏省电力有限公司电力科学研究院 | Multipoint-access area load inertia time constant identification method |
CN115483707A (en) * | 2022-10-12 | 2022-12-16 | 东南大学溧阳研究院 | Novel power system frequency situation prediction method considering photovoltaic frequency modulation |
Family Cites Families (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111276973A (en) * | 2020-03-09 | 2020-06-12 | 国网江苏省电力有限公司 | Method for evaluating inertia requirement of power system considering wind power fluctuation |
CN114638074A (en) * | 2021-12-31 | 2022-06-17 | 国网辽宁省电力有限公司电力科学研究院 | Inertia evaluation method based on quantum derivation algorithm |
CN114792055A (en) * | 2022-05-05 | 2022-07-26 | 山西大学 | Asynchronous motor equivalent inertia evaluation method based on transient reactance post-potential |
CN114841005A (en) * | 2022-05-12 | 2022-08-02 | 南京理工大学 | Method for evaluating equivalent inertia of asynchronous motor in inertia response stage |
CN115622149A (en) * | 2022-09-06 | 2023-01-17 | 华中科技大学 | System frequency response modeling method and system for double-fed fan participating in primary frequency modulation |
-
2023
- 2023-04-23 CN CN202310439751.2A patent/CN116191478B/en active Active
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113972654A (en) * | 2021-11-01 | 2022-01-25 | 国网江苏省电力有限公司电力科学研究院 | Multipoint-access area load inertia time constant identification method |
CN115483707A (en) * | 2022-10-12 | 2022-12-16 | 东南大学溧阳研究院 | Novel power system frequency situation prediction method considering photovoltaic frequency modulation |
Also Published As
Publication number | Publication date |
---|---|
CN116191478A (en) | 2023-05-30 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN103887815B (en) | Based on wind energy turbine set parameter identification and the Dynamic Equivalence of service data | |
Nadour et al. | Comparative analysis between PI & backstepping control strategies of DFIG driven by wind turbine | |
Leithead et al. | Control of variable speed wind turbines: Design task | |
Zhong et al. | A novel frequency regulation strategy for a PV system based on the curtailment power-current curve tracking algorithm | |
CN110492479B (en) | Method for identifying rotational inertia and damping of distributed grid-connected equipment | |
Mensou et al. | An efficient nonlinear Backstepping controller approach of a wind power generation system based on a DFIG | |
CN112542855B (en) | Modeling and simulation method for phasor model of double-fed wind power generation system | |
CN110829487A (en) | Dynamic frequency prediction method for power system | |
CN111725840B (en) | Parameter identification method for direct-drive wind generating set controller | |
CN107947228B (en) | Stochastic stability analysis method for power system containing wind power based on Markov theory | |
CN111756038B (en) | New energy power system equal frequency difference inertia estimation method considering frequency modulation characteristics | |
Li et al. | Research on clustering equivalent modeling of large-scale photovoltaic power plants | |
CN111311021A (en) | Theoretical power prediction method, device, equipment and storage medium for wind power plant | |
CN116738636A (en) | Multi-machine equivalent method of doubly-fed wind power station considering impedance characteristics and synchronization mechanism | |
CN115578016A (en) | Online evaluation method for frequency modulation capability of wind power plant with incomplete model | |
Weijie et al. | Investigating instability of the wind turbine simulator with the conventional inertia emulation scheme | |
CN110829491A (en) | Grid-connected photovoltaic power generation system parameter identification method based on transient disturbance | |
Beniss et al. | Improvement of Power Quality Injected into the Grid by Using a FOSMC-DPC for Doubly Fed Induction Generator. | |
CN104617578A (en) | Method for acquiring available power transmission capability of power system with wind power plant | |
CN116191478B (en) | Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors | |
Yan et al. | Transient modelling of doubly‐fed induction generator based wind turbine on full operation condition and rapid starting period based on low voltage ride‐through testing | |
Jin et al. | Maximum power control of wind turbines with practical prescribed time stability based on wind estimation | |
CN109657380A (en) | A kind of double-fed fan motor field Dynamic Equivalence based on Extended Kalman filter | |
Bouguerra | Comparative study between PI, FLC, SMC and Fuzzy sliding mode controllers of DFIG wind turbine | |
Manjeera et al. | Design and Implementation of Fuzzy logic-2DOF controller for Emulation of wind turbine System |
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 |