Offshore wind energy short-time forecast wave effect correction method and system

Abstract

The invention discloses an offshore wind energy short-time forecast wave effect correction method and system. The method comprises the following steps: S1, generating a three-dimensional wind field initial field of a wind power plant and wave spectrum characteristic parameters of an air-sea interface based on meteorological data; s2, a fan is initialized based on the three-dimensional wind field initial field of the wind power plant, and fan parameters and a fan indicator function are obtained; s3, on the basis of the fan parameters and the fan indicator function, fan grid point thrust is obtained through calculation; s4, based on the wind power plant three-dimensional wind field initial field and the wave spectrum characteristic parameters, total moving surface drag force is obtained through calculation; s5, based on the wind power plant three-dimensional wind field initial field, the fan grid point thrust and the total moving surface drag force, short-time wind energy data are obtained through calculation; s6, outputting the short-time wind energy data based on a preset format; according to the method, the influence of sea surface waves on wind energy forecasting can be considered during short-time wind energy forecasting, and the accuracy of short-time wind energy forecasting is improved.

Description

Technical Field

[0001] The invention belongs to the field of computer processing technology, and in particular relates to a method and system for correcting wave effects of offshore wind energy short-term forecasts. Background Art

[0002] Offshore wind power is becoming a key development direction for wind power construction due to its great potential and advantages such as not occupying land resources. However, compared with onshore wind power, the average power and short-term power fluctuations of offshore wind power are regulated by the phenomenon of sea surface waves, which brings new challenges to the stability and safety of power system operation.

[0003] High-resolution wind energy simulation based on computational fluid dynamics is one of the main means of hourly short-term wind energy forecasting. The wind field data and wind power data it outputs can be used to predict offshore wind power grid disturbances, improve wind power grid connection efficiency, and reduce grid operation costs, which is of great significance for the stable operation of the power system and improving economic benefits. In recent years, with the rapid development of computer computing power, the ultra-high-resolution large eddy simulation technology in the branch of computational fluid dynamics technology is gradually being applied to short-term wind energy forecasting due to its advantages of accurately simulating atmospheric boundary layer turbulence, wind turbine wakes and their interactions. However, the high computational complexity of large eddy simulation is still the main reason that restricts its large-scale calculation and accurate analysis of boundary wave processes. The current short-term wind energy forecasting method cannot accurately analyze the impact of sea surface waves on wind energy forecasting, resulting in low accuracy of short-term wind energy forecasting. Summary of the invention

[0004] To solve the above problems, the present invention proposes a method and system for correcting wave effects in offshore wind energy short-term forecasts. This method can take into account the impact of sea surface waves on wind energy forecasts when forecasting short-term wind energy, thereby improving the accuracy of short-term wind energy forecasts.

[0005] To achieve the above object, the present invention adopts the following technical solutions:

[0006] A method for correcting wave effects in offshore wind energy short-term forecasting comprises the following steps:

[0007] S1. Generate the initial three-dimensional wind field of the wind farm and the wave spectrum characteristic parameters of the sea-air interface based on meteorological data;

[0008] S2. Initializing the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters and wind turbine indicator functions;

[0009] S3, based on the fan parameters and the fan indicator function, calculate the fan grid thrust;

[0010] S4. Based on the three-dimensional wind field initial field and wave spectrum characteristic parameters of the wind farm, the total moving surface drag force is calculated;

[0011] S5. Calculate the short-term wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the grid thrust of the wind turbine and the total moving surface drag force;

[0012] S6. Output short-term wind energy data based on a preset format.

[0013] Preferably, the specific process of step S1 is:

[0014] S11, obtaining meteorological data, where the meteorological data includes wind field data, wind energy data and wave data;

[0015] S12, simulating the three-dimensional wind farm initial field of the wind farm based on the wind farm data, and then generating a gridded time series based on the meteorological data, the size of the simulation area and the corresponding grid spacing of the simulation area as the wind energy value for simulation;

[0016] S13. Obtain the phase velocity based on the wave data and generate the wave spectrum characteristic parameters of the sea-air interface. The calculation formula of the phase velocity is:

[0017]

[0018] c p =(c px ,c py ),

[0019] In the formula, x is the horizontal abscissa; y is the horizontal ordinate; t is time; c px is the phase velocity in the x direction; is the phase velocity in the y direction; η is the wave height; represents the square of the magnitude of the horizontal gradient of wave height; c p is the phase velocity.

[0020] Preferably, the specific process of step S2 is:

[0021] S21, initializing the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters; wherein the wind turbine parameters include the wind turbine position coordinates (x wt ,y wt ), fan thrust coefficient C T , fan brake disc diameter D, fan brake disc thickness s, fan brake disc axis height H and fan direction θ wt ;

[0022] S22, calculating a fan indicator function according to the fan parameters; wherein the fan indicator function includes a normalized indicator function and the fan smoothing indicator function;

[0023] The calculation formula of the normalized indicator function is:

[0024]

[0025] In the formula, I wt (x, y, z) is the normalized fan indicator function; H(t') is the Heaviside function; is the axial local coordinate of the fan brake disc; s is the thickness of the fan brake disc; D is the diameter of the fan brake disc; is the radial distance centered on the axis of the fan brake disc; V is the spatial volume of the fan brake disc;

[0026] The calculation formula of the Heaviside function is:

[0027]

[0028] Where H(t') is the Heaviside function; t' is the function independent variable;

[0029] Among them, the calculation formula of the fan brake disc coordinates is:

[0030]

[0031] In the formula, The local coordinate of the fan brake disc axial direction; The local coordinates of the wind turbine brake disc span; is the vertical local coordinate of the wind turbine brake disc; z is the vertical coordinate; θ wt is the direction of the fan; x wt is the horizontal coordinate of the fan position; wt is the vertical coordinate of the fan position;

[0032] The calculation formula of the radial distance centered on the axis of the fan brake disc is:

[0033]

[0034] In the formula, is the radial distance centered on the axis of the fan brake disc; The local coordinate of the fan brake disc axial direction; The local coordinates of the wind turbine brake disc span;

[0035] Among them, the calculation formula for the space volume of the fan brake disc is:

[0036] V=πsD 2 / 4,

[0037] Where, V is the volume of the fan brake disc; s is the thickness of the fan brake disc; D is the diameter of the fan brake disc;

[0038] The calculation formula of the fan smoothing indicator function is:

[0039] R wt(x,y,z)=∫∫∫G(x-x')G(y-y')(z-z')I wt (x',y',z')dx'dy'dz',

[0040] In the formula, Rn wt (x, y, z) is the fan smoothing indicator function; G(t') is the integral Gaussian kernel function; x', y' and z' are the integral independent variables;

[0041] Among them, the calculation formula of the integral Gaussian kernel function is:

[0042]

[0043] Where G(t') is the integrated Gaussian kernel function; Δ g is the filter width; a is the adjustable coefficient; Δx is the grid spacing in the x direction; Δy is the grid spacing in the y direction; Δz is the grid spacing in the z direction.

[0044] Preferably, the specific process of step S3 is:

[0045] S31, calculate the average wind speed of the fan space (u) according to the fan position coordinates, current wind speed and fan smoothing indicator function d ,v d ), the calculation formula for the average wind speed in the fan space is:

[0046] u d =∫∫∫Rn wt (x,y,z)u(x,y,z)dxdydz,

[0047] v d =∫∫∫Rn wt (x,y,z)v(x,y,z)dxdydz,

[0048] In the formula, u d is the x-direction component of the average wind speed in the wind turbine space; u(x,y,z) is the x-direction component of the wind speed in the wind field; v d is the y-direction component of the average wind speed in the wind turbine space; v(x,y,z) is the y-direction component of the wind speed in the wind field;

[0049] The average wind speed in the fan space (u d ,v d ) and the spatial and temporal average wind speed of the wind turbine at the previous moment Perform e-index time filtering and averaging to calculate the current wind turbine's temporal and spatial average wind speed The calculation formula for the current wind turbine time-space average wind speed is:

[0050]

[0051] In the formula, is the x-direction component of the current spatiotemporal average wind speed of the wind turbine; Δt is the time integration step of the Navier-Stokes equation; τ is the adjustable relaxation time scale; is the x-direction component of the spatial and temporal average wind speed of the wind turbine at the previous moment; is the y-direction component of the current spatial and temporal average wind speed of the wind turbine; is the y-direction component of the wind turbine's temporal and spatial average wind speed at the previous moment;

[0052] S32. Calculate the total thrust of the fan brake disc. The calculation formula for the total thrust of the fan brake disc is:

[0053]

[0054] In the formula, F wt is the total thrust of the fan brake disc; D is the diameter of the fan brake disc; ρ is the air density; C T is the fan thrust coefficient;

[0055] S33, distribute the total thrust of the wind turbine brake disc to each simulation area in the initial field of the three-dimensional wind farm of the wind farm, and calculate the grid thrust f of the wind turbine wt =(f x (x,y,z),f y (x,y,z)), the calculation formula of the fan grid thrust is:

[0056] f x (x,y,z)=F wt R wt cos(θ wt ),

[0057] f y (x,y,z)=F wt R wt sin(θ wt ),

[0058] In the formula, f x (x, y, z) is the x-direction component of the fan grid thrust; F wt is the total thrust of the fan brake disc; θ wt is the direction of the fan; f y (x,y,z) is the y-direction component of the fan grid thrust.

[0059] Preferably, the specific process of step S4 is:

[0060] S41. In the initial field of the three-dimensional wind field of the wind farm, take the wind field (u, v) at a height of 2.5Δz, perform low-pass filtering on the 2Δ scale, and obtain the air-sea interface filtered wind field in,

[0061] S42, obtaining the Euler velocity field, sub-grid effective wave height, wave height and wave phase velocity of the wave in the wave spectrum characteristic parameters;

[0062] S43. Calculate the tangential friction of waves relative to the atmosphere based on the filtered wind field, Euler velocity field and sub-grid effective wave height at the sea-air interface. The calculation formula is:

[0063]

[0064]

[0065] In the formula, τ eqw is the tangential friction of the wave relative to the atmosphere; is τ eqw The x-direction component of is the subgrid effective wave height η rms The relevant wall model friction factor; u w is the x-direction component of the wave velocity; The x-direction component of the wind field is filtered for the air-sea interface; The y-direction component of the wind field is filtered for the air-sea interface; is τ eqw The y-direction component of v w is the y-direction component of the wave velocity;

[0066] S44. Calculate the shape stress of waves relative to the atmosphere based on the filtered wind field, wave height and wave phase velocity at the sea-air interface. The calculation formula is:

[0067]

[0068] In the formula, τ wpm is the shape stress of the wave relative to the atmosphere; is τ wpm The x-direction component of The air-sea interface filters the wind field; c p is the wave phase velocity; is the wave height gradient direction; is the unit vector in the x direction; is τ wpm y-direction component; is the unit vector in the y direction; c px is the wave phase velocity in the x direction; c py is the wave phase velocity in the y direction;

[0069] S45. The total moving surface drag force is calculated by summing the tangential friction force of the wave relative to the atmosphere and the shape stress of the wave relative to the atmosphere. The calculation formula is:

[0070] τ MOSD =τ eqw +τwpm ,

[0071] In the formula, τ MOSD is the total moving surface drag force; τ eqw is the tangential friction of the wave relative to the atmosphere; τ wpm is the shape stress of the wave relative to the atmosphere.

[0072] Preferably, the specific process of step S5 is:

[0073] S51, obtaining the background pressure field of the initial field of the three-dimensional wind field of the wind farm;

[0074] S52. Calculate the short-term wind speed of the wind farm based on the background pressure field, wind turbine grid thrust and total moving surface drag force using the Navier-Stokes equation:

[0075]

[0076] In the formula, u c is the short-term wind speed of the wind farm; p is the disturbed air pressure field; p b is the background pressure field; Ω is the angular velocity of the earth's rotation; f wt is the fan drag force;

[0077] S53. Calculate the short-term wind power of a single wind turbine according to the current average wind speed of the wind turbine in time and space, the diameter of the wind turbine brake disk and the wind energy utilization function. The calculation formula is:

[0078]

[0079] Where P wt is the short-term wind power of a single wind turbine; C wt is the wind energy utilization function; D is the diameter of the wind turbine brake disc;

[0080] S54. Use the short-term wind speed of the wind farm and the short-term wind power of a single wind turbine as short-term wind energy data.

[0081] An offshore wind energy short-term forecast wave effect correction system is used to implement the offshore wind energy short-term forecast wave effect correction method. The offshore wind energy short-term forecast wave effect correction system includes a meteorological data initialization module, a wind turbine parameter initialization module, a wind turbine drag force module, a wave drag force module, a numerical calculation module, and a data output module, wherein:

[0082] The meteorological data initialization module is used to generate the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters of the sea-air interface based on the meteorological data;

[0083] The wind turbine parameter initialization module is used to initialize the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters and wind turbine indicator functions;

[0084] The fan drag force module is used to calculate the fan grid thrust based on the fan parameters and the fan indicator function;

[0085] The wave drag force module is used to calculate the total moving surface drag force based on the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters;

[0086] The numerical calculation module is used to calculate the short-term wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the grid thrust of the wind turbine and the total moving surface drag force;

[0087] The data output module is used to output short-term wind energy data based on a preset format.

[0088] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the method for correcting the wave effect of offshore wind energy short-term forecast.

[0089] A computer device comprises a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein the processor implements the method for correcting the wave effect of offshore wind energy short-term forecast when executing the computer program.

[0090] After adopting the above technical solution, the present invention has the following beneficial effects:

[0091] 1. The present invention can generate a three-dimensional wind field initial field of a wind farm and wave spectrum characteristic parameters of the sea-air interface according to meteorological data, and can calculate the wind turbine grid thrust based on the wind turbine parameters and wind turbine indicator function initialized from the three-dimensional wind field initial field of the wind farm, and can calculate the total moving surface drag force based on the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters, and further can calculate the short-term wind energy data based on the three-dimensional wind field initial field of the wind farm, the wind turbine grid thrust and the total moving surface drag force, so that the influence of sea surface waves on the wind energy forecast can be taken into account when forecasting short-term wind energy, thereby improving the accuracy of short-term wind energy forecast.

[0092] 2. The present invention can initialize the wind turbines in the initial field of the three-dimensional wind farm of the wind farm, and then initialize the wind turbine parameters such as the wind turbine position coordinates, wind turbine thrust coefficient, and wind turbine brake disk diameter, so that the obtained wind turbine parameters are more accurate, and the accuracy of the wind turbine indicator function obtained based on the wind turbine parameters is also improved.

[0093] 3. The present invention can obtain the total thrust of the brake disk of the wind machine based on the calculated current spatiotemporal average wind speed of the wind machine combined with the wind machine parameters, and can distribute the obtained total thrust of the brake disk of the wind machine to each simulation area in the three-dimensional wind field initial field of the wind farm, and determine the wind turbine grid point thrust in each simulation area, thereby making the distributed wind turbine grid point thrust more accurate.

[0094] 4. The present invention can obtain the temporal and spatial average wind speed of the fan at the previous moment, and can calculate the temporal and spatial average wind speed of the fan at the previous moment with the calculated spatial average wind speed of the fan to obtain the current temporal and spatial average wind speed of the fan. The current temporal and spatial average wind speed of the fan calculated based on the temporal and spatial average wind speed of the fan at the previous moment is more reasonable.

[0095] 5. The present invention can calculate the background pressure field, wind turbine grid thrust and total moving surface drag force through the NS equation to obtain the short-term wind speed of the wind farm; it can also calculate the short-term wind power of a single wind turbine based on the current wind turbine spatiotemporal average wind speed, wind turbine brake disk diameter and wind energy utilization function, so that the short-term wind energy data contains multiple types of data. BRIEF DESCRIPTION OF THE DRAWINGS

[0096] Figure 1 It is a schematic flow chart of the method for correcting the wave effect of offshore wind energy short-term forecast according to the present invention;

[0097] Figure 2 It is a structural schematic diagram of the offshore wind energy short-term forecast wave effect correction system of the present invention. DETAILED DESCRIPTION

[0098] In order to make the purpose, technical solution and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with the embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0099] like Figure 1 to Figure 2 As shown, a method for correcting wave effects in short-term offshore wind energy forecasting includes the following steps:

[0100] S1. Generate the initial three-dimensional wind field of the wind farm and the wave spectrum characteristic parameters of the sea-air interface based on meteorological data;

[0101] S2. Initializing the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters and wind turbine indicator functions;

[0102] S3, based on the fan parameters and the fan indicator function, calculate the fan grid thrust;

[0103] S4. Based on the three-dimensional wind field initial field and wave spectrum characteristic parameters of the wind farm, the total moving surface drag force is calculated;

[0104] S5. Calculate the short-term wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the grid thrust of the wind turbine and the total moving surface drag force;

[0105] S6. Output short-term wind energy data based on a preset format.

[0106] The specific process of step S1 is:

[0107] S11, obtaining meteorological data, where the meteorological data includes wind field data, wind energy data and wave data;

[0108] S12, simulating the three-dimensional wind farm initial field of the wind farm based on the wind farm data, and then generating a gridded time series based on the meteorological data, the size of the simulation area and the corresponding grid spacing of the simulation area as the wind energy value for simulation;

[0109] S13. Obtain the phase velocity based on the wave data and generate the wave spectrum characteristic parameters of the sea-air interface. The calculation formula of the phase velocity is:

[0110]

[0111] In the formula, x is the horizontal abscissa; y is the horizontal ordinate; t is time; c px is the phase velocity in the x direction; is the phase velocity in the y direction; η is the wave height; represents the square of the magnitude of the horizontal gradient of wave height; c p is the phase velocity.

[0112] The specific process of step S2 is:

[0113] S21, initializing the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters; wherein the wind turbine parameters include the wind turbine position coordinates (x wt ,y wt ), fan thrust coefficient C T , fan brake disc diameter D, fan brake disc thickness s, fan brake disc axis height H and fan direction θ wt ;

[0114] S22, calculating a fan indicator function according to the fan parameters; wherein the fan indicator function includes a normalized indicator function and the fan smoothing indicator function;

[0115] The calculation formula of the normalized indicator function is:

[0116]

[0117] In the formula, I wt (x, y, z) is the normalized fan indicator function; H(t') is the Heaviside function; is the axial local coordinate of the fan brake disc; s is the thickness of the fan brake disc; D is the diameter of the fan brake disc; is the radial distance centered on the axis of the fan brake disc; V is the spatial volume of the fan brake disc;

[0118] The calculation formula of the Heaviside function is:

[0119]

[0120] Where H(t') is the Heaviside function; t' is the function independent variable;

[0121] Among them, the calculation formula of the fan brake disc coordinates is:

[0122]

[0123]

[0124] In the formula, The local coordinate of the fan brake disc axial direction; The local coordinates of the wind turbine brake disc span; is the vertical local coordinate of the wind turbine brake disc; z is the vertical coordinate; θ wt is the direction of the fan; x wt is the horizontal coordinate of the fan position; wt is the vertical coordinate of the fan position;

[0125] The calculation formula of the radial distance centered on the axis of the fan brake disc is:

[0126]

[0127] In the formula, is the radial distance centered on the axis of the fan brake disc; The local coordinate of the fan brake disc axial direction; The local coordinates of the wind turbine brake disc span;

[0128] Among them, the calculation formula for the space volume of the fan brake disc is:

[0129] V=πsD 2 / 4,

[0130] Where, V is the volume of the fan brake disc; s is the thickness of the fan brake disc; D is the diameter of the fan brake disc;

[0131] The calculation formula of the fan smoothing indicator function is:

[0132] R wt (x,y,z)=∫∫∫G(x-x')G(y-y')(z-z')I wt (x',y',z')dx'dy'dz',

[0133] In the formula, Rn wt (x, y, z) is the fan smoothing indicator function; G(t') is the integral Gaussian kernel function; x', y' and z' are the integral independent variables;

[0134] Among them, the calculation formula of the integral Gaussian kernel function is:

[0135]

[0136] Where G(t') is the integrated Gaussian kernel function; Δ g is the filter width; a is the adjustable coefficient; Δx is the grid spacing in the x direction; Δy is the grid spacing in the y direction; Δz is the grid spacing in the z direction.

[0137] The specific process of step S3 is:

[0138] S31, calculate the average wind speed of the fan space (u) according to the fan position coordinates, current wind speed and fan smoothing indicator function d ,v d ), the calculation formula for the average wind speed in the fan space is:

[0139] u d =∫∫∫Rn wt (x,y,z)u(x,y,z)dxdydz,

[0140] v d =∫∫∫Rn wt (x,y,z)v(x,y,z)dxdydz,

[0141] In the formula, u d is the x-direction component of the average wind speed in the wind turbine space; u(x,y,z) is the x-direction component of the wind speed in the wind field; v d is the y-direction component of the average wind speed in the wind turbine space; v(x,y,z) is the y-direction component of the wind speed in the wind field;

[0142] The average wind speed in the fan space (u d ,v d ) and the spatial and temporal average wind speed of the wind turbine at the previous moment Perform e-index time filtering and averaging to calculate the current wind turbine's temporal and spatial average wind speed The calculation formula for the current wind turbine time-space average wind speed is:

[0143]

[0144] In the formula, is the x-direction component of the current spatiotemporal average wind speed of the wind turbine; Δt is the time integration step of the Navier-Stokes equation; τ is the adjustable relaxation time scale; is the x-direction component of the spatial and temporal average wind speed of the wind turbine at the previous moment; is the y-direction component of the current spatial and temporal average wind speed of the wind turbine; is the y-direction component of the wind turbine's temporal and spatial average wind speed at the previous moment;

[0145] S32. Calculate the total thrust of the fan brake disc. The calculation formula for the total thrust of the fan brake disc is:

[0146]

[0147] In the formula, F wt is the total thrust of the fan brake disc; D is the diameter of the fan brake disc; ρ is the air density; C T is the fan thrust coefficient;

[0148] S33, distribute the total thrust of the wind turbine brake disc to each simulation area in the initial field of the three-dimensional wind farm of the wind farm, and calculate the grid thrust f of the wind turbine wt =(f x (x,y,z),f y (x,y,z)), the calculation formula of the fan grid thrust is:

[0149] f x (x,y,z)=F wt R wt cos(θ wt ),

[0150] f y (x,y,z)=F wt R wt sin(θ wt ),

[0151] In the formula, f x (x, y, z) is the x-direction component of the fan grid thrust; F wt is the total thrust of the fan brake disc; θ wt is the direction of the fan; f y (x,y,z) is the y-direction component of the fan grid thrust.

[0152] The specific process of step S4 is:

[0153] S41. In the initial field of the three-dimensional wind field of the wind farm, take the wind field (u, v) at a height of 2.5Δz, perform low-pass filtering on the 2Δ scale, and obtain the air-sea interface filtered wind field in,

[0154] S42, obtaining the Euler velocity field, sub-grid effective wave height, wave height and wave phase velocity of the wave in the wave spectrum characteristic parameters;

[0155] S43. Calculate the tangential friction of waves relative to the atmosphere based on the filtered wind field, Euler velocity field and sub-grid effective wave height at the sea-air interface. The calculation formula is:

[0156]

[0157] In the formula, τ eqw is the tangential friction of the wave relative to the atmosphere; is τ eqw The x-direction component of is the subgrid effective wave height η rms The relevant wall model friction factor; u w is the x-direction component of the wave velocity; The x-direction component of the wind field is filtered for the air-sea interface; The y-direction component of the wind field is filtered for the air-sea interface; is τ eqw The y-direction component of b w is the y-direction component of the wave velocity;

[0158] S44. Calculate the shape stress of waves relative to the atmosphere based on the filtered wind field, wave height and wave phase velocity at the sea-air interface. The calculation formula is:

[0159]

[0160] In the formula, τ wpm is the shape stress of the wave relative to the atmosphere; is τ wpm The x-direction component of The air-sea interface filters the wind field; c p is the wave phase velocity; is the wave height gradient direction; is the unit vector in the x direction; is τ wpm y-direction component; is the unit vector in the y direction; c px is the wave phase velocity in the x direction; c py is the wave phase velocity in the y direction;

[0161] S45. The total moving surface drag force is calculated by summing the tangential friction force of the wave relative to the atmosphere and the shape stress of the wave relative to the atmosphere. The calculation formula is:

[0162] τ MOSD =τ eqw +τ wpm ,

[0163] In the formula, τ MOSD is the total moving surface drag force; τ eqw is the tangential friction of the wave relative to the atmosphere; τ wpm is the shape stress of the wave relative to the atmosphere.

[0164] The specific process of step S5 is:

[0165] S51, obtaining the background pressure field of the initial field of the three-dimensional wind field of the wind farm;

[0166] S52. Calculate the short-term wind speed of the wind farm based on the background pressure field, wind turbine grid thrust and total moving surface drag force using the Navier-Stokes equation:

[0167]

[0168] In the formula, u c is the short-term wind speed of the wind farm; p is the disturbed air pressure field; p b is the background pressure field; Ω is the angular velocity of the earth's rotation; f wt is the fan drag force;

[0169] S53. Calculate the short-term wind power of a single wind turbine according to the current average wind speed of the wind turbine in time and space, the diameter of the wind turbine brake disk and the wind energy utilization function. The calculation formula is:

[0170]

[0171] Where P wt is the short-term wind power of a single wind turbine; C wt is the wind energy utilization function; D is the diameter of the wind turbine brake disc;

[0172] S54. Use the short-term wind speed of the wind farm and the short-term wind power of a single wind turbine as short-term wind energy data.

[0173] A short-term forecast wave effect correction system for offshore wind energy is used to implement the short-term forecast wave effect correction method for offshore wind energy. The short-term forecast wave effect correction system for offshore wind energy includes a meteorological data initialization module 1, a wind turbine parameter initialization module 2, a wind turbine drag force module 3, a wave drag force module 4, a numerical calculation module 5, and a data output module 6, wherein:

[0174] The meteorological data initialization module 1 is used to generate the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters of the sea-air interface based on the meteorological data;

[0175] The wind turbine parameter initialization module 2 is used to initialize the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters and wind turbine indicator functions;

[0176] The fan drag force module 3 is used to calculate the fan grid thrust based on the fan parameters and the fan indicator function;

[0177] The wave drag force module 4 is used to calculate the total moving surface drag force based on the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters;

[0178] The numerical calculation module 5 is used to calculate the short-term wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the grid thrust of the wind turbine and the total moving surface drag force;

[0179] The data output module 6 is used to output short-term wind energy data based on a preset format.

[0180] A computer-readable storage medium stores a computer program, which, when executed by a processor, implements the method for correcting the wave effect of offshore wind energy short-term forecast.

[0181] A computer device comprises a processor, a memory and a computer program stored in the memory and configured to be executed by the processor, wherein the processor implements the method for correcting the wave effect of offshore wind energy short-term forecast when executing the computer program.

[0182] The above is only a preferred specific embodiment of the present invention, but the protection scope of the present invention is not limited thereto. Any changes or substitutions that can be easily thought of by a person skilled in the art within the technical scope disclosed by the present invention should be included in the protection scope of the present invention. Therefore, the protection scope of the present invention should be based on the protection scope of the claims.

Claims

1. A method for correcting wave effects in short-term offshore wind energy forecast, characterized in that: The following steps are involved: S1. Generate the initial three-dimensional wind field of the wind farm and the wave spectrum characteristic parameters of the sea-air interface based on meteorological data; S2. Initializing the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters and wind turbine indicator functions; S3, based on the fan parameters and the fan indicator function, calculate the fan grid thrust; S4. Based on the three-dimensional wind field initial field and wave spectrum characteristic parameters of the wind farm, the total moving surface drag force is calculated; S5. Calculate the short-term wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the grid thrust of the wind turbine and the total moving surface drag force; S6. Output short-term wind energy data based on a preset format.

2. A method for correcting wave effects in offshore wind energy short-term forecast according to claim 1, characterized in that: The specific process of step S1 is: S11, obtaining meteorological data, the meteorological data including wind field data, wind energy data and wave data; S12, simulating the three-dimensional wind farm initial field of the wind farm based on the wind farm data, and then generating a gridded time series based on the meteorological data, the size of the simulation area and the corresponding grid spacing of the simulation area as the wind energy value for simulation; S13. Obtain the phase velocity based on the wave data and generate the wave spectrum characteristic parameters of the sea-air interface. The calculation formula of the phase velocity is: c p =(c px ,c py ), In the formula, x is the horizontal abscissa; y is the horizontal ordinate; t is time; c px is the phase velocity in the x direction; c py is the phase velocity in the y direction; η is the wave height; represents the square of the horizontal gradient magnitude of the wave height; c p is the phase velocity.

3. A method for correcting wave effects in short-term forecast of offshore wind energy according to claim 2, characterized in that: The specific process of step S2 is: S21, initializing the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters; wherein the wind turbine parameters include the wind turbine position coordinates (x wt ,y wt ), fan thrust coefficient C T , fan brake disc diameter D, fan brake disc thickness s, fan brake disc axis height H and fan direction θ wt ; S22, calculating a fan indicator function according to the fan parameters; wherein the fan indicator function includes a normalized indicator function and the fan smoothing indicator function; The calculation formula of the normalized indicator function is: In the formula, I wt (x, y, z) is the normalized fan indicator function; H(t') is the Heaviside function; is the axial local coordinate of the fan brake disc; s is the thickness of the fan brake disc; D is the diameter of the fan brake disc; is the radial distance centered on the axis of the fan brake disc; V is the spatial volume of the fan brake disc; The calculation formula of the Heaviside function is: Where H(t') is the Heaviside function; t' is the function independent variable; Among them, the calculation formula of the fan brake disc coordinates is: In the formula, The local coordinate of the fan brake disc axial direction; The local coordinates of the wind turbine brake disc span; is the vertical local coordinate of the wind turbine brake disc; z is the vertical coordinate; θ wt is the direction of the fan; x wt is the horizontal coordinate of the fan position; wt is the vertical coordinate of the fan position; The calculation formula of the radial distance centered on the axis of the fan brake disc is: In the formula, is the radial distance centered on the axis of the fan brake disc; The local coordinate of the fan brake disc axial direction; The local coordinates of the wind turbine brake disc span; Among them, the calculation formula for the space volume of the fan brake disc is: V=πsD 2 / 4, Where, V is the volume of the fan brake disc; s is the thickness of the fan brake disc; D is the diameter of the fan brake disc; The calculation formula of the fan smoothing indicator function is: Rn wt (x,y,z)=∫∫∫G(x-x')G(y-y')(z-z')I wt (x',y',z')dx'dy'dz', In the formula, Rn wt (x, y, z) is the fan smoothing indicator function; G(t') is the integral Gaussian kernel function; x', y' and z' are the integral independent variables; Among them, the calculation formula of the integral Gaussian kernel function is: Where G(t') is the integrated Gaussian kernel function; Δ g is the filter width; a is the adjustable coefficient; Δx is the grid spacing in the x direction; Δy is the grid spacing in the y direction; Δz is the grid spacing in the z direction.

4. A method for correcting wave effects in short-term offshore wind energy forecasting as claimed in claim 3, characterized in that: The specific process of step S3 is: S31, calculate the average wind speed of the fan space (u) according to the fan position coordinates, current wind speed and fan smoothing indicator function d ,v d ), the calculation formula for the average wind speed in the fan space is: u d =∫∫∫Rn wt (x,y,z)u(x,y,z)dxdydz, v d =∫∫∫Rn wt (x,y,z)v(x,y,z)dxdydz, In the formula, u d is the x-direction component of the average wind speed in the wind turbine space; u(x,y,z) is the x-direction component of the wind speed in the wind field; v d is the y-direction component of the average wind speed in the wind turbine space; v(x,y,z) is the y-direction component of the wind speed in the wind field; The average wind speed in the fan space (u d ,v d ) and the average wind speed of the wind turbine in time and space at the previous moment Perform e-index time filtering and averaging to calculate the current wind turbine's temporal and spatial average wind speed The calculation formula for the current wind turbine time-space average wind speed is: In the formula, is the x-direction component of the current spatiotemporal average wind speed of the wind turbine; Δt is the time integration step of the Navier-Stokes equation; τ is the adjustable relaxation time scale; is the x-direction component of the spatial and temporal average wind speed of the wind turbine at the previous moment; is the y-direction component of the current spatial and temporal average wind speed of the wind turbine; is the y-direction component of the wind turbine's temporal and spatial average wind speed at the previous moment; S32. Calculate the total thrust of the fan brake disc. The calculation formula for the total thrust of the fan brake disc is: In the formula, F wt is the total thrust of the fan brake disc; D is the diameter of the fan brake disc; ρ is the air density; C T is the fan thrust coefficient; S33, distribute the total thrust of the wind turbine brake disc to each simulation area in the initial field of the three-dimensional wind farm of the wind farm, and calculate the grid thrust f of the wind turbine wt =(f x (x,y,z),f y (x,y,z)), the calculation formula of the fan grid thrust is: f x (x,y,z)=F wt R wt cos(θ wt ), f y (x,y,z)=F wt R wt sin(θ wt ), In the formula, f x (x, y, z) is the x-direction component of the fan grid thrust; F wt is the total thrust of the fan brake disc; θ wt is the direction of the fan; f y (x,y,z) is the y-direction component of the fan grid thrust.

5. A method for correcting wave effects in short-term offshore wind energy forecasting as claimed in claim 4, characterized in that: The specific process of step S4 is: S41. In the initial field of the three-dimensional wind field of the wind farm, take the wind field (u, v) at a height of 2.5Δz, perform low-pass filtering on the 2Δ scale, and obtain the air-sea interface filtered wind field in, S42, obtaining the Euler velocity field, sub-grid effective wave height, wave height and wave phase velocity of the wave in the wave spectrum characteristic parameters; S43. Calculate the tangential friction of waves relative to the atmosphere based on the filtered wind field, Euler velocity field and sub-grid effective wave height at the sea-air interface. The calculation formula is: In the formula, τ eqw is the tangential friction of the wave relative to the atmosphere; is τ eqw The x-direction component of is the subgrid effective wave height η rms The relevant wall model friction factor; u w is the x-direction component of the wave velocity; The x-direction component of the wind field is filtered for the air-sea interface; The y-direction component of the wind field is filtered for the air-sea interface; is τ eqw The y-direction component of v w is the y-direction component of the wave velocity; S44. Calculate the shape stress of waves relative to the atmosphere based on the filtered wind field, wave height and wave phase velocity at the sea-air interface. The calculation formula is: In the formula, τ wpm is the shape stress of the wave relative to the atmosphere; is τ wpm The x-direction component of The air-sea interface filters the wind field; c p is the wave phase velocity; is the wave height gradient direction; is the unit vector in the x direction; is τ wpm y-direction component; is the unit vector in the y direction; c px is the wave phase velocity in the x direction; c py is the wave phase velocity in the y direction; S45. The total moving surface drag force is calculated by summing the tangential friction force of the wave relative to the atmosphere and the shape stress of the wave relative to the atmosphere. The calculation formula is: t MOSD =t eqw +t wpm , In the formula, τ MOSD is the total moving surface drag force; τ eqw is the tangential friction of the wave relative to the atmosphere; τ wpm is the shape stress of the wave relative to the atmosphere.

6. A method for correcting wave effects in short-term offshore wind energy forecasting as claimed in claim 5, characterized in that: The specific process of step S5 is: S51, obtaining the background pressure field of the initial field of the three-dimensional wind field of the wind farm; S52. Calculate the short-term wind speed of the wind farm based on the background pressure field, wind turbine grid thrust and total moving surface drag force using the Navier-Stokes equation: In the formula, u c is the short-term wind speed of the wind farm; p is the disturbed air pressure field; p b is the background pressure field; Ω is the angular velocity of the earth's rotation; f wt is the fan drag force; S53. Calculate the short-term wind power of a single wind turbine according to the current average wind speed of the wind turbine in time and space, the diameter of the wind turbine brake disk and the wind energy utilization function. The calculation formula is: Where P wt is the short-term wind power of a single wind turbine; C wt is the wind energy utilization function; D is the diameter of the wind turbine brake disc; S54. Use the short-term wind speed of the wind farm and the short-term wind power of a single wind turbine as short-term wind energy data.

7. An offshore wind energy short-term forecast wave effect correction system, characterized in that: The method for correcting the wave effect of offshore wind energy short-term forecast according to any one of claims 1 to 6 is used to implement the method, wherein the offshore wind energy short-term forecast wave effect correction system comprises a meteorological data initialization module, a wind turbine parameter initialization module, a wind turbine drag force module, a wave drag force module, a numerical calculation module, and a data output module, wherein: The meteorological data initialization module is used to generate the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters of the sea-air interface based on the meteorological data; The wind turbine parameter initialization module is used to initialize the wind turbine based on the three-dimensional wind farm initial field of the wind farm to obtain wind turbine parameters and wind turbine indicator functions; The fan drag force module is used to calculate the fan grid thrust based on the fan parameters and the fan indicator function; The wave drag force module is used to calculate the total moving surface drag force based on the three-dimensional wind field initial field of the wind farm and the wave spectrum characteristic parameters; The numerical calculation module is used to calculate the short-term wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the grid thrust of the wind turbine and the total moving surface drag force; The data output module is used to output short-term wind energy data based on a preset format.

8. A computer-readable storage medium having a computer program stored thereon, characterized in that: When the program is executed by a processor, the method for correcting wave effects in offshore wind energy short-term forecasting as described in any one of claims 1 to 6 is implemented.

9. A computer device, characterized in that: The method comprises a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein when the processor executes the computer program, the method for correcting the wave effect of offshore wind energy short-term forecast as described in any one of claims 1 to 6 is implemented.

Images