A wind power distribution forecasting method and system applicable to complex terrain

Abstract

The present invention discloses a wind power distribution prediction method and system applicable to complex terrain. The method includes: establishing the primary distribution characteristics of the wind field from the spatial distribution of low-altitude wind direction and wind speed; using a fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area; combining the terrain grid point data and the spatio-temporally non-uniform wind profile in the network cell to optimize and adjust the primary distribution characteristics of the wind field to form the refined distribution characteristics of the wind field; calculating the total turbulent kinetic energy of the full spectrum from the simulated atmospheric turbulent kinetic energy of the target area using the sampling frequency limit parameter; substituting the refined distribution characteristics of the wind field and the total turbulent kinetic energy of the full spectrum into the mesoscale meteorological model simulation process of the target area to obtain the wind power distribution characteristics. The present invention can more accurately simulate the wind power distribution, better adapt to complex terrain conditions, and improve the accuracy and reliability of the wind power distribution simulation.

Description

Technical Field

[0001] The present invention belongs to the technical field of atmospheric pollution diffusion simulation, and particularly relates to a wind power distribution forecasting method and system applicable to complex terrain. Background Art

[0002] The simulation of wind fields and turbulence fields mainly relies on mesoscale numerical models, which have the advantages of covering a large area and providing continuous meteorological and environmental data. However, with the development of urbanization and the increase in terrain complexity, the limitations of mesoscale models in resolving local spatial heterogeneity and non-uniformity have become increasingly prominent, and the simulation accuracy of wind energy forecasting and pollution diffusion in complex terrain and urban areas has significantly decreased.

[0003] Although high-resolution large eddy simulation (LES) and computational fluid dynamics (CFD) models can simulate local flow fields in detail, their simulation ranges are limited and the computational costs are extremely high, making it difficult to meet the needs of large-scale forecasting. At the same time, the surface changes brought about by rapid urbanization are difficult to update input data in a timely manner, thus also limiting the application of traditional models in dynamic environments. Summary of the Invention

[0004] The present invention provides a wind power distribution forecasting method and system applicable to complex terrain, which solves the problem of insufficient forecasting accuracy of existing mesoscale meteorological models under complex terrain and ground object conditions.

[0005] To achieve the above object, the technical solution adopted by the present invention is:

[0006] The first aspect of the present invention provides a wind power distribution forecasting method applicable to complex terrain, including:

[0007] Dividing the target area into a number of grid cells, and generating terrain grid point data according to the terrain distribution characteristics and ground object spatial distribution characteristics in each grid cell;

[0008] Using a mesoscale meteorological model to simulate the spatial distribution of low-altitude wind direction, wind speed, and atmospheric turbulent kinetic energy in the target area, and establishing the primary distribution characteristics of the wind field from the spatial distribution of low-altitude wind direction and wind speed;

[0009] Obtaining the historical wind speed observation data of the target area, using a fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area; combining the terrain grid point data and the spatio-temporally non-uniform wind profile in the grid cell to optimize and adjust the primary distribution characteristics of the wind field to form the refined distribution characteristics of the wind field;

[0010] Calculating a sampling frequency limit parameter based on the historical wind speed observation data of the target area, and calculating the total turbulent kinetic energy of the full spectrum from the simulated atmospheric turbulent kinetic energy of the target area using the sampling frequency limit parameter;

[0011] Substitute the refined distribution characteristics of the wind field and the total turbulent kinetic energy in the full spectral band into the mesoscale meteorological model simulation process of the target area to obtain the wind power distribution characteristics.

[0012] Furthermore, use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area, specifically including:

[0013] When the wind in the network cell flows vertically in, use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the first wind profile function at the wind inlet and outlet positions and the second wind profile function at the intermediate position; the vertical inflow wind profile is composed of the first wind profile function and the second wind profile function;

[0014] When the wind in the network cell flows in parallel, use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the third wind profile function at the wind inlet position, the fourth wind profile function at the wind outlet position, and the fifth wind profile function at the intermediate position; the parallel inflow wind profile is composed of the third wind profile function, the fourth wind profile function, and the fifth wind profile function;

[0015] The spatio-temporally non-uniform wind profile of the target area is composed of the vertical inflow wind profile and the parallel inflow wind profile.

[0016] Furthermore, use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the first wind profile function at the wind inlet and outlet positions, specifically including:

[0017] Construct a three-dimensional coordinate in the target area, and set the coordinate axis directions to be axis direction, axis direction, and axis direction;

[0018] Use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the leeward side axis direction component of the first wind profile at the wind inlet and outlet positions, and the expression formula is:

[0019]

[0020]

[0021] In the formula, is the height variable, is the height of the first segmentation point, is the height of the second segmentation point, is the height of the third segmentation point, is the leeward side axis direction component function of the first wind profile; , , and is the leeward side of the first wind profile fitting parameters in the axial direction component function; 、 、 is the leeward side of the first wind profile height constant term in the axial direction component function; and is the leeward side of the first wind profile wind speed constant term in the axial direction component function; is an exponential function; is the initial wind speed; is the ratio of the height variable to the surface height; is the height ratio to the surface height; is the height ratio to the surface height; is the height ratio to the surface height;

[0022] The leeward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using a fluid dynamics algorithm axial direction component, and the expression formula is:

[0023]

[0024]

[0025] In the formula, 、 、 are the fitting parameters of the axial direction component on the leeward side of the first wind profile ; and are the height constant terms of the axial direction component on the leeward side of the first wind profile ; is the leeward side of the first wind profile wind speed constant term of the axial direction component; is the leeward side of the first wind profile axial direction component function;

[0026] The leeward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using a fluid dynamics algorithm axial direction component, and the expression formula is:

[0027]

[0028]

[0029] In the formula, On the leeward side of the first wind profile Axial direction component function; 、 、 and On the leeward side of the first wind profile Fitting parameters of the axial direction component; 、 、 and On the leeward side of the first wind profile Height constant term of the axial direction component; On the leeward side of the first wind profile Wind speed constant term of the axial direction component;

[0030] The axial direction component on the leeward side of the first wind profile Axial direction component, the axial direction component on the leeward side of the first wind profile Axial direction component and the axial direction component on the leeward side of the first wind profile are fitted to form the function of the leeward side of the first wind profile;

[0031] The axial direction component of the windward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm, and the expression formula is: In the formula,

[0032]

[0033]

[0034] In the formula, is the axial direction component function of the windward side of the first wind profile ; 、 、 and are the fitting parameters in the axial direction component function of the windward side of the first wind profile; ; 、 and are the height constant terms in the axial direction component function of the windward side of the first wind profile; ; is the wind speed constant term in the axial direction component function of the windward side of the first wind profile; ;

[0035] The axial direction component of the windward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm, and the expression formula is: In the formula,

[0036]

[0037]

[0038] In the formula, is the w-axis direction component function on the windward side of the first wind profile; Axis direction component function; , , and are the fitting parameters in the w-axis direction component function on the windward side of the first wind profile; Axis direction component function fitting parameters; , and are the height constant terms in the w-axis direction component function on the windward side of the first wind profile; Height constant term; and are the wind speed constant terms in the w-axis direction component function on the windward side of the first wind profile; Wind speed constant term;

[0039] The w-axis direction component on the windward side of the first wind profile is obtained by fitting the historical wind speed observation data using the fluid dynamics algorithm. The expression formula is:

[0040]

[0041]

[0042] In the formula, is the w-axis direction component function on the windward side of the first wind profile; Axis direction component function; and are the fitting parameters of the w-axis direction component function on the windward side of the first wind profile; Axis direction component function fitting parameters; and are the height constant terms of the w-axis direction component function on the windward side of the first wind profile; Height constant term; and are the wind speed constant terms of the w-axis direction component function on the windward side of the first wind profile; Wind speed constant term;

[0043] The w-axis direction component on the windward side of the first wind profile, the w-axis direction component on the windward side of the first wind profile, and the w-axis direction component on the windward side of the first wind profile are fitted to form the function on the windward side of the first wind profile; Axis direction component, Axis direction component, and Axis direction component on the windward side of the first wind profile are fitted to form the function on the windward side of the first wind profile;

[0044] The first wind profile function is composed of the function on the leeward side of the first wind profile and the function on the windward side of the first wind profile.

[0045] Furthermore, the second wind profile function at the intermediate position is obtained by fitting the historical wind speed observation data using the fluid dynamics algorithm. Specifically, it includes:

[0046] The second wind profile leeward side at the middle position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm axial direction component, and the expression formula is:

[0047]

[0048]

[0049] In the formula, is the axial direction component function of the second wind profile leeward side ; , , and are the fitting parameters in the axial direction component function of the second wind profile leeward side ; , and are the height constant terms in the axial direction component function of the second wind profile leeward side ; and are the wind speed constant terms in the axial direction component function of the second wind profile leeward side ;

[0050] The second wind profile leeward side at the middle position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm axial direction component, and the expression formula is:

[0051]

[0052]

[0053] In the formula, , , are the fitting parameters of the axial direction component of the second wind profile leeward side ; and are the height constant terms of the axial direction component of the second wind profile leeward side ; is the wind speed constant term of the axial direction component of the second wind profile leeward side ; is the axial direction component function of the second wind profile leeward side ;

[0054] The second wind profile leeward side at the middle position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm axial direction component, and the expression formula is:

[0055]

[0056]

[0057] In the formula, 、 、 and are the fitting parameters of the leeward side axial direction component of the second wind profile; 、 、 and are the height constant terms of the leeward side axial direction component of the second wind profile; is the wind speed constant term of the leeward side axial direction component of the second wind profile; is the leeward side of the second wind profile axial direction component function;

[0058] The leeward side of the second wind profile axial direction component, the leeward side of the second wind profile axial direction component, and the leeward side of the second wind profile axial direction component are fitted to form the leeward side function of the second wind profile;

[0059] Using the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the axial direction component of the windward side of the second wind profile at the intermediate position, and the expression formula is:

[0060]

[0061]

[0062] In the formula, 、 、 and are the fitting parameters of the windward side axial direction component of the second wind profile; 、 and are the height constant terms of the windward side axial direction component of the second wind profile; and are the wind speed constant terms of the windward side axial direction component of the second wind profile; is the windward side of the second wind profile axial direction component function;

[0063] Using the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the windward side of the second wind profile at the intermediate position The axial direction component, and its expression formula is:

[0064]

[0065]

[0066] In the formula, , , and are the fitting parameters of the axial direction component on the windward side of the second wind profile; The height constant term of the axial direction component on the windward side of the second wind profile; , and are the height constant terms of the axial direction component on the windward side of the second wind profile; The height constant term of the axial direction component on the windward side of the second wind profile; and are the wind speed constant terms of the axial direction component on the windward side of the second wind profile; The wind speed constant term of the axial direction component on the windward side of the second wind profile; is the function of the axial direction component on the windward side of the second wind profile; The function of the axial direction component on the windward side of the second wind profile;

[0067] The axial direction component on the windward side of the second wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, and its expression formula is: The axial direction component, and its expression formula is:

[0068]

[0069]

[0070] In the formula, , and are the fitting parameters of the axial direction component on the windward side of the second wind profile; The height constant term of the axial direction component on the windward side of the second wind profile; , and are the height constant terms of the axial direction component on the windward side of the second wind profile; The height constant term of the axial direction component on the windward side of the second wind profile; and are the wind speed constant terms of the axial direction component on the windward side of the second wind profile; The wind speed constant term of the axial direction component on the windward side of the second wind profile; is the function of the axial direction component on the windward side of the second wind profile; The function of the axial direction component on the windward side of the second wind profile;

[0071] The axial direction component on the windward side of the second wind profile, the axial direction component on the windward side of the second wind profile, and the axial direction component on the windward side of the second wind profile are fitted to form the function of the windward side of the second wind profile; The axial direction component on the windward side of the second wind profile, The axial direction component on the windward side of the second wind profile, and The axial direction component on the windward side of the second wind profile are fitted to form the function of the windward side of the second wind profile;

[0072] The second wind profile function is composed of the second leeward side function of the wind profile and the second windward side function of the wind profile.

[0073] Furthermore, the third wind profile function at the air inlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, specifically including:

[0074] The third leeward side of the wind profile at the air inlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0075]

[0076]

[0077] In the formula, 、 and are the fitting parameters of the axial direction component of the third leeward side of the wind profile ; and are the height constant terms of the axial direction component of the third leeward side of the wind profile ; is the wind speed constant term of the axial direction component of the third leeward side of the wind profile ; is the function of the axial direction component of the third leeward side of the wind profile ;

[0078] The third windward side of the wind profile at the air inlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0079]

[0080]

[0081] In the formula, 、 and are the fitting parameters of the axial direction component of the third windward side of the wind profile ; and are the height constant terms of the axial direction component of the third windward side of the wind profile ; is the wind speed constant term of the axial direction component of the third windward side of the wind profile ; is the function of the axial direction component of the third windward side of the wind profile ;

[0082] The third wind profile at the wind inlet position is obtained by fitting the historical wind speed observation data using a fluid mechanics algorithm. The axial component is expressed as:

[0083]

[0084]

[0085] In the formula, , and The third wind profile Fitting parameters of the axial components; and The third wind profile The height constant term of the axial component; The third wind profile Axis direction component function;

[0086] The third wind profile at the wind inlet position is obtained by fitting the historical wind speed observation data using a fluid mechanics algorithm. The axial component is expressed as:

[0087]

[0088]

[0089] In the formula, and The third wind profile Fitting parameters of the axial components; and The third wind profile The height constant term of the axial component; The third wind profile The wind speed constant term of the axial component; The third wind profile Axis direction component function;

[0090] From the leeward side of the third wind profile Axial component function, windward side of the third wind profile Axial component function, third wind profile Axial component function and third wind profile The axial component function is fitted to form the third wind profile function.

[0091] Furthermore, the fourth wind profile function of the wind outlet position is obtained by fitting the historical wind speed observation data using a fluid mechanics algorithm, specifically including:

[0092] The fourth wind profile leeward side of the air outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0093]

[0094]

[0095] In the formula, and are the fitting parameters of the axial direction component of the fourth wind profile leeward side ; and are the height constant terms of the axial direction component of the fourth wind profile leeward side ; is the wind speed constant term of the axial direction component of the fourth wind profile leeward side ; is the axial direction component function of the fourth wind profile leeward side ;

[0096] The fourth wind profile windward side of the air outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0097]

[0098]

[0099] In the formula, and are the fitting parameters of the axial direction component of the fourth wind profile windward side ; and are the height constant terms of the axial direction component of the fourth wind profile windward side ; is the wind speed constant term of the axial direction component of the fourth wind profile windward side ; is the axial direction component function of the fourth wind profile windward side ;

[0100] The fourth wind profile of the air outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0101]

[0102]

[0103] In the formula, 、 and is the fitting parameter of the axial direction component of the fourth wind profile ; and is the height constant term of the axial direction component of the fourth wind profile ; is the axial direction component function of the fourth wind profile ;

[0104] The fourth wind profile at the wind outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0105]

[0106]

[0107] In the formula, and is the fitting parameter of the axial direction component of the fourth wind profile ; and is the height constant term of the axial direction component of the fourth wind profile ; and is the wind speed constant term of the axial direction component of the fourth wind profile ; is the axial direction component function of the fourth wind profile ;

[0108] The fourth wind profile function is composed of the axial direction component function on the leeward side of the fourth wind profile the axial direction component function on the windward side of the fourth wind profile the axial direction component function of the fourth wind profile and the axial direction component function of the fourth wind profile ;

[0109] Furthermore, the fifth wind profile function at the intermediate position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, specifically including:

[0110] The axial direction component on the leeward side of the fifth wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm and the expression formula is:

[0111]

[0112]

[0113] In the formula, and On the leeward side of the fifth wind profile Fitting parameters for the axial direction component; and On the leeward side of the fifth wind profile Height constant term for the axial direction component of the fifth wind profile; On the leeward side of the fifth wind profile Wind speed constant term for the axial direction component of the fifth wind profile; On the leeward side of the fifth wind profile Axial direction component function of the fifth wind profile;

[0114] The axial direction component of the fifth wind profile on the windward side at the intermediate position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, and the expression formula is: The axial direction component, and the expression formula is:

[0115]

[0116]

[0117] In the formula, and On the windward side of the fifth wind profile Fitting parameters for the axial direction component; and On the windward side of the fifth wind profile Height constant term for the axial direction component of the fifth wind profile; On the windward side of the fifth wind profile Wind speed constant term for the axial direction component of the fifth wind profile; On the windward side of the fifth wind profile Axial direction component function of the fifth wind profile;

[0118] The axial direction component of the fifth wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, and the expression formula is: The axial direction component, and the expression formula is:

[0119]

[0120]

[0121] In the formula, 、 and For the axial direction component of the fifth wind profile Fitting parameters; and For the fifth wind profile Height constant term for the axial direction component; For the fifth wind profile Axial direction component function;

[0122] The historical wind speed observation data is fitted using a hydrodynamic algorithm to obtain the fifth wind profile at the intermediate position. The axial direction component, and the expression formula is:

[0123]

[0124]

[0125] In the formula, and are the fitting parameters of the axial direction component of the fifth wind profile; are the height constant terms of the axial direction component of the fifth wind profile; and are the wind speed constant terms of the axial direction component of the fifth wind profile; is the function of the axial direction component of the fifth wind profile; and are the functions of the axial direction component on the leeward side of the fifth wind profile, the functions of the axial direction component on the windward side of the fifth wind profile, the functions of the axial direction component of the fifth wind profile, and the functions of the axial direction component of the fifth wind profile. The functions of the axial direction components on the leeward side, windward side, and at different positions of the fifth wind profile are fitted to form the fifth wind profile function. is the function of the axial direction component of the fifth wind profile; is the fifth wind profile axis direction component function;

[0126] It is composed of the function of the axial direction component on the leeward side of the fifth wind profile, the function of the axial direction component on the windward side of the fifth wind profile, the function of the axial direction component of the fifth wind profile, and the function of the axial direction component of the fifth wind profile. The function of the axial direction component on the leeward side, windward side, and at different positions of the fifth wind profile is fitted to form the fifth wind profile function. axis direction component function, the function of the axial direction component on the windward side of the fifth wind profile axis direction component function, the fifth wind profile axis direction component function and the fifth wind profile axis direction component function fit to form the fifth wind profile function.

[0127] Furthermore, the sampling frequency limit parameter is calculated based on the historical wind speed observation data of the target area, and the total turbulent kinetic energy of the full spectrum is calculated from the simulated atmospheric turbulent kinetic energy of the target area using the sampling frequency limit parameter. Specifically, it includes:

[0128] The sampling frequency limit parameter is calculated based on the historical wind speed observation data of the target area, and the expression formula is:

[0129]

[0130]

[0131]

[0132] In the formula, is the historical turbulent kinetic energy of the observation spectrum; is the total historical turbulent kinetic energy of the full spectrum; is the turbulent energy spectrum in the historical wind speed observation data; is the set constant term; is the wave number corresponding to the response frequency of the observation instrument; is the wave number corresponding to the turbulent dissipation scale; k is the variable wave number; is the wave number corresponding to the largest-scale turbulent eddies in the atmosphere; is the sampling frequency limit parameter;

[0133] The total turbulent kinetic energy of the entire spectral band is calculated from the turbulent kinetic energy of the atmosphere simulated in the target area using the sampling frequency limit parameter, and the expression formula is:

[0134]

[0135] In the formula, is the turbulent kinetic energy of the atmosphere simulated in the target area; is the total turbulent kinetic energy of the entire spectral band in the target area.

[0136] The second aspect of the present invention provides an air pollution diffusion prediction method applicable to complex terrains, including:

[0137] Adopting the wind power distribution prediction method described in the first aspect to generate wind power distribution characteristics;

[0138] Obtaining the background pollution concentration and pollutant emission data of the target area, and combining the wind power distribution to calculate the pollution diffusion speed and the spatio-temporal distribution characteristics of the local pollutant concentration.

[0139] The third aspect of the present invention provides a wind power distribution prediction system applicable to complex terrains, including:

[0140] A surface division unit that divides the target area into a number of grid cells, and generates surface grid point data according to the terrain distribution characteristics and the spatial distribution characteristics of ground objects in each grid cell;

[0141] A simulation unit that uses a mesoscale meteorological model to simulate the spatial distribution of low-altitude wind direction, wind speed, and turbulent kinetic energy of the atmosphere in the target area, and establishes a primary wind field distribution characteristic from the spatial distribution of low-altitude wind direction and wind speed;

[0142] A wind field optimization unit that obtains the historical wind speed observation data of the target area, and uses a fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area; in the grid cell, combines the surface grid point data and the spatio-temporally non-uniform wind profile to optimize and adjust the primary wind field distribution characteristic to form a refined wind field distribution characteristic;

[0143] A turbulence optimization unit that calculates the sampling frequency limit parameter according to the historical wind speed observation data of the target area, and calculates the total turbulent kinetic energy of the entire spectral band from the turbulent kinetic energy of the atmosphere simulated in the target area using the sampling frequency limit parameter;

[0144] An output unit that substitutes the refined wind field distribution characteristic and the total turbulent kinetic energy of the entire spectral band into the mesoscale meteorological model simulation process of the target area to obtain the wind power distribution characteristic.

[0145] Advantages of the present invention compared with the prior art:

[0146] In the present invention, by dividing the target area into several network cells and generating ground table point data in combination with the terrain and ground feature characteristics in each cell, the spatial distribution characteristics of complex terrain and ground features can be more accurately reflected. By combining the ground table point data and the spatio-temporally non-uniform wind profile in the network cell to optimize and adjust the primary distribution characteristics of the wind field to form the refined distribution characteristics of the wind field, the influence of complex terrain on the wind field is simulated, thereby improving the accuracy of wind field simulation.

[0147] In the present invention, the sampling frequency limit parameter is calculated based on the historical wind speed observation data of the target area, and the total turbulent kinetic energy of the full spectrum is calculated from the turbulent kinetic energy of the simulated atmospheric turbulence in the target area using the sampling frequency limit parameter. By extending from the turbulent kinetic energy of the observed spectrum to the total turbulent kinetic energy of the full spectrum, the distribution of turbulent energy can be considered more comprehensively, and the accuracy of turbulent kinetic energy simulation is significantly improved.

[0148] Substituting the optimized refined distribution characteristics of the wind field and the total turbulent kinetic energy of the full spectrum into the mesoscale meteorological model in the present invention can more accurately simulate the wind power distribution, better adapt to complex terrain conditions, and improve the accuracy and reliability of wind power distribution simulation. Description of the Drawings

[0149] Figure 1 is a flowchart of a wind power distribution prediction method applicable to complex terrain provided in Embodiment 1 of the present invention;

[0150] Figure 2 is a wind profile characteristic diagram under complex terrain and ground conditions provided in Embodiment 1 of the present invention;

[0151] Figure 3 is a comparison diagram of the simulated value and the observed value of the surface sensible heat flux provided in Embodiment 1 of the present invention;

[0152] Figure 4 is a wind energy power generation distribution diagram provided in Embodiment 1 of the present invention;

[0153] Figure 5 is a time series diagram of the simulated value of the atmospheric turbulence diffusion coefficient provided in Embodiment 2 of the present invention;

[0154] Figure 6 is PM provided in Embodiment 2 of the present invention 2.5 comparison diagram of the simulated value and the observation of the concentration. Detailed Embodiments

[0155] The present invention will be further described below in conjunction with the accompanying drawings. The following embodiments are only used to more clearly illustrate the technical solutions of the present invention, and cannot be used to limit the protection scope of the present invention.

[0156] Embodiment 1

[0157] As Figure 1 shown, this embodiment provides a wind power distribution prediction method applicable to complex terrain, including:

[0158] Dividing the target area into a number of network cells, and generating ground grid point data according to the terrain distribution characteristics and the spatial distribution characteristics of ground objects in each network cell;

[0159] Using a mesoscale meteorological model to simulate the spatial distribution of low-altitude wind direction, wind speed and atmospheric turbulent kinetic energy in the target area, and establishing the primary distribution characteristics of the wind field from the spatial distribution of low-altitude wind direction and wind speed;

[0160] Obtaining the historical wind speed observation data of the target area, and using the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area, specifically including:

[0161] As Figure 2 shown, when the wind in the network cell flows vertically in, using the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the first wind profile function at the wind inlet and outlet positions and the second wind profile function at the intermediate position; The vertical inflow wind profile is composed of the first wind profile function and the second wind profile function;

[0162] When the wind in the network cell flows in parallel, using the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the third wind profile function at the wind inlet position, the fourth wind profile function at the wind outlet position and the fifth wind profile function at the intermediate position; The parallel inflow wind profile is composed of the third wind profile function, the fourth wind profile function and the fifth wind profile function;

[0163] The spatio-temporally non-uniform wind profile of the target area is composed of the vertical inflow wind profile and the parallel inflow wind profile. This embodiment can more accurately reflect the wind field structure under complex terrain and ground object conditions, and provide a more scientific basis for wind energy resource assessment and wind farm site selection. By optimizing and adjusting the wind field distribution characteristics, the spatio-temporal changes of wind speed can be better captured, thereby improving the accuracy of wind power prediction.

[0164] Among them, using the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the first wind profile function at the wind inlet and outlet positions, specifically including:

[0165] Constructing a three-dimensional coordinate in the target area, and setting the axis directions to be divided into axis direction, axis direction and Axis direction; set the first segmentation point, the second segmentation point, and the third segmentation point at different heights, and divide the atmospheric environment into the near-ground section, the street canyon interior section, the building roof section, and the free atmosphere section above the roof through the first segmentation point, the second segmentation point, and the third segmentation point;

[0166] Use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the leeward side of the first wind profile at the inlet and outlet positions Axis direction component, and the expression formula is:

[0167]

[0168]

[0169] In the formula, is the height variable, is the height of the first segmentation point, is the height of the second segmentation point, is the height of the third segmentation point, is the leeward side of the first wind profile Axis direction component function; , , and are the fitting parameters in the leeward side of the first wind profile Axis direction component function; , , are the height constant terms in the leeward side of the first wind profile Axis direction component function; and are the wind speed constant terms in the leeward side of the first wind profile Axis direction component function; is the exponential function; is the initial wind speed; is the ratio of the height variable to the surface height; is the height Ratio to the surface height; is the height Ratio to the surface height; is the height Ratio to the surface height; where the surface height is the height of the terrain and the spatial distribution of ground objects;

[0170] Use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the leeward side of the first wind profile at the inlet and outlet positions Axis direction component, and the expression formula is:

[0171]

[0172]

[0173] In the formula, , , are the fitting parameters of the leeward side of the first wind profile in the axis direction component; and are the height constant terms of the leeward side of the first wind profile in the axis direction component; is the wind speed constant term of the leeward side of the first wind profile in the axis direction component; is the function of the leeward side of the first wind profile in the axis direction component;

[0174] The leeward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using the fluid dynamics algorithm for the axis direction component, and the expression formula is:

[0175]

[0176]

[0177] In the formula, is the function of the leeward side of the first wind profile in the axis direction component; , , and are the fitting parameters of the leeward side of the first wind profile in the axis direction component; , , and are the height constant terms of the leeward side of the first wind profile in the axis direction component; is the wind speed constant term of the leeward side of the first wind profile in the axis direction component;

[0178] The axis direction component of the leeward side of the first wind profile, the axis direction component of the leeward side of the first wind profile, and the axis direction component of the leeward side of the first wind profile are fitted to form the function of the leeward side of the first wind profile;

[0179] The windward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using the fluid dynamics algorithm for the axis direction component, and the expression formula is:

[0180]

[0181]

[0182] In the formula, is the function of the axial direction component on the windward side of the first wind profile; is the axial direction component function; , , and are the fitting parameters in the function of the axial direction component on the windward side of the first wind profile; is the fitting parameter in the function of the axial direction component on the windward side of the first wind profile; , and are the height constant terms in the function of the axial direction component on the windward side of the first wind profile; is the height constant term in the function of the axial direction component on the windward side of the first wind profile; is the wind speed constant term in the function of the axial direction component on the windward side of the first wind profile; is the wind speed constant term in the function of the axial direction component on the windward side of the first wind profile;

[0183] The axial direction component on the windward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using the fluid dynamics algorithm. The expression formula is: The expression formula is:

[0184]

[0185]

[0186] In the formula, is the function of the axial direction component on the windward side of the first wind profile; is the axial direction component function; , , and are the fitting parameters in the function of the axial direction component on the windward side of the first wind profile; is the fitting parameter in the function of the axial direction component on the windward side of the first wind profile; , and are the height constant terms in the function of the axial direction component on the windward side of the first wind profile; is the height constant term in the function of the axial direction component on the windward side of the first wind profile; and are the wind speed constant terms in the function of the axial direction component on the windward side of the first wind profile; is the wind speed constant term in the function of the axial direction component on the windward side of the first wind profile;

[0187] The axial direction component on the windward side of the first wind profile at the inlet and outlet positions is obtained by fitting the historical wind speed observation data using the fluid dynamics algorithm. The expression formula is:

[0188]

[0189]

[0190] In the formula, The windward side of the first wind profile Axis component function; and The windward side of the first wind profile Fitting parameters of the axial component function; and The windward side of the first wind profile The height constant term of the axial component function; and The windward side of the first wind profile The wind speed constant term of the axial component function;

[0191] The first wind profile is on the windward side Axial component, windward side of the first wind profile Axial component and windward side of the first wind profile The axial direction component is fitted to form the windward side function of the first wind profile;

[0192] The first wind profile function is composed of the first wind profile leeward side function and the first wind profile windward side function.

[0193] Among them, the fluid mechanics algorithm is used to fit the historical wind speed observation data to obtain the second wind profile function at the middle position, specifically including:

[0194] The fluid dynamics algorithm is used to fit the historical wind speed observation data to obtain the leeward side of the second wind profile at the middle position The axial component is expressed as:

[0195]

[0196]

[0197] In the formula, The leeward side of the second wind profile Axis component function; , , and The leeward side of the second wind profile Fitting parameters in the axial component function; , and The leeward side of the second wind profile The height constant term in the axial component function; and The leeward side of the second wind profile The wind speed constant term in the axial component function;

[0198] The historical wind speed observation data is fitted using a hydrodynamic algorithm to obtain the leeward side of the second wind profile at the intermediate position axial direction component, and the expression formula is:

[0199]

[0200]

[0201] In the formula, 、 、 are the fitting parameters of the axial direction component of the leeward side of the second wind profile; and and are the height constant terms of the axial direction component of the leeward side of the second wind profile; axial direction component, and the expression formula is: is the wind speed constant term of the axial direction component of the leeward side of the second wind profile; axial direction component; is the function of the axial direction component of the leeward side of the second wind profile; axial direction component, and the expression formula is:

[0202] The historical wind speed observation data is fitted using a hydrodynamic algorithm to obtain the leeward side of the second wind profile at the intermediate position axial direction component, and the expression formula is:

[0203]

[0204]

[0205] In the formula, 、 、 and are the fitting parameters of the axial direction component of the leeward side of the second wind profile; axial direction component, and the expression formula is: 、 、 and are the height constant terms of the axial direction component of the leeward side of the second wind profile; axial direction component, and the expression formula is: is the wind speed constant term of the axial direction component of the leeward side of the second wind profile; axial direction component; is the function of the axial direction component of the leeward side of the second wind profile; axial direction component;

[0206] The axial direction component of the leeward side of the second wind profile, the axial direction component of the leeward side of the second wind profile, and the axial direction component of the leeward side of the second wind profile are fitted to form the function of the leeward side of the second wind profile; axial direction component, the axial direction component of the leeward side of the second wind profile axial direction component, and the axial direction component of the leeward side of the second wind profile axial direction component are fitted to form the function of the leeward side of the second wind profile;

[0207] Obtain the second wind profile on the windward side at the intermediate position by fitting the historical wind speed observation data using the fluid dynamics algorithm The axial direction component, and the expression formula is:

[0208]

[0209]

[0210] In the formula, , , and are the fitting parameters of the axial direction component of the second wind profile on the windward side ; , and are the height constant terms of the axial direction component of the second wind profile on the windward side ; and are the wind speed constant terms of the axial direction component of the second wind profile on the windward side ; is the function of the axial direction component of the second wind profile on the windward side ;

[0211] Obtain the second wind profile on the windward side at the intermediate position by fitting the historical wind speed observation data using the fluid dynamics algorithm The axial direction component, and the expression formula is:

[0212]

[0213]

[0214] In the formula, , , and are the fitting parameters of the axial direction component of the second wind profile on the windward side ; , and are the height constant terms of the axial direction component of the second wind profile on the windward side ; and are the wind speed constant terms of the axial direction component of the second wind profile on the windward side ; is the function of the axial direction component of the second wind profile on the windward side ;

[0215] Obtain the second wind profile on the windward side at the intermediate position by fitting the historical wind speed observation data using the fluid dynamics algorithm The axial direction component, and the expression formula is:

[0216]

[0217]

[0218] In the formula, 、 and are the fitting parameters of the axial direction component on the windward side of the second wind profile; 、 、 and are the height constant terms of the axial direction component on the windward side of the second wind profile; and are the wind speed constant terms of the axial direction component on the windward side of the second wind profile; is the axial direction component function on the windward side of the second wind profile;

[0219] The axial direction component on the windward side of the second wind profile, the axial direction component on the windward side of the second wind profile, and the axial direction component on the windward side of the second wind profile are fitted to form the function on the windward side of the second wind profile;

[0220]

[0221]

[0222]

[0223]

[0224]

[0225] In the formula, 、 and are the fitting parameters of the axial direction component on the leeward side of the third wind profile; and are the height constant terms of the axial direction component on the leeward side of the third wind profile; is the wind speed constant term of the axial direction component on the leeward side of the third wind profile; ​​​On the leeward side of the third wind profile Axial direction component function;

[0226] Using the hydrodynamics algorithm to fit the historical wind speed observation data to obtain the windward side of the third wind profile at the incoming wind position Axial direction component, the expression formula is:

[0227]

[0228]

[0229] In the formula, 、 and Are the fitting parameters of the axial direction component on the windward side of the third wind profile And and Are the height constant terms of the axial direction component on the windward side of the third wind profile Axial direction component height constant term; Is the axial direction component on the windward side of the third wind profile Axial direction component wind speed constant term; Is the axial direction component on the windward side of the third wind profile Axial direction component function;

[0230] Using the hydrodynamics algorithm to fit the historical wind speed observation data to obtain the third wind profile at the incoming wind position Axial direction component, the expression formula is:

[0231]

[0232]

[0233] In the formula, 、 and Are the fitting parameters of the axial direction component of the third wind profile Axial direction component fitting parameter; and Are the height constant terms of the axial direction component of the third wind profile Axial direction component height constant term; Is the axial direction component of the third wind profile Axial direction component function;

[0234] Using the hydrodynamics algorithm to fit the historical wind speed observation data to obtain the third wind profile at the incoming wind position Axial direction component, the expression formula is:

[0235]

[0236]

[0237] In the formula, and are the fitting parameters of the axial direction component of the third wind profile; axis direction component of the third wind profile; and are the height constant terms of the axial direction component of the third wind profile; axis direction component of the third wind profile; is the wind speed constant term of the axial direction component of the third wind profile; axis direction component of the third wind profile; is the function of the axial direction component of the third wind profile; axis direction component of the third wind profile;

[0238] The third wind profile function is composed of the axial direction component function on the leeward side of the third wind profile, the axial direction component function on the windward side of the third wind profile, the axial direction component function of the third wind profile, and the axial direction component function of the third wind profile. axis direction component of the third wind profile; axis direction component of the third wind profile; axis direction component of the third wind profile; axis direction component of the third wind profile;

[0239] Among them, the fourth wind profile function at the air outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, specifically including:

[0240] The axial direction component on the leeward side of the fourth wind profile at the air outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, and the expression formula is: axis direction component of the fourth wind profile;

[0241]

[0242]

[0243] In the formula, and are the fitting parameters of the axial direction component on the leeward side of the fourth wind profile; axis direction component of the fourth wind profile; and are the height constant terms of the axial direction component on the leeward side of the fourth wind profile; axis direction component of the fourth wind profile; is the wind speed constant term of the axial direction component on the leeward side of the fourth wind profile; axis direction component of the fourth wind profile; is the function of the axial direction component on the leeward side of the fourth wind profile; axis direction component of the fourth wind profile;

[0244] The axial direction component on the windward side of the fourth wind profile at the air outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm, and the expression formula is: axis direction component of the fourth wind profile;

[0245]

[0246]

[0247] In the formula, and are the fitting parameters of the axial direction component on the windward side of the fourth wind profile; The height constant term of the axial direction component on the windward side of the fourth wind profile; and are the height constant terms of the axial direction component on the windward side of the fourth wind profile; The height constant term of the axial direction component on the windward side of the fourth wind profile; are the height constant terms of the axial direction component on the windward side of the fourth wind profile; The wind speed constant term of the axial direction component on the windward side of the fourth wind profile; are the height constant terms of the axial direction component on the windward side of the fourth wind profile; The function of the axial direction component on the windward side of the fourth wind profile;

[0248] The fourth wind profile of the outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0249]

[0250]

[0251] In the formula, , and are the fitting parameters of the axial direction component of the fourth wind profile; The fitting parameters of the axial direction component of the fourth wind profile; and are the height constant terms of the axial direction component of the fourth wind profile; The height constant term of the axial direction component of the fourth wind profile; are the height constant terms of the axial direction component of the fourth wind profile; The function of the axial direction component of the fourth wind profile;

[0252] The fourth wind profile of the outlet position is obtained by fitting the historical wind speed observation data using the fluid mechanics algorithm The axial direction component, and the expression formula is:

[0253]

[0254]

[0255] In the formula, and are the fitting parameters of the axial direction component of the fourth wind profile; The fitting parameters of the axial direction component of the fourth wind profile; and are the height constant terms of the axial direction component of the fourth wind profile; The height constant term of the axial direction component of the fourth wind profile; and are the height constant terms of the axial direction component of the fourth wind profile; The wind speed constant term of the axial direction component of the fourth wind profile; are the height constant terms of the axial direction component of the fourth wind profile; Axial direction component function;

[0256] From the leeward side of the fourth wind profile Axial direction component function, windward side of the fourth wind profile Axial direction component function, fourth wind profile Axial direction component function and the fourth wind profile The axial direction component function is fitted to form the fourth wind profile function.

[0257] Among them, the fifth wind profile function at the intermediate position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm, specifically including:

[0258] The leeward side of the fifth wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm Axial direction component, the expression formula is:

[0259]

[0260]

[0261] In the formula, and are the fitting parameters of the axial direction component on the leeward side of the fifth wind profile ; and are the height constant terms of the axial direction component on the leeward side of the fifth wind profile ; is the wind speed constant term of the axial direction component on the leeward side of the fifth wind profile ; is the axial direction component function on the leeward side of the fifth wind profile ;

[0262] The windward side of the fifth wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm Axial direction component, the expression formula is:

[0263]

[0264]

[0265] In the formula, and are the fitting parameters of the axial direction component on the windward side of the fifth wind profile ; and are the height constant terms of the axial direction component on the windward side of the fifth wind profile ; is the wind speed constant term of the axial direction component on the windward side of the fifth wind profile ; For the windward side of the fifth wind profile Axial direction component function;

[0266] The fifth wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm Axial direction component, and the expression formula is:

[0267]

[0268]

[0269] In the formula, 、 and are the fitting parameters of the axial direction component of the fifth wind profile ; and are the height constant terms of the axial direction component of the fifth wind profile ; is the axial direction component function of the fifth wind profile ;

[0270] The fifth wind profile at the intermediate position is obtained by fitting the historical wind speed observation data using the hydrodynamics algorithm Axial direction component, and the expression formula is:

[0271]

[0272]

[0273] In the formula, and are the fitting parameters of the axial direction component of the fifth wind profile ; and are the height constant terms of the axial direction component of the fifth wind profile ; and are the wind speed constant terms of the axial direction component of the fifth wind profile ; is the axial direction component function of the fifth wind profile ;

[0274] From the axial direction component function on the leeward side of the fifth wind profile 、the axial direction component function on the windward side of the fifth wind profile 、the axial direction component function of the fifth wind profile and the axial direction component function of the fifth wind profile are fitted to form the fifth wind profile function

[0275] Optimize and adjust the primary distribution characteristics of the wind field by combining the grid point data in the network cell and the spatio-temporal non-uniform wind profile to form the refined distribution characteristics of the wind field. In this embodiment, the mesoscale meteorological model simulation process in the prior art is used as the control group, and the improved simulation process disclosed in this embodiment is used as the experimental group. As Figure 3 shown, by comparing the simulated sensible heat flux values of the control group and the experimental group with the observed values, the rationality of the wind profile function fitted in this embodiment can be verified.

[0276] In this embodiment, by fitting the historical wind speed observation data using the hydrodynamic algorithm, the spatio-temporal non-uniform wind profile of the target area can be accurately obtained. The influence of complex terrain and ground objects on the wind field is fully considered, so that the wind profile can more realistically reflect the variation characteristics of the actual wind speed with height and time. Combining the grid point data in the network cell to optimize and adjust the primary distribution characteristics of the wind field obtained by the mesoscale meteorological model simulation can form a more refined wind field distribution characteristic, significantly improving the accuracy and reliability of the wind field simulation.

[0277] Calculate the sampling frequency limit parameter based on the historical wind speed observation data of the target area, and calculate the total turbulent kinetic energy of the full spectrum from the simulated atmospheric turbulent kinetic energy of the target area using the sampling frequency limit parameter, specifically including:

[0278] Calculate the sampling frequency limit parameter based on the historical wind speed observation data of the target area, and the expression formula is:

[0279]

[0280]

[0281]

[0282] In the formula, is the historical turbulent kinetic energy of the observed spectrum segment; is the total historical turbulent kinetic energy of the full spectrum; is the turbulent energy spectrum in the historical wind speed observation data; is the set constant term; is the wave number corresponding to the response frequency of the observation instrument, and its frequency is usually between 0.025 - 2 Hz; is the wave number corresponding to the turbulent dissipation scale, and its frequency is about 1000 Hz; is the variable wave number; is the wave number corresponding to the largest scale turbulent eddy in the atmosphere; is the sampling frequency limit parameter;

[0283] Calculate the total turbulent kinetic energy of the full spectrum from the simulated atmospheric turbulent kinetic energy of the target area using the sampling frequency limit parameter, and the expression formula is:

[0284]

[0285] In the formula, is the atmospheric turbulent kinetic energy simulated for the target area; is the total turbulent kinetic energy of the full spectral band in the target area.

[0286] Traditional methods can only capture low-frequency turbulent kinetic energy, ignoring the contribution of high-frequency turbulence. By introducing a sampling frequency limit parameter in this embodiment, this deficiency can be effectively compensated, thus more comprehensively reflecting the distribution characteristics of turbulent energy. This embodiment not only improves the simulation ability of turbulent kinetic energy under complex terrain and surface conditions, but also provides a more accurate turbulent flow field input for wind power prediction and atmospheric pollution diffusion simulation.

[0287] Substitute the refined distribution characteristics of the wind field and the total turbulent kinetic energy of the full spectral band into the mesoscale meteorological model simulation process of the target area to obtain the wind power distribution characteristics. The wind power distribution characteristics consist of two parts: the wind kinetic energy contributed by the mean wind speed and the wind kinetic energy contributed by the turbulent pulsation.

[0288] Calculate the wind power contributed by the mean wind speed based on the mean wind speed under complex terrain. The expression formula is:

[0289]

[0290] In the formula, is the wind power contributed by the mean wind speed; is the air density, is the mean wind speed, and A is the area swept by the wind.

[0291] The wind power contributed by the turbulent pulsation is determined by the total turbulent kinetic energy of the full spectral band. The expression formula is:

[0292]

[0293] In the formula, is the turbulent pulsation wind speed; is the wind power contributed by the turbulent pulsation;

[0294] Superimpose the wind power contributed by the turbulent pulsation and the wind power contributed by the mean wind speed to obtain the wind power distribution characteristics.

[0295] This embodiment can more accurately predict the wind power distribution. For example, Figure 4 as shown in the wind power generation power distribution map, it provides a more scientific basis for wind energy resource assessment and wind farm site selection, optimizes the wind farm site selection and layout, thereby improving the wind energy utilization efficiency and reducing the development cost. In addition, this embodiment provides more accurate technical support for atmospheric pollution diffusion simulation and air quality forecasting.

[0296] Example 2

[0297] This example provides an air pollution diffusion prediction method applicable to complex terrains, including:

[0298] Using the wind power distribution prediction method described in Example 1 to generate wind power distribution characteristics;

[0299] Obtain the background pollution concentration and pollutant emission data of the target area, and calculate the pollution diffusion speed and the spatio-temporal distribution characteristics of local pollutant concentrations in combination with the wind power distribution, which are used as the experimental group. There is existing technology for calculating the pollution diffusion speed and the spatio-temporal distribution characteristics of local pollutant concentrations from the background pollution concentration, pollutant emission data, and wind power distribution, and this example will not repeat it here.

[0300] Use the original mesoscale meteorological model as the control group to simulate the spatial distribution of low-altitude wind direction, wind speed, and turbulent kinetic energy within the target area; and calculate the atmospheric turbulent diffusion coefficient and pollutant concentration to generate preliminary wind field, turbulent field, and pollutant concentration distribution results.

[0301] Combining the refined wind field and turbulent field simulated by the experimental group and the control group, recalculate the atmospheric turbulent diffusion coefficient to obtain the atmospheric turbulent diffusion coefficient as shown in Figure 5 The calculation formula is as follows:

[0302]

[0303] In the formula, is the atmospheric turbulent diffusion coefficient, is the total turbulent kinetic energy of the full spectral band of the target area, is the characteristic scale used to calculate the turbulent diffusion coefficient, is a constant with a value of 0.4.

[0304] Combining the refined anthropogenic emission data, based on the improved turbulent diffusion coefficient, re-simulate the spatio-temporal distribution of pollutant concentrations under complex terrain conditions and generate air quality prediction results. Process the anthropogenic emission dataset with a resolution of 1 km into anthropogenic source data and input it into the improved mesoscale meteorological prediction model to simulate the atmospheric pollution diffusion process, and obtain the simulated value of PM Figure 6 concentration as shown in 2.5 and compare it with the observation results to prove that this example is more suitable for simulating the air pollution diffusion characteristics under complex terrain and ground features.

[0305] Example 3

[0306] This example provides a wind power distribution prediction system applicable to complex terrains. The wind power distribution prediction system can execute the wind power distribution prediction method described in Example 1. The wind power distribution prediction system includes:

[0307] Surface division unit, which divides the target area into a grid to form a number of network cells, and generates surface grid point data according to the topographic distribution characteristics and the spatial distribution characteristics of ground objects in each network cell;

[0308] Simulation unit, which uses a mesoscale meteorological model to simulate the spatial distribution of low-altitude wind direction, wind speed and atmospheric turbulent kinetic energy in the target area, and establishes the primary distribution characteristics of the wind field from the spatial distribution of low-altitude wind direction and wind speed;

[0309] Wind field optimization unit, which obtains the historical wind speed observation data of the target area, and uses the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area; in the network cell, it combines the surface grid point data and the spatio-temporally non-uniform wind profile to optimize and adjust the primary distribution characteristics of the wind field to form the refined distribution characteristics of the wind field;

[0310] Turbulence optimization unit, which calculates the sampling frequency limit parameter according to the historical wind speed observation data of the target area, and uses the sampling frequency limit parameter to calculate the total turbulent kinetic energy of the full spectrum from the simulated atmospheric turbulent kinetic energy of the target area;

[0311] Output unit, which substitutes the refined distribution characteristics of the wind field and the total turbulent kinetic energy of the full spectrum into the mesoscale meteorological model simulation process of the target area to obtain the wind power distribution characteristics.

[0312] Those skilled in the art should understand that the embodiments of the present application can be provided as a method, a system, or a computer program product. Therefore, the present application can adopt the form of a complete hardware embodiment, a complete software embodiment, or an embodiment combining software and hardware aspects. Moreover, the present application can adopt the form of a computer program product implemented on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0313] The present application is described with reference to the flowcharts and / or block diagrams of methods, devices (systems), and computer program products according to the embodiments of the present application. It should be understood that each process and / or block in the flowchart and / or block diagram can be realized by computer program instructions, and the combination of the processes and / or blocks in the flowchart and / or block diagram can also be realized by computer program instructions. These computer program instructions can be provided to the processor of a general-purpose computer, a special-purpose computer, an embedded processor, or other programmable data processing devices to generate a machine, so that the instructions executed by the processor of the computer or other programmable data processing devices generate a device for realizing the functions specified in Figure 1 one process or multiple processes and / or blocks Figure 1 one block or multiple blocks.

[0314] These computer program instructions can also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to operate in a particular manner, such that the instructions stored in the computer-readable memory produce a manufacture including an instruction device that implements the functions specified in one or more of the processes and / or blocks Figure 1 one or more of the processes and / or blocks Figure 1 specified in the function.

[0315] These computer program instructions can also be loaded onto a computer or other programmable data processing apparatus to cause a series of operational steps to be performed on the computer or other programmable apparatus to produce a computer-implemented process, whereby the instructions executed on the computer or other programmable apparatus provide steps for implementing the functions specified in one or more of the processes and / or blocks Figure 1 one or more of the processes and / or blocks Figure 1 specified in the function.

[0316] The above is only the preferred embodiment of the present invention. It should be pointed out that for those of ordinary skill in the art, without departing from the technical principle of the present invention, several improvements and modifications can be made, and these improvements and modifications should also be regarded as the protection scope of the present invention.

Claims

1. A wind power distribution forecasting method applicable to complex terrain, characterized in that, Including: Dividing the target area into a number of network cells by grid division, and generating ground grid point data according to the terrain distribution characteristics and the spatial distribution characteristics of ground objects in each network cell; Simulating the spatial distribution of low-altitude wind direction, wind speed and atmospheric turbulent kinetic energy in the target area by using a mesoscale meteorological model, and establishing the primary distribution characteristics of the wind field from the spatial distribution of low-altitude wind direction and wind speed; Obtaining the historical wind speed observation data of the target area, and fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the spatio-temporally non-uniform wind profile of the target area, specifically including: When the wind in the network cell flows vertically in, fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the first wind profile function at the wind inlet and outlet positions and the second wind profile function at the middle position; The vertical inflow wind profile is composed of the first wind profile function and the second wind profile function; When the wind in the network cell flows in parallel, fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the third wind profile function at the wind inlet position, the fourth wind profile function at the wind outlet position, and the fifth wind profile function at the middle position; The parallel inflow wind profile is composed of the third wind profile function, the fourth wind profile function and the fifth wind profile function; The spatio-temporally non-uniform wind profile of the target area is composed of the vertical inflow wind profile and the parallel inflow wind profile; Combining the ground grid point data and the spatio-temporally non-uniform wind profile in the network cell to optimize and adjust the primary distribution characteristics of the wind field to form the refined distribution characteristics of the wind field; Calculating the sampling frequency limit parameter according to the historical wind speed observation data of the target area, and calculating the total turbulent kinetic energy of the full spectrum by using the sampling frequency limit parameter from the atmospheric turbulent kinetic energy simulated in the target area; Substituting the refined distribution characteristics of the wind field and the total turbulent kinetic energy of the full spectrum into the simulation process of the mesoscale meteorological model of the target area to obtain the wind power distribution characteristics.

2. The wind power distribution prediction method according to claim 1, characterized in that Fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the first wind profile function at the wind inlet and outlet positions, specifically including: Build a three-dimensional coordinate in the target area and set the axis directions to be divided into axis direction, axis direction, and axis direction; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the first wind profile on the leeward side of the incoming and outgoing wind positions, the axial direction component, and the expression formula is: Axial direction component, the expression formula is: ； ； In the formula, is the height variable, is the height of the first segmentation point, is the height of the second segmentation point, is the height of the third segmentation point, is the leeward side of the first wind profile axial direction component function; , , and are the fitting parameters in the leeward side of the first wind profile axial direction component function; , , are the height constant terms in the leeward side of the first wind profile axial direction component function; and are the wind speed constant terms in the leeward side of the first wind profile axial direction component function; is the exponential function; is the initial wind speed; is the ratio of the height variable to the surface height; is the ratio of the height to the surface height; is the ratio of the height to the surface height; is the ratio of the height to the surface height; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the first wind profile on the leeward side of the incoming and outgoing wind positions, the axial direction component, and the expression formula is: Axial direction component, the expression formula is: ； ； In the formula, , , are the fitting parameters of the leeward side of the first wind profile in the axis direction component; and are the height constant terms of the leeward side of the first wind profile in the axis direction component; is the wind speed constant term of the leeward side of the first wind profile in the axis direction component; is the function of the leeward side of the first wind profile in the axis direction component; The first wind profile on the leeward side of the incoming and outgoing wind positions is obtained by fitting the historical wind speed observation data using a hydrodynamics algorithm The axial direction component, and the expression formula is as follows: ； ； In the formula, is the leeward side of the first wind profile axial direction component function; , , and are the fitting parameters of the axial direction component on the leeward side of the first wind profile; , , and are the height constant terms of the axial direction component on the leeward side of the first wind profile; is the wind speed constant term of the axial direction component on the leeward side of the first wind profile; On the leeward side of the first wind profile Axial direction component, on the leeward side of the first wind profile Axial direction component and on the leeward side of the first wind profile The axial direction components are fitted to form a function of the leeward side of the first wind profile; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the first wind profile on the windward side of the incoming and outgoing wind positions, the axial direction component, and the expression formula is: The axial direction component, and the expression formula is: ； ； In the formula, is the first wind profile upwind side axial direction component function; , , and are the fitting parameters in the first wind profile upwind side axial direction component function; , and are the height constant terms in the first wind profile upwind side axial direction component function; is the first wind profile upwind side axial direction component function of the wind speed constant term; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the first wind profile on the windward side of the incoming and outgoing wind positions, the axial direction component, and the expression formula is: Axial direction component, the expression formula is: ； ； In the formula, is the axial direction component function on the windward side of the first wind profile; The axial direction component function; , , and are the fitting parameters in the axial direction component function on the windward side of the first wind profile; The fitting parameters in the axial direction component function; , and are the height constant terms in the axial direction component function on the windward side of the first wind profile; The height constant terms in the axial direction component function; and are the wind speed constant terms in the axial direction component function on the windward side of the first wind profile; The wind speed constant terms in the axial direction component function; Fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the w-axis direction component on the windward side of the first wind profile at the wind inlet and outlet positions, and the expression formula is: ； ； In the formula, is the axial direction component function on the windward side of the first wind profile; is the axial direction component function; and are the fitting parameters of the axial direction component function on the windward side of the first wind profile; is the fitting parameter of the axial direction component function; and are the height constant terms of the axial direction component function on the windward side of the first wind profile; is the height constant term of the axial direction component function; and are the wind speed constant terms of the axial direction component function on the windward side of the first wind profile; is the wind speed constant term of the axial direction component function; On the windward side of the first wind profile Axial direction component, on the windward side of the first wind profile Axial direction component and on the windward side of the first wind profile The axial direction components are fitted to form a function of the windward side of the first wind profile; The first wind profile function is composed of the first wind profile leeward side function and the first wind profile windward side function.

3. The wind power distribution forecasting method according to claim 2, wherein Fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the second wind profile function at the middle position, specifically including: Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the second wind profile on the leeward side of the middle position The axial direction component, and the expression formula is: ； ； In the formula, is the leeward side of the second wind profile axial direction component function; , , and are the fitting parameters in the leeward side axial direction component function of the second wind profile; , and are the height constant terms in the leeward side axial direction component function of the second wind profile; and are the wind speed constant terms in the leeward side axial direction component function of the second wind profile; The historical wind speed observation data is fitted using a hydrodynamic algorithm to obtain the second wind profile on the leeward side at the intermediate position, and the axial direction component is expressed by the formula: The axial direction component is expressed by the formula: ； ； In the formula, , , are the fitting parameters of the leeward side of the second wind profile in the axis direction component; and are the height constant terms of the axis direction component of the leeward side of the second wind profile; is the wind speed constant term of the axis direction component of the leeward side of the second wind profile; is the axis direction component function of the leeward side of the second wind profile; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the second wind profile on the leeward side of the middle position The axial direction component, and the expression formula is: ； ； In the formula, , , and are the fitting parameters of the leeward side axial direction component of the second wind profile; , , and are the height constant terms of the leeward side axial direction component of the second wind profile; is the wind speed constant term of the leeward side axial direction component of the second wind profile; is the function of the leeward side axial direction component of the second wind profile; On the leeward side of the second wind profile Axial direction component, on the leeward side of the second wind profile Axial direction component and on the leeward side of the second wind profile The axial direction components are fitted to form a function of the leeward side of the second wind profile; The second wind profile on the windward side at the intermediate position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial direction component, and the expression formula is: ； ； In the formula, , , and are the fitting parameters of the axial direction component on the windward side of the second wind profile; , and are the height constant terms of the axial direction component on the windward side of the second wind profile; and are the wind speed constant terms of the axial direction component on the windward side of the second wind profile; is the axial direction component function on the windward side of the second wind profile; ​​​​ The historical wind speed observation data is fitted using a hydrodynamic algorithm to obtain the second wind profile on the windward side at the intermediate position The axial direction component, and the expression formula is as follows: ； ； In the formula, , , and are the fitting parameters of the wind profile's upwind side axis direction component; , and are the height constant terms of the wind profile's upwind side axis direction component; and are the wind speed constant terms of the wind profile's upwind side axis direction component; is the function of the wind profile's upwind side axis direction component; The historical wind speed observation data is fitted using a hydrodynamic algorithm to obtain the second wind profile on the windward side at the intermediate position axial direction component, and the expression formula is as follows: ； ； In the formula, , and are the fitting parameters of the axial direction component on the windward side of the second wind profile; are the height constant terms of the axial direction component on the windward side of the second wind profile; , and are the height constant terms of the axial direction component on the windward side of the second wind profile; are the height constant terms of the axial direction component on the windward side of the second wind profile; and are the wind speed constant terms of the axial direction component on the windward side of the second wind profile; are the wind speed constant terms of the axial direction component on the windward side of the second wind profile; is the function of the axial direction component on the windward side of the second wind profile; is the function of the axial direction component on the windward side of the second wind profile; On the windward side of the second wind profile Axial direction component, on the windward side of the second wind profile Axial direction component and on the windward side of the second wind profile The axial direction components are fitted to form a function of the windward side of the second wind profile; The second wind profile function is composed of the second wind profile leeward side function and the second wind profile windward side function.

4. The wind power distribution prediction method according to claim 2, characterized in that Fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the third wind profile function at the wind inlet position, specifically including: The third wind profile on the leeward side of the incoming wind position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial direction component, and the expression formula is as follows: ； ； In the formula, , and are fitting parameters of the leeward side of the third wind profile in the axis direction component; and are the height constant terms of the leeward side of the third wind profile in the axis direction component; is the wind speed constant term of the leeward side of the third wind profile in the axis direction component; is the function of the leeward side of the third wind profile in the axis direction component; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the third wind profile on the windward side of the incoming wind position The axial direction component, and the expression formula is as follows: ； ； In the formula, , and are the fitting parameters of the axial direction component on the windward side of the third wind profile; The height constant term of the axial direction component on the windward side of the third wind profile; and are the height constant terms of the axial direction component on the windward side of the third wind profile; The height constant term of the axial direction component on the windward side of the third wind profile; is the wind speed constant term of the axial direction component on the windward side of the third wind profile; The wind speed constant term of the axial direction component on the windward side of the third wind profile; is the function of the axial direction component on the windward side of the third wind profile; The function of the axial direction component on the windward side of the third wind profile; Using a hydrodynamics algorithm to fit the historical wind speed observation data to obtain the third wind profile at the air inlet position The axial direction component, and the expression formula is: ； ； In the formula, , and are fitting parameters of the axial direction component of the third wind profile ; and are the height constant terms of the axial direction component of the third wind profile ; is the function of the axial direction component of the third wind profile . The third wind profile of the incoming wind position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial direction component, and the expression formula is: ； ； In the formula, and are the fitting parameters of the axial direction component of the third wind profile; axis direction component; and are the height constant terms of the axial direction component of the third wind profile; axis direction component; is the wind speed constant term of the axial direction component of the third wind profile; axis direction component; is the function of the axial direction component of the third wind profile; axis direction component function; On the leeward side of the third wind profile Axial direction component function, on the windward side of the third wind profile Axial direction component function, the third wind profile Axial direction component function and the third wind profile The axial direction component functions are fitted to form the third wind profile function.

5. The wind power distribution prediction method according to claim 2, characterized in that Fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the fourth wind profile function at the wind outlet position, specifically including: Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the fourth wind profile on the leeward side of the wind outlet position The axial direction component, and the expression formula is: ； ； In the formula, and are the fitting parameters of the leeward side of the fourth wind profile in the axis direction component; and are the height constant terms of the leeward side of the fourth wind profile in the axis direction component; is the wind speed constant term of the leeward side of the fourth wind profile in the axis direction component; is the function of the leeward side of the fourth wind profile in the axis direction component; Using a hydrodynamics algorithm to fit the historical wind speed observation data to obtain the fourth wind profile on the windward side of the air outlet position The axial direction component, and the expression formula is as follows: ； ； In the formula, and are the fitting parameters of the axial direction component on the windward side of the fourth wind profile; is the height constant term of the axial direction component on the windward side of the fourth wind profile; and are the height constant terms of the axial direction component on the windward side of the fourth wind profile; is the height constant term of the axial direction component on the windward side of the fourth wind profile; is the wind speed constant term of the axial direction component on the windward side of the fourth wind profile; is the wind speed constant term of the axial direction component on the windward side of the fourth wind profile; is the function of the axial direction component on the windward side of the fourth wind profile; is the function of the axial direction component on the windward side of the fourth wind profile; The fourth wind profile of the air outlet position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial direction component, and the expression formula is as follows: ； ； In the formula, , and are the fitting parameters of the axial direction component of the fourth wind profile ; and are the height constant terms of the axial direction component of the fourth wind profile ; is the function of the axial direction component of the fourth wind profile . The fourth wind profile of the air outlet position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial direction component, and the expression formula is as follows: ； ； In the formula, and are the fitting parameters of the axial component of the fourth wind profile; is the height constant term of the axial component of the fourth wind profile; and are the height constant terms of the axial component of the fourth wind profile; is the wind speed constant term of the axial component of the fourth wind profile; and are the wind speed constant terms of the axial component of the fourth wind profile; is the function of the axial component of the fourth wind profile; is the function of the axial component of the fourth wind profile; is the function of the axial component of the fourth wind profile; On the leeward side of the fourth wind profile Axial direction component function, on the windward side of the fourth wind profile Axial direction component function, the fourth wind profile Axial direction component function and the fourth wind profile The axial direction component function fitting constitutes the fourth wind profile function.

6. The wind power distribution forecasting method according to claim 2, characterized in that, Fitting the historical wind speed observation data by using a fluid mechanics algorithm to obtain the fifth wind profile function at the middle position, specifically including: The fifth wind profile on the leeward side of the intermediate position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial direction component, and the expression formula is as follows: ； ； In the formula, and are the fitting parameters of the leeward side of the fifth wind profile in the axis direction component; and are the height constant terms of the leeward side of the fifth wind profile in the axis direction component; is the wind speed constant term of the leeward side of the fifth wind profile in the axis direction component; is the function of the leeward side of the fifth wind profile in the axis direction component; The fifth wind profile on the windward side of the intermediate position is obtained by fitting the historical wind speed observation data using a hydrodynamic algorithm The axial component is expressed by the formula: ； ； In the formula, and are the fitting parameters of the axial direction component on the windward side of the fifth wind profile; is the height constant term of the axial direction component on the windward side of the fifth wind profile; and are the height constant terms of the axial direction component on the windward side of the fifth wind profile; is the height constant term of the axial direction component on the windward side of the fifth wind profile; is the wind speed constant term of the axial direction component on the windward side of the fifth wind profile; is the wind speed constant term of the axial direction component on the windward side of the fifth wind profile; is the function of the axial direction component on the windward side of the fifth wind profile; is the function of the axial direction component on the windward side of the fifth wind profile; Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the fifth wind profile at the intermediate position The axial direction component, and the expression formula is: ； ； In the formula, , and are the fitting parameters of the axial direction component of the fifth wind profile ; and are the height constant terms of the axial direction component of the fifth wind profile ; is the function of the axial direction component of the fifth wind profile . Using a hydrodynamic algorithm to fit the historical wind speed observation data to obtain the fifth wind profile at the intermediate position The axial direction component, and the expression formula is: ； ； In the formula, and are the fitting parameters of the axial direction component of the fifth wind profile ; and are the height constant terms of the axial direction component of the fifth wind profile ; and are the wind speed constant terms of the axial direction component of the fifth wind profile ; is the function of the axial direction component of the fifth wind profile . On the leeward side of the fifth wind profile Axial direction component function, on the windward side of the fifth wind profile Axial direction component function, the fifth wind profile Axial direction component function and the fifth wind profile The axial direction component functions are fitted to form the fifth wind profile function.

7. The wind power distribution prediction method according to claim 1, characterized in that Calculating the sampling frequency limit parameter according to the historical wind speed observation data of the target area, and calculating the total turbulent kinetic energy of the full spectrum by using the sampling frequency limit parameter from the atmospheric turbulent kinetic energy simulated in the target area, specifically including: The sampling frequency limit parameter is calculated based on the historical wind speed observation data of the target area, and the expression formula is: ； ； ； In the formula, is the historical turbulent kinetic energy in the observed spectral band; is the historical total turbulent kinetic energy in the full spectral band; is the turbulent energy spectrum in the historical observed data of wind speed; is the set constant term; is the wave number corresponding to the response frequency of the observation instrument; is the wave number corresponding to the turbulent dissipation scale; k is the variable wave number; is the wave number corresponding to the largest-scale turbulent eddy in the atmosphere; is the sampling frequency limit parameter; The total turbulent kinetic energy of the full spectral band is calculated from the simulated atmospheric turbulent kinetic energy of the target area by using the sampling frequency limit parameter, and the expression formula is: ； In the formula, is the atmospheric turbulent kinetic energy simulated for the target area; is the total turbulent kinetic energy of the full spectral band of the target area.

8. A method for predicting the diffusion of air pollution applicable to complex terrain, characterized in that The wind power distribution characteristics are generated by using the wind power distribution prediction method described in any one of claims 1 to 7. The air pollution diffusion prediction method includes: Obtain the background pollution concentration and pollutant emission data of the target area, and calculate the pollution diffusion speed and the spatio-temporal distribution characteristics of the local pollutant concentration in combination with the wind power distribution.

9. A wind power distribution forecasting system applicable to complex terrain, characterized in that, It includes: Surface division unit: The target area is divided into grids to form a number of network cells, and the surface grid point data is generated according to the terrain distribution characteristics and the spatial distribution characteristics of ground objects in each network cell; Simulation unit: In the target area, the spatial distributions of the low-altitude wind direction, wind speed and atmospheric turbulent kinetic energy are simulated by using the mesoscale meteorological model, and the primary distribution characteristics of the wind field are established from the spatial distributions of the low-altitude wind direction and wind speed; Wind field optimization unit: Obtain the historical wind speed observation data of the target area, and use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area; In the network cell, the primary distribution characteristics of the wind field are optimized and adjusted in combination with the surface grid point data and the spatio-temporally non-uniform wind profile to form the refined distribution characteristics of the wind field; Turbulence optimization unit: The sampling frequency limit parameter is calculated based on the historical wind speed observation data of the target area, and the total turbulent kinetic energy of the full spectral band is calculated from the simulated atmospheric turbulent kinetic energy of the target area by using the sampling frequency limit parameter; Output unit: Substitute the refined distribution characteristics of the wind field and the total turbulent kinetic energy of the full spectral band into the mesoscale meteorological model simulation process of the target area to obtain the wind power distribution characteristics; The wind field optimization unit uses the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the spatio-temporally non-uniform wind profile of the target area, which specifically includes: When the wind in the network cell flows in vertically, use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the first wind profile function at the wind inlet and outlet positions and the second wind profile function at the middle position; The vertical inflow wind profile is composed of the first wind profile function and the second wind profile function; When the wind in the network cell flows in parallel, use the fluid mechanics algorithm to fit the historical wind speed observation data to obtain the third wind profile function at the wind inlet position, the fourth wind profile function at the wind outlet position and the fifth wind profile function at the middle position; The parallel inflow wind profile is composed of the third wind profile function, the fourth wind profile function and the fifth wind profile function; The spatio-temporally non-uniform wind profile of the target area is composed of the vertical inflow wind profile and the parallel inflow wind profile.

Images