Forecasting and Optimization of Wind Speed over the Gobi Grassland Wind Farm in Western Inner Mongolia

: Wind power, as one of the primary clean energies, is an important way to achieve the goals of carbon peak and carbon neutrality. Therefore, high-resolution measurement and accurate forecasting of wind speed are very important in the organization and dispatching of the wind farm. In this study, several methodologies, including the mesoscale WRF (Weather Research and Forecasting(WRF) model, mathematical statistics algorithms, and machine learning algorithms, were adopted to systematically explore the predictability and optimization of wind speed of a Gobi grassland wind farm located in western Inner Mongolia. Results show that the rear-row turbines were signiﬁcantly affected by upwind turbine wakes. The output power of upwind-group turbines was 591 KW with an average wind speed of 7.66 m/s, followed by 532 KW and 7.02 m/s in the middle group and 519 KW and 6.92 m/s in the downwind group. The higher the wind speed was, the more signiﬁcantly the wake effect was presented. Intercomparison between observations and WRF simulations showed an average deviation of 3.73 m/s. Two postprocessing methods of bilinear interpolation and nearest replacement could effectively reduce the errors by 34.85% and 36.19%, respectively, with average deviations of 2.43 m/s and 2.38 m/s. A cycle correction algorithm named Average Variance–Trend (AVT) can further optimize the errors to 2.14 m/s and 2.13 m/s. In another aspect, the categorical boosting (CatBoost) artiﬁcial intelligence algorithm also showed a great performance in improving the accuracy of WRF outputs, and the four-day average deviation of 26–29 September decreased from 3.21 m/s to around 2.50 m/s. However, because of the inﬂuence of large-scale circulations, there still exist large errors in the results of various correction algorithms. It is therefore suggested through the investigation that data assimilation of the northwest and Mongolian plateau, boundary layer parameterization scheme optimization, and embedding of high-resolution topographic data could have great potential for obtaining more accurate forecasting products.


Introduction
Wind energy, as one of the most promising clean energy sources, has been maintaining a strong momentum of development all around the world [1][2][3][4][5]. China is a country with abundant wind resources. According to the national wind energy resource assessment data in 2020, the exploitable wind energy at 70 m in mainland China was about 5 billion kilowatts, which provides a good prospect for the development and utilization of wind power in China [6][7][8]. Especially under the major national strategic goals of carbon peak and carbon neutrality, China's power system is changing from traditional fossil energy to a new-energy-based structure in order to achieve the goal of low-carbon economic development [9][10][11].
Wind speed forecasting is considered one of the most effective and economical means to increase the peak load regulation capacity of the power grid, the capacity of the power grid to accommodate wind power, and the safety of power system operations [12][13][14]. Through more accurate prediction of the regional wind speed, the dispatching department can reasonably arrange the power generation plan and reduce the reserve capacity of the power system, thereby strengthening the ability of risk resistance of the power grid and promoting the development of renewable energy [15][16][17][18][19][20][21].
The Urat Gobi grassland, located in western Inner Mongolia, has a characteristic continental arid climate, far away from the sea and lying deep in the mainland. This region has four distinctive seasons with a short summer, a dry and windy spring, a mild and cool autumn, and a long-lasting winter from November to March. There are less frequent rain and snow in winter and spring, whereas concentrated rainfall usually occurs in summer with a large temperature difference between daytime and nighttime. In addition, the region is rich in wind energy resources, second only to the southeast coast of China. The number of hours with a wind speed larger than 8 m/s in a year is 2382 h, accounting for 27.2% of the whole year. The superior wind energy reserves provide favorable conditions for the development of wind power in the Inner Mongolia Autonomous Region.
This study takes a Gobi grassland wind farm located in western inner Mongolia as the research object. By using one-month observations (September 2021) of the wind speed made by the sensors mounted on turbines, in combination with a mesoscale meteorological model, multiple mathematical statistics algorithms, and the CatBoost machine learning algorithm, we systemically explored the predictability and optimization of wind speed of the Gobi grassland wind farm.

Observational Site
The Gobi grassland wind farm (41.55 • N, 106.99 • E) was built north of the Bayannur city of Inner Mongolia in 2009, which is located in the Urat Gobi grassland at the northern foot of Yinshan Mountains, as shown in Figure 1a. This wind Farm has 68 turbines with a rated power of 1500 KW. A meteorological wind sensor was mounted behind the turbine hub to measure the wind at hub height. Among these 68 turbines, six numbered 08, 09, 24, 33, 37, and 45 were out of service, marked as "×" in Figure 1d. There were also some failures in Turbine 60, but it could still generate power.

Weather Research and Forecasting (WRF) Model
The Weather Research and Forecasting (WRF) model is a mesoscale numerical weather prediction framework designed for operational weather forecasting and atmospheric scientific research [22][23][24][25]. Here we adopted version 4.2, which is a fully compressible and nonhydrostatic mode with an Arakawa C-grid in the horizontal direction and high-terrain following coordinates in a vertical direction by default. In general, the WRF model can not only be applied to the case simulation of realistic weather but is also widely processed in the forecast of wind speed and power.
The simulation region was centered on (41.5 • N, 107.0 • E), and the three nested grid domains were designed with a spatial resolution of 27 km, 9 km, and 3 km, respectively, as shown in Figure 2a. The vertically layered encryption was 55 layers, and the top layer was 50 hPa. The parameterization schemes selected in this simulation are summarized in Table 1, in which the cumulus cloud parameterization scheme is disabled for Domain 03. For a detailed scheme introduction, please refer to   [23].

Weather Research and Forecasting (WRF) Model
The Weather Research and Forecasting (WRF) model is a mesoscale nume weather prediction framework designed for operational weather forecasting and at pheric scientific research [22][23][24][25]. Here we adopted version 4.2, which is a fully comp ible and nonhydrostatic mode with an Arakawa C-grid in the horizontal direction high-terrain following coordinates in a vertical direction by default. In general, the model can not only be applied to the case simulation of realistic weather but is also w processed in the forecast of wind speed and power.
The simulation region was centered on (41.5° N, 107.0° E), and the three nested domains were designed with a spatial resolution of 27 km, 9 km, and 3 km, respecti as shown in Figure 2a. The vertically layered encryption was 55 layers, and the top was 50 hPa. The parameterization schemes selected in this simulation are summariz Table 1, in which the cumulus cloud parameterization scheme is disabled for Domai For a detailed scheme introduction, please refer to   [23].  A previous study has pointed out that the meteorological variables such as 10 m wind, sea-land terrain, surface net radiation, surface pressure, 10 m zonal wind component, 10 m meridional wind component, sensible heat flux, total column water vapor, 2 m dew point temperature, mean sea level pressure, and total cloud cover all have essentia impacts on wind speed forecasting [26]. In addition to the influence of large-scale circulation, the boundary layer wind field is also affected by topographic relief, vegetation covers  A previous study has pointed out that the meteorological variables such as 10 m wind, sea-land terrain, surface net radiation, surface pressure, 10 m zonal wind component, 10 m meridional wind component, sensible heat flux, total column water vapor, 2 m dew point temperature, mean sea level pressure, and total cloud cover all have essential impacts on wind speed forecasting [26]. In addition to the influence of large-scale circulation, the boundary layer wind field is also affected by topographic relief, vegetation covers, as well as land use types. Simple smoothing operation of coarse-resolution model over terrains will produce potential errors in the boundary layer wind field, especially near the ground surface [27][28][29][30]. What is more, directly increasing the resolution of the model may lead to the complicated "gray zone" of boundary layer parameterization schemes, thereby generating more problems. In order to obtain a better result, the uses of mathematical statistics methods and artificial intelligence methods to learn the impact of terrain and other factors have been widely proposed for the postprocessing correction of wind speed.

Mathematical Statistics Algorithms Bilinear Interpolation and Nearest Replacement
The output of WRF model simulations was stored in grid mode according to netCDF format. Since the grid position cannot precisely correspond to the turbine position, it is necessary to use the spatial interpolation method or replacement method to extrapolate and obtain accurate turbine position data. In this paper, bilinear interpolation and nearest replacement methodologies were adopted to process the data.
Bilinear interpolation is to estimate f (x,y) with the four nearest grid points ( Figure 2b). By performing linear interpolation in the X direction first and then linear interpolation in the Y direction, we can obtain the final result of f (x,y). This interpolation output contains effective information on the four nearest grid points and will, therefore, be more representative and meaningful. More importantly, the computational efficiency is also reduced. The specific formula is as follows: The nearest replacement directly assigns the value of f (x,y) equal to the nearest grid point value, which has the advantage of higher computational efficiency compared with bilinear interpolation.

Average Variance-Trend (AVT) Method
The Average Variance-Trend (AVT) method was used to correct the wind speed calculated by bilinear interpolation and nearest replacement. This methodology, first proposed by Zhang et al. (2017) [31], de-trends both observational and simulated data series in the first step. Then, the two primary AVT parameters of average and variance of the two de-trending series are calculated. Finally, by correcting the average and variance of the de-trending simulated wind speed with an empirical equation [32], the corrected wind speed series can be obtained, which has the same average and variance as the observations. For a detailed introduction and specific equations of the AVT method, please refer to Zhang et al. (2018; [32,33]. Zhang et al. (2017) [31] applied the AVT method to three wind farms (Heiyazi wind farm, Machangshan wind farm, and Nanqiu wind farm) located in Gansu province. They found that the corrections based on AVT have satisfactory performance and result in an overall reduction in the forecast error rate of 30%. Recently, Han et al. (2022) [34] used this method to correct the prediction results of solar radiation forecast and reduce the bias effectively.

Machine Learning Algorithm
CatBoost, a categorical boosting algorithm, is a Gradient Boosting Decision Tree (GBDT) framework proposed on the basis of the Symmetrical Decision Tree algorithm, with characteristics of fewer parameters, class-based variables, and high accuracy. CatBoost solves the problems of gradient bias and prediction shift, thereby reducing the occurrence of overfitting and improving the accuracy and generalization ability of the algorithm [35,36]. Its advantages are: (1) It automatically treats categorical features in a special way. First, we make some statistics on categorical features and calculate the frequency of a certain category. Then we add the hyperparameter to generate new numerical features. With CatBoost, one no longer has to manually process categorical features; (2) CatBoost also uses composite category features to take advantage of the connections between features, which greatly enriches feature dimensions; (3) The method of Ordered Boost is used to avoid the deviation of gradient estimation and solve the problem of prediction offset; (4) The base model of CatBoost uses symmetric trees, and the method of calculating leaf value is also different from the traditional boosting algorithm. CatBoost has optimized calculating the average value and adopted other algorithms, all of which can prevent model overfitting.

Wind Speed Distribution and Their Wake Effect
Analysis of the wind data observed on the turbines in September 2021 in the Gobi grassland wind farm shows that southerly wind dominates this area throughout the month, primarily distributed around 175-185 • . According to this characteristic, we grouped these turbines into three categories along north-south wind direction, i.e., upwind group (Turbines 10, 16, 31, 32, 42, 44, and 46), middle group (Turbines 06, 22, 23, 27, 28, 57, and 64), and downwind group (Turbines 02, 03, 50, 51, 52, 53, and 60), shown in Figure 3a. By averaging the generated powers and wind speeds of the three groups, it is observed that the average power of the upwind group is 591 KW, with an average wind speed of 7.66 m/s, followed by 532 KW of the middle group (an average wind speed of 7.02 m/s), and 519 KW of downwind group (an average wind speed of 6.92 m/s). The results indicated that the wake effects of the front turbines would significantly affect the generating capacity of rear-row turbines, with an average decrease of 12%. 7.02 m/s), and 519KW of downwind group (an average wind speed of 6.92 m/s). The sults indicated that the wake effects of the front turbines would significantly affect generating capacity of rear-row turbines, with an average decrease of 12%. The power and wind speed in the middle and downwind groups were 365 KW and 5 m/s, 118 KW and 3.65 m/s, respectively. By comparison, it is found that the wind wa effect significantly reduced the wind speed of the rear-row turbines, but the power ge eration of downwind turbines was higher than that of the middle group and upwi group, and the reasons for this phenomenon were not clear. In Figure 4b, the proporti of the conditions with the southerly wind of 5-10 m/s was 53.7%, and the downwi power and average wind speed were 523 KW and 6.97 m/s, respectively. The power a average wind speed of the middle and upwind groups were 537 KW and 7.03 m/s, 5 KW and 7.40 m/s, respectively. Results for the two scenarios indicated that in low wi speed conditions, the wake effect was not significant, and the powers generated by bo front and rear groups were similar. In Figure 4c, the proportion of south wind 10-15 m was 16.7%, downwind power and wind speed were 688 KW and 8.17 m/s, smaller th 711 KW and 8.28 m/s in the middle group and 1283 KW and 11.83 m/s in upwind grou It is observed that the increase in wind speed made the upwind wake effect more p nounced and the power generation in the rear row decreased by ~ 46%. In Figure 4d, proportion of south wind >15 m/s was 2.0%. Upwind power and wind speed were 14 KW and 16.16 m/s, respectively, higher than those of 1234 KW and 11.83 m/s in the mid group and 1234 KW and 11.77 m/s in the downwind group. At higher wind conditio the wake effect perturbed by upwind turbines almost affected the whole wind farm, a the rear turbines could not reach full capacities.  According to different wind speed regimes, the turbines are divided into four types: 0-5 m /s, 5-10 m/s, 10-15 m/s, and >15 m/s, as shown in Figure 4 and Table 2. In Figure 4a, the occurrences of southerly 0-5 m/s wind accounted for 27.6%, with an average of generating power and wind speed in upwind direction of 368 KW and 5.90 m/s, respectively. The power and wind speed in the middle and downwind groups were 365 KW and 5.87 m/s, 118 KW and 3.65 m/s, respectively. By comparison, it is found that the wind wake effect significantly reduced the wind speed of the rear-row turbines, but the power generation of downwind turbines was higher than that of the middle group and upwind group, and the reasons for this phenomenon were not clear. In Figure 4b, the proportion of the conditions with the southerly wind of 5-10 m/s was 53.7%, and the downwind power and average wind speed were 523 KW and 6.97 m/s, respectively. The power and average wind speed of the middle and upwind groups were 537 KW and 7.03 m/s, 524 KW and 7.40 m/s, respectively. Results for the two scenarios indicated that in low wind speed conditions, the wake effect was not significant, and the powers generated by both front and rear groups were similar. In Figure 4c, the proportion of south wind 10-15 m/s was 16.7%, downwind power and wind speed were 688 KW and 8.17 m/s, smaller than 711 KW and 8.28 m/s in the middle group and 1283 KW and 11.83 m/s in upwind groups. It is observed that the increase in wind speed made the upwind wake effect more pronounced and the power generation in the rear row decreased by~46%. In Figure 4d, the proportion of south wind >15 m/s was 2.0%. Upwind power and wind speed were 1494 KW and 16.16 m/s, respectively, higher than those of 1234 KW and 11.83 m/s in the middle group and 1234 KW and 11.77 m/s in the downwind group. At higher wind conditions, the wake effect perturbed by upwind turbines almost affected the whole wind farm, and the rear turbines could not reach full capacities.   Figure 5 shows the deviation distribution of hourly WRF wind speed correction after bilinear interpolation and nearest replacement. Specifically, we used the observed and simulated data of turbines 06, 43, 50, and 64 for comparison. Intercomparison between WRF simulations and observations for the whole of September shows a root mean square error (RMSE) of 3.73 m/s. By postprocessing the WRF simulations with the two mathe matical algorithms, we found that the two methods all exhibited great performance in presenting the trend of wind speed. The error of bilinear interpolation was reduced by 34.85%, and the average deviation was 2.43 m/s; in comparison, the nearest replacemen method shows a much better result with a deviation of 2.38 m/s, decreasing by 36.19% Results indicated that the implication of the adjacent turbine agreed better with observations than the multipoint interpolation. A possible explanation was that the adjacent wind turbines disturbed the local flow field and introduced more unknown errors, although this methodology introduced data from four grid points.   Figure 5 shows the deviation distribution of hourly WRF wind speed correction after bilinear interpolation and nearest replacement. Specifically, we used the observed and simulated data of turbines 06, 43, 50, and 64 for comparison. Intercomparison between WRF simulations and observations for the whole of September shows a root mean square error (RMSE) of 3.73 m/s. By postprocessing the WRF simulations with the two mathematical algorithms, we found that the two methods all exhibited great performance in presenting the trend of wind speed. The error of bilinear interpolation was reduced by 34.85%, and the average deviation was 2.43 m/s; in comparison, the nearest replacement method shows a much better result with a deviation of 2.38 m/s, decreasing by 36.19%. Results indicated that the implication of the adjacent turbine agreed better with observations than the multipoint interpolation. A possible explanation was that the adjacent wind turbines disturbed the local flow field and introduced more unknown errors, although this methodology introduced data from four grid points. Considering the better performance of the nearest replacement, an AVT method was used to further correct its results (shown in Figure 6). Before the correction, there were three high deviation periods: 1:00-3:00, 11:00-15:00, and 20:00-24:00 after the evening. After processing the AVT corrections, the overall absolute deviation was significantly reduced. In Figure 7, it is observed that when the simulated wind speed increased to 8 m/s, the deviation decreased gradually as well. Whereas when the simulated wind speed continued to increase from 10 m/s, the deviation inversely started to increase [23]. Therefore, it is evident that the correction algorithm was able to make a great contribution in reducing the prediction deviation for the wind speed simulations smaller than 10 m/s, but for those beyond 10 m/s, the corrections were not so significant. In terms of real wind speed distribution, the samples with wind speed below 10 m/s accounted for as high as 86% of the total samples, much higher than that of high wind speed cases. It had a negative deviation of 41% and a positive deviation of 45%, showing no significant difference. Samples higher than 10 m/s accounted for only 14% of the total samples, including 10% negative bias and 4% positive bias, and the number of samples with negative bias played a dominant role. The cumulative wind power for the cases smaller than 10 m/s was 1.82 million KW, while the cumulative powers of the samples beyond 10 m/s was 1.24 million KW, indicating that the accuracy of high wind speed prediction was very important and needed to be further improved. Considering the better performance of the nearest replacement, an AVT method was used to further correct its results (shown in Figure 6). Before the correction, there were three high deviation periods: 1:00-3:00, 11:00-15:00, and 20:00-24:00 after the evening. After processing the AVT corrections, the overall absolute deviation was significantly reduced. In Figure 7, it is observed that when the simulated wind speed increased to 8 m/s, the deviation decreased gradually as well. Whereas when the simulated wind speed continued to increase from 10 m/s, the deviation inversely started to increase [23]. Therefore, it is evident that the correction algorithm was able to make a great contribution in reducing the prediction deviation for the wind speed simulations smaller than 10 m/s, but for those beyond 10 m/s, the corrections were not so significant. In terms of real wind speed distribution, the samples with wind speed below 10 m/s accounted for as high as 86% of the total samples, much higher than that of high wind speed cases. It had a negative deviation of 41% and a positive deviation of 45%, showing no significant difference. Samples higher than 10 m/s accounted for only 14% of the total samples, including 10% negative bias and 4% positive bias, and the number of samples with negative bias played a dominant role. The cumulative wind power for the cases smaller than 10 m/s was 1.82 million KW, while the cumulative powers of the samples beyond 10 m/s was 1.24 million KW, indicating that the accuracy of high wind speed prediction was very important and needed to be further improved.  In the AVT algorithm, the correction effect from the first day to the first 21 days was the best choice. In the experiment of this paper, historical data from 1-3 days and 12-15 days were used as cycles, which can achieve relatively good results. After rolling correction, the mean absolute error of wind speed after 1-3 days of cycle correction was obviously reduced. The margin of error was reduced by 10%. The error of the bilinear interpolation scheme decreased from 2.43 m/s to 2.14 m/s. The nearest replacement scheme error decreased from 2.38 m/s to 2.13 m/s.

Wind Speed Correction Based on Mathematical Statistics Algorithms
The results show that the 15-minute deviation in simulation results can be corrected to the correct range at four different turbines, and the deviation of some periods can be significantly reduced. In the mid-latitude region, the near-stratum wind field mainly considers the large-scale dynamic effect, while Inner Mongolia is located in the region under the common influence of Mongolian cyclones, westerlies, and western Pacific subtropical  In the AVT algorithm, the correction effect from the first day to the first 21 days wa the best choice. In the experiment of this paper, historical data from 1-3 days and 12-1 days were used as cycles, which can achieve relatively good results. After rolling correc tion, the mean absolute error of wind speed after 1-3 days of cycle correction was obv ously reduced. The margin of error was reduced by 10%. The error of the bilinear inter polation scheme decreased from 2.43 m/s to 2.14 m/s. The nearest replacement schem error decreased from 2.38 m/s to 2.13 m/s.
The results show that the 15-minute deviation in simulation results can be correcte to the correct range at four different turbines, and the deviation of some periods can b In the AVT algorithm, the correction effect from the first day to the first 21 days was the best choice. In the experiment of this paper, historical data from 1-3 days and 12-15 days were used as cycles, which can achieve relatively good results. After rolling correction, the mean absolute error of wind speed after 1-3 days of cycle correction was obviously reduced. The margin of error was reduced by 10%. The error of the bilinear interpolation scheme decreased from 2.43 m/s to 2.14 m/s. The nearest replacement scheme error decreased from 2.38 m/s to 2.13 m/s. The results show that the 15-minute deviation in simulation results can be corrected to the correct range at four different turbines, and the deviation of some periods can be significantly reduced. In the mid-latitude region, the near-stratum wind field mainly considers the large-scale dynamic effect, while Inner Mongolia is located in the region under the common influence of Mongolian cyclones, westerlies, and western Pacific subtropical high. Once there is an obvious deviation in the forecast of the location of the local weather system, the wind speed result is bound to have a large deviation [37].

Wind Speed Correction Based on CatBoost Artificial Intelligence Algorithm
Two modes of cumulative learning method (CUM) and cyclic learning method (CYC) of the CatBoost artificial intelligence algorithm were used to correct the WRF simulation results, as shown in Figure 8. The CUM methodology means that all the datasets from before the prediction date serve as the training data for CatBoost, while in the CYC methodology, only the datasets with designated cyclic days prior to the prediction date are used for CatBoost training. In this study, for CUM methodology (Table 3), the data from 1 to 25 September were used for sample training, and the model prediction of the WRF-simulated wind field was based on 26 September data. Then, sample training was carried out using data from 1 to 26 September, and the prediction was made for 27 September. By parity of reasoning, the average RMSE of the wind field every 15 min is shown in Table 3. The correction has a better improvement on WRF forecast results. The four-day average RMSE can be decreased from 3.21 m/s to 2.55 m/s. For CYC methodology (Table 3), the 10day data were designated for sample cyclic training; for example, 16-25 September data were used for sample training, and 26 September was used for model prediction; sample training was conducted from 17 to 26 September, and the model prediction was made on 27 September, and so on. The four-day average RMSE decreased from 3.21 m/s to 2.45 m/s. It is observed that the performance of the circular design is better than the sum correction. After testing the cycle of about 1-30 days, it was found that the cycle of about 10 days had the best performance, which is similar to the cycling scale of atmospheric circulation. Compared with the AVT output, it is observed that the four-day average of 2.55 m/s and 2.45 m/s of CUM_RMSE and CYC_RMSE showed a better performance than the 2.66 m/s of AVT_RMSE.

Summary, Discussions, and Outlooks
In this study, multiple approaches, including the mesoscale WRF model, mathematical statistics algorithms, and CatBoost machine learning algorithm, were adopted to explore the predictability of wind speed in the Gobi grassland wind farm in the Urat area of western Inner Mongolia. We found that the wake effect of front wind turbines significantly reduced the power generated by the rear-row turbines. The higher the wind speed is, the more pronounced the wake effect is on the rear-row group. The AVT postprocessing algorithm combined with bilinear interpolation and nearest replacement effectively optimized the WRF simulations, and CatBoost artificial intelligence algorithm also showed a good perfor-mance in reducing the errors of WRF output. However, under the influence of large-scale circulation, there are still relatively large deviations in the multi correction algorithms.
On the basis of these optimization results, we put forward the following suggestions as the main outlooks of this study: (1) Inner Mongolia is located in the region jointly affected by Mongolian cyclones, westerlies, and western Pacific Subtropical high, so the data assimilation of the northwest and Mongolian plateau should be strengthened in order to improve the simulation accuracy of wind speed; (2) In the boundary layer scheme, there are still major problems in parameterizing explicit variables of small-scale processes, and the results obtained by spatial interpolation or replacement methods not only have large deviations but may even cover the real values. Therefore, it is suggested that an optimized parameterization scheme should be proposed and constructed on the basis of the measurement of boundary layer structure over the underlying surface type and boundary layer structure in different seasons in the Gobi grassland landscape; (3) Local heterogeneous terrain has a pronounced influence on wind speed, and the terrain at wind power sites may have a potential impact on simulation results. It is therefore suggested to introduce high-resolution topographic data within the 30 km domain of the wind farm to improve the simulation precision of the near-surface wind field.