A method and system for correcting short-term offshore wind energy prediction wave effects

Abstract

The application discloses a kind of offshore wind energy short-time forecast wave effect correction method and system, the method includes the following steps: S1, initial field of three-dimensional wind field of wind farm is generated based on meteorological data and wave spectrum characteristic parameter of sea-air interface;S2, based on three-dimensional wind field of wind farm initial field, initialization is carried out to wind turbine, and wind turbine parameter and wind turbine indication function are obtained;S3, based on wind turbine parameter and wind turbine indication function, wind turbine grid point thrust is calculated;S4, based on three-dimensional wind field of wind farm initial field and wave spectrum characteristic parameter, total moving surface drag force is calculated;S5, based on three-dimensional wind field of wind farm initial field, wind turbine grid point thrust and total moving surface drag force, short-time wind energy data are calculated;S6, short-time wind energy data are output based on preset format;The method can consider the influence of sea surface wave on wind energy forecast when short-time wind energy is forecasted, and the accuracy of short-time wind energy forecast is improved.

Description

Technical Field

[0001] This invention belongs to the field of computer processing technology, specifically relating to a method and system for correcting wave effects in short-term forecasts of offshore wind energy. Background Technology

[0002] Offshore wind power is becoming a key development direction for wind power construction due to its large wind energy potential and the fact that it does not occupy land resources. However, compared with onshore wind power, the average power and short-term power fluctuations of offshore wind power are affected by sea surface waves, which brings new challenges to the stability and security of power system operation.

[0003] High-resolution wind energy simulation based on computational fluid dynamics (CFD) is one of the main methods for hourly short-term wind energy forecasting. Its output wind field and wind power data can be used to predict offshore wind power grid disturbances, improve wind power grid connection efficiency, and reduce grid operating costs, which is of great significance for the stable operation of the power system and improving economic efficiency. In recent years, with the rapid development of computing power, large eddy simulation (LES), a branch of CFD technology with ultra-high resolution, is gradually being applied to short-term wind energy forecasting due to its advantages such as accurately simulating atmospheric boundary layer turbulence, wind turbine wakes, and their interactions. However, the high computational cost of LES remains the main reason limiting its large-scale computation and accurate analysis of boundary wave processes. Current short-term wind energy forecasting methods cannot accurately analyze the impact of sea surface waves on wind energy forecasting, resulting in low accuracy in short-term wind energy forecasts. Summary of the Invention

[0004] To address the aforementioned issues, this invention proposes a method and system for correcting wave effects in short-term offshore wind energy forecasting. This method takes into account the impact of sea surface waves on wind energy forecasting, thereby improving the accuracy of short-term wind energy forecasting.

[0005] To achieve the above objectives, the present invention adopts the following technical solution:

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

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

[0008] S2. Initialize the wind turbine based on the initial three-dimensional wind field of the wind farm to obtain the wind turbine parameters and wind turbine indication function;

[0009] S3. Based on the fan parameters and fan indication function, the fan grid thrust is calculated;

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

[0011] S5. Short-term wind energy data are calculated based on the three-dimensional initial wind field of the wind farm, the wind turbine grid thrust 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 as follows:

[0014] S11. Obtain meteorological data, including wind field data, wind energy data, and wave data;

[0015] S12. Based on wind field data, the initial three-dimensional wind field of the wind farm is simulated. The initial three-dimensional wind field of the wind farm is then used to generate a gridded time series as wind energy values ​​for simulation based on meteorological data, the size of the simulation area and the corresponding grid spacing of the simulation area.

[0016] S13. Obtain the phase velocity from the wave data, and generate the wave spectrum characteristic parameters of the air-sea interface. The formula for calculating 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; c represents the square of the magnitude of the horizontal gradient of the wave height. p The phase velocity is given.

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

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

[0022] S22. Calculate the fan indicator function based on the fan parameters; wherein, the fan indicator function includes a normalized indicator function and the fan smoothing indicator function;

[0023] The formula for calculating the normalized indicator function is:

[0024]

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

[0026] The formula for calculating the Herveside function is as follows:

[0027]

[0028] In the formula, H(t') is the Herveside function; t' is the independent variable of the function;

[0029] The formula for calculating the wind turbine brake disc coordinates is as follows:

[0030]

[0031] In the formula, For the local coordinates of the wind turbine brake disc axis; The wind turbine's braking disk extends to local coordinates; Z represents the vertical local coordinates of the wind turbine's brake disc; Z represents the vertical coordinates; θ represents the vertical coordinates. wt For the orientation of the fan; x wt y is the x-coordinate of the wind turbine's location; wt The vertical coordinate represents the location of the wind turbine.

[0032] The formula for calculating the radial distance centered on the axis of the wind turbine brake disc is as follows:

[0033]

[0034] In the formula, The radial distance centered on the axis of the wind turbine brake disc; For the local coordinates of the wind turbine brake disc axis; The wind turbine's braking disk extends to local coordinates;

[0035] The formula for calculating the spatial volume of the wind turbine brake disc is as follows:

[0036] V=πsD 2 / 4,

[0037] In the formula, V is the spatial volume of the wind turbine brake disc; s is the thickness of the wind turbine brake disc; and D is the diameter of the wind turbine brake disc.

[0038] The formula for calculating the wind turbine smoothing indicator function is:

[0039] Rn 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 wind turbine smoothing indicator function; G(t') is the integral Gaussian kernel function; x',y' and z' are the independent variables of the integral;

[0041] The formula for calculating the integral Gaussian kernel function is as follows:

[0042]

[0043] In the formula, G(t') is the integral Gaussian kernel function; Δ g Δx is the filter width; Δy 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 as follows:

[0045] S31. Calculate the spatial average wind speed (u) of the wind turbine based on the turbine's location coordinates, current wind speed, and the turbine smoothing indicator function. d ,v d The formula for calculating 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 The x-axis component of the average wind speed in the wind turbine space; u(x,y,z) is the x-axis component of the wind speed in the wind field; v d Let v(x,y,z) be the y-component of the average wind speed in the wind turbine space; v(x,y,z) be the y-component of the wind speed in the wind field.

[0049] The average wind speed in the fan space (u) d ,v d The spatiotemporal average wind speed of the wind turbine at the previous moment By performing an e-exponential time-filtered average, the current spatiotemporal average wind speed of the wind turbine can be calculated. The formula for calculating the current spatiotemporal average wind speed of the wind turbine is:

[0050]

[0051] In the formula, denoted as 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 equations; τ is the adjustable relaxation time scale. The x-direction component of the spatiotemporal average wind speed of the wind turbine at the previous moment; The y-direction component of the current spatiotemporal average wind speed of the wind turbine; The y-direction component of the spatiotemporal average wind speed of the wind turbine at the previous moment;

[0052] S32. Calculate the total thrust of the wind turbine brake disc. The formula for calculating the total thrust of the wind turbine brake disc is:

[0053]

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

[0055] S33. Distribute the total thrust of the wind turbine's brake disc to each simulated region in the initial three-dimensional wind field of the wind farm, and calculate the wind turbine grid thrust f. wt =(f x (x,y,z),f y (x,y,z)), the formula for calculating the wind turbine 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) represents the x-axis component of the wind turbine grid thrust; F wt For the total thrust of the wind turbine brake disc; θ wt For the orientation of the fan; f y (x,y,z) represents the y-direction component of the wind turbine grid thrust.

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

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

[0061] S42. Obtain the Euler velocity field, subgrid effective wave height, wave height and wave phase velocity of the wave from the wave spectrum characteristic parameters;

[0062] S43. Based on the filtered wind field at the air-sea interface, the Eulerian velocity field, and the effective wave height of the subgrid, calculate the tangential frictional force of the waves relative to the atmosphere. The calculation formula is as follows:

[0063]

[0064]

[0065] In the formula, τ eqw This refers to the tangential frictional force of the wave relative to the atmosphere. For τ eqw The x-direction component; To match the effective wave height η of the subgrid rms Related wall model friction factor; u w The wave velocity is represented by the x-axis component. The x-direction component of the filtered wind field at the air-sea interface; The y-direction component of the filtered wind field at the air-sea interface; For τ eqw The y-direction component; v w The wave velocity component in the y-direction;

[0066] S44. Based on the air-sea interface filtered wind field, wave height, and wave phase velocity, calculate the wave shape stress relative to the atmosphere. The calculation formula is as follows:

[0067]

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

[0069] S45. Summing the tangential frictional force of the wave relative to the atmosphere and the shape stress of the wave relative to the atmosphere, calculate the total drag force of the moving surface. The calculation formula is:

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

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

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

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

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

[0075]

[0076] In the formula, u c p represents the short-term wind speed at the wind farm; p represents the disturbed air pressure field; p b The background pressure field is represented by Ω; the Earth's rotational angular velocity is represented by f. wt For the drag force of the wind turbine;

[0077] S53. Based on the current spatiotemporal average wind speed of the wind turbine, the diameter of the wind turbine's brake disc, and the wind energy utilization function, calculate the short-term wind power of a single wind turbine. The calculation formula is as follows:

[0078]

[0079] In the formula, P wt C represents the short-term wind power of a single wind turbine unit. wt D 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] A short-term forecast wave effect correction system for offshore wind energy is provided to implement the aforementioned short-term forecast wave effect correction method for offshore wind energy. The 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 initial field of the wind farm and the wave spectrum characteristic parameters of the air-sea interface based on the meteorological data.

[0083] The wind turbine parameter initialization module is used to initialize the wind turbine based on the initial three-dimensional wind field of the wind farm, and obtain the wind turbine parameters and the wind turbine indication function.

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

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

[0086] The numerical calculation module is used to calculate short-term wind energy data based on the three-dimensional initial wind field of the wind farm, the wind turbine grid thrust, 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 having a computer program stored thereon, which, when executed by a processor, implements the aforementioned method for correcting wave effects in short-term forecasts of marine wind energy.

[0089] A computer device includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor, when executing the computer program, implements the aforementioned method for correcting wave effects in short-term forecasts of marine wind energy.

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

[0091] 1. This invention can generate a three-dimensional initial wind field of a wind farm and wave spectrum characteristic parameters of the air-sea interface based on meteorological data. It can also calculate the wind turbine grid thrust based on the wind turbine parameters and wind turbine indicator function obtained from the initialization of the three-dimensional initial wind field of the wind farm, and calculate the total moving surface drag force based on the three-dimensional initial wind field of the wind farm and the wave spectrum characteristic parameters. Furthermore, it can calculate short-term wind energy data based on the three-dimensional initial wind field of the wind farm, the wind turbine grid thrust, and the total moving surface drag force. This allows the influence of sea surface waves on wind energy forecasting to be considered when forecasting short-term wind energy, thus improving the accuracy of short-term wind energy forecasting.

[0092] 2. This invention can initialize the wind turbines in the initial three-dimensional wind field of a wind farm, and then initialize wind turbine parameters such as wind turbine position coordinates, wind turbine thrust coefficient, and wind turbine brake disc diameter, thereby making the obtained wind turbine parameters more accurate and improving the accuracy of the wind turbine indication function obtained based on the wind turbine parameters.

[0093] 3. This invention can calculate the total thrust of the wind turbine brake disc by combining the calculated spatiotemporal average wind speed of the current wind turbine with the wind turbine parameters. Furthermore, the obtained total thrust of the wind turbine brake disc can be allocated to each simulated region in the initial field of the three-dimensional wind farm to determine the wind turbine grid thrust of each simulated region, thereby making the allocated wind turbine grid thrust more accurate.

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

[0095] 5. This invention can calculate the background pressure field, wind turbine grid thrust, and total moving surface drag force using the Navier-Stokes equations 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 spatiotemporal average wind speed of the wind turbine, the diameter of the wind turbine brake disc, and the wind energy utilization function, thus enabling the short-term wind energy data to contain multiple types of data. Attached Figure Description

[0096] Figure 1 This is a flowchart illustrating the wave effect correction method for short-term offshore wind energy forecasting of the present invention.

[0097] Figure 2 This is a schematic diagram of the structure of the short-term forecast wave effect correction system for marine wind energy of the present invention. Detailed Implementation

[0098] To make the objectives, technical solutions, and advantages of this invention clearer, the invention will be further described in detail below with reference to embodiments. It should be understood that the specific embodiments described herein are merely illustrative and not intended to limit the invention.

[0099] like Figures 1 to 2 As shown, a method for correcting wave effects in short-term offshore wind energy forecasts 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 air-sea interface based on meteorological data;

[0101] S2. Initialize the wind turbine based on the initial three-dimensional wind field of the wind farm to obtain the wind turbine parameters and wind turbine indication function;

[0102] S3. Based on the fan parameters and fan indication function, the fan grid thrust is calculated;

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

[0104] S5. Short-term wind energy data are calculated based on the three-dimensional initial wind field of the wind farm, the wind turbine grid thrust 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 as follows:

[0107] S11. Obtain meteorological data, including wind field data, wind energy data, and wave data;

[0108] S12. Based on wind field data, the initial three-dimensional wind field of the wind farm is simulated. The initial three-dimensional wind field of the wind farm is then used to generate a gridded time series as wind energy values ​​for simulation based on meteorological data, the size of the simulation area and the corresponding grid spacing of the simulation area.

[0109] S13. Obtain the phase velocity from the wave data, and generate the wave spectrum characteristic parameters of the air-sea interface. The formula for calculating 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; c represents the square of the magnitude of the horizontal gradient of the wave height. p The phase velocity is given.

[0112] The specific process of step S2 is as follows:

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

[0114] S22. Calculate the fan indicator function based on the fan parameters; wherein, the fan indicator function includes a normalized indicator function and the fan smoothing indicator function;

[0115] The formula for calculating the normalized indicator function is:

[0116]

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

[0118] The formula for calculating the Herveside function is as follows:

[0119]

[0120] In the formula, H(t') is the Herveside function; t' is the independent variable of the function;

[0121] The formula for calculating the wind turbine brake disc coordinates is as follows:

[0122]

[0123]

[0124] In the formula, For the local coordinates of the wind turbine brake disc axis; The wind turbine's braking disk extends to local coordinates; Z represents the vertical local coordinates of the wind turbine's brake disc; Z represents the vertical coordinates; θ represents the vertical coordinates. wt For the orientation of the fan; x wt y is the x-coordinate of the wind turbine's location; wt The vertical coordinate represents the location of the wind turbine.

[0125] The formula for calculating the radial distance centered on the axis of the wind turbine brake disc is as follows:

[0126]

[0127] In the formula, The radial distance centered on the axis of the wind turbine brake disc; For the local coordinates of the wind turbine brake disc axis; The wind turbine's braking disk extends to local coordinates;

[0128] The formula for calculating the spatial volume of the wind turbine brake disc is as follows:

[0129] V=πsD 2 / 4,

[0130] In the formula, V is the spatial volume of the wind turbine brake disc; s is the thickness of the wind turbine brake disc; and D is the diameter of the wind turbine brake disc.

[0131] The formula for calculating the wind turbine smoothing indicator function is:

[0132] Rn 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 wind turbine smoothing indicator function; G(t') is the integral Gaussian kernel function; x',y' and z' are the independent variables of the integral;

[0134] The formula for calculating the integral Gaussian kernel function is as follows:

[0135]

[0136] In the formula, G(t') is the integral Gaussian kernel function; Δ g Δx is the filter width; Δy 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 as follows:

[0138] S31. Calculate the spatial average wind speed (u) of the wind turbine based on the turbine's location coordinates, current wind speed, and the turbine smoothing indicator function. d ,v d The formula for calculating 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 The x-axis component of the average wind speed in the wind turbine space; u(x,y,z) is the x-axis component of the wind speed in the wind field; v d Let v(x,y,z) be the y-component of the average wind speed in the wind turbine space; v(x,y,z) be the y-component of the wind speed in the wind field.

[0142] The average wind speed in the fan space (u) d ,v d The spatiotemporal average wind speed of the wind turbine at the previous moment By performing an e-exponential time-filtered average, the current spatiotemporal average wind speed of the wind turbine can be calculated. The formula for calculating the current spatiotemporal average wind speed of the wind turbine is:

[0143]

[0144] In the formula, denoted as 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 equations; τ is the adjustable relaxation time scale. The x-direction component of the spatiotemporal average wind speed of the wind turbine at the previous moment; The y-direction component of the current spatiotemporal average wind speed of the wind turbine; The y-direction component of the spatiotemporal average wind speed of the wind turbine at the previous moment;

[0145] S32. Calculate the total thrust of the wind turbine brake disc. The formula for calculating the total thrust of the wind turbine brake disc is:

[0146]

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

[0148] S33. Distribute the total thrust of the wind turbine's brake disc to each simulated region in the initial three-dimensional wind field of the wind farm, and calculate the wind turbine grid thrust f. wt =(f x (x,y,z),f y (x,y,z)), the formula for calculating the wind turbine 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) represents the x-axis component of the wind turbine grid thrust; F wt For the total thrust of the wind turbine brake disc; θ wt For the orientation of the fan; f y (x,y,z) represents the y-direction component of the wind turbine grid thrust.

[0152] The specific process of step S4 is as follows:

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

[0154] S42. Obtain the Euler velocity field, subgrid effective wave height, wave height and wave phase velocity of the wave from the wave spectrum characteristic parameters;

[0155] S43. Based on the filtered wind field at the air-sea interface, the Eulerian velocity field, and the effective wave height of the subgrid, calculate the tangential frictional force of the waves relative to the atmosphere. The calculation formula is as follows:

[0156]

[0157] In the formula, τ eqw This refers to the tangential frictional force of the wave relative to the atmosphere. For τ eqw The x-direction component; To match the effective wave height η of the subgrid rms Related wall model friction factor; u w The wave velocity is represented by the x-axis component. The x-direction component of the filtered wind field at the air-sea interface; The y-direction component of the filtered wind field at the air-sea interface; For τ eqw The y-direction component; b w The wave velocity is represented by the y-direction component.

[0158] S44. Based on the air-sea interface filtered wind field, wave height, and wave phase velocity, calculate the wave shape stress relative to the atmosphere. The calculation formula is as follows:

[0159]

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

[0161] S45. Summing the tangential frictional force of the wave relative to the atmosphere and the shape stress of the wave relative to the atmosphere, calculate the total drag force of the moving surface. The calculation formula is:

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

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

[0164] The specific process of step S5 is as follows:

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

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

[0167]

[0168] In the formula, u c p represents the short-term wind speed at the wind farm; p represents the disturbed air pressure field; p b The background pressure field is represented by Ω; the Earth's rotational angular velocity is represented by f. wt For the drag force of the wind turbine;

[0169] S53. Based on the current spatiotemporal average wind speed of the wind turbine, the diameter of the wind turbine's brake disc, and the wind energy utilization function, calculate the short-term wind power of a single wind turbine. The calculation formula is as follows:

[0170]

[0171] In the formula, P wt C represents the short-term wind power of a single wind turbine unit. wt D 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 provided to implement the aforementioned short-term forecast wave effect correction method for offshore wind energy. The system 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 initial field of the wind farm and the wave spectrum characteristic parameters of the air-sea 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 initial field of the wind farm to obtain wind turbine parameters and wind turbine indication function;

[0176] The wind turbine drag force module 3 is used to calculate the wind turbine grid thrust based on the wind turbine parameters and the wind turbine indication 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 and wave spectrum characteristic parameters of the wind farm.

[0178] The numerical calculation module 5 is used to calculate short-time wind energy data based on the three-dimensional initial field of the wind farm, the wind turbine grid thrust 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 having a computer program stored thereon, which, when executed by a processor, implements the aforementioned method for correcting wave effects in short-term forecasts of marine wind energy.

[0181] A computer device includes a processor, a memory, and a computer program stored in the memory and configured to be executed by the processor, wherein the processor, when executing the computer program, implements the aforementioned method for correcting wave effects in short-term forecasts of marine wind energy.

[0182] The above description is merely a preferred embodiment of the present invention, but the scope of protection of the present invention is not limited thereto. Any variations or substitutions that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be included within the scope of protection of the present invention. Therefore, the scope of protection of the present invention should be determined by the scope of the claims.

Claims

1. A method for offshore wind energy short-term forecast wave effect correction, characterized in that, The method comprises the following steps: S1, generating an initial field of a three-dimensional wind field of a wind farm and a wave spectrum characteristic parameter of a sea-air interface based on meteorological data; S2, initializing the wind turbine based on the initial field of the three-dimensional wind field of the wind farm to obtain wind turbine parameters and a wind turbine indicator function; S3, calculating the wind turbine grid point thrust based on the wind turbine parameters and the wind turbine indicator function; S4, calculating the total moving surface drag force based on the initial field of the three-dimensional wind field of the wind farm and the wave spectrum characteristic parameter; S5, calculating short-time wind energy data based on the initial field of the three-dimensional wind field of the wind farm, the wind turbine grid point thrust and the total moving surface drag force; S6, outputting the short-time wind energy data in a preset format.

2. A short-term offshore wind energy forecast wave effect correction method according to claim 1, characterized in that, The specific process of step S1 is as follows: S11, obtaining meteorological data, which comprises wind field data, wind energy data and wave data; S12, simulating the initial field of the three-dimensional wind field of the wind farm based on the wind field data, and generating a grid time sequence based on the meteorological data, the size of the simulation area and the corresponding grid distance of the simulation area as wind energy numerical simulation; S13, obtaining the phase velocity based on the wave data to generate the wave spectrum characteristic parameter of the sea-air interface, and the calculation formula of the phase velocity is as follows: c p = (c px ,c py ), where x is the horizontal abscissa; y is the horizontal ordinate; t is time; c px is the x-direction phase velocity; c py is the y-direction phase velocity; η is the wave height; represents the square of the order of magnitude of the horizontal gradient of wave height; c p is the phase velocity.

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

4. A short-term offshore wind energy forecast wave effect correction method according to claim 3, characterized in that, The specific process of step S3 is as follows: S31、According to the fan position coordinates, the current wind speed and the fan smoothing indication function, calculate the fan space average wind speed (u d ,v d ), and the calculation formula of the fan space average wind speed 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, where u d is the fan space-averaged wind speed x-direction component; u(x, y, z) is the wind farm wind speed x-direction component; v d is the fan space-averaged wind speed y-direction component; v(x, y, z) is the wind farm wind speed y-direction component; The fan space average wind speed (u d ,v d ) is subjected to e exponential time filtering average with the previous time fan space-time average wind speed to obtain the current fan space-time average wind speed The calculation formula of the current fan space-time average wind speed is: wherein is the current time-averaged wind speed in the x-direction of the fan; Δt is the time integration step of the Navier-Stokes equation; τ is the adjustable relaxation time scale; is the previous time-averaged wind speed in the x-direction of the fan; is the current time-averaged wind speed in the y-direction of the fan; is the previous time-averaged wind speed in the y-direction of the fan; S32, calculating the total thrust of the wind turbine brake disc, and the calculation formula of the total thrust of the wind turbine brake disc is as follows: F = 0.5 * D * D * ρ * C wt F = 0.5 * D * D * ρ * C T F = 0.5 * D * D * ρ * C S33, distributing the total thrust of the wind turbine brake disc to each simulation area in the initial field of the three-dimensional wind field of the wind farm, calculating the wind turbine grid point thrust f wt = (f x (x, y, z), f y (x, y, z), the calculation formula of the wind turbine grid point thrust is: f x (x,y,z) = F wt R wt cos(θ wt ), f y (x,y,z) = F wt R wt sin(θ wt ), wherein f x (x, y, z) is the wind turbine grid point thrust x component; F wt is the total wind turbine brake disc thrust; θ wt is the wind turbine heading; f y (x, y, z) is the wind turbine grid point thrust y component.

5. A short-term offshore wind energy forecast wave effect correction method according to claim 4, characterized in that, The specific process of step S4 is as follows: S41, in the initial three-dimensional wind field of the wind farm, the wind field (u, v) at the height of 2.5Δz is taken, low-pass filtering is performed on the scale of 2Δ, and the filtered wind field of the sea-air interface is obtained wherein, S42, obtaining the Euler velocity field of the wave, the sub-grid effective wave height, the wave height and the wave phase velocity in the wave spectrum characteristic parameter; S43, calculating the tangential friction of the wave relative to the atmosphere based on the filtered wind field, the Euler velocity field and the sub-grid effective wave height of the sea-air interface, and the calculation formula is as follows: where τ eqw is the tangential friction force of the wave relative to the atmosphere; is the x-direction component of τ eqw ; is the wall friction factor associated with the subgrid significant wave height η rms ; u w is the x-direction component of the wave velocity; is the x-direction component of the filtered wind field at the sea-air interface; is the y-direction component of the filtered wind field at the sea-air interface; is the y-direction component of τ eqw ; v w is the y-direction component of the wave velocity; S44, calculating the shape stress of the wave relative to the atmosphere based on the filtered wind field, the wave height and the wave phase velocity of the sea-air interface, and the calculation formula is as follows: where τ wpm is the shape stress of the wave relative to the atmosphere; is the x-component of τ wpm ; is the filtered wind field at the sea-air interface;c p is the phase velocity of the wave; is the direction of the wave height gradient; is the unit vector in the x-direction; is the y-component of τ wpm ; is the unit vector in the y-direction;c px is the x-direction phase velocity of the wave;c py is the y-direction phase velocity of the wave; S45, summing the tangential friction of the wave relative to the atmosphere and the shape stress of the wave relative to the atmosphere to calculate the total moving surface drag force, and the calculation formula is as follows: τ MOSD = τ eqw + τ wpm , where τ MOSD is the total moving surface drag force; τ eqw is the tangential friction of the wave relative to the atmosphere; and τ wpm is the form stress of the wave relative to the atmosphere.

6. A short-term offshore wind energy forecast wave effect correction method according to claim 5, characterized in that, The specific process of step S5 is as follows: S51, obtaining the background pressure field of the initial field of the three-dimensional wind field of the wind farm; S52, calculating the short-time wind speed of the wind farm based on the Navier-Stokes equation and the background pressure field, the wind turbine grid point thrust and the total moving surface drag force: wherein u c is the short-term wind farm wind speed; p is the perturbation pressure field; p b is the background pressure field; Ω is the earth rotation angular velocity; f wt is the wind turbine drag force; S53, calculating the short-time wind power of the wind turbine based on the current wind turbine space-time average wind speed, the diameter of the wind turbine brake disc and the wind energy utilization function, and the calculation formula is as follows: wherein P wt is the short-time wind power of the wind turbine; C wt is the wind energy utilization function; and D is the diameter of the brake disc of the wind turbine. S54, taking the short-time wind speed of the short-time wind farm and the short-time wind power of the single wind turbine as the short-time wind energy data.

7. A short-term offshore wind energy forecast wave effect correction system, characterized by, The offshore wind energy short-time prediction 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 configured to generate an initial wind field of a wind farm and wave spectrum characteristic parameters of a sea-air interface based on meteorological data. The wind turbine parameter initialization module is configured to initialize the wind turbine based on the initial wind field of the wind farm to obtain wind turbine parameters and a wind turbine indication function. The wind turbine drag force module is configured to calculate a wind turbine grid point thrust based on the wind turbine parameters and the wind turbine indication function. The wave drag force module is configured to calculate a total moving surface drag force based on the initial wind field of the wind farm and the wave spectrum characteristic parameters. The numerical calculation module is configured to calculate short-time wind energy data based on the initial wind field of the wind farm, the wind turbine grid point thrust, and the total moving surface drag force. The data output module is configured to output the short-time wind energy data in a preset format.

8. A computer-readable storage medium having stored thereon a computer program, characterized in that, The program is executed by a processor to implement the offshore wind energy short-time prediction wave effect correction method according to any one of claims 1-6.

9. A computer device, comprising: The offshore wind energy short-time prediction wave effect correction method according to any one of claims 1-6 is implemented by a processor executing a computer program stored in a memory.

Images