WIND FIELD DYNAMIC DOWNSCALING METHOD BASED ON AERODYNAMIC PARAMETERS OF SIMPLIFIED TERRAIN

Abstract

A wind field dynamic downscaling method based on aerodynamic parameters of simplified terrain, the method comprises steps of: numerically simulating the simplified terrain based on computational fluid dynamics to obtain the aerodynamic parameters of the simplified terrain; redistributing the wind speed at the corner point of a mesoscale grid within the downscaling grid based on terrain elevation data, land use type data and the aerodynamic parameters, to implement the wind field downscaling calculation. It is based on the aerodynamic parameters of the two-dimensional simplified terrain, and a new wind field dynamic downscaling scheme is designated by adding the high-resolution terrain elevation data and the land use type data.

Claims

1. A wind field dynamic downscaling method based on aerodynamic parameters of simplified terrain, wherein the method comprises steps of: numerically simulating the simplified terrain based on computational fluid dynamics to obtain the aerodynamic parameters of the simplified terrain; redistributing the wind speed at the corner point of a mesoscale grid within the downscaling grid based on terrain elevation data, land use type data and the aerodynamic parameters, to implement the wind field downscaling calculation.

2. The method according to claim 1, wherein the simplified terrain is a two-dimensional simplified terrain including two basic formations, namely, two-dimensional wind speed inflow upslope and two-dimensional wind speed inflow downslope.

3. The method according to claim 1, wherein after numerically simulating the simplified terrain based on computational fluid dynamics to obtain the aerodynamic parameters of the simplified terrain, the method further comprises: performing a wind tunnel test to verify the aerodynamic parameters of the simplified terrain obtained by numerical simulation.

4. The method according to claim 1, wherein before redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, the method further comprises: acquiring the mesoscale wind field of a region based on an area-limited numerical weather prediction mode.

5. The method according to claim 4, wherein said aerodynamic parameters of the simplified terrain comprise an average wind speed ratio at a midpoint of the slope, and said average wind speed ratio satisfies: $R_{\mathrm{ix}}=\frac{U_{\mathrm{ix}}}{U_f},R_{\mathrm{iz}}=\frac{U_{\mathrm{iz}}}{U_f},$ wherein R.sub.ix and R.sub.iz are the average wind speed ratios of the measured point i in the downwind direction and the vertical direction respectively, U.sub.ix and U.sub.iz are the average wind speeds of the measured point i in the downwind direction and the vertical direction respectively, and U.sub.f is the average wind speed of the uniform inflow.

6. The method according to claim 5, wherein redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, comprises: assigning the wind speed at the corner point of the mesoscale grid to the downscaling grid; calculating the wind speed at the midpoint of the mesoscale grid based on the downscaling grid with the assignment; calculating the wind speed at the center point of the mesoscale grid based on the wind speed at the midpoint of the mesoscale grid; completing a first downscaling calculation based on the wind speed at the center point of the mesoscale grid.

7. The method according to claim 6, wherein the wind speed at the center point of the mesoscale grid comprises a zonal wind speed and a meridional wind speed at the center point of the mesoscale grid; the zonal wind speed at the center point of the mesoscale grid satisfies: $u_{m+1,n+1}=\{\begin{array}{c}{u}_{m,n+1}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(u_{m,n+1}>0\right)\\ u_{m+2,n+1}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(u_{m+2,n+1}<0\right)\end{array},$ where u.sub.m+1,n+1 represents the zonal wind speed at the center point of the mesoscale grid, u.sub.m,n+1 and u.sub.m+2,n+1 represent the zonal wind speeds at the midpoints of the two opposite sides in the latitudinal direction of the mesoscale grid respectively, and Ri.sub.m+1,n+1 represents the average wind speed ratio at the center point of the mesoscale grid; the meridional wind speed at the center point of the mesoscale grid satisfies: $v_{m+1,n+1}=\{\begin{array}{c}{v}_{m+1,n}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(v_{m+1,n}>0\right)\\ v_{m+1,n+2}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(v_{m+1,n+2}<0\right)\end{array},$ where v.sub.m+1,n+1 represents the meridional wind speed at the center point of the mesoscale grid, and v.sub.m+1,n and v.sub.m+1,n+2 represent the meridional wind speeds at the midpoints of the two opposite sides in the longitudinal direction of the mesoscale grid respectively.

8. The method according to claim 7, wherein after completing the first downscaling calculation based on the wind speed at the center point of the mesoscale grid, redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, further comprises: repeating the assignment and performing the downscaling calculation based on the downscaling grid after the first downscaling calculation until the resolution of the downscaling grid meets the requirement, and calculating the average wind speed of the downscaled wind field based on the downscaling grid with the resolution meeting the requirements.

9. The method according to claim 8, wherein after the resolution meets the downscaling requirement, redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, further comprises: according to the form of logarithmic law of the vertical wind profile, the change of the average wind speed of the downscaled wind field along the vertical direction being expressed as: $v_z=\frac{u_{*}}{}\left[\ln .Math.\frac{z}{z_0}-{}_m\left(\frac{z}{L_m}\right)\right],$ where v.sub.z is the average wind speed at height z, u* is the friction speed, is the von Krmn constant, z.sub.0 is the surface roughness length, .sub.m is the stability correction function for the logarithmic profile of the average wind speed, and L.sub.m is Obukhov length.

10. The method according to claim 9, wherein the height z ranges from 0 to 150 m.

Description

BRIEF DESCRIPTION OF THE DRAWINGS

[0035] In order to make the objects, technical solutions and beneficial effects of the present invention clearer, the present invention is illustrated by providing the following drawings:

[0036] FIG. 1 is a flow chart of an embodiment of the present invention;

[0037] FIG. 2 is a schematic view of cross-section of a model according to an embodiment of the present invention;

[0038] FIG. 3 is a wind tunnel aerodynamic plane schematic view according to an embodiment of the present invention;

[0039] FIG. 4 is a schematic view of a latitudinal-direction downscaling calculation process according to an embodiment of the present invention;

[0040] FIG. 5 is a schematic view of a longitudinal-direction downscaling calculation process according to an embodiment of the present invention.

DETAILED DESCRIPTION

[0041] Hereinafter, preferred embodiments of the present invention will be described in detail in combination with the accompanying drawings.

[0042] The present invention provides a wind field dynamic downscaling method based on aerodynamic parameters of simplified terrain. As shown in FIG. 1, the method comprises the following steps:

[0043] Numerically simulating the simplified terrain based on computational fluid dynamics to obtain the aerodynamic parameters of the simplified terrain;

[0044] Redistributing the wind speed at the corner point of a mesoscale grid within the downscaling grid based on terrain elevation data, land use type data and the aerodynamic parameters, to implement the wind field downscaling calculation.

[0045] In this embodiment, numerically simulating the simplified terrain based on computational fluid dynamics to obtain the aerodynamic parameters of the simplified terrain comprises that the CFD numerical simulation is completed by using the commercial software FLUENT. The coupling system of pressure and speed in CFD calculation is SIMPLEC algorithm, the governing equation is solved by a segregated method, and the pressure difference mode is PRESTO. The turbulence model uses the Realizable k-c model, the convection term of the governing equation adopts the second-order upwind format, and the residual value is taken as 510.sup.4 according to the computational convergence criterion.

[0046] The numerical wind field of CFD adopts the uniform inflow boundary condition, the inflow adopts the speed inlet boundary condition, and the outflow adopts the pressure outlet boundary condition. The top surface uses a symmetrical boundary condition, which is equivalent to a free-slip wall surface. The bottom surface of the calculation domain and the surface of the research object adopt non-slip wall surfaces.

[0047] CFD numerical simulations were carried out for both upslope and downslope conditions, in which 12 kinds of slope angles such as 5, 10, 15, 20, 25, 30, 35, 40, 45, 50, 55, 60 are considered.

[0048] In this embodiment, as shown in FIG. 2, the simplified terrain in the embodiment specifically refers to a two-dimensional simplified terrain, in which the slope a is used as a characteristic parameter and the two basic formations of two-dimensional wind speed inflow upslope and downslope are considered.

[0049] In this embodiment, before redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, the method further comprises: acquiring the mesoscale wind field of a region based on an area-limited numerical weather prediction mode.

[0050] As to the specific Terrain Elevation Data, this embodiment adopts the data obtained in the Shuttle Radar Topography Mission (SRTM) conducted in 2000 by the National Aeronautics and Space Administration (NASA), which is referred to as SRTM data.

[0051] The land use type data adopts the global 30-meter ground-surface coverage data (GlobeLand30) released by the National Geomatic Center of China in 2010. The images classified in the GlobeLand30 are multi-spectral images with a horizontal resolution of 30 m, including United States Land Resources Satellite (Landsat) TM5, ETM+ multispectral images and China Environmental Disaster Reduction Satellite (HJ-1) multispectral images. The GlobeLand30 data consists of 10 types of cultivated land, forests, grasslands, shrubs, water areas, wetlands, tundra, man-made cover, bare land, glaciers and permanent snow.

[0052] Acquiring the mesoscale wind field of a region based on an area-limited numerical weather prediction mode comprises, firstly, acquiring global numerical forecast data by a global numerical weather prediction mode; then acquiring the numerical forecast data of the region by the area-limited numerical weather prediction mode, including a near-surface wind field, i.e., the mesoscale wind field having horizontal resolution of kilometer level.

[0053] In this embodiment, after numerically simulating the simplified terrain based on computational fluid dynamics to obtain the aerodynamic parameters of the simplified terrain, the method further comprises: performing a wind tunnel test to verify the aerodynamic parameters of the simplified terrain obtained by numerical simulation.

[0054] As shown in FIG. 3, it is a wind tunnel aerodynamic plane schematic view according to an embodiment of the present invention. The wind tunnel used in this embodiment is a series-wound double-test-section back-flow/direct-flow large-scale multi-functional boundary layer wind tunnel, wherein its low-speed test section is 4.4 m wide, 3.0 m high and 24.0 m long, and the maximum wind speed is greater than 30.0 m/s; its high-speed test section is 2.2 m wide, 2 m high and 5.0 m long, and the maximum wind speed is greater than 80.0 m/s, as shown in FIG. 3.

[0055] In this embodiment, the wind tunnel test is a pressure-measuring test for rigid models. In this embodiment, four types of models M1, M2, M3, and M4 are experimentally fabricated. As shown in FIG. 2, the cross-sectional dimensions of four models M1, M2, M3, and M4 are shown in Table 1:

TABLE-US-00001 TABLE 1 Parameters of Models M1~M4 Slope Characteristic Model Blocking Model No. angle Dimension L/mm Height H/mm Rate/% M1 15 500 133.97 4.47 M2 30 400 230.94 7.70 M3 45 200 200 6.67 M4 60 150 259.81 8.66
The specific test parameters of the wind tunnel test in this embodiment are shown in Table 2.

TABLE-US-00002 TABLE 2 Technical Parameters of Wind Tunnel Test Classification of Parameters Name of Parameters Indexes Meteorology, Inflow Conditions Uniform Flow Field Geomorphology and Building Environment Model Parameters Dimensions of Test 4.4 m width 3 m Section height 24 m length Form of Model Rigid Model Number of Models 5 Test Parameters Wind Speed of Free 7 m/s Inflow Pressure-Measuring 330 Hz Sampling frequency Pressure-Measuring 60 s Sampling Time Wind Speed Sampling 625 Hz Frequency Wind Speed Sampling 60 s Time

[0056] Optionally, in this embodiment, the aerodynamic parameter of the two-dimensional simplified terrain acquired by the CFD numerical simulation and wind tunnel test is the average wind speed ratio of the midpoint of the slope, and the average wind speed ratio satisfies:

[00005]$R_{\mathrm{ix}}=\frac{U_{\mathrm{ix}}}{U_f},R_{\mathrm{iz}}=\frac{U_{\mathrm{iz}}}{U_f},$

[0057] Wherein R.sub.ix and R.sub.iz are the average wind speed ratios of the measured point i in the downwind direction and the vertical direction respectively, U.sub.ix and U.sub.iz are the average wind speeds of the measured point i in the downwind direction and the vertical direction respectively, and U.sub.f is the average wind speed of the uniform inflow.

[0058] Optionally, in this embodiment, redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, comprises:

[0059] Assigning the wind speed at the corner point of the mesoscale grid to the downscaling grid;

[0060] Calculating the wind speed at the midpoint of the mesoscale grid based on the downscaling grid with the assignment;

[0061] Calculating the wind speed at the center point of the mesoscale grid based on the wind speed at the midpoint of the mesoscale grid;

[0062] Completing a first downscaling calculation based on the wind speed at the center point of the mesoscale grid.

[0063] Specifically, the downscaling calculation includes calculations in two directions, namely zonal wind speed and meridional wind speed, corresponding to the grid division diagrams of FIGS. 4 and 5 respectively, which are further described below based on the contents of FIGS. 4 and 5.

[0064] As shown in FIG. 4, for zonal wind speed, the calculation steps are described as follows:

[0065] 1) firstly, assigning the wind speed U at the four corner points of the mesoscale grid indicated by the thick solid lines in FIG. 4 to the downscaling grid u indicated by the thin solid lines in FIG. 4:

u.sub.m,n=U.sub.i,j

u.sub.m+2,n=U.sub.i+1,j

u.sub.m,n+2=U.sub.i,j+1

u.sub.m+2,n+2=U.sub.i+1,j+1.

[0066] 2) finding the zonal wind speed u at the midpoints of the four sides of the mesoscale grid indicated by the thick solid lines in FIG. 4:

[00006]$.Math.u_{m,n+1}=\frac{u_{m,n}+u_{m,n+2}}{2}=\frac{U_{i,j}+U_{i,j+1}}{2}$$.Math.u_{m+2,n+1}=\frac{u_{m+2,n}+u_{m+2,n+2}}{2}=\frac{U_{i+1,j}+U_{i+1,j+1}}{2}$$.Math.u_{m+1,n}=\{\begin{array}{c}{u}_{m,n}.Math.{\mathrm{Ri}}_{m+1,n}=U_{i,j}.Math.{\mathrm{Ri}}_{m+1,n}\left(U_{i,j}>0\right)\\ u_{m+2,n}.Math.{\mathrm{Ri}}_{m+1,n}=U_{i+1,j}.Math.{\mathrm{Ri}}_{m+1,n}\left(U_{i+1,j}<0\right)\end{array}.Math..Math.u_{m+1,n+2}=\{\begin{array}{c}{u}_{m,n+2}.Math.{\mathrm{Ri}}_{m+1,n+2}=U_{i,j+1}.Math.{\mathrm{Ri}}_{m+1,n+2}\left(U_{i,j+1}>0\right)\\ u_{m+2,n+2}.Math.{\mathrm{Ri}}_{m+1,n+2}=U_{i+1,j+1}.Math.{\mathrm{Ri}}_{m+1,n+2}\left(U_{i+1,j+1}<0\right)\end{array}.$

[0067] 3) finding the zonal wind speed u at the center point of the mesoscale grid:

[00007]$u_{m+1,n+1}=\{\begin{array}{c}{u}_{m,n+1}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(u_{m,n+1}>0\right)\\ u_{m+2,n+1}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(u_{m+2,n+1}<0\right)\end{array}.$

[0068] As shown in FIG. 5, for the meridional wind speed, the calculation steps are described as follows:

[0069] 1) firstly, assigning the wind speed Vat the four corner points of the mesoscale grid indicated by the thick solid lines in FIG. 5 to the downscaling grid v indicated by the thin solid lines in FIG. 5:

v.sub.m,n=V.sub.i,j

v.sub.m+2,n=V.sub.i+1,j

v.sub.m,n+2=V.sub.i,j+1

v.sub.m+2,n+2=V.sub.i+1,j+1.

[0070] 2) finding the meridional wind speed at the midpoints of the four sides of the mesoscale grid indicated by thick solid lines:

[00008]$v_{m+1,n}=\frac{v_{m,n}+v_{m+2,n}}{2}=\frac{V_{i,j}+V_{i+1,j}}{2}$$v_{m+1,n+2}=\frac{v_{m,n+2}+v_{m+2,n+2}}{2}=\frac{V_{i,j+1}+V_{i+1,j+1}}{2}$$v_{m,n+1}=\{\begin{array}{c}{v}_{m,n}.Math.{\mathrm{Ri}}_{m,n+1}=V_{i,j}.Math.{\mathrm{Ri}}_{m,n+1}\left(V_{i,j}>0\right)\\ v_{m,n+2}.Math.{\mathrm{Ri}}_{m,n+1}=V_{i,j+1}.Math.{\mathrm{Ri}}_{m,n+1}\left(V_{i,j+1}<0\right)\end{array}.Math..Math.v_{m+2,n+1}=\{\begin{array}{c}{v}_{m+2,n}.Math.{\mathrm{Ri}}_{m+2,n+1}=V_{i+1,j}.Math.{\mathrm{Ri}}_{m+2,n+1}\left(V_{i+1,j}>0\right)\\ v_{m+2,n+2}.Math.{\mathrm{Ri}}_{m,n+1}=V_{i+1,j+1}.Math.{\mathrm{Ri}}_{m+2,n+1}\left(V_{i+1,j+1}<0\right)\end{array}.$

[0071] 3) finding the meridional wind speed v at the center point of the mesoscale grid:

[00009]$v_{m+1,n+1}=\{\begin{array}{c}{v}_{m+1,n}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(v_{m+1,n}>0\right)\\ v_{m+1,n+2}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(v_{m+1,n+2}<0\right)\end{array}.$

[0072] That is, the wind speed at the center point of the mesoscale grid includes the zonal wind speed at the center point of the mesoscale grid and the meridional wind speed at the center point of the mesoscale grid.

[0073] The zonal wind speed at the center point of the mesoscale grid satisfies:

[00010]$u_{m+1,n+1}=\{\begin{array}{c}{u}_{m,n+1}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(u_{m,n+1}>0\right)\\ u_{m+2,n+1}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(u_{m+2,n+1}<0\right)\end{array},$

[0074] Where u.sub.m+1,n+1 represents the zonal wind speed at the center point of the mesoscale grid, and u.sub.m,n+1 and u.sub.m+2,n+1 represent the zonal wind speeds at the midpoints of the two opposite sides in the latitudinal direction of the mesoscale grid respectively, and Ri.sub.m+1,n+1 represents the average wind speed ratio at the center point of the mesoscale grid;

[0075] The meridional wind speed at the center point of the mesoscale grid satisfies:

[00011]$v_{m+1,n+1}=\{\begin{array}{c}{v}_{m+1,n}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(v_{m+1,n}>0\right)\\ v_{m+1,n+2}.Math.{\mathrm{Ri}}_{m+1,n+1}\left(v_{m+1,n+2}<0\right)\end{array},$

[0076] Where v.sub.m+1,n+1 represents the meridional wind speed at the center point of the mesoscale grid, and v.sub.m+1,n and v.sub.m+1,n+2 represent the meridional wind speeds at the midpoints of the two opposite sides in the longitudinal direction of the mesoscale grid respectively.

[0077] Optionally, after completing the first downscaling calculation based on the wind speed at the center point of the mesoscale grid, redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, further comprises:

[0078] Repeating the assignment and performing the downscaling calculation based on the downscaling grid after the first downscaling calculation until the resolution of the downscaling grid meets the requirement, and calculating the average wind speed of the downscaled wind field based on the downscaling grid with the resolution meeting the requirements.

[0079] The first downscaling calculation was completed. To further downscale, the downscaling calculation steps 1) 3) for the zonal wind speed and the meridional wind speed are repeated within the four downscaling grids indicated by thin solid lines in FIGS. 4 and 5, and the second downscaling calculation is performed to find the zonal wind speed u and the meridional wind speed v at various points inside the downscaling grids. It is repeated like this until the horizontal resolution meets the downscaling requirement.

[0080] According to the foregoing equations, in the present embodiment, the average wind speed calculation of the downscaled wind field with a height of 10 m is taken as an example. U and V are the average zonal wind speed and the average meridional wind speed at a height of 10 m from the ground surface in the mesoscale mode, which may be derived from the output results of the area-limited numerical weather mode; u, v are the average zonal wind speed and the average meridional wind speed at a height of 10 m from the ground surface in the downscaled wind field; Ri is the average wind speed ratio obtained by the CFD numerical simulation. Assuming that the horizontal resolution of the mesoscale grid is x of kilometer level, the horizontal resolution is x/2 after the first downscaling is completed; the horizontal resolution is x/4 after the second downscaling is completed; its horizontal resolution is x/8 after the third downscaling is completed, and so on.

[0081] Based on the foregoing scheme, the average wind speed at a height of 10 m from the ground surface in the downscaled wind field can be obtained:

[00012]$v.Math..Math.{10}_{m+p,n+q}=\sqrt{u_{m+p,n+q}^2+v_{m+p,n+q}^2};$$p=0,1,\mathrm{.Math.}.Math.$$q=0,1,\mathrm{.Math.}$$.Math..Math.{10}_{m+p,n+q}=\arctan \left(\frac{v_{m+p,n+q}}{u_{m+p,n+q}}\right);$$p=0,1,\mathrm{.Math.}.Math.$$q=0,1,\mathrm{.Math.}$

[0082] Where v10 and 10 are the average total wind speed and average wind direction at a height of 10 m from the ground surface in the downscaled wind field, respectively.

[0083] Optionally, after the grid resolution meets the downscaling requirement, redistributing the wind speed at the corner point of the mesoscale grid within the downscaling grid based on the terrain elevation data, the land use type data and the aerodynamic parameters to implement the wind field downscaling calculation, further comprises:

[0084] According to the form of logarithmic law of the vertical wind profile, the change of the average wind speed of the downscaled wind field along the vertical direction being expressed as:

[00013]$v_z=\frac{u_{*}}{}\left[\ln .Math.\frac{z}{z_0}-{}_m\left(\frac{z}{L_m}\right)\right]$

[0085] Where v.sub.z is the average wind speed at height z, the height z ranges from 0 to 150 m, u* is the friction speed, is the von Kmn constant, =0.40, z.sub.0 is the surface roughness length (obtained by land use type data), .sub.m is the stability correction function for the logarithmic profile of the average wind speed, and L.sub.m is Obukhov length.

[0086] According to the average wind speed and average wind direction at a height of 10 m from the ground surface obtained based on the above calculation, the average wind speeds at other heights are expressed as:

[00014]$v_z=v.Math..Math.10.Math.\frac{\ln .Math.\frac{z}{z_0}-{}_m\left(\frac{z}{L_m}\right)}{\ln .Math.\frac{10}{z_0}-{}_m\left(\frac{10}{L_m}\right)}.$

[0087] The method of the present invention is further described below in conjunction with specific embodiments.

[0088] In the current area-limited model named Shanghai Meteorological Bureau-WRF ADAS Rapid Refresh System (SMB-WARR), the horizontal grid spacing is 3 km and the number of grid points is 852792.

[0089] After one downscaling calculation is performed using the downscaling method of the present invention, the horizontal grid spacing is 1.5 km, the number of grid points is 17031583 (about 2.7 million), and the calculation time for completing the 24-hour forecast is about 1.5 minutes;

[0090] After two downscaling calculations are performed using the downscaling method of the present invention, the horizontal grid spacing is 750 m, the number of grid points is 34053165 (about 10.78 million), and the calculation time for completing the 24-hour forecast is about 4 minutes;

[0091] After three downscaling calculations are performed, the horizontal grid spacing obtained is 375 m, the number of grid points is 68096329 (about 43.09 million), and the calculation time for completing the 24-hour forecast is about 12 minutes, as shown in Table 3.

TABLE-US-00003 TABLE 3 Statistics Table of Calculation Efficiency of the Downscaling Method in the Present Invention Number of Calculation Times of STIDM Horizontal Grid Time Downscaling Spacing (m) Number of Grids (minutes) 1 1500 1703 1583 1.5 (About 2.7 million) 2 750 3405 3165 4 (About 10.78 million) 3 375 6809 6329 12 (About 43.09 million)

[0092] Based on the above technical solutions, the method of the present invention has the following advantages:

[0093] 1. Downscaled area is large. Since the downscaling method of the present invention is to perform the downscaling calculation based on the mesoscale grid, its calculable area is very large and its horizontal range can reach thousands of kilometers, satisfying the need for large-area wind field forecast of the systems such as power grids, railways, and highways; the downscaled area in the method of the invention is larger than other existing downscaling techniques, e.g., CFD software (FLUENT or OpenFOAM) having the calculation horizontal range generally within several kilometers, and software such as CALMET and WT having the horizontal range generally within hundreds of kilometers.

[0094] 2. Computational efficiency is high. The aerodynamic parameters of the simplified terrain, i.e., the average wind speed ratios, are made off-line, and a dynamic downscaling method based on a mesoscale grid is used, so the calculation efficiency is very high; the calculation of this method only takes a few minutes in contrast with the traditional dynamic downscaling schemes, for example, CFD software (FLUENT or OpenFOAM) having the calculation time of usually several hours or even days, and for example, the CALMET and WT tools having the calculation time of several hours.

[0095] 3. The resolution of wind field is high. The method of the present invention will realize the refined wind field forecast under complex terrain conditions, and will provide important technical guarantee and business supports for the construction of the typhoon wind disaster early warning system within large areas such as transmission lines, railways and highways, and it is conducive to the development of disaster prevention and mitigation in coastal areas during typhoon landing; the horizontal resolution of the downscaled wind field obtained by this method is up to 100 m level, which will make it possible to accurately calculate the distribution of wind energy in any geographical location; in future life, if humans increase the utilization of wind power resources, the method of the present invention will provide technical support for efficient and accurate utilization of wind power resources, and contribute to clean energy and sustainable development of human beings.

[0096] It shall be noted that, in the text, the term comprises, or includes, or any other variants thereof, covers non-exclusive inclusion, such that a process, a method, an article, or a device comprising a series of elements includes not only those elements but also other elements that are not explicitly listed, or inherent elements of such a process, method, article, or device. Without further limitation, an element defined by the phrase comprising a . . . does not exclude the existence of additional identical elements in the process, method, article or device that comprises the element.

[0097] The serial numbers of the embodiments of the present invention described above are merely for depiction, and do not mean that the embodiment is good or bad.

[0098] Through the description of the above embodiments, those skilled in the art can clearly realize that the method in the above embodiments can be implemented by means of software plus a necessary general hardware platform, and of course, it can also be implemented by hardware, but in many cases, the former is a better embodiment. Based on such understanding, the technical solution of the present invention essentially, or the part that contributes to the prior art, may be embodied in the form of a software product, and the computer software product is stored in a storage medium such as ROM/RAM, disk, optical disc, including a number of instructions for causing a terminal, which may be a cell phone, a computer, a server, an air conditioner, or a network device, etc., to perform the methods described in various embodiments of the present invention.