CN116345541A - Wind power penetration power limit calculation method, device and storage medium - Google Patents

Wind power penetration power limit calculation method, device and storage medium Download PDF

Info

Publication number
CN116345541A
CN116345541A CN202310147603.3A CN202310147603A CN116345541A CN 116345541 A CN116345541 A CN 116345541A CN 202310147603 A CN202310147603 A CN 202310147603A CN 116345541 A CN116345541 A CN 116345541A
Authority
CN
China
Prior art keywords
power
offshore oilfield
power grid
oilfield group
fan
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
CN202310147603.3A
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.)
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Original Assignee
Beijing Research Center of CNOOC China Ltd
CNOOC China Ltd
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Beijing Research Center of CNOOC China Ltd, CNOOC China Ltd filed Critical Beijing Research Center of CNOOC China Ltd
Priority to CN202310147603.3A priority Critical patent/CN116345541A/en
Publication of CN116345541A publication Critical patent/CN116345541A/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
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/381—Dispersed generators
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/15—Correlation function computation including computation of convolution operations
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
    • G06F17/10—Complex mathematical operations
    • G06F17/18—Complex mathematical operations for evaluating statistical data, e.g. average values, frequency distributions, probability functions, regression analysis
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/24—Arrangements for preventing or reducing oscillations of power in networks
    • H02J3/241—The oscillation concerning frequency
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J3/00—Circuit arrangements for ac mains or ac distribution networks
    • H02J3/38—Arrangements for parallely feeding a single network by two or more generators, converters or transformers
    • H02J3/46—Controlling of the sharing of output between the generators, converters, or transformers
    • H02J3/466—Scheduling the operation of the generators, e.g. connecting or disconnecting generators to meet a given demand
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2203/00—Indexing scheme relating to details of circuit arrangements for AC mains or AC distribution networks
    • H02J2203/20—Simulating, e g planning, reliability check, modelling or computer assisted design [CAD]
    • H—ELECTRICITY
    • H02—GENERATION; CONVERSION OR DISTRIBUTION OF ELECTRIC POWER
    • H02J—CIRCUIT ARRANGEMENTS OR SYSTEMS FOR SUPPLYING OR DISTRIBUTING ELECTRIC POWER; SYSTEMS FOR STORING ELECTRIC ENERGY
    • H02J2300/00—Systems for supplying or distributing electric power characterised by decentralized, dispersed, or local generation
    • H02J2300/20—The dispersed energy generation being of renewable origin
    • H02J2300/28—The renewable source being wind energy
    • 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
    • Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
    • Y02E10/00—Energy generation through renewable energy sources
    • Y02E10/70—Wind energy
    • Y02E10/76—Power conversion electric or electronic aspects

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • General Physics & Mathematics (AREA)
  • Data Mining & Analysis (AREA)
  • Mathematical Physics (AREA)
  • Computational Mathematics (AREA)
  • Mathematical Analysis (AREA)
  • Mathematical Optimization (AREA)
  • Pure & Applied Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Power Engineering (AREA)
  • Databases & Information Systems (AREA)
  • Software Systems (AREA)
  • General Engineering & Computer Science (AREA)
  • Algebra (AREA)
  • Computing Systems (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Bioinformatics & Cheminformatics (AREA)
  • Bioinformatics & Computational Biology (AREA)
  • Evolutionary Biology (AREA)
  • Operations Research (AREA)
  • Probability & Statistics with Applications (AREA)
  • Supply And Distribution Of Alternating Current (AREA)

Abstract

The invention relates to a wind power penetration power limit calculation method, a wind power penetration power limit calculation device and a storage medium, wherein the wind power penetration power limit calculation method comprises the following steps: calculating to obtain a wind power penetration power limit under transient frequency safety constraint according to a preset first disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function; calculating to obtain a wind power penetration power limit under a steady-state frequency safety constraint according to a preset second disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function; and selecting a smaller value of wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as the wind power transmission power limit of the power grid. The method can comprehensively consider the influence of low inertia and fluctuation of the fan on the transient state and steady state frequency safety of the power grid, quickly solve the wind power transmission power limit of the offshore oilfield group power grid, and better meet the safety and stability requirements of the offshore oilfield group power grid.

Description

Wind power penetration power limit calculation method, device and storage medium
Technical Field
The invention relates to the field of wind power grid connection control, in particular to a wind power penetration power limit calculation method, a wind power penetration power limit calculation device and a storage medium.
Background
In the "two carbon" background, offshore oilfield cluster grids will gradually increase wind power access capacity. However, the low inertia and volatility of wind power are extremely easy to threaten the safe and stable operation of the offshore oilfield group power grid, and the wind power penetration power limit of the offshore oilfield group power grid is necessary to be researched and determined. The wind power penetration power limit is defined as the ratio of the maximum wind power capacity accessible by the power grid to the maximum load of the power grid under the condition of meeting safe and stable operation.
Because the offshore oilfield group power grid has small power generation capacity and is often operated in an island state, the influence of wind power access on the frequency stability of the offshore oilfield group power grid is very remarkable. In particular, the frequency of the small power grid is easy to change greatly due to the fluctuation of the output of a fan or the faults of a submarine cable, and the frequency of the small power grid is also challenged by the transient frequency safety problem besides the steady-state frequency safety problem. However, the wind power penetration power limit calculation method comprehensively considering transient state and steady state frequency safety constraint is aimed at a land power grid, the calculation model takes a reheat steam turbine as a prime mover-speed regulator, and is not applicable to an offshore oilfield group power grid taking a gas turbine as a main body; the method for calculating the wind power penetration power limit of the offshore power grid by considering the response speed of the gas turbine mainly takes steady-state frequency safety as constraint according to the expected wind power output or a preset scheduling period.
Disclosure of Invention
Aiming at the problems, the invention aims to provide a wind power transmission power limit calculation method, a wind power transmission power limit calculation device and a storage medium, which make up the defect of considering transient problems in wind power transmission power limit calculation of an offshore power grid and better meet the requirements of safe and stable operation of offshore oilfield group power grids.
In order to achieve the above purpose, the present invention adopts the following technical scheme:
a wind power penetration power limit calculation method is applied to an offshore oilfield group power grid and comprises the following steps:
calculating to obtain a wind power penetration power limit under transient frequency safety constraint according to a preset first disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function;
calculating to obtain a wind power penetration power limit under a steady-state frequency safety constraint according to a preset second disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function;
and selecting a smaller value of wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as the wind power transmission power limit of the offshore oilfield group power grid.
Further, the equivalent frequency response transfer function of the offshore oilfield group power grid is expressed as follows:
Figure BDA0004089650220000021
wherein DeltaP d For disturbance power ΔP d With the increase of the power generation side output or the decrease of the load side power positive, Δf is any disturbance power ΔP d Corresponding frequency response results, wherein H is an inertial time constant of the offshore oilfield group power grid, D is an offshore oilfield group power grid damping coefficient, R is an offshore oilfield group power grid difference adjustment coefficient, X is a gas turbine feeding control parameter, T F Adjusting time constant, T, for gas turbine valve position CD Is the gas turbine compressor exhaust time constant.
Further, the constructing the fan output fluctuation probability density function includes:
calculating fan output fluctuation data delta P w (t) the calculation formula is:
Figure BDA0004089650220000022
wherein P is w (t) is the total output of the wind power plant at the moment t, and the time span of the output data of the wind power plant is set to be 0-t 0 The value range of t is [0, t 0 -Δt]Δt is the fluctuation time interval in the process of calculating fluctuation data, and n is the number of fan access stations; fan output fluctuation data Δp w The value range of (C) is [ -1,1]；
The output fluctuation data delta P of the fan is processed w Dividing the interval from-1 to 1 into N fluctuation intervals, and dividing the data quantity in each fluctuation interval by the fan output fluctuation data delta P w Obtaining probability density of each fluctuation interval and establishing windOutput fluctuation data delta P w Probability density histogram of (a);
fitting the probability density histogram into a continuous curve to obtain a fan output fluctuation probability density function f p (ΔP w )。
Further, the first disturbance power expression corresponding to the most serious fault of the offshore oilfield group power grid includes:
if the most serious fault of the offshore oilfield group power grid is that the power grid is cut off in an emergency due to the fault of the sending-out system, when the offshore oilfield group power grid is cut off in an emergency due to the fault of the sending-out system, the first disturbance power delta P d1 The expression of (2) is:
Figure BDA0004089650220000031
wherein n is the number of fan access stations, P n Rated capacity of single fan, P B For the active power base number of the offshore oilfield group power grid, the calculation formula is as follows:
Figure BDA0004089650220000032
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure BDA0004089650220000033
rated active sum of all generator sets of offshore oilfield group power grid is +.>
Figure BDA0004089650220000034
Is the sum of the loads of all platforms of the offshore oilfield group power grid, and k is +.>
Figure BDA0004089650220000035
And->
Figure BDA0004089650220000036
Is a ratio of (2).
Further, the establishing a second disturbance power expression according to the fan output fluctuation probability density function includes:
calculating different fan output fluctuation data delta P according to the fan output fluctuation probability density function w Corresponding second disturbance power DeltaP d2 The calculation formula is as follows:
Figure BDA0004089650220000037
wherein n is the number of fan access stations, P n Rated capacity of single fan, P B For the active power base number of the offshore oilfield group power grid, the calculation formula is as follows:
Figure BDA0004089650220000038
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure BDA0004089650220000039
rated active sum of all generator sets of offshore oilfield group power grid is +.>
Figure BDA00040896502200000310
Is the sum of the loads of all platforms of the offshore oilfield group power grid, and k is +.>
Figure BDA00040896502200000311
And->
Figure BDA00040896502200000312
Is a ratio of (2).
Further, the calculating the wind power penetration power limit under the transient frequency safety constraint according to the first disturbance power expression and the offshore oilfield group power grid equivalent frequency response transfer function includes:
assuming that the number of fan access stations is 1, calculating and inputting first disturbance power delta P d1 The frequency response delta f (t) of the offshore oilfield group power grid;
calculating the maximum frequency change of the offshore oilfield group power grid, and judging whether the maximum frequency response change of the offshore oilfield group power grid meets the transient frequency safety constraint max|delta f (t) |is less than or equal to the transient frequency safety constraint max|delta f (t) |Δf T Wherein Δf T As the maximum variation allowable value of transient frequency, deltaf T The value is 3Hz;
if the offshore oilfield group power grid meets the transient frequency safety constraint, adding 1 to the number of fan access stations, and calculating the corresponding input first disturbance power delta P again d1 Judging whether the transient frequency safety constraint is met or not again according to the frequency response delta f (t) of the offshore oilfield group power grid, and stopping circulation until the offshore oilfield group power grid does not meet the transient frequency safety constraint;
the number of fan access stations at the end of the cycle is reduced by 1 to be recorded as the final number of fan access stations, and the formula is passed
Figure BDA0004089650220000041
And calculating to obtain the wind power penetration power limit under the transient frequency safety constraint.
Further, the calculating the wind power penetration power limit under the steady-state frequency safety constraint according to the second disturbance power expression and the offshore oilfield group power grid equivalent frequency response transfer function includes:
assuming that the number of fan access stations is 1, calculating different fan output fluctuation data delta P according to a fan output fluctuation probability density function w Corresponding to the input disturbance power delta P d2 The frequency response delta f (t) of the offshore oilfield group power grid;
according to the steady-state frequency deviation requirement max|Δf (t) |is less than or equal to Δf of the offshore oilfield group power grid S Determining maximum allowable fan output fluctuation delta P of offshore oilfield group power grid w max Wherein Δf S As the maximum variation allowable value of transient frequency, deltaf S The value is 0.5Hz;
calculating accumulated probability Pr (devia) of allowable fan output fluctuation of the offshore oilfield group power grid, judging whether Pr (devia) meets steady-state frequency safety constraint Pr (devia) is more than or equal to beta, and adopting a calculation formula:
Figure BDA0004089650220000042
wherein, beta is the minimum cumulative probability of the allowable fan output fluctuation of the offshore oilfield group power grid, and the value of beta is 0.8;
if the offshore oilfield group power grid meets the steady-state frequency safety constraint, adding 1 to the number of fan access stations, and calculating delta P again d2 And delta f (t) to obtain the maximum allowable fan output fluctuation delta P of the offshore oilfield group power grid w max And accumulating the probability Pr (devia) and re-judging whether the accumulated probability Pr (devia) meets the steady-state frequency safety constraint or not, and stopping the loop until the accumulated probability Pr (devia) does not meet the steady-state frequency safety constraint;
the number of fan access stations at the end of the cycle is reduced by 1 to be recorded as the final number of fan access stations, and the formula is passed
Figure BDA0004089650220000051
And calculating to obtain the wind power penetration power limit under the steady-state frequency safety constraint.
To achieve the above objective, another aspect of the present application further provides a wind power penetration power limit calculating device, applied to an offshore oilfield group power grid, including:
the first calculation module is used for calculating the wind power penetration power limit under the transient frequency safety constraint according to a preset first disturbance power expression and the equivalent frequency response transfer function of the offshore oilfield group power grid;
the second calculation module is used for calculating the wind power penetration power limit under the steady-state frequency safety constraint according to a preset second disturbance power expression and the equivalent frequency response transfer function of the offshore oilfield group power grid;
and the comparison module is used for selecting smaller values in wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as wind power transmission power limits of the offshore oilfield group power grid.
To achieve the above object, another aspect of the present application further provides a computer readable storage medium, including a computer program for implementing the above method.
To achieve the above object, another aspect of the present application further provides a computing device, including a memory and a processor, where the memory is configured to store a computer program, and the computer program is configured to implement the above method when executed by the processor.
Due to the adoption of the technical scheme, the invention has the following advantages:
according to the method and device for calculating the wind power transmission power limit of the offshore oilfield group power grid and the storage medium, provided by the invention, the influence of low inertia and fluctuation of a fan on transient state and steady state frequency safety of the offshore oilfield group power grid can be comprehensively considered, the wind power transmission power limit of the offshore oilfield group power grid can be rapidly solved, and the obtained result can better meet the safe and stable requirements of the offshore oilfield group power grid.
Drawings
Various other advantages and benefits will become apparent to those of ordinary skill in the art upon reading the following detailed description of the preferred embodiments. The drawings are only for purposes of illustrating the preferred embodiments and are not to be construed as limiting the invention. Like parts are designated with like reference numerals throughout the drawings. In the drawings:
fig. 1 is a flowchart of a method for calculating a power limit of a wind turbine in an embodiment of the invention.
Detailed Description
Exemplary embodiments of the present invention will be described in more detail below with reference to the accompanying drawings. While exemplary embodiments of the present invention are shown in the drawings, it should be understood that the present invention may be embodied in various forms and should not be limited to the embodiments set forth herein. Rather, these embodiments are provided so that this disclosure will be thorough and complete, and will fully convey the scope of the invention to those skilled in the art.
Aiming at the characteristics of small power generation capacity of the offshore oilfield group power grid and easy occurrence of large fluctuation of frequency after wind power access, the application provides a wind power penetration power limit simplified calculation method which takes the maximum frequency change under the most serious fault as a transient frequency safety constraint index and takes the cumulative probability density of the allowable fan output fluctuation as a steady frequency safety constraint index.
Referring to fig. 1, the present embodiment provides a wind power penetration power limit calculation method applied to an offshore oilfield group power grid, including:
s101: constructing an equivalent frequency response transfer function and a fan output fluctuation probability density function of an offshore oilfield group power grid;
in this embodiment, the offshore oilfield group power grid mainly uses a single-shaft gas turbine (SG) gas turbine as a prime mover-speed regulator, and for this feature, an offshore oilfield group power grid equivalent frequency response transfer function is established, where the offshore oilfield group power grid equivalent frequency response transfer function is expressed as:
Figure BDA0004089650220000071
wherein DeltaP d For disturbance power ΔP d With the increase of the power generation side output or the decrease of the load side power positive, Δf is any disturbance power ΔP d Corresponding frequency response results, wherein H is an inertial time constant of the offshore oilfield group power grid, D is an offshore oilfield group power grid damping coefficient, R is an offshore oilfield group power grid difference adjustment coefficient, X is a gas turbine feeding control parameter, T F Adjusting time constant, T, for gas turbine valve position CD Is the gas turbine compressor exhaust time constant.
In this embodiment, since the fan output fluctuation data of each region is not easy to directly obtain, it needs to be separated from the actual wind speed and the fan output data, in order to definitely threat the fan output fluctuation situation of the offshore oilfield group power grid steady-state frequency safety and the mathematical representation thereof, a first-order differential transformation method is used to determine the fluctuation time sampling interval Δt, and then it needs to calculate the fan output fluctuation data Δp w (t) the calculation formula is:
Figure BDA0004089650220000072
wherein P is w (t) is the total output of the wind power plant at the moment t, and the time span of the output data of the wind power plant is set to be 0-t 0 The value range of t is [0, t 0 -Δt]Δt is the fluctuation time interval in the process of calculating fluctuation data, and n is the number of fan access stations; fan output fluctuation data Δp w The value range of (C) is [ -1,1]；
The output fluctuation data delta P of the fan is processed w Dividing the interval from-1 to 1 into N fluctuation intervals, and dividing the data quantity in each fluctuation interval by the fan output fluctuation data delta P w Obtaining probability density of each fluctuation interval and establishing fan output fluctuation data delta P w Probability density histogram of (a);
fitting the probability density histogram into a continuous curve to obtain a fan output fluctuation probability density function f p (ΔP w )。
S102: calculating to obtain a wind power penetration power limit under transient frequency safety constraint according to the first disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function; calculating to obtain a wind power penetration power limit under a steady-state frequency safety constraint according to the second disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function;
in this embodiment, the most serious fault and the mathematical representation thereof are used for definitely threatening the transient frequency safety of the offshore oilfield group power grid. The largest gas turbine generator set cut-out, wind farm cut-out, etc. can generally be selected as the most severe fault. Assuming that the most serious fault is the whole emergency cut-out of the offshore wind farm due to the fault of the delivery system, describing the most serious fault and the disturbance power expression thereof by using a step function due to the short cutting-out action time, the first disturbance power delta P can be obtained d1 The expression of (2) is:
Figure BDA0004089650220000081
wherein n is the number of fan access stations, P n Rated capacity of single fan, P B For offshore oilfield group powerThe net active power base number is calculated by the following formula:
Figure BDA0004089650220000082
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure BDA0004089650220000083
rated active sum of all generator sets of offshore oilfield group power grid is +.>
Figure BDA0004089650220000084
Is the sum of the loads of all platforms of the offshore oilfield group power grid, and k is +.>
Figure BDA0004089650220000085
And->
Figure BDA0004089650220000086
Is a ratio of (2).
In the present embodiment, different fan output fluctuation data Δp is calculated according to the fan output fluctuation probability density function w Corresponding second disturbance power DeltaP d2 The calculation formula is as follows:
Figure BDA0004089650220000087
wherein n is the number of fan access stations, P n Rated capacity of single fan, P B For the active power base number of the offshore oilfield group power grid, the calculation formula is as follows:
Figure BDA0004089650220000088
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure BDA0004089650220000089
rated active sum of all generator sets of offshore oilfield group power grid is +.>
Figure BDA00040896502200000810
Is the sum of the loads of all platforms of the offshore oilfield group power grid, and k is +.>
Figure BDA00040896502200000811
And->
Figure BDA00040896502200000812
Is a ratio of (2).
S103: calculating to obtain a wind power penetration power limit under transient frequency safety constraint according to the first disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function; calculating to obtain a wind power penetration power limit under a steady-state frequency safety constraint according to the second disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function;
in this embodiment, based on the most serious fault and the disturbance power expression thereof determined in step S102, the wind power penetration power limit under the transient frequency safety constraint is calculated with the maximum frequency variation under the most serious fault as the constraint.
Assuming that the number of fan access stations is 1, calculating and inputting first disturbance power delta P d1 The frequency response delta f (t) of the offshore oilfield group power grid;
calculating the maximum frequency change of the offshore oilfield group power grid, and judging whether the maximum frequency response change of the offshore oilfield group power grid meets the transient frequency safety constraint max|delta f (t) |delta f or not T Wherein Δf T As the maximum variation allowable value of transient frequency, deltaf T The value is 3Hz, namely the lowest frequency of the offshore oilfield group power grid representing the transient process is not lower than 47Hz and the highest frequency is not higher than 53Hz.
If the offshore oilfield group power grid meets the transient frequency safety constraint, adding 1 to the number of fan access stations, and calculating the corresponding input first disturbance power delta P again d1 Judging whether the transient frequency safety constraint is met or not again according to the frequency response delta f (t) of the offshore oilfield group power grid, and stopping circulation until the offshore oilfield group power grid does not meet the transient frequency safety constraint;
the number of fan access stations at the end of the cycle is reduced by 1 to be recorded as the final number of fan access stations, and the formula is passed
Figure BDA0004089650220000091
And calculating to obtain the wind power penetration power limit under the transient frequency safety constraint.
In this embodiment, based on the fan output fluctuation probability density function and the corresponding disturbance power expression determined in step S102, the wind power penetration power limit under the steady-state frequency security constraint is calculated with the cumulative probability of allowable fan output fluctuation as the constraint.
Assuming that the number of fan access stations is 1, calculating different fan output fluctuation data delta P according to a fan output fluctuation probability density function w Corresponding to the input disturbance power delta P d2 The frequency response delta f (t) of the offshore oilfield group power grid;
according to the steady-state frequency deviation requirement max|Δf (t) |is less than or equal to Δf of the offshore oilfield group power grid S Determining maximum allowable fan output fluctuation delta P of offshore oilfield group power grid w max Wherein Δf S As the maximum variation allowable value of transient frequency, deltaf S The value is 0.5Hz;
calculating accumulated probability Pr (devia) of allowable fan output fluctuation of the offshore oilfield group power grid, judging whether Pr (devia) meets steady-state frequency safety constraint Pr (devia) is more than or equal to beta, and adopting a calculation formula:
Figure BDA0004089650220000101
wherein, beta is the minimum cumulative probability of the allowable fan output fluctuation of the offshore oilfield group power grid, and the value of beta is 0.8;
if the offshore oilfield group power grid meets the steady-state frequency safety constraint, adding 1 to the number of fan access stations, and calculating delta P again d2 And delta f (t) to obtain the maximum allowable fan output fluctuation delta P of the offshore oilfield group power grid w max And accumulating the probability Pr (devia) and re-judging whether the accumulated probability Pr (devia) meets the steady-state frequency safety constraint until the accumulated probability Pr (devia) does not meet the steady-state frequency safety constraintTerminating the cycle at bundle time;
the number of fan access stations at the end of the cycle is reduced by 1 to be recorded as the final number of fan access stations, and the formula is passed
Figure BDA0004089650220000102
And calculating to obtain the wind power penetration power limit under the steady-state frequency safety constraint.
S104: and selecting a smaller value of wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as the wind power transmission power limit of the offshore oilfield group power grid.
In this embodiment, in order to make the final calculation result meet each frequency safety constraint, comparing the wind power penetration power limit result under the transient frequency safety constraint with the wind power penetration power limit result under the steady frequency safety constraint, and selecting the lowest value of the two, to determine the final wind power penetration power limit result of the offshore oilfield group power grid.
The application also provides a wind power penetration power limit calculation device, is applied to offshore oilfield group electric wire netting, includes:
the first calculation module is used for calculating the wind power penetration power limit under the transient frequency safety constraint according to a preset first disturbance power expression and the equivalent frequency response transfer function of the offshore oilfield group power grid;
the second calculation module is used for calculating the wind power penetration power limit under the steady-state frequency safety constraint according to a preset second disturbance power expression and the equivalent frequency response transfer function of the offshore oilfield group power grid;
and the comparison module is used for selecting smaller values in wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as wind power transmission power limits of the offshore oilfield group power grid.
The apparatus provided in this embodiment is used to execute the foregoing method embodiment, and specific flow and details refer to the foregoing embodiment and are not repeated herein.
The present application also provides a computer readable storage medium comprising a computer program for implementing the above method.
The present application also provides a computing device comprising a memory and a processor, the memory for storing a computer program which, when executed by the processor, implements the method described above.
From the above, the technical scheme provided by the application comprehensively considers transient and steady-state frequency safety constraints, overcomes the defect of transient problem consideration of wind power transmission power limit calculation of the offshore power grid, can comprehensively consider the influence of low inertia and fluctuation of a fan on transient and steady-state frequency safety of the offshore oilfield group power grid and rapidly solve the wind power transmission power limit of the power grid, is suitable for offshore oilfield group power grids and other small power grids with frequency which is easy to change greatly, and can better meet the safe and steady operation requirements of the power grid. The provided maximum frequency change constraint index under the most serious faults and the cumulative probability constraint index allowing the fan output to fluctuate can effectively embody the transient state and steady state safety characteristics of the power grid and simplify the limit solving flow. And the equivalent frequency response transfer function of the power grid is established according to the characteristics of the prime mover and the speed regulator of the offshore oilfield group power grid, so that the wind power penetration power limit can be rapidly solved. The wind power penetration power limit calculation method can be expressed as wind power maximum access capacity, wind power bearing capacity of a power grid and the like; in addition, other frequency response calculation methods can be adopted in various wind power penetration power limit calculation methods, so that the power grid frequency response calculation under each disturbance power in the transient state and steady state frequency safety constraint indexes provided by the invention is realized, and the same function is realized typically like electromagnetic/electromechanical transient state simulation or equivalent frequency response calculation; the method provided by the application is not only suitable for offshore oilfield group power grids, but also can be applied to offshore oilfield group power grids and other small power grids taking single-shaft and split-shaft gas turbines as main bodies.
Finally, it should be noted that: the above embodiments are only for illustrating the technical solution of the present invention, and are not limiting; 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. The wind power penetration power limit calculation method is applied to an offshore oilfield group power grid and is characterized by comprising the following steps of:
calculating to obtain a wind power penetration power limit under transient frequency safety constraint according to a preset first disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function;
calculating to obtain a wind power penetration power limit under a steady-state frequency safety constraint according to a preset second disturbance power expression and an offshore oilfield group power grid equivalent frequency response transfer function;
and selecting a smaller value of wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as the wind power transmission power limit of the offshore oilfield group power grid.
2. The method of claim 1, wherein the offshore oilfield group grid equivalent frequency response transfer function is expressed as:
Figure FDA0004089650210000011
wherein DeltaP d For disturbance power ΔP d With the increase of the power generation side output or the decrease of the load side power positive, Δf is any disturbance power ΔP d Corresponding frequency response results, wherein H is an inertial time constant of the offshore oilfield group power grid, D is an offshore oilfield group power grid damping coefficient, R is an offshore oilfield group power grid difference adjustment coefficient, X is a gas turbine feeding control parameter, T F Adjusting time constant, T, for gas turbine valve position CD Is the gas turbine compressor exhaust time constant.
3. The method of claim 1, wherein the fan output fluctuation probability density function comprises:
calculating fan output fluctuation data delta P w (t) the calculation formula is:
Figure FDA0004089650210000012
wherein P is w (t) is the total output of the wind power plant at the moment t, and the time span of the output data of the wind power plant is set to be 0-t 0 The value range of t is [0, t 0 -Δt]Δt is the fluctuation time interval in the process of calculating fluctuation data, and n is the number of fan access stations; fan output fluctuation data Δp w The value range of (C) is [ -1,1]；
The output fluctuation data delta P of the fan is processed w Dividing the interval from-1 to 1 into N fluctuation intervals, and dividing the data quantity in each fluctuation interval by the fan output fluctuation data delta P w Obtaining probability density of each fluctuation interval and establishing fan output fluctuation data delta P w Probability density histogram of (a);
fitting the probability density histogram into a continuous curve to obtain a fan output fluctuation probability density function f p (ΔP w )。
4. The method of claim 1, wherein the first disturbance power expression comprises:
if the most serious fault of the offshore oilfield group power grid is that the power grid is cut off in an emergency due to the fault of the sending-out system, when the offshore oilfield group power grid is cut off in an emergency due to the fault of the sending-out system, the first disturbance power delta P d1 The expression of (2) is:
Figure FDA0004089650210000021
wherein n is the number of fan access stations, P n Rated capacity of single fan, P B For the active power base number of the offshore oilfield group power grid, the calculation formula is as follows:
Figure FDA0004089650210000022
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure FDA0004089650210000023
rated active sum of all generator sets of offshore oilfield group power grid is +.>
Figure FDA0004089650210000024
Is the sum of the loads of all platforms of the offshore oilfield group power grid, and k is +.>
Figure FDA0004089650210000025
And->
Figure FDA0004089650210000026
Is a ratio of (2).
5. A method according to claim 3, wherein the second disturbance power expression comprises:
calculating different fan output fluctuation data delta P according to the fan output fluctuation probability density function w Corresponding second disturbance power DeltaP d2 The calculation formula is as follows:
Figure FDA0004089650210000027
wherein n is the number of fan access stations, P n Rated capacity of single fan, P B For the active power base number of the offshore oilfield group power grid, the calculation formula is as follows:
Figure FDA0004089650210000031
wherein, the liquid crystal display device comprises a liquid crystal display device,
Figure FDA0004089650210000032
rated active sum of all generator sets of offshore oilfield group power grid is +.>
Figure FDA0004089650210000033
Is the sum of the loads of all platforms of the offshore oilfield group power grid, and k is +.>
Figure FDA0004089650210000034
And->
Figure FDA0004089650210000035
Is a ratio of (2).
6. The method of claim 4, wherein calculating the wind power penetration power limit under the transient frequency security constraint according to the preset first disturbance power expression and the offshore oilfield group grid equivalent frequency response transfer function comprises:
assuming that the number of fan access stations is 1, calculating and inputting first disturbance power delta P d1 The frequency response delta f (t) of the offshore oilfield group power grid;
calculating the maximum frequency change of the offshore oilfield group power grid, and judging whether the maximum frequency response change of the offshore oilfield group power grid meets the transient frequency safety constraint max|delta f (t) |delta f or not T Wherein Δf T As the maximum variation allowable value of transient frequency, deltaf T The value is 3Hz;
if the offshore oilfield group power grid meets the transient frequency safety constraint, adding 1 to the number of fan access stations, and calculating the corresponding input first disturbance power delta P again d1 Judging whether the transient frequency safety constraint is met or not again according to the frequency response delta f (t) of the offshore oilfield group power grid, and stopping circulation until the offshore oilfield group power grid does not meet the transient frequency safety constraint;
the number of fan access stations at the end of the cycle is reduced by 1 to be recorded as the final number of fan access stations, and the formula is passed
Figure FDA0004089650210000036
And calculating to obtain the wind power penetration power limit under the transient frequency safety constraint.
7. The method of claim 5, wherein calculating the wind power penetration power limit under the steady-state frequency security constraint from the preset second disturbance power expression and the offshore oilfield group grid equivalent frequency response transfer function comprises:
assuming that the number of fan access stations is 1, calculating different fan output fluctuation data delta P according to a fan output fluctuation probability density function w Corresponding to the input disturbance power delta P d2 The frequency response delta f (t) of the offshore oilfield group power grid;
according to the steady-state frequency deviation requirement max|Δf (t) |is less than or equal to Δf of the offshore oilfield group power grid S Determining maximum allowable fan output fluctuation delta P of offshore oilfield group power grid w max Wherein Δf S As the maximum variation allowable value of transient frequency, deltaf S The value is 0.5Hz;
calculating accumulated probability Pr (devia) of allowable fan output fluctuation of the offshore oilfield group power grid, judging whether Pr (devia) meets steady-state frequency safety constraint Pr (devia) is more than or equal to beta, and adopting a calculation formula:
Figure FDA0004089650210000041
wherein, beta is the minimum cumulative probability of the allowable fan output fluctuation of the offshore oilfield group power grid, and the value of beta is 0.8;
if the offshore oilfield group power grid meets the steady-state frequency safety constraint, adding 1 to the number of fan access stations, and calculating delta P again d2 And delta f (t) to obtain the maximum allowable fan output fluctuation delta P of the offshore oilfield group power grid w max And accumulating the probability Pr (devia) and re-judging whether the accumulated probability Pr (devia) satisfies the steady-state frequency safety constraint until the accumulated probability Pr (devia) does not satisfy the steady-state frequency safety constraintTerminating the cycle while constrained;
the number of fan access stations at the end of the cycle is reduced by 1 to be recorded as the final number of fan access stations, and the formula is passed
Figure FDA0004089650210000042
And calculating to obtain the wind power penetration power limit under the steady-state frequency safety constraint.
8. A wind power penetration power limit calculation device applied to an offshore oilfield group power grid, comprising:
the first calculation module is used for calculating the wind power penetration power limit under the transient frequency safety constraint according to a preset first disturbance power expression and the equivalent frequency response transfer function of the offshore oilfield group power grid;
the second calculation module is used for calculating the wind power penetration power limit under the steady-state frequency safety constraint according to a preset second disturbance power expression and the equivalent frequency response transfer function of the offshore oilfield group power grid;
and the comparison module is used for selecting smaller values in wind power transmission power limits under transient frequency safety constraint and steady frequency safety constraint as wind power transmission power limits of the offshore oilfield group power grid.
9. A computer readable storage medium, characterized in that the computer readable storage medium comprises a computer program for implementing the method of any of claims 1 to 7.
10. A computing device comprising a memory and a processor, the memory for storing a computer program that, when executed by the processor, implements the method of any of claims 1 to 7.
CN202310147603.3A 2023-02-20 2023-02-20 Wind power penetration power limit calculation method, device and storage medium Pending CN116345541A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202310147603.3A CN116345541A (en) 2023-02-20 2023-02-20 Wind power penetration power limit calculation method, device and storage medium

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202310147603.3A CN116345541A (en) 2023-02-20 2023-02-20 Wind power penetration power limit calculation method, device and storage medium

Publications (1)

Publication Number Publication Date
CN116345541A true CN116345541A (en) 2023-06-27

Family

ID=86879858

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202310147603.3A Pending CN116345541A (en) 2023-02-20 2023-02-20 Wind power penetration power limit calculation method, device and storage medium

Country Status (1)

Country Link
CN (1) CN116345541A (en)

Similar Documents

Publication Publication Date Title
Ahmed et al. Modified TID controller for load frequency control of a two-area interconnected diverse-unit power system
WO2021164112A1 (en) Frequency control method and system during using wind farm as black-start power source by means of optimal configuration of energy storage
Sharma et al. Frequency excursion mitigation strategy using a novel COA optimised fuzzy controller in wind integrated power systems
Wang et al. Dynamic load frequency control for high-penetration wind power considering wind turbine fatigue load
Chaine et al. Performance of CSA optimized controllers of DFIGs and AGC to improve frequency regulation of a wind integrated hydrothermal power system
CN112366726B (en) Primary frequency modulation coefficient optimization method for thermal power generating unit and related equipment
Taher et al. Optimal gain scheduling controller design of a pitch-controlled VS-WECS using DE optimization algorithm
CN114844064A (en) Self-adaptive variable parameter frequency modulation method and device for double-fed variable-speed pumped storage unit
CN110513158B (en) Feed-forward multistage speed regulation method for steam turbine
Kazemi et al. Adaptive fractional‐order control of power system frequency in the presence of wind turbine
Meisami‐Azad et al. Anti‐windup linear parameter‐varying control of pitch actuators in wind turbines
US20230344242A1 (en) Method for automatic adjustment of power grid operation mode base on reinforcement learning
CN116345541A (en) Wind power penetration power limit calculation method, device and storage medium
Lee et al. Operation scheme for a wind farm to mitigate output power variation
Xu et al. Parameter optimization of robust non-fragile fractional order PID controller for pump turbine governing system
CN108063461B (en) Wind power-containing electric power system active power scheduling method considering small interference stability risk
CN115842376A (en) Method, device and medium for evaluating equivalent inertia trend and safety state of power system
CN115940148A (en) Minimum inertia requirement evaluation method and device, electronic equipment and storage medium
CN115207941A (en) Method for evaluating inertia level of power system
Rezaei et al. Robust fault tolerant pitch control of wind turbines
Sina et al. An optimal controller for load frequency control in multi-area deregulated power system
CN110854852B (en) Configuration method of key parameters of AGC main station in high-proportion hydropower area
CN113991647A (en) Frequency response capacity planning-oriented random production simulation method for power system
CN110233486B (en) Energy storage auxiliary fuzzy PID control method based on active prediction SOC
CN112134291A (en) Reactive power voltage regulation control method for large wind power plant

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