CN113113908A - Time domain analysis method and system suitable for frequency response of modern large power grid - Google Patents

Time domain analysis method and system suitable for frequency response of modern large power grid Download PDF

Info

Publication number
CN113113908A
CN113113908A CN202110287946.0A CN202110287946A CN113113908A CN 113113908 A CN113113908 A CN 113113908A CN 202110287946 A CN202110287946 A CN 202110287946A CN 113113908 A CN113113908 A CN 113113908A
Authority
CN
China
Prior art keywords
frequency response
power grid
large power
time domain
model
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Pending
Application number
CN202110287946.0A
Other languages
Chinese (zh)
Inventor
韩洪兴
鞠平
熊浩
秦川
付伟
江叶峰
周航
张振华
梁文腾
罗凯明
刘林
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
State Grid Jiangsu Electric Power Co Ltd
Hohai University HHU
Original Assignee
State Grid Jiangsu Electric Power Co Ltd
Hohai University HHU
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 Jiangsu Electric Power Co Ltd, Hohai University HHU filed Critical State Grid Jiangsu Electric Power Co Ltd
Priority to CN202110287946.0A priority Critical patent/CN113113908A/en
Publication of CN113113908A publication Critical patent/CN113113908A/en
Pending legal-status Critical Current

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00—Computer-aided design [CAD]
    • G06F30/20—Design optimisation, verification or simulation
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00—Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/06—Electricity, gas or water supply
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • 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]

Abstract

The invention provides a time domain analysis method and a time domain analysis system suitable for frequency response of a modern large power grid, wherein a general frequency response model suitable for the large power grid containing thermal power, hydropower and new energy power generation is constructed on the basis of a classical system frequency response model; according to the dynamic and steady state data of the actually measured large power grid, all parameters of the universal frequency response model are uniquely determined; and finally, obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation. The method is applicable to modern power grids containing various types of power supplies such as thermal power, hydropower, new energy power generation and the like, can improve the accuracy of frequency response calculation of a large power grid system, can quickly and accurately calculate the dynamic process of system frequency response when the power grid generates power disturbances of different sizes, has a simple and quick analysis process, and has wide application prospects and engineering practical values in the field of large power grid frequency safety analysis and control.

Description

Time domain analysis method and system suitable for frequency response of modern large power grid
Technical Field
The invention relates to the technical field of smart power grids, in particular to a frequency safety analysis and control technology of a modern large power grid containing thermal power, hydropower and new energy power generation, and particularly relates to a time domain analysis method and a time domain analysis system suitable for frequency response of the modern large power grid containing thermal power, hydropower and new energy power generation.
Background
The frequency is a core index of the operation of a large power grid and has particularly important influence on the power grid and users. The large grid frequency depends on the real-time balance of the power source and the load. When the water-fire power supply is mainly used in the tradition, the regulating capability is strong, and the frequency safety can be fully guaranteed. Due to uneven distribution of primary energy of a power grid in China, the proportion of long-distance alternating current and direct current power transmission and large-scale new energy power generation is rapidly improved. However, when power is transmitted across the region in a long distance, faults occur, and new energy resources are random and intermittent, so that the controllability of a power supply is reduced, the power supply and demand instantaneous balance difficulty is greatly increased, the frequency safety problem is increasingly severe, and the problem that the safe operation of a large power grid needs to be solved urgently is formed. In order to deal with the frequency safety risk, the primary task is to quickly and accurately calculate the frequency response of the large power grid caused by power shortage, so that corresponding frequency modulation measures can be made.
At present, the method for determining the system frequency response mainly comprises an all-state time domain simulation method, a linearization model analysis method, an artificial intelligence method, a single-machine equivalent model method and the like. Although the principle is clear and the application is common, the full-state time domain simulation method is limited by the accuracy of a full-network model, and the prediction precision of the full-state time domain simulation method on the dynamic frequency is difficult to improve; the linear model analysis method is to carry out partial linearization on the model on the basis of the full-state model method and also depends on the accuracy of the full-network model; the accuracy of the artificial intelligence method depends on a large amount of measured data, and the artificial intelligence method is difficult to popularize and apply in an actual power grid at present; the single-machine equivalent model is simple and is widely applied to the fields of setting of low-frequency load shedding, system frequency safety evaluation and the like, but the classical system frequency response model is only suitable for a thermal power generation system, and the application of the traditional system frequency response model to the frequency response of a modern large power grid containing thermal power, wind power and new energy power generation is limited.
Disclosure of Invention
Aiming at the defect that the rapidity and the accuracy of the existing frequency response calculation method are difficult to coordinate, the invention provides a time domain analysis method and a time domain analysis system which are suitable for the frequency response of a modern large power grid, not only can be suitable for the modern power grid containing various power supplies such as thermal power, hydropower and new energy power generation, but also can quickly and accurately calculate the dynamic process of the system frequency response when the power grid generates different power disturbances, thereby being beneficial to carrying out frequency control on the power grid fluctuation.
In order to achieve the above object, a first aspect of the present invention provides a time domain analysis method for frequency response of a modern large power grid including thermal power generation, hydroelectric power generation and new energy power generation, including the following steps:
step 1, constructing a general frequency response model suitable for a large power grid containing thermal power, hydropower and new energy power generation on the basis of a classical system frequency response model;
step 2, according to the dynamic and steady state data of the measured large power grid, uniquely determining all parameters of the universal frequency response model;
and 3, obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation.
On the basis of the first aspect of the present invention, a time domain analysis system of a modern large power grid frequency response is further provided, which includes:
the module is used for constructing a general frequency response model suitable for a large power grid containing thermal power, hydropower and new energy power generation on the basis of a classical system frequency response model;
a module for uniquely determining all parameters of the universal frequency response model according to the dynamic and steady-state data of the actually measured large power grid;
and the module is used for obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation.
On the basis of the first aspect of the present invention, a time domain analysis system of a modern large power grid frequency response is further provided, which includes:
one or more processors;
a memory storing instructions that are operable, when executed by the one or more processors, to cause the one or more processors to perform operations comprising processes to perform the aforementioned time domain resolution method of modern large grid frequency response.
Compared with the prior art, the technical scheme of the invention can obtain the following beneficial effects:
(1) the analysis method for the frequency response of the large power grid can be generally used for the large power grid containing various types of power supplies such as thermal power, hydropower and new energy power generation, and compared with the traditional calculation process of the classical model SFR, the accuracy of the frequency response calculation of the large power grid system can be improved through the processing of the method, because the traditional SFR model aims at the power grid system of the thermal power generation; the traditional frequency response obtained based on detailed model simulation has low calculation speed, too many parameters (such as a generator, a power grid, load and the like) and difficult parameter determination, and the analytic process of the large power grid frequency response of the invention has few determined parameters and accurate calculation;
(2) according to the analytic calculation method for the frequency response of the large power grid, provided by the invention, the dynamic frequency response process of the power grid under different power disturbances can be quickly (calculation time is nearly zero) and accurately calculated only after the model parameters and the analytic formula are determined in the whole calculation process, the analytic process is simple and quick, and the analytic calculation method has wide application prospects and engineering practical values in the field of frequency safety analysis and control of the large power grid.
It should be understood that all combinations of the foregoing concepts and additional concepts described in greater detail below can be considered as part of the inventive subject matter of this disclosure unless such concepts are mutually inconsistent. In addition, all combinations of claimed subject matter are considered a part of the presently disclosed subject matter.
The foregoing and other aspects, embodiments and features of the present teachings can be more fully understood from the following description taken in conjunction with the accompanying drawings. Additional aspects of the present invention, such as features and/or advantages of exemplary embodiments, will be apparent from the description which follows, or may be learned by practice of specific embodiments in accordance with the teachings of the present invention.
Drawings
The drawings are not intended to be drawn to scale. In the drawings, each identical or nearly identical component that is illustrated in various figures may be represented by a like numeral. For purposes of clarity, not every component may be labeled in every drawing. Embodiments of various aspects of the present invention will now be described, by way of example, with reference to the accompanying drawings, in which:
FIG. 1 is a flow chart of an implementation of the time domain analysis method applicable to the frequency response of a modern large power grid according to the present invention;
fig. 2 is a structural diagram of an IEEE 10 machine 39 node calculation system according to an embodiment of the present invention, in which generator sets corresponding to buses 31, 32, 33, and 37 are hydroelectric generator sets, generator sets corresponding to buses 30 and 34 are doubly-fed wind turbine generator sets, and generator sets corresponding to buses 35, 36, 38, and 39 are thermal generator sets.
FIG. 3 is a graph comparing the actual frequency response curve with the calculation result of the frequency response analytic formula of the present invention when the sample system has a 5% power shortage.
FIG. 4 is a graph comparing the actual frequency response curve with the frequency response analytic formula of the present invention when the sample system has a 2.5% power deficit.
FIG. 5 is a graph comparing the actual frequency response curve with the frequency response analytical formula calculation results of the present invention when the example system has a 7.5% power deficit.
Fig. 6 is a comparison graph of the measured frequency response curve of the east China power grid and the calculation result of the frequency response analytic formula of the invention.
Detailed Description
In order to better understand the technical content of the present invention, specific embodiments are described below with reference to the accompanying drawings.
In this disclosure, aspects of the present invention are described with reference to the accompanying drawings, in which a number of illustrative embodiments are shown. Embodiments of the present disclosure are not necessarily intended to include all aspects of the invention. It should be appreciated that the various concepts and embodiments described above, as well as those described in greater detail below, may be implemented in any of numerous ways, as the disclosed concepts and embodiments are not limited to any one implementation. In addition, some aspects of the present disclosure may be used alone, or in any suitable combination with other aspects of the present disclosure.
With reference to the drawings, the time domain analysis method suitable for the frequency response of the modern large power grid containing thermal power, hydropower and new energy power generation according to the embodiment of the invention comprises the following steps: step 1, constructing a general frequency response model suitable for a large power grid containing thermal power, hydropower and new energy power generation on the basis of a classical system frequency response model; step 2, according to the dynamic and steady state data of the measured large power grid, uniquely determining all parameters of the universal frequency response model; and 3, obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation.
Specifically, in step 1, on the basis of a classical system frequency response model, a general frequency response model suitable for a large power grid is established, including:
step 11, on the basis of a classical system frequency response model SFR, adopting a standard second-order transfer function Gm(s) the system replaces a classical system frequency response model SFR prime mover and a speed regulator equivalent model so as to be suitable for a large power grid containing thermal power, hydropower and new energy power generation; where the standard second order transfer function is expressed as follows:
Figure BDA0002981263140000041
in the formula, a0、a1、b0、b1Is the coefficient of the standard second order transfer function, Δ PmThe shortage of power of the power grid (the frequency of the power grid fluctuates, mainly because of the power imbalance of the power grid, and the imbalance of power generation and power utilization); delta f is the frequency response of the power grid, and s is a Laplace operator;
step 12, obtaining a general frequency response model G(s) of the large power grid, namely a power grid frequency response delta f and a power grid power disturbance delta PdThe general transfer function between:
Figure BDA0002981263140000042
in the formula, A0、A1、A2、A3、B0、B1The coefficients of a general frequency response model G(s) of the large power grid, namely the coefficients simplified aiming at the model G(s), are convenient for subsequent calculation and processing, H is the equivalent inertia of the power grid, namely the equivalent inertia of the large power grid, D is the equivalent damping coefficient of the power grid, K is the equivalent damping coefficient of the power gridGThe equivalent frequency modulation coefficient of the generator.
Because the classical system frequency response model SFR is a frequency response model aiming at the traditional thermal power grid, a speed regulator model special for thermal power is arranged in the traditional SFR model, the model structure and each parameter have practical physical significance, and the model structure and each parameter aiming at a certain thermal power prime motor and a speed regulator cannot be applied to the modern large power grid system containing thermal power, hydropower and new energy power generation. Therefore, the invention adopts a standard second-order transfer function to replace a prime motor and a speed regulator equivalent model in a classical SFR model aiming at the characteristics of a modern large power grid system so as to be suitable for power grids containing thermal power, hydropower and new energy power generation (photovoltaic or wind power).
In step 2, according to the dynamic and steady state data of the measured large power grid, all parameters of the universal frequency response model are uniquely determined, including:
step 21: according to the steady state data of the measured large power grid, namely the steady state power disturbance delta P in the measured large power gridd∞And steady state frequency deviation Δ f∞Obtaining the coefficient A in the general frequency response model G(s) of the large power grid3；
Figure BDA0002981263140000043
Step 22: according to the dynamic data of the modern measured large power grid, namely the dynamic data of power disturbance and frequency response, the rest coefficients theta in the general frequency response model G(s) of the large power grid are estimated by adopting a least square method, and the following are obtained:
Figure BDA0002981263140000051
θ＝[A0，A1，A2，B0，B1]T
in the formula, subscript c represents the frequency dynamic response data calculated by using the general frequency response model g(s), subscript a represents the actually measured power disturbance and frequency response data, and N represents the number of points of the actually measured power disturbance and frequency response data.
In the foregoing step 3, obtaining a time domain analysis of the modern large grid frequency response through inverse laplace transform includes:
step 31: to power shortage of electric network delta Pd(t) carrying out Laplace transformation to obtain complex frequency domain expression of frequency response of the large power grid system:
Figure BDA0002981263140000052
Figure BDA0002981263140000053
after Laplace transformation is carried out, a transfer function equation of the third order is substituted, and complex frequency domain expression of system frequency response is obtained:
Figure BDA0002981263140000054
in the formula, epsilon (t) is a unit step function;
step 32: performing inverse laplacian transform on Δ f(s) to obtain a time domain analysis formula Δ f (t) of the system frequency response:
Figure BDA0002981263140000055
in the formula, C0、C1、C2、T1、T2、ω2Theta is a parameter expressed by a complex frequency domain of the frequency response of the large power grid system and is obtained by calculation of inverse Laplace transform, and t is time.
In the specific implementation, in step 32, the transfer function characteristic equation A is considered0s3+A1s2+A2s+A30 root. One solid root and a pair of conjugated multiple roots are selected from the three roots:
Figure BDA0002981263140000056
the expansion is as follows:
Figure BDA0002981263140000061
if plural number K1X + yj, the coefficients are as follows:
Figure BDA0002981263140000062
by adopting inverse Laplace transform, a time domain analytic solution Δ f (t) of a frequency response can be obtained, and the system frequency response comprises three terms, namely a constant term, a monotonic attenuation term and an oscillation attenuation term:
Figure BDA0002981263140000063
the invention is illustrated below with reference to specific examples. The system used in the embodiment is shown in fig. 2, which is an IEEE 10 machine 39 node calculation system. The generator sets corresponding to the buses 31, 32, 33 and 37 are hydroelectric generator sets, the generator sets corresponding to the buses 30 and 34 are doubly-fed wind turbine generator sets, and the generator sets corresponding to the buses 35, 36, 38 and 39 are thermal power generator sets. The total load power of the system is 1104MW and the steady-state frequency is 50 Hz.
Setting 5% load power disturbance in a calculation system to cause the system to generate power shortage, obtaining parameters of a system general frequency response model G(s) according to the steps 1 and 2 of the method, wherein the parameters are shown in a table 1, and obtaining a system frequency response analytic formula according to the step C, wherein the analytic formula is as follows:
Δf(t)＝ΔPd[0.037040-0.052558e-t/0.581767+0.106958e-t/15.586034cos(0.138612t-1.412009)]ε(t)
TABLE 1 IEEE39 node system general frequency response model parameter estimation results
Figure BDA0002981263140000064
When a 5% power disturbance occurs, the example system actual frequency response is compared with the analytic formula calculation result, as shown in fig. 3. Therefore, the calculation result of the analytic formula is highly consistent with the actual frequency response result.
In order to verify the adaptability of the analytic formula to different power disturbances, 2.5% and 7.5% of load power disturbances are respectively set in the example system, and the comparison graph of the actual frequency response curve and the analytic formula calculation result is shown in fig. 4 and 5. Therefore, under different power disturbance sizes, the calculation result of the analytic expression is highly consistent with the actual frequency response result.
Further, the system frequency response of 20-month 20-day bingjin dc single-pole blocking fault in 2015 of east china power grid is taken as an example. In 2015, the east China power grid has 241.8GW of thermal power generating units, 20.18GW of hydroelectric power generating units, 14.01GW of nuclear power generating units, 9.08GW of wind power generating units and 3.77GW of photovoltaic power generation units, and has 7 direct current lines for transmitting 31.76GW of electric power to the east China power grid. 10/20/03: 05:14 in 2015, the guest gold direct current line has a single-pole blocking fault, so that the power shortage is 3700 MW. Since the accident occurs at midnight, the total pre-accident load is only around 160GW, i.e. the power shortage caused by the accident is about 2.313% of the total pre-accident load. In this event, the system frequency drops from 50.01Hz to 49.77Hz and then returns to 49.87 Hz.
Parameters of a general frequency response model G(s) of the east China Power grid obtained according to the steps 1 and 2 of the method are shown in Table 2, and a system frequency response analytic formula obtained according to the step C is as follows:
Δf(t)＝ΔPd[0.117622-0.061684e-t/3.943777+0.274610e-t/10.695187cos(0.141211t-1.843771)]ε(t)
TABLE 2 estimation results of the frequency response model for the east China Power grid
Figure BDA0002981263140000071
A comparison graph of the measured frequency response curve of the east China power grid and the calculation result of the analytic formula is shown in FIG. 6, which shows that the method can quickly and accurately calculate the dynamic process of the system frequency response when the power disturbance occurs to the power grid.
According to another embodiment of the present invention, in combination with the example shown in fig. 1, there is also provided a time domain analysis system suitable for a frequency response of a modern large power grid, including:
the module is used for constructing a general frequency response model suitable for a large power grid on the basis of a classical system frequency response model;
a module for uniquely determining all parameters of the universal frequency response model according to the dynamic and steady-state data of the actually measured large power grid;
and the module is used for obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation.
The specific implementation of each module may be implemented according to the exemplary implementation process of the above embodiment, and is not described herein again.
According to another embodiment of the present invention, in conjunction with the example shown in fig. 1, there is also provided a time domain resolving system suitable for a modern large grid frequency response, for example implemented in the form of a server or a server array, including:
one or more processors;
a memory storing instructions that are operable, when executed by the one or more processors, to cause the one or more processors to perform operations comprising performing a process of time domain resolution of a large grid frequency response of any of the foregoing embodiments, in particular the implementation of the embodiment of fig. 1.
Although the present invention has been described with reference to the preferred embodiments, it is not intended to be limited thereto. Those skilled in the art can make various changes and modifications without departing from the spirit and scope of the invention. Therefore, the protection scope of the present invention should be determined by the appended claims.

Claims (9)

1. A time domain analysis method suitable for modern large power grid frequency response is characterized by comprising the following steps:
step 1, constructing a general frequency response model suitable for a large power grid containing thermal power, hydropower and new energy power generation on the basis of a classical system frequency response model;
step 2, according to the dynamic and steady state data of the measured large power grid, uniquely determining all parameters of the universal frequency response model;
and 3, obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation.
2. The time domain analysis method suitable for the frequency response of the modern large power grid according to claim 1, wherein the building of the general frequency response model suitable for the large power grid is based on a classical system frequency response model, and comprises the following steps:
step 11, on the basis of a classical system frequency response model SFR, adopting a standard second-order transfer function Gm(s) replacing a classical system frequency response model SFR prime mover and a speed regulator equivalent model so as to be suitable for a large power grid; where the standard second order transfer function is expressed as follows:
Figure FDA0002981263130000011
in the formula, a0、a1、b0、b1Is the coefficient of the standard second order transfer function, Δ PmThe shortage of power of the power grid; delta f is the frequency response of the power grid, and s is a Laplace operator;
step 12, obtaining a general frequency response model G(s) of the large power grid, namely power grid frequency response delta f and power grid power disturbance delta PdThe general transfer function between:
Figure FDA0002981263130000012
in the formula, A0、A1、A2、A3、B0、B1Coefficients of general frequency response model G(s) for large power gridH is the equivalent inertia of the power grid, D is the equivalent damping coefficient of the power grid, KGThe equivalent frequency modulation coefficient of the generator.
3. The time domain analysis method suitable for the frequency response of the modern large power grid according to claim 2, wherein all the parameters of the universal frequency response model are uniquely determined according to the dynamic and steady state data of the measured large power grid, and the method comprises the following steps:
step 21: according to the steady state data of the measured large power grid, namely the steady state power disturbance delta P thereind∞And steady state frequency deviation Δ f∞，
Obtaining a coefficient A in the general frequency response model G(s) of the large power grid3；
Figure FDA0002981263130000013
Step 22: according to the dynamic data of the large power grid, namely the dynamic data of power disturbance and frequency response, which are actually measured, the rest coefficients theta in the general frequency response model G(s) of the large power grid are estimated by adopting a least square method, and the following are obtained:
Figure FDA0002981263130000021
θ＝[A0，A1，A2，B0，B1]T
in the formula, subscript c represents the frequency dynamic response data calculated by using the general frequency response model g(s), subscript a represents the actually measured power disturbance and frequency response data, and N represents the total amount of the actually measured power disturbance and frequency response data.
4. The time domain analysis method suitable for the modern large power grid frequency response according to claim 3, wherein the obtaining of the time domain analysis of the large power grid frequency response through inverse Laplace transform comprises:
step 31: for large power grid power shortage delta Pd(t) carrying out Laplace transformation to obtain complex frequency domain expression of frequency response of the large power grid system:
Figure FDA0002981263130000022
Figure FDA0002981263130000023
Figure FDA0002981263130000024
in the formula, epsilon (t) is a unit step function;
step 32: inverse laplacian transform is performed on Δ f(s) to obtain a time domain analysis formula Δ f (t) of the system frequency response:
Figure FDA0002981263130000025
in the formula, C0、C1、C2、T1、T2、ω2Theta is a parameter expressed by a complex frequency domain of the frequency response of the large power grid system and is obtained by calculation of inverse Laplace transform, and t is time.
5. A time domain analysis system adapted for modern large grid frequency response, comprising:
the module is used for constructing a general frequency response model suitable for a large power grid containing thermal power, hydropower and new energy power generation on the basis of a classical system frequency response model;
a module for uniquely determining all parameters of the universal frequency response model according to the dynamic and steady-state data of the actually measured large power grid;
and the module is used for obtaining time domain analysis of the frequency response of the large power grid through inverse Laplace transformation.
6. Time domain analysis system adapted to modern large power grid frequency response according to claim 5, wherein the module for building a generic frequency response model adapted to said large power grid based on a classical system frequency response model is arranged to use a standard second order transfer function G based on a classical system frequency response model SFRm(s) replacing a classical system frequency response model SFR prime mover and speed regulator equivalent model, and then obtaining a general frequency response model G(s) of the large power grid, namely a power grid frequency response delta f and a power grid power disturbance delta PdGeneral transfer function g(s) between:
Figure FDA0002981263130000031
in the formula, A0、A1、A2、A3、B0、B1Coefficient of general frequency response model G(s) of large power grid, H is equivalent inertia of power grid, D is equivalent damping coefficient of power grid, KGThe equivalent frequency modulation coefficient of the generator;
where the standard second order transfer function is expressed as follows:
Figure FDA0002981263130000032
in the formula, a0、a1、b0、b1Is the coefficient of the standard second order transfer function, Δ PmThe shortage of power of the power grid; and delta f is the frequency response of the power grid, and s is a Laplace operator.
7. The time domain analysis system for the frequency response of a modern large power grid according to claim 6, wherein the module for uniquely determining all parameters of the universal frequency response model based on measured large power grid dynamic and steady state data is configured to determine model parameters as follows:
according to the steady state data of the measured large power grid, namely the steady state power disturbance delta P thereind∞And steady state frequency deviation Δ f∞Obtaining the coefficient A in the general frequency response model G(s) of the large power grid3；
Figure FDA0002981263130000033
According to the dynamic data of the modern measured large power grid, namely the dynamic data of power disturbance and frequency response, the rest coefficients theta in the general frequency response model G(s) of the large power grid are estimated by adopting a least square method, and the following are obtained:
Figure FDA0002981263130000034
θ＝[A0，A1，A2，B0，B1]T
in the formula, subscript c represents the frequency dynamic response data calculated by using the general frequency response model g(s), subscript a represents the actually measured power disturbance and frequency response data, and N represents the total amount of the actually measured power disturbance and frequency response data.
8. Time domain resolution system applicable to modern large grid frequency response according to claim 7, characterized in that said module for obtaining the time domain resolution of the large grid frequency response by inverse laplace transform is arranged to the grid power deficit Δ Pd(t) carrying out Laplace transformation to obtain complex frequency domain expression of frequency response of the large power grid system:
Figure FDA0002981263130000041
Figure FDA0002981263130000042
Figure FDA0002981263130000043
in the formula, epsilon (t) is a unit step function;
then, by inverse laplacian transform on Δ f(s), a time domain analysis formula Δ f (t) of the system frequency response is obtained:
Figure FDA0002981263130000044
in the formula, C0、C1、C2、T1、T2、ω2Theta is a parameter expressed by a complex frequency domain of the frequency response of the large power grid system and is obtained by calculation of inverse Laplace transform, and t is time.
9. A time domain analysis system adapted for modern large grid frequency response, comprising:
one or more processors;
a memory storing instructions that are operable, when executed by the one or more processors, to cause the one or more processors to perform operations comprising performing a process of any one of claims 1-4.
CN202110287946.0A 2021-03-17 2021-03-17 Time domain analysis method and system suitable for frequency response of modern large power grid Pending CN113113908A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110287946.0A CN113113908A (en) 2021-03-17 2021-03-17 Time domain analysis method and system suitable for frequency response of modern large power grid

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110287946.0A CN113113908A (en) 2021-03-17 2021-03-17 Time domain analysis method and system suitable for frequency response of modern large power grid

Publications (1)

Publication Number Publication Date
CN113113908A true CN113113908A (en) 2021-07-13

Family

ID=76711931

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110287946.0A Pending CN113113908A (en) 2021-03-17 2021-03-17 Time domain analysis method and system suitable for frequency response of modern large power grid

Country Status (1)

Country Link
CN (1) CN113113908A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114204611A (en) * 2021-12-08 2022-03-18 华北电力大学 Frequency response analysis calculation method suitable for all damping states

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
李东辉等: ""电力系统频率响应的改进模型与参数估计"", 《电力工程技术》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114204611A (en) * 2021-12-08 2022-03-18 华北电力大学 Frequency response analysis calculation method suitable for all damping states

Similar Documents

Publication Publication Date Title
Bhadu et al. Robust linear quadratic Gaussian‐based discrete mode wide area power system damping controller
Azizipanah-Abarghooee et al. A new approach to the online estimation of the loss of generation size in power systems
CN106548410B (en) Method for evaluating voltage unbalance probability of power distribution network containing distributed power supply
Ju et al. Analytic assessment of the power system frequency security
CN111293686A (en) ARMAX system identification-based real-time evaluation method for inertia of power system
Han et al. Nonlinear adaptive robust control strategy of doubly fed induction generator based on virtual synchronous generator
CN111130135B (en) Power system inertia calculation method suitable for high-proportion new energy access
Zhao et al. Distributed continuation power flow method for integrated transmission and active distribution network
Fanglei et al. Estimating maximum penetration level of renewable energy based on frequency stability constrains in power grid
CN109149566B (en) Modeling method of simulation model for predicting lowest frequency point under high-power deficiency
CN113113908A (en) Time domain analysis method and system suitable for frequency response of modern large power grid
Zhang et al. Dynamic state estimation of a grid-connected converter of a renewable generation system using adaptive cubature Kalman filtering
Joseph et al. Unbiased minimum variance filter-based generator state estimation using PMU measurements for unknown generator input
Rao et al. Wideband impedance online identification of wind farms based on combined data-driven and knowledge-driven
Tang et al. Research on forced oscillations disturbance source locating through an energy approach
Vicente et al. Stochastic simulations and stability to determine maximum wind power penetration of an island network
CN113076628B (en) Analysis method and system suitable for frequency safety index of modern large power grid
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
Ayodele Comparative assessment of SVC and TCSC controllers on the small signal stability margin of a power system incorporating intermittent wind power generation
CN110048428B (en) Power system probability load flow calculation method based on probability conservation principle
Cai et al. Location and identification method for low frequency oscillation source considering control devices of generator
Sjödin The Price of Synchrony: Evaluating Transient Power Losses inRenewable Energy IntegratedPower Networks
Huang et al. Adaptive Cubature Kalman Filter Based Dynamic State Estimation for Grid-connected Photovoltaic System
Yang et al. Simplified Prediction Model of Frequency Nadir for Power Systems Penetrated With Renewable Energy
CN116191478B (en) Equivalent inertia evaluation and frequency response modeling method for multiple asynchronous motors

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
RJ01 Rejection of invention patent application after publication

Application publication date: 20210713

RJ01 Rejection of invention patent application after publication