CN116008656A - Power grid frequency change rate measurement window determining method, device and storage medium - Google Patents
Power grid frequency change rate measurement window determining method, device and storage medium Download PDFInfo
- Publication number
- CN116008656A CN116008656A CN202211625125.4A CN202211625125A CN116008656A CN 116008656 A CN116008656 A CN 116008656A CN 202211625125 A CN202211625125 A CN 202211625125A CN 116008656 A CN116008656 A CN 116008656A
- Authority
- CN
- China
- Prior art keywords
- change rate
- frequency change
- target
- rate measurement
- value
- 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
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y04—INFORMATION OR COMMUNICATION TECHNOLOGIES HAVING AN IMPACT ON OTHER TECHNOLOGY AREAS
- Y04S—SYSTEMS INTEGRATING TECHNOLOGIES RELATED TO POWER NETWORK OPERATION, COMMUNICATION OR INFORMATION TECHNOLOGIES FOR IMPROVING THE ELECTRICAL POWER GENERATION, TRANSMISSION, DISTRIBUTION, MANAGEMENT OR USAGE, i.e. SMART GRIDS
- Y04S10/00—Systems supporting electrical power generation, transmission or distribution
- Y04S10/50—Systems or methods supporting the power network operation or management, involving a certain degree of interaction with the load-side end user applications
Landscapes
- Supply And Distribution Of Alternating Current (AREA)
Abstract
The invention relates to the technical field of frequency change rate measurement of power systems, and discloses a method and a device for determining a power grid frequency change rate measurement window and a storage medium. The method comprises the steps of determining the total kinetic energy of a system and typical N-2 faults in a target operation mode, determining unbalanced power of the system after each typical N-2 fault occurs, and further calculating a theoretical value of a system frequency change rate corresponding to the corresponding target fault; determining a target frequency change rate measurement point; aiming at measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values, constructing a target constraint condition, taking a minimized frequency change rate measurement window as a target function, and constructing a frequency change rate measurement window optimization model; and carrying out optimization solving on the model to obtain an optimal frequency change rate measurement window and outputting the optimal frequency change rate measurement window. The method can solve the optimal value of the power grid frequency change rate measurement window under high-proportion new energy permeation, and has the advantages of wide applicability and high flexibility.
Description
Technical Field
The present invention relates to the field of power system frequency change rate measurement technologies, and in particular, to a method and apparatus for determining a power grid frequency change rate measurement window, and a storage medium.
Background
The system frequency change rate is a measurable index which is closely related to the inertia of the power system, accurate measurement of the index is helpful for determining the system inertia and disturbance rejection capability under high-proportion new energy permeation, and a judgment basis is possibly provided for the stability control of the regional power grid in the future. The method for measuring the frequency change rate mainly relates to a filter, a measuring algorithm and a measuring window, wherein the selection of the measuring window has great influence on the measurement of the frequency change rate, and is a key factor for accurately measuring the frequency change rate.
Currently, a measurement window with a fixed value is adopted for measuring the frequency change rate, and the fixed value is generally 500ms. For the traditional power system mainly comprising synchronous machines, the system has larger inertia and strong disturbance rejection capability, so that the system frequency change rate is smaller under the conventional large disturbance, and the frequency change rate can be accurately measured by adopting a mode of measuring the frequency change rate by a fixed-value measuring window.
However, along with the rapid development of the novel power system in the high-proportion new energy permeation scene, the synchronous machine gradually exits from the power system, the system inertia of the power grid is reduced, the disturbance rejection capability is weakened, the system frequency is changed rapidly in the large disturbance scene, the frequency change rate is increased, and the scheme of measuring the frequency change rate by adopting a fixed value measuring window is difficult to be suitable for the novel power system which is rapidly developed. The excessively long measurement window can cause delay of frequency change rate measurement and reduction of absolute value of measurement data, and meanwhile, the dynamic response characteristic of an actual local power grid system is difficult to reflect, and related measurement, configuration and response of part of new energy inverters are possibly difficult to match with a power system, so that new problems are caused. The too short frequency change rate measurement window can cause that the deviation of each sampling point of the power grid to the measured frequency change rate is larger, and meanwhile, the measured frequency change rate error is larger possibly due to the influence of waveform distortion and the like.
Disclosure of Invention
The invention provides a method, a device and a storage medium for determining a power grid frequency change rate measurement window, which can be suitable for calculation of frequency change rate measurement windows of power systems in different areas and different running modes, play an important role in inertia estimation, risk assessment and stability control decision of a novel power system under high-proportion new energy permeation in the future, and solve the technical problem that the existing scheme for measuring the frequency change rate by adopting a fixed-value measurement window is difficult to be suitable for the power system under the high-proportion new energy permeation.
The first aspect of the invention provides a method for determining a power grid frequency change rate measurement window, which comprises the following steps:
determining the total system kinetic energy of a target power grid in a target operation mode; the target operation mode is an operation mode corresponding to the minimum system rotational inertia of the target power grid;
determining a typical N-2 fault of the target power grid in the target operation mode to obtain a target fault set;
determining unbalanced power of the system after each target fault occurs; the target faults are faults in the target fault set;
calculating the system frequency change rate corresponding to the corresponding target faults according to the unbalanced power and the total kinetic energy of the system to obtain a system frequency change rate theoretical value corresponding to each target fault;
Determining each frequency change rate measuring point distributed in a key area of a target power grid to obtain a target frequency change rate measuring point set;
constructing target constraint conditions aiming at measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the measurement parameter value comprises a frequency change rate measurement value maximum value and a frequency change rate measurement deviation in a target fault occurrence period; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the target constraint includes: the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the maximum value of all frequency change rate measurement deviations is not more than a frequency change rate measurement deviation fixed value, and the maximum value of each frequency change rate measurement value is not more than a corresponding system frequency change rate theoretical value;
taking the minimum frequency change rate measurement window as an objective function, and establishing a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition;
And carrying out optimization solution on the frequency change rate measurement window optimization model to obtain and output an optimal frequency change rate measurement window.
According to an implementation manner of the first aspect of the present invention, the determining a total system kinetic energy of the target power grid in the target operation manner includes:
determining a target operation mode of the target power grid by taking the minimum starting number of the synchronous units, the maximum total output of the new energy units and the minimum inertia constant of the synchronous units as constraint conditions of the operation mode;
and calculating the total system kinetic energy of the target power grid in a target operation mode.
According to an implementation manner of the first aspect of the present invention, the calculating the total system kinetic energy of the target power grid in the target operation manner includes:
calculating the total system kinetic energy of the target power grid in a target operation mode according to the following formula:
E sys ＝E SG +E IM +E V(VS) +E V(CS) +E Load
wherein E is sys Representing the total system kinetic energy of the target power grid in a target operation mode, E SG Representing the total kinetic energy of the synchronous units, E IM Representing the total kinetic energy of the asynchronous induction motor, E V(VS) Energy form representing virtual inertia of voltage source type, E V(CS) Energy form representing virtual inertia of current source, E Load The energy form representing the static load voltage equivalent inertia.
According to one implementation manner of the first aspect of the present invention, the calculating the system frequency change rate corresponding to the corresponding target fault according to the unbalanced power and the total kinetic energy of the system to obtain the theoretical value of the system frequency change rate corresponding to each target fault includes:
calculating the system frequency change rate corresponding to each target fault according to the following formula:
in the formula, roCoF theory (i) For the system frequency change rate corresponding to the target fault i, f 0 For the nominal frequency of the system, ΔP i Unbalanced power of system after occurrence of target fault i, E sys And the total kinetic energy of the system of the target power grid in the target operation mode is obtained.
According to one implementation manner of the first aspect of the present invention, the target constraint condition uses a target fault occurrence time as a lower bound of the target fault occurrence period, and uses a sum of the target fault occurrence time and a preset fault time factor as an upper bound of the target fault occurrence period.
According to one implementation manner of the first aspect of the present invention, the establishing a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition includes:
the frequency change rate measurement window optimization model is established as follows:
s.t.
T 0 ≤T≤T max
max{[max(RoCoF i,j,T (t))-min(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
max{RoCoF i,j,T (t),t∈[t 0 ,t 0 +τ]}≤RoCoF theory (i)
Wherein T represents the frequency change rate measurement window value, T 0 Representing the initial value of a preset window, T max Representing the maximum allowable value of a preset window omega F Representing the target fault set, Ω M Representing a set of target frequency rate of change measurement points, rocofs i,j,T (T) represents a time T frequency change rate measurement value, T, of the frequency change rate measurement point j corresponding to the target fault i in the occurrence period of the target fault when the frequency change rate measurement window value is T 0 For the occurrence time of the target fault, tau is a preset fault time factor, rocofs theory (i) And measuring a deviation constant value for the system frequency change rate corresponding to the target fault i and the frequency change rate.
According to one implementation manner of the first aspect of the present invention, the performing optimization solution on the frequency change rate measurement window optimization model includes:
and solving the frequency change rate measurement window optimization model by gradually increasing the frequency change rate measurement window value.
A second aspect of the present invention provides a power grid frequency change rate measurement window determining apparatus, including:
the first determining module is used for determining the total system kinetic energy of the target power grid in a target running mode; the target operation mode is an operation mode corresponding to the minimum system rotational inertia of the target power grid;
The second determining module is used for determining a typical N-2 fault of the target power grid in the target operation mode to obtain a target fault set;
a third determining module, configured to determine an unbalanced power of the system after each target fault occurs; the target faults are faults in the target fault set;
the calculation module is used for calculating the system frequency change rate corresponding to the corresponding target faults according to the unbalanced power and the total kinetic energy of the system to obtain the system frequency change rate theoretical value corresponding to each target fault;
the fourth determining module is used for determining each frequency change rate measuring point spread over the key area of the target power grid to obtain a target frequency change rate measuring point set;
the first construction module is used for constructing target constraint conditions aiming at measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the measurement parameter value comprises a frequency change rate measurement value maximum value and a frequency change rate measurement deviation in a target fault occurrence period; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the target constraint includes: the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the maximum value of all frequency change rate measurement deviations is not more than a frequency change rate measurement deviation fixed value, and the maximum value of each frequency change rate measurement value is not more than a corresponding system frequency change rate theoretical value;
The second construction module is used for taking the minimum frequency change rate measurement window as an objective function, and establishing a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition;
and the solving module is used for carrying out optimization solving on the frequency change rate measurement window optimization model to obtain and output an optimal frequency change rate measurement window.
According to one implementation manner of the second aspect of the present invention, the first determining module includes:
the determining unit is used for determining a target operation mode of the target power grid by taking the minimum starting number of the synchronous units, the maximum total output of the new energy units and the minimum inertia constant of the synchronous units as constraint conditions of the operation mode;
and the first calculation unit is used for calculating the total system kinetic energy of the target power grid in a target operation mode.
According to one possible implementation manner of the second aspect of the present invention, the first computing unit is specifically configured to:
calculating the total system kinetic energy of the target power grid in a target operation mode according to the following formula:
E sys ＝E SG +E IM +E V(VS) +E V(CS) +E Load
wherein E is sys Representing the total system kinetic energy of the target power grid in a target operation mode, E SG Representing the total kinetic energy of the synchronous units, E IM Representing the total kinetic energy of the asynchronous induction motor, E V(VS) Energy form representing virtual inertia of voltage source type, E V(CS) Energy form representing virtual inertia of current source, E Load The energy form representing the static load voltage equivalent inertia.
According to one manner in which the second aspect of the present invention can be implemented, the computing module includes:
the second calculating unit is used for calculating the system frequency change rate corresponding to each target fault according to the following formula:
in the formula, roCoF theory (i) For the system frequency change rate corresponding to the target fault i, f 0 For the nominal frequency of the system, ΔP i Unbalanced power of system after occurrence of target fault i, E sys And the total kinetic energy of the system of the target power grid in the target operation mode is obtained.
According to one implementation manner of the second aspect of the present invention, the target constraint condition uses a target fault occurrence time as a lower bound of the target fault occurrence period, and uses a sum of the target fault occurrence time and a preset fault time factor as an upper bound of the target fault occurrence period.
According to one manner of implementation of the second aspect of the present invention, the second building block includes:
the model building unit is used for building a frequency change rate measurement window optimization model as follows:
T 0 ≤T≤T max
max{[max(RoCoF i,j,T (t))-min(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
max{RoCoF i,j,T (t),t∈[t 0 ,t 0 +τ]}≤RoCoF theory (i)
Wherein T represents the frequency change rate measurement window value, T 0 Representing the initial value of a preset window, T max Representing the maximum allowable value of a preset window omega F Representing the target fault set, Ω M Representing a set of target frequency rate of change measurement points, rocofs i,j,T (T) represents a time T frequency change rate measurement value, T, of the frequency change rate measurement point j corresponding to the target fault i in the occurrence period of the target fault when the frequency change rate measurement window value is T 0 For the purpose ofMarking fault occurrence time, wherein tau is a preset fault time factor, and RoCoF theory (i) And measuring a deviation constant value for the system frequency change rate corresponding to the target fault i and the frequency change rate.
According to one manner of implementation of the second aspect of the present invention, the solving module includes:
and the solving unit is used for solving the frequency change rate measurement window optimization model by gradually increasing the frequency change rate measurement window value.
A third aspect of the present invention provides a power grid frequency change rate measurement window determining apparatus, including:
a memory for storing instructions; the instructions are used for realizing the method for determining the power grid frequency change rate measurement window in the mode which can be realized by any one of the above;
and the processor is used for executing the instructions in the memory.
A fourth aspect of the present invention is a computer-readable storage medium having stored thereon a computer program which, when executed by a processor, implements the grid frequency change rate measurement window determination method according to any one of the modes that can be implemented as described above.
From the above technical scheme, the invention has the following advantages:
the method takes an operation mode corresponding to the minimum system rotational inertia of a target power grid as a target operation mode, determines the total kinetic energy of the system and typical N-2 faults in the target operation mode, determines unbalanced power of the system after each typical N-2 fault occurs, and further calculates the system frequency change rate corresponding to the corresponding target fault to obtain a system frequency change rate theoretical value corresponding to each target fault; determining each frequency change rate measuring point distributed in a key area of a target power grid as a target frequency change rate measuring point; aiming at measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values, constructing a target constraint condition, taking a minimized frequency change rate measurement window as a target function, and constructing a frequency change rate measurement window optimization model; carrying out optimization solution on the frequency change rate measurement window optimization model to obtain and output an optimal frequency change rate measurement window; according to the method, the optimal value of the frequency change rate measurement window of the power grid under high-proportion new energy permeation can be solved, corresponding parameters can be flexibly adjusted according to actual requirements, the solved frequency change rate measurement window value is ensured to be suitable for a system, when parameters such as relevant stability control and power electronic control of the system change, the updated frequency change rate measurement window value can be calculated only by changing constraint values, and therefore the method is suitable for calculating the frequency change rate measurement window of power systems in different areas and in different operation modes, and plays an important role in inertia estimation, risk assessment and stability control decision of the novel power system under future high-proportion new energy permeation.
Drawings
In order to more clearly illustrate the embodiments of the invention or the technical solutions of the prior art, the drawings which are used in the description of the embodiments or the prior art will be briefly described, it being obvious that the drawings in the description below are only some embodiments of the invention, and that other drawings can be obtained from these drawings without inventive faculty for a person skilled in the art.
FIG. 1 is a flowchart of a method for determining a power grid frequency change rate measurement window according to an alternative embodiment of the present invention;
fig. 2 is a block diagram of structural connection of a device for determining a frequency change rate measurement window of a power grid according to an alternative embodiment of the present invention.
Reference numerals:
1-a first determination module; 2-a second determination module; 3-a third determination module; 4-a calculation module; 5-a fourth determination module; 6-a first building block; 7-a second building block; 8-solving module.
Detailed Description
The embodiment of the invention provides a method, a device and a storage medium for determining a power grid frequency change rate measurement window, which are used for solving the technical problem that the existing scheme for measuring the frequency change rate by adopting a fixed-value measurement window is difficult to be suitable for a power system under high-proportion new energy permeation.
In order to make the objects, features and advantages of the present invention more comprehensible, the technical solutions in the embodiments of the present invention are described in detail below with reference to the accompanying drawings, and it is apparent that the embodiments described below are only some embodiments of the present invention, but not all embodiments of the present invention. All other embodiments, which can be made by those skilled in the art based on the embodiments of the invention without making any inventive effort, are intended to be within the scope of the invention.
The invention provides a method for determining a network frequency change rate measurement window.
Referring to fig. 1, fig. 1 shows a flowchart of a method for determining a power grid frequency change rate measurement window according to an embodiment of the present invention.
The method for determining the network frequency change rate measurement window provided by the embodiment of the invention comprises the steps S1-S8.
Step S1, determining the total system kinetic energy of a target power grid in a target operation mode; the target operation mode is an operation mode corresponding to the condition that the system moment of inertia of the target power grid is minimum.
When determining the total system kinetic energy of the target power grid in the target operation mode, the target operation mode needs to be determined. In the embodiment of the application, the operation mode corresponding to the time when the system moment of inertia of the target power grid is minimum is taken as the target operation mode.
The determination of the typical limit operation mode of the target power grid needs to comprehensively consider 3 constraint conditions:
(1) The starting quantity of the synchronous units is as small as possible;
(2) The total output of the new energy unit is as high as possible;
(3) The synchronous machine set selects the machine set with small opening inertia (small inertia constant) as far as possible to start.
By considering the above three conditions in combination, a typical limit operation mode of the power grid can be established, in which the system inertia is small, the system frequency is most sensitive to disturbance, and therefore, the boundary (i.e., the maximum value) of the frequency change rate measurement window can be determined based on the operation mode.
Therefore, as a specific implementation mode, the target operation mode of the target power grid is determined by taking the minimum starting number of the synchronous units, the maximum total output of the new energy units and the minimum inertia constant of the synchronous units as constraint conditions of the operation mode.
In one possible implementation, the total system kinetic energy of the target grid in the target operating mode is calculated according to the following formula:
E sys ＝E SG +E IM +E V(VS) +E V(CS) +E Load
wherein E is sys Representing the total system kinetic energy of the target power grid in a target operation mode, E SG Representing the total kinetic energy of the synchronous units, E IM Representing the total kinetic energy of the asynchronous induction motor, E V(VS) Energy form representing virtual inertia of voltage source type, E V(CS) Energy form representing virtual inertia of current source, E Load The energy form representing the static load voltage equivalent inertia.
And S2, determining a typical N-2 fault of the target power grid in the target operation mode, and obtaining a target fault set.
If there are P typical N-2 faults of the target power grid in the target operation mode, the target fault set may be expressed as:
Ω F ＝{i|i＝1,2,…,P}
in omega F Representing the target set of faults.
It should be noted that the type of the typical N-2 fault may be selected from the existing N-2 fault types according to actual situations, which is not limited in this embodiment.
Step S3, determining unbalanced power of the system after each target fault occurs; the target faults are faults in the target fault set.
And S4, calculating the system frequency change rate corresponding to the corresponding target faults according to the unbalanced power and the total kinetic energy of the system, and obtaining the theoretical value of the system frequency change rate corresponding to each target fault.
In one implementation manner, the calculating the system frequency change rate corresponding to the corresponding target fault according to the unbalanced power and the total kinetic energy of the system to obtain the theoretical value of the system frequency change rate corresponding to each target fault includes:
Calculating the system frequency change rate corresponding to each target fault according to the following formula:
in the formula, roCoF theory (i) For the system frequency change rate corresponding to the target fault i, f 0 For the nominal frequency of the system, ΔP i Unbalanced power of system after occurrence of target fault i, E sys And the total kinetic energy of the system of the target power grid in the target operation mode is obtained.
In this embodiment, the system frequency change rate corresponding to each target fault is calculated according to the above formula, and the method is simple and convenient.
And S5, determining each frequency change rate measuring point spread over the key area of the target power grid to obtain a target frequency change rate measuring point set.
The target power grid key area can be set according to actual conditions. If Q frequency change rate measurement points are arranged in the target power grid critical area, the target frequency change rate measurement point set may be expressed as:
Ω M ＝{j|j＝1,2,…,Q}
in omega M Representing the set of target frequency change rate measurement points.
And S6, constructing target constraint conditions according to measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values.
The target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the measurement parameter value comprises a frequency change rate measurement value maximum value and a frequency change rate measurement deviation in a target fault occurrence period; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the target constraint includes: the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the maximum value of all frequency change rate measurement deviations is not more than a frequency change rate measurement deviation fixed value, and the maximum value of each frequency change rate measurement value is not more than a corresponding system frequency change rate theoretical value.
In the present embodiment, three constraint conditions are set as target preset conditions. The first constraint condition is that the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the second constraint condition is that the maximum value of all frequency change rate measurement deviations is not larger than a frequency change rate measurement deviation fixed value, and the third constraint condition is that the maximum value of each frequency change rate measurement value is not larger than a corresponding system frequency change rate theoretical value.
As a specific embodiment, the frequency change rate measurement deviation in the target failure occurrence period is a difference between the maximum value of the frequency change rate measurement value and the minimum value of the frequency change rate measurement value in the target failure occurrence period. Based on this embodiment, the corresponding second constraint is as follows:
max{[max(RoCoF i,j,T (t))-min(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
in omega F Representing the target fault set, Ω M Representing a set of target frequency rate of change measurement points, rocofs i,j,T (T) represents a time T frequency change rate measurement value, T, of the frequency change rate measurement point j corresponding to the target fault i in the occurrence period of the target fault when the frequency change rate measurement window value is T 0 For the occurrence time of the target fault, τ is a preset fault time factor, and e is a frequency change rate measurement deviation fixed value.
In other embodiments, an average of the maximum value of the frequency change rate measurement value and the next-largest value of the frequency change rate measurement value in the target failure occurrence period may be calculated as the first calculation value, an average of the minimum value of the frequency change rate measurement value and the next-smallest value of the frequency change rate measurement value in the target failure occurrence period may be calculated as the second calculation value, and a difference between the first calculation value and the second calculation value may be calculated as the frequency change rate measurement deviation in the target failure occurrence period. Based on this embodiment, the corresponding second constraint is as follows:
max{[max'(RoCoF i,j,T (t))-min'(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
wherein:
wherein max -1 (RoCoF i,j,T (t)) is the next largest value of the frequency change rate measurement, min -1 (RoCoF i,j,T (t)) is the frequency change rate measurement value in the target failure occurrence period.
In one implementation manner, the target constraint condition uses a target fault occurrence time as a lower bound of the target fault occurrence period, and uses a sum of the target fault occurrence time and a preset fault time factor as an upper bound of the target fault occurrence period.
And S7, taking the minimum frequency change rate measurement window as an objective function, and establishing a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition.
In one manner that can be implemented, the frequency change rate measurement window optimization model is built as follows:
s.t.
T 0 ≤T≤T max
max{[max(RoCoF i,j,T (t))-min(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
max{RoCoF i,j,T (t),t∈[t 0 ,t 0 +τ]}≤RoCoF theory (i)
wherein T represents the frequency change rate measurement window value, T 0 Representing the initial value of a preset window, T max Representing the maximum allowable value of a preset window omega F Representing the target fault set, Ω M Representing a set of target frequency rate of change measurement points, rocofs i,j,T (T) represents a time T frequency change rate measurement value, T, of the frequency change rate measurement point j corresponding to the target fault i in the occurrence period of the target fault when the frequency change rate measurement window value is T 0 For the occurrence time of the target fault, tau is a preset fault time factor, rocofs theory (i) And measuring a deviation constant value for the system frequency change rate corresponding to the target fault i and the frequency change rate.
The frequency change rate measurement deviation fixed value epsilon is determined by the measurement precision required by the power grid stability control; the preset fault time factor tau may be determined by the grid fault response characteristics; presetting a window initial value T 0 Can be determined by the system original measurement window value; presetting a maximum allowable value T of a window max The method can be determined by the constant value of various stability control systems of the power grid and the response time of power electronic equipment, and the like, so that the measurement delay caused by the measurement window is ensured not to influence the stability control strategy action of the system and the action of the power electronic equipment participating in the system control.
And S8, carrying out optimization solution on the frequency change rate measurement window optimization model to obtain an optimal frequency change rate measurement window and outputting the optimal frequency change rate measurement window.
In one implementation manner, the optimizing the frequency change rate measurement window optimization model includes:
and solving the frequency change rate measurement window optimization model by gradually increasing the frequency change rate measurement window value.
In this embodiment, the calculation meaning of the frequency change rate measurement window optimization model is that all typical N-2 faults and all frequency change rate measurement points in a fault set are calculated through simulation, the frequency change rate measurement deviation at each moment in a period of time after the fault is increased gradually, and when a proper frequency change rate measurement window value is selected so that all frequency change rate measurement deviations are smaller than a given frequency change rate measurement deviation fixed value, the frequency change rate measurement window value can be determined to be a required value.
In consideration of the calculation amount problem, when the frequency change rate measurement window is solved, the time of each increase can be flexibly adjusted. As a way of being able to be achieved, each time the frequency change rate measurement window value is set to increase by at least 10ms while the frequency change rate measurement window value is gradually increased.
The embodiment of the invention can have unexpected technical effects:
corresponding parameters can be flexibly adjusted according to actual requirements, and the solved frequency change rate measurement window value is ensured to be suitable for a system; when parameters such as system related stability control and power electronic control change, the updated frequency change rate measurement window value can be calculated only by changing the constraint fixed value, the flexibility is high, the method is applicable to calculation of frequency change rate measurement windows of power systems in different areas and different running modes, and plays an important role in inertia estimation, risk assessment and stability control decision of a novel power system under high-proportion new energy penetration in the future.
The invention also provides a device for determining the power grid frequency change rate measurement window, which can be used for executing the power grid frequency change rate measurement window determining method according to any one of the embodiments of the invention.
Referring to fig. 2, fig. 2 is a block diagram illustrating structural connection of a device for determining a frequency change rate measurement window of a power grid according to an embodiment of the present invention.
The device for determining the power grid frequency change rate measurement window provided by the embodiment of the invention comprises the following components:
the first determining module 1 is used for determining the total system kinetic energy of the target power grid in a target operation mode; the target operation mode is an operation mode corresponding to the minimum system rotational inertia of the target power grid;
The second determining module 2 is used for determining a typical N-2 fault of the target power grid in the target operation mode to obtain a target fault set;
a third determining module 3, configured to determine an unbalanced power of the system after each target fault occurs; the target faults are faults in the target fault set;
the calculation module 4 is used for calculating the system frequency change rate corresponding to the corresponding target faults according to the unbalanced power and the total kinetic energy of the system to obtain the system frequency change rate theoretical value corresponding to each target fault;
a fourth determining module 5, configured to determine each frequency change rate measurement point distributed over the target power grid key area, so as to obtain a target frequency change rate measurement point set;
the first construction module 6 is configured to construct a target constraint condition for measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the measurement parameter value comprises a frequency change rate measurement value maximum value and a frequency change rate measurement deviation in a target fault occurrence period; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the target constraint includes: the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the maximum value of all frequency change rate measurement deviations is not more than a frequency change rate measurement deviation fixed value, and the maximum value of each frequency change rate measurement value is not more than a corresponding system frequency change rate theoretical value;
A second construction module 7, configured to set up a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition by using the minimized frequency change rate measurement window as an objective function;
and the solving module 8 is used for carrying out optimization solving on the frequency change rate measurement window optimization model to obtain and output an optimal frequency change rate measurement window.
In one possible implementation, the first determining module 1 includes:
the determining unit is used for determining a target operation mode of the target power grid by taking the minimum starting number of the synchronous units, the maximum total output of the new energy units and the minimum inertia constant of the synchronous units as constraint conditions of the operation mode;
and the first calculation unit is used for calculating the total system kinetic energy of the target power grid in a target operation mode.
In one implementation, the first computing unit is specifically configured to:
calculating the total system kinetic energy of the target power grid in a target operation mode according to the following formula:
E sys ＝E SG +E IM +E V ( VS )+E V ( CS )+E Load
wherein E is sys Representing the total system kinetic energy of the target power grid in a target operation mode, E SG Representing the total kinetic energy of the synchronous units, E IM Representing the total kinetic energy of the asynchronous induction motor, E V(VS) Energy form representing virtual inertia of voltage source type, E V(CS) Energy form representing virtual inertia of current source, E Load The energy form representing the static load voltage equivalent inertia.
In one possible implementation, the computing module 4 comprises:
the second calculating unit is used for calculating the system frequency change rate corresponding to each target fault according to the following formula:
in the formula, roCoF theory (i) For the target fault iThe rate of change of system frequency, f 0 For the nominal frequency of the system, ΔP i Unbalanced power of system after occurrence of target fault i, E sys And the total kinetic energy of the system of the target power grid in the target operation mode is obtained.
In one implementation manner, the target constraint condition uses a target fault occurrence time as a lower bound of the target fault occurrence period, and uses a sum of the target fault occurrence time and a preset fault time factor as an upper bound of the target fault occurrence period.
In one possible implementation, the second building block 7 comprises:
the model building unit is used for building a frequency change rate measurement window optimization model as follows:
s.t.
T 0 ≤T≤T max
max{[max(RoCoF i,j,T (t))-min(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
max{RoCoF i,j,T (t),t∈[t 0 ,t 0 +τ]}≤RoCoF theory (i)
wherein T represents the frequency change rate measurement window value, T 0 Representing the initial value of a preset window, T max Representing the maximum allowable value of a preset window omega F Representing the target fault set, Ω M Representing a set of target frequency rate of change measurement points, rocofs i,j,T (T) represents a time T frequency change rate measurement value, T, of the frequency change rate measurement point j corresponding to the target fault i in the occurrence period of the target fault when the frequency change rate measurement window value is T 0 For the occurrence time of the target fault, tau is a preset fault time factor, rocofs theory (i) And measuring a deviation constant value for the system frequency change rate corresponding to the target fault i and the frequency change rate.
In one possible implementation, the solving module 8 includes:
and the solving unit is used for solving the frequency change rate measurement window optimization model by gradually increasing the frequency change rate measurement window value.
The invention also provides a device for determining the power grid frequency change rate measurement window, which comprises the following steps:
a memory for storing instructions; the instructions are used for implementing the method for determining the power grid frequency change rate measurement window according to any one of the embodiments;
and the processor is used for executing the instructions in the memory.
The invention also provides a computer readable storage medium, wherein a computer program is stored on the computer readable storage medium, and the computer program is executed by a processor to realize the method for determining the power grid frequency change rate measurement window according to any one of the embodiments.
It will be clearly understood by those skilled in the art that, for convenience and brevity of description, specific working procedures of the above-described apparatus and module may refer to corresponding procedures in the foregoing method embodiments, and specific beneficial effects of the above-described apparatus and module may refer to corresponding beneficial effects in the foregoing method embodiments, which are not described herein again.
In the several embodiments provided in this application, it should be understood that the disclosed apparatus and method may be implemented in other ways. For example, the apparatus embodiments described above are merely illustrative, and for example, the division of the modules is merely a logical function division, and there may be additional divisions when actually implemented, for example, multiple modules or components may be combined or integrated into another apparatus, or some features may be omitted or not performed.
The modules described as separate components may or may not be physically separate, and components shown as modules may or may not be physical modules, i.e., may be located in one place, or may be distributed over a plurality of network modules. Some or all of the modules may be selected according to actual needs to achieve the purpose of the solution of this embodiment.
In addition, each functional module in each embodiment of the present invention may be integrated into one processing module, or each module may exist alone physically, or two or more modules may be integrated into one module. The integrated modules may be implemented in hardware or in software functional modules.
The integrated modules, if implemented in the form of software functional modules and sold or used as a stand-alone product, may be stored in a computer readable storage medium. Based on such understanding, the technical solution of the present invention may be embodied essentially or in part or all of the technical solution or in part in the form of a software product stored in a storage medium, including instructions for causing a computer device (which may be a personal computer, a server, or a network device, etc.) to perform all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a random access Memory (RAM, random Access Memory), a magnetic disk, or an optical disk, or other various media capable of storing program codes.
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 technical features thereof can be replaced by equivalents; such modifications and substitutions do not depart from the spirit and scope of the technical solutions of the embodiments of the present invention.
Claims (10)
1. A method for determining a frequency change rate measurement window of a power grid, comprising:
determining the total system kinetic energy of a target power grid in a target operation mode; the target operation mode is an operation mode corresponding to the minimum system rotational inertia of the target power grid;
determining a typical N-2 fault of the target power grid in the target operation mode to obtain a target fault set;
determining unbalanced power of the system after each target fault occurs; the target faults are faults in the target fault set;
calculating the system frequency change rate corresponding to the corresponding target faults according to the unbalanced power and the total kinetic energy of the system to obtain a system frequency change rate theoretical value corresponding to each target fault;
Determining each frequency change rate measuring point distributed in a key area of a target power grid to obtain a target frequency change rate measuring point set;
constructing target constraint conditions aiming at measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the measurement parameter value comprises a frequency change rate measurement value maximum value and a frequency change rate measurement deviation in a target fault occurrence period; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the target constraint includes: the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the maximum value of all frequency change rate measurement deviations is not more than a frequency change rate measurement deviation fixed value, and the maximum value of each frequency change rate measurement value is not more than a corresponding system frequency change rate theoretical value;
taking the minimum frequency change rate measurement window as an objective function, and establishing a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition;
And carrying out optimization solution on the frequency change rate measurement window optimization model to obtain and output an optimal frequency change rate measurement window.
2. The method for determining a power grid frequency change rate measurement window according to claim 1, wherein determining the total system kinetic energy of the target power grid in the target operation mode comprises:
determining a target operation mode of the target power grid by taking the minimum starting number of the synchronous units, the maximum total output of the new energy units and the minimum inertia constant of the synchronous units as constraint conditions of the operation mode;
and calculating the total system kinetic energy of the target power grid in a target operation mode.
3. The method for determining a power grid frequency change rate measurement window according to claim 2, wherein the calculating the total system kinetic energy of the target power grid in the target operation mode includes:
calculating the total system kinetic energy of the target power grid in a target operation mode according to the following formula:
E sys ＝E SG +E IM +E V(VS) +E V(CS) +E Load
wherein E is sys Representing the total system kinetic energy of the target power grid in a target operation mode, E SG Representing the total kinetic energy of the synchronous units, E IM Representing the total kinetic energy of the asynchronous induction motor, E V(VS) Energy form representing virtual inertia of voltage source type, E V(CS) Energy form representing virtual inertia of current source, E Load The energy form representing the static load voltage equivalent inertia.
4. The method for determining a power grid frequency change rate measurement window according to claim 1, wherein the calculating the system frequency change rate corresponding to the corresponding target fault according to the unbalanced power and the total kinetic energy of the system to obtain the theoretical value of the system frequency change rate corresponding to each target fault includes:
calculating the system frequency change rate corresponding to each target fault according to the following formula:
in the formula, roCoF theory (i) For the system frequency change rate corresponding to the target fault i, f 0 For the nominal frequency of the system, ΔP i Unbalanced power of system after occurrence of target fault i, E sys And the total kinetic energy of the system of the target power grid in the target operation mode is obtained.
5. The method according to claim 1, wherein the target constraint condition uses a target fault occurrence time as a lower bound of the target fault occurrence period and uses a sum of the target fault occurrence time and a preset fault time factor as an upper bound of the target fault occurrence period.
6. The grid frequency rate of change measurement window determination method of claim 5, wherein said establishing a frequency rate of change measurement window optimization model from said objective function and said objective constraint comprises:
The frequency change rate measurement window optimization model is established as follows:
s.t.
T 0 ≤T≤T max
max{[max(RoCoF i,j,T (t))-min(RoCoF i,j,T (t)),t∈[t 0 ,t 0 +τ]],i∈Ω F ,j∈Ω M }≤∈
max{RoCoF i,j,T (t),t∈[t 0 ,t 0 +τ]}≤RoCoF theory (i)
wherein T represents the frequency change rate measurement window value, T 0 Representing the initial value of a preset window, T max Representing the maximum allowable value of a preset window omega F Representing the target fault set, Ω M Representing a set of target frequency rate of change measurement points, rocofs i,j,T (t) represents the target fault i pairThe corresponding frequency change rate measurement point j is a frequency change rate measurement value T at time T in the occurrence period of the target fault when the frequency change rate measurement window value is T, and T 0 For the occurrence time of the target fault, tau is a preset fault time factor, rocofs theory (i) And measuring a deviation constant value for the system frequency change rate corresponding to the target fault i and the frequency change rate.
7. The method for determining a frequency change rate measurement window of a power grid according to claim 6, wherein the optimizing the frequency change rate measurement window optimization model includes:
and solving the frequency change rate measurement window optimization model by gradually increasing the frequency change rate measurement window value.
8. A grid frequency rate of change measurement window determination apparatus, comprising:
the first determining module is used for determining the total system kinetic energy of the target power grid in a target running mode; the target operation mode is an operation mode corresponding to the minimum system rotational inertia of the target power grid;
The second determining module is used for determining a typical N-2 fault of the target power grid in the target operation mode to obtain a target fault set;
a third determining module, configured to determine an unbalanced power of the system after each target fault occurs; the target faults are faults in the target fault set;
the calculation module is used for calculating the system frequency change rate corresponding to the corresponding target faults according to the unbalanced power and the total kinetic energy of the system to obtain the system frequency change rate theoretical value corresponding to each target fault;
the fourth determining module is used for determining each frequency change rate measuring point spread over the key area of the target power grid to obtain a target frequency change rate measuring point set;
the first construction module is used for constructing target constraint conditions aiming at measurement parameter values of different target frequency change rate measurement points corresponding to each target fault under different frequency change rate measurement window values; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the measurement parameter value comprises a frequency change rate measurement value maximum value and a frequency change rate measurement deviation in a target fault occurrence period; the target frequency change rate measurement points are frequency change rate measurement points in the target frequency change rate measurement point set; the target constraint includes: the frequency change rate measurement window value is between a preset window initial value and a preset window maximum allowable value, the maximum value of all frequency change rate measurement deviations is not more than a frequency change rate measurement deviation fixed value, and the maximum value of each frequency change rate measurement value is not more than a corresponding system frequency change rate theoretical value;
The second construction module is used for taking the minimum frequency change rate measurement window as an objective function, and establishing a frequency change rate measurement window optimization model according to the objective function and the objective constraint condition;
and the solving module is used for carrying out optimization solving on the frequency change rate measurement window optimization model to obtain and output an optimal frequency change rate measurement window.
9. A grid frequency rate of change measurement window determination apparatus, comprising:
a memory for storing instructions; wherein the instructions are for implementing a grid frequency change rate measurement window determination method as claimed in any one of claims 1 to 7;
and the processor is used for executing the instructions in the memory.
10. A computer readable storage medium, characterized in that the computer readable storage medium has stored thereon a computer program which, when executed by a processor, implements the grid frequency change rate measurement window determination method according to any of claims 1-7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211625125.4A CN116008656A (en) | 2022-12-16 | 2022-12-16 | Power grid frequency change rate measurement window determining method, device and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211625125.4A CN116008656A (en) | 2022-12-16 | 2022-12-16 | Power grid frequency change rate measurement window determining method, device and storage medium |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116008656A true CN116008656A (en) | 2023-04-25 |
Family
ID=86037361
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211625125.4A Pending CN116008656A (en) | 2022-12-16 | 2022-12-16 | Power grid frequency change rate measurement window determining method, device and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116008656A (en) |
-
2022
- 2022-12-16 CN CN202211625125.4A patent/CN116008656A/en active Pending
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Wall et al. | Estimation of generator inertia available during a disturbance | |
| US11208984B2 (en) | Method and apparatus for controlling pitch of wind turbine in extreme turbulence wind conditions | |
| CN111444596B (en) | Method and device for detecting frequency safety after power disturbance of power system | |
| CN112434936A (en) | Power system inertia security domain evaluation method and system, electronic equipment and readable storage medium | |
| CN109005152B (en) | Method and system for evaluating attack hazards of source-network-load system | |
| CN116008656A (en) | Power grid frequency change rate measurement window determining method, device and storage medium | |
| CN117394436A (en) | Photovoltaic bearing capacity assessment method and system based on random scene analysis | |
| CN112421650B (en) | Method and system for evaluating power grid inertia | |
| CN116596333A (en) | Energy storage power station grade determining method and device based on index evaluation system | |
| CN113471972B (en) | Method and system for monitoring inertia of power system | |
| CN112366733B (en) | Method, device, system and medium for optimizing PID (proportion integration differentiation) parameters of speed regulator of hydroelectric generating set | |
| CN110707729A (en) | Method and device for analyzing multi-region primary frequency modulation condition and storage medium | |
| CN107122254A (en) | A kind of computer repairs control method and system, restorative procedure and system | |
| JP4001994B2 (en) | Assessing accident stability of power systems | |
| CN114123236A (en) | Frequency response control method and device for power electronic equipment | |
| CN114977191A (en) | Static voltage stability analysis method and device and computer storage medium | |
| CN113839416A (en) | Wind power maximum access proportion determining method and device based on frequency response model | |
| CN109088413B (en) | Dynamic reactive power planning and site selection method and device for direct-current receiving-end power grid | |
| CN116760067A (en) | Safety inertia determination method and device for power grid system and computer equipment | |
| JP2022016046A (en) | Electric power system stabilization system | |
| CN110943480A (en) | Power system frequency modulation method and device, computer equipment and storage medium | |
| CN110417041B (en) | Method and system for supporting direct-current high-power failure | |
| CN112737049B (en) | Battery equalization control method and device and battery module | |
| CN117335493B (en) | Regional power grid-oriented inertia evaluation method for each control region | |
| CN111007396B (en) | Method and device for detecting accuracy of PSS2A/B model and storage medium |
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 |