CN114243805B - Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting - Google Patents
Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting Download PDFInfo
- Publication number
- CN114243805B CN114243805B CN202111593772.7A CN202111593772A CN114243805B CN 114243805 B CN114243805 B CN 114243805B CN 202111593772 A CN202111593772 A CN 202111593772A CN 114243805 B CN114243805 B CN 114243805B
- Authority
- CN
- China
- Prior art keywords
- synchronous machine
- frequency response
- inertia
- power
- amplitude limiting
- 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/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
- H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
- H02J3/48—Controlling the sharing of the in-phase component
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- 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/40—Synchronising a generator for connection to a network or to another generator
-
- 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)
- Mathematical Physics (AREA)
- Theoretical Computer Science (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Power Engineering (AREA)
- Pure & Applied Mathematics (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Mathematical Optimization (AREA)
- Mathematical Analysis (AREA)
- Control Of Eletrric Generators (AREA)
Abstract
The invention discloses a synchronous machine system frequency response analysis calculation method considering amplitude limiting of a speed regulator, which comprises the following steps: acquiring all frequency modulation parameters of the whole network synchronous machine, including spare capacity, inertia, difference adjustment coefficient, reheating time constant and the like; the amplitude limiting time of full generation of mechanical power is calculated through the speed regulator parameters of the synchronous machine, the frequency response process of the system is segmented according to the amplitude limiting time, the synchronous machines of all the segments participating in the frequency response are aggregated into a single-machine equivalent model according to capacity, the non-zero initial state of the frequency response process of all the segments is calculated, and the transfer function of the inertia center frequency is subjected to inverse Laplace transform through a non-zero initial state analysis method to obtain a frequency response segmented time domain analysis model of the whole network inertia center. The invention can analyze and calculate the frequency response of the inertia center of the whole network more accurately on the premise of considering the amplitude limiting of the speed regulator, and simultaneously, the calculation effect of the invention is further close to the real power network, thereby better making day-ahead scheduling plan service for the power network.
Description
Technical Field
The invention belongs to the field of frequency modulation of power systems, and particularly relates to a synchronous machine system frequency response analysis calculation method considering amplitude limiting of a speed regulator.
Background
With the gradual exploitation of fossil energy, the demand for energy is gradually increased, and the relationship between human beings and the nature is increasingly tense. Under the background, constructing a novel power system mainly based on new energy, and considering new energy grid connection and related technologies thereof become an important direction for future power technology development.
The new energy occupation ratio of China increases year by year, but due to the characteristics of the fluctuation and weak inertia of the new energy, the new energy has poor supporting effect on the system frequency after the regional power grid generates load disturbance, and therefore the synchronous machine still plays a main role in frequency modulation.
Although the frequency response of the synchronous machine system is studied more deeply at present, the influence of reserve capacity cannot be considered by the mainstream frequency response calculation method at present, the modeling of a nonlinear amplitude limiting link of a speed regulator is still to be perfected, and how to better establish day-ahead scheduling plan service with strong disturbance rejection capability for a power grid is still to be studied.
Object of the Invention
The invention aims to solve the problems and provides a synchronous machine system frequency response analysis and calculation method considering the amplitude limit of a speed regulator, which can accurately analyze and calculate the frequency response of the inertia center of the whole network on the premise of considering the amplitude limit of the speed regulator and simultaneously ensure that the calculation effect of the frequency response is further close to the real power grid, thereby better establishing day-ahead scheduling plan service for the power grid.
Disclosure of Invention
The invention provides a synchronous machine system frequency response analysis and calculation method considering speed regulator amplitude limiting, which comprises the following steps:
step 1, acquiring key frequency modulation parameters of a whole network synchronous machine, including spare capacity P mg Inertia piece, difference adjustment coefficient R and reheating time constant T R High pressure turbine power fraction F H Damping coefficient xi and mechanical gain coefficient K m Disturbance power P Step ;
Step 2, calculating a transfer function of the output mechanical power of each synchronizer according to the key frequency modulation parameters of each synchronizer in the whole-network synchronizers obtained in the step 1, and judging the amplitude limiting time of each synchronizer;
and step 3: according to the amplitude limiting time of every synchronous machine making system frequency response processSegmenting each segment of the synchronous machine participating in frequency response, and setting a mechanical gain coefficient K according to the installed capacity and the total system capacity of the synchronous machine m Using the mechanical gain coefficient K of the adjusted synchronous machine mg Setting the reserve capacity P mg Inertia piece, difference adjustment coefficient R and reheating time constant T R High pressure turbine power fraction F H Finally, aggregating the frequency modulation parameters to obtain a single-machine equivalent model of each period of time;
and 4, step 4: and calculating the non-zero initial state of the transfer function of each time frequency response segment, and performing inverse Laplace transform on the transfer function of the inertia center frequency by using a non-zero initial state analysis method to obtain a frequency response segmented time domain analysis model of the whole network inertia center.
Preferably, in step 1, the key frequency modulation parameters are obtained by a human-computer interface according to specific operation requirements of the whole network synchronous machine.
Preferably, the governor output mechanical power transfer function of each synchronous machine in step 2 is calculated by equation (1):
wherein L is i 、R、F H 、T R 、H、K m 、P Step The power division coefficient, the difference adjustment coefficient, the power fraction of the high-pressure turbine, the reheating time constant, the inertia, the mechanical gain coefficient and the unbalanced power of the single-machine equivalent model are respectively, xi is a damping ratio, omega is n Is the natural frequency of the second order model;
calculating the time of reaching the amplitude limit value of the mechanical power; calculating the fractional power coefficient according to the speed regulator parameters of each synchronous machine, wherein the calculation formula is shown as formulas (2) to (4):
preferably, the mechanical gain coefficient K of each synchronous machine after being set in the step 3 mg Calculated from equation (5):
wherein S is g The installed capacity of each synchronous machine;
the equivalent frequency modulation parameters of the single-machine equivalent model are calculated by the following formulas (6) to (10):
wherein R is g 、F Hg 、T Rg 、H g Respectively the adjusted difference coefficient, the power fraction of the high-pressure turbine, the reheating time constant, the inertia and lambda of the synchronous machine g Are intermediate variables that assist in the computation.
Preferably, the frequency response segmented time domain analytic model of the whole network inertia center obtained in the step 4 is expressed as shown in equations (11) to (15):
wherein, ω is 1 、ω 2 、A 1 、A 2 、A 7 、A 8 、K 1 、K 5 The intermediate variables of the formula are conveniently written in columns, and no specific physical significance is realized.
Drawings
Fig. 1 is a topology diagram of a 3-machine 9-node synchronous machine system according to an embodiment of the present invention.
Fig. 2 is a simulation experiment diagram of the synchronous machine G1 for verifying the correctness of the mechanical power calculated by the algorithm.
Fig. 3 is a simulation experiment diagram of the synchronous machine G2 for verifying the correctness of the algorithm for calculating the mechanical power.
FIG. 4 is a simulation experiment diagram of a synchronous machine G3 for verifying the correctness of the mechanical power calculated by the algorithm.
Fig. 5 is a comparison graph of the effectiveness of the verification algorithm.
Detailed Description
The technical solutions in the embodiments of the present invention are 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 of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
The invention provides a synchronous machine system frequency response analysis and calculation method considering amplitude limiting of a speed regulator, which comprises the following specific steps of:
step 1: and acquiring key frequency modulation parameters.
According to specific operation requirements, key frequency modulation parameters are obtained through a human-computer interface, such as: spare capacity P mg Inertia H, difference adjustment coefficient R and reheating time constant T R High pressure turbine power fraction F H Damping coefficient xi and mechanical gain coefficient K m Disturbance power P Step And entering the step 2 after the completion.
Step 2: clipping time calculation
Calculating the transfer function of the output mechanical power of each synchronous machine according to the frequency modulation parameters of the speed regulators of each synchronous machine, and judging the amplitude limiting time of the speed regulators of each synchronous machine, wherein the transfer function of the output mechanical power of the speed regulators of each synchronous machine is calculated by the formula (1):
wherein L is i 、R、F H 、T R 、H、K m 、P Step The sub-power coefficient, the difference adjustment coefficient, the high-pressure turbine power fraction, the reheating time constant, the inertia, the mechanical gain coefficient and the unbalanced power of the single-machine equivalent model are respectively, xi is a damping ratio, and omega is n Is the natural frequency of the second order model;
calculating the time for the mechanical power to reach the amplitude limit value; calculating the sub-power coefficient according to the speed regulator parameters of each synchronous machine, wherein the calculation formula is shown as formulas (2) to (4):
and (4) entering the step 3 after the completion.
And step 3: multi-machine polymerization;
setting mechanical gain coefficient K of each synchronous machine mg Calculated from equation (5):
wherein S is g The installed capacity of each synchronous machine;
according to the set mechanical gain coefficient of the synchronous machine, setting inertia H, difference adjustment coefficient R and reheating time constant T R High pressure turbine power fraction F H And finally, aggregating the key frequency modulation parameters to obtain a single-machine equivalent model, wherein the equivalent frequency modulation parameters of the single-machine equivalent model are calculated by the following formulas (6) to (10):
wherein R is g 、F Hg 、T Rg 、H g Respectively is the difference adjustment coefficient, the high-pressure turbine power fraction, the reheating time constant, the inertia and lambda of the set synchronous machine g Are intermediate variables that assist in the computation.
And 4, step 4: non-zero initial state analysis.
Performing inverse Laplace transform on the transfer function of each frequency response segment by using a non-zero initial state analysis method to obtain a frequency response segmented time domain analysis model of the whole network inertia center, which is expressed as formulas (11) to (15):
wherein, ω is 1 、ω 2 、A 1 、A 2 、A 7 、A 8 、K 1 、K 5 The intermediate variables of the formula are conveniently written in columns, and no specific physical significance is realized.
Examples
The process of the invention is illustrated below by means of a specific example. Fig. 1 is a topological diagram of a 3-machine 9-node synchronous machine system, which includes three synchronous machines G1, G2, and G3 and 9 nodes as shown in the figure. The system voltage class in the region is 220kV, the total capacity of the synchronous machine is 402MW, the total load is 350MW, and the key frequency modulation parameters of the system are as follows:
the key frequency modulation parameter is brought into the mechanical power output power, as shown in figures 2-4, the frequency response of the system after disturbance obtained by the existing method and the method provided by the invention is shown in figure 5, and the comparison shows that the model can calculate the frequency response of the system after disturbance on the basis of considering the amplitude limiting of the speed regulator.
The invention has the following beneficial effects:
the method can carry out multi-machine aggregation according to the installed capacity of the synchronous machine and the total system capacity on the basis of the known key frequency modulation parameters of the power grid, consider the influence of the spare capacity in a single-machine equivalent model, and calculate and obtain a frequency response time domain analysis model of the whole-network inertia center by using a non-zero initial state analysis method, wherein the model can calculate and consider the power grid frequency response after the amplitude limiting disturbance of the speed regulator.
Claims (2)
1. A synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting is characterized by comprising the following steps:
step 1, acquiring key frequency modulation parameters of a whole network synchronous machine, including spare capacity P mg Inertia H, difference adjustment coefficient R and reheating time constant T R High pressure turbine power fraction F H Damping coefficient xi and mechanical gain coefficient K m Disturbance power P Step ;
Step 2, calculating a transfer function of the output mechanical power of each synchronizer according to the key frequency modulation parameters of each synchronizer in the whole-network synchronizer obtained in the step 1, and judging the amplitude limiting time of each synchronizer; the transfer function of each synchronous machine output mechanical power is calculated by the following formula (1):
wherein L is i 、R、F H 、T R 、H、K m 、P Step The sub-power coefficient, the difference adjustment coefficient, the high-pressure turbine power fraction, the reheating time constant, the inertia, the mechanical gain coefficient and the unbalanced power of the single-machine equivalent model are respectively, xi is a damping ratio, and omega is n Is the natural frequency of the second order model;
calculating the time for the mechanical power to reach the amplitude limit value; calculating the sub-power coefficient according to the speed regulator parameters of each synchronous machine, wherein the calculation formula is shown as formulas (2) to (4):
and 3, step 3: according to the amplitude limiting time of each synchronous machine, the frequency response process of the system is segmented, and each segment of synchronous machine participating in the frequency response sets a mechanical gain coefficient K according to the installed capacity and the total system capacity of the synchronous machine m Using the adjusted mechanical gain coefficient K of the synchronous machine mg Setting the reserve capacity P mg Inertia H, difference adjustment coefficient R and reheating time constant T R High pressure turbine power fraction F H Finally, the frequency modulation parameters are aggregated to obtainA single machine equivalent model of each period of time; the set mechanical gain coefficient K of the synchronous machine mg Calculated from equation (5):
wherein S is g The installed capacity of each synchronous machine;
the equivalent frequency modulation parameters of the single-machine equivalent model are calculated by the following formulas (6) to (10):
wherein R is g 、FH g 、T Rg 、H g Respectively is the difference adjustment coefficient, the high-pressure turbine power fraction, the reheating time constant, the inertia and lambda of the set synchronous machine g An intermediate variable for secondary calculations;
and 4, step 4: calculating a non-zero initial state of a transfer function of each time frequency response, and performing inverse Laplace transform on the transfer function of the inertia center frequency by using a non-zero initial state analysis method to obtain a frequency response segmented time domain analysis model of the whole network inertia center, wherein the frequency response segmented time domain analysis model is expressed as formulas (11) to (15):
wherein, ω is 1 、ω 2 、A 1 、A 2 、A 7 、A 8 、K 1 、K 5 In order to conveniently column-write intermediate variables of the formula, no specific physical significance exists.
2. The method for analyzing and calculating the frequency response of the synchronous machine system considering the amplitude limit of the speed regulator according to claim 1, wherein the key frequency modulation parameters are acquired by a human-computer interface in step 1 according to specific operation requirements of the full-network synchronous machine.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111593772.7A CN114243805B (en) | 2021-12-20 | 2021-12-20 | Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202111593772.7A CN114243805B (en) | 2021-12-20 | 2021-12-20 | Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting |
Publications (2)
Publication Number | Publication Date |
---|---|
CN114243805A CN114243805A (en) | 2022-03-25 |
CN114243805B true CN114243805B (en) | 2023-01-10 |
Family
ID=80762236
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202111593772.7A Active CN114243805B (en) | 2021-12-20 | 2021-12-20 | Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN114243805B (en) |
Families Citing this family (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115085287B (en) * | 2022-07-20 | 2023-04-18 | 华北电力大学 | Modeling method considering frequency stability constraint of amplitude limiting link of speed regulator |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014101515A1 (en) * | 2012-12-24 | 2014-07-03 | 国家电网公司 | Method for designing automatic generation control model under grid connection of intermittent energy |
CN110750882A (en) * | 2019-09-29 | 2020-02-04 | 武汉大学 | Wind power ratio limit value analytical calculation method considering frequency constraint |
CN112821424A (en) * | 2021-01-29 | 2021-05-18 | 国网辽宁省电力有限公司大连供电公司 | Power system frequency response analysis method based on data-model fusion drive |
CN113346482A (en) * | 2021-05-19 | 2021-09-03 | 电子科技大学 | Method for predicting frequency space-time distribution of wide area power system based on SFR model |
Family Cites Families (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106910142B (en) * | 2017-02-17 | 2020-09-22 | 三峡大学 | Method for calculating frequency characteristic of power system with wind power active power-frequency coupling effect |
CN113746134A (en) * | 2021-08-20 | 2021-12-03 | 华北电力大学(保定) | Feasible domain calculation method for inertia and primary frequency modulation control parameters of photovoltaic unit |
-
2021
- 2021-12-20 CN CN202111593772.7A patent/CN114243805B/en active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2014101515A1 (en) * | 2012-12-24 | 2014-07-03 | 国家电网公司 | Method for designing automatic generation control model under grid connection of intermittent energy |
CN110750882A (en) * | 2019-09-29 | 2020-02-04 | 武汉大学 | Wind power ratio limit value analytical calculation method considering frequency constraint |
CN112821424A (en) * | 2021-01-29 | 2021-05-18 | 国网辽宁省电力有限公司大连供电公司 | Power system frequency response analysis method based on data-model fusion drive |
CN113346482A (en) * | 2021-05-19 | 2021-09-03 | 电子科技大学 | Method for predicting frequency space-time distribution of wide area power system based on SFR model |
Non-Patent Citations (5)
Title |
---|
Adaptive parameter control strategy of VSG for improving system transient stability;W. Fan;《2017 IEEE 3rd International Future Energy Electronics Conference and ECCE Asia (IFEEC 2017 - ECCE Asia)》;20170727;全文 * |
基于改进SFR模型的含风电虚拟惯性/一次调频电力系统频率响应特性;李世春等;《可再生能源》;20190516(第05期);全文 * |
孤网发电机调速系统参数定量调整理论分析;李飞等;《电力系统保护与控制》;20161101(第21期);全文 * |
电力系统稳定计算中火电厂调速系统模型及其应用分析;宋新立等;《电网技术》;20081205(第23期);全文 * |
考虑控制量限幅的发电机非线性分散鲁棒汽门控制;常乃超等;《继电器》;20031225(第12期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN114243805A (en) | 2022-03-25 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109586269B (en) | Direct-current micro-grid virtual inertia control method and system considering parameter self-optimization | |
CN115085287B (en) | Modeling method considering frequency stability constraint of amplitude limiting link of speed regulator | |
CN103474986B (en) | A kind of long-time scale power system frequency fluctuation simulation method | |
CN104866665A (en) | Hybrid simulation method including power electronic equipment based on interface equivalence and interaction | |
CN101969199A (en) | Fault loss estimation method for risk assessment of transient rotor angle stability | |
CN105162164A (en) | Method of building low-order dynamic frequency response model with wind power integration system | |
CN108879667B (en) | Power grid closed-loop control power flow simulation method | |
CN106202793B (en) | A kind of mains frequency stability simulation method based on primary frequency modulation limitation | |
CN109449937B (en) | Method for determining new energy bearing capacity of power grid through system frequency stability constraint | |
CN106169115B (en) | Method for estimating lowest frequency of two-area power grid by using analytic-numerical method | |
CN106786529B (en) | A kind of distribution static security analysis method | |
CN103236026A (en) | Optimizing method of high-permeability throughput type power system planning scheme | |
CN114243805B (en) | Synchronous machine system frequency response analysis calculation method considering speed regulator amplitude limiting | |
CN102163845B (en) | Optimal configuration method of distributed generations (DG) based on power moment algorithm | |
CN114123344A (en) | Power system inertia evaluation method and device based on adaptive recursive least squares | |
CN114398761B (en) | Synchronous machine system frequency response analysis calculation method considering speed regulator reinforced dead zone | |
CN103618339B (en) | Total regulating power distribution method in automatic generation control | |
CN114204611B (en) | Frequency response analysis calculation method suitable for all damping states | |
CN106532785A (en) | Load frequency control method of considering novel cloud model | |
CN109149566A (en) | A kind of modeling method of the simulation model of the high-power minimum point prediction of missing lower frequency | |
CN104240151A (en) | Transient stability optimal correcting and control system and method for power system | |
CN114221396B (en) | Frequency response analytic calculation method considering general dead zone of speed regulator | |
CN104376195A (en) | Method for verifying transient state model of photovoltaic power station | |
CN113991647A (en) | Frequency response capacity planning-oriented random production simulation method for power system | |
CN114298478A (en) | Small disturbance stability identification method and system for wind power grid-connected 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 |