CN118135138A - Terrain modeling method for complex terrain wind field numerical simulation - Google Patents
Terrain modeling method for complex terrain wind field numerical simulation Download PDFInfo
- Publication number
- CN118135138A CN118135138A CN202410571925.5A CN202410571925A CN118135138A CN 118135138 A CN118135138 A CN 118135138A CN 202410571925 A CN202410571925 A CN 202410571925A CN 118135138 A CN118135138 A CN 118135138A
- Authority
- CN
- China
- Prior art keywords
- terrain
- mesoscale
- microscale
- topography
- fusion
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 38
- 238000004088 simulation Methods 0.000 title claims abstract description 33
- 238000012876 topography Methods 0.000 claims abstract description 66
- 230000004927 fusion Effects 0.000 claims abstract description 51
- 238000004364 calculation method Methods 0.000 claims abstract description 31
- 230000008878 coupling Effects 0.000 claims abstract description 13
- 238000010168 coupling process Methods 0.000 claims abstract description 13
- 238000005859 coupling reaction Methods 0.000 claims abstract description 13
- 230000008569 process Effects 0.000 claims abstract description 6
- 230000006870 function Effects 0.000 claims description 25
- 230000007704 transition Effects 0.000 claims description 10
- 238000009826 distribution Methods 0.000 claims description 7
- 238000011835 investigation Methods 0.000 claims description 2
- 238000000605 extraction Methods 0.000 abstract description 4
- 238000010586 diagram Methods 0.000 description 13
- 238000007500 overflow downdraw method Methods 0.000 description 9
- 238000004590 computer program Methods 0.000 description 7
- 230000005540 biological transmission Effects 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 3
- 238000013075 data extraction Methods 0.000 description 3
- 230000014759 maintenance of location Effects 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000003860 storage Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 238000010276 construction Methods 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000010248 power generation Methods 0.000 description 2
- 241000208340 Araliaceae Species 0.000 description 1
- 235000005035 Panax pseudoginseng ssp. pseudoginseng Nutrition 0.000 description 1
- 235000003140 Panax quinquefolius Nutrition 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000052 comparative effect Effects 0.000 description 1
- 238000005520 cutting process Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 235000008434 ginseng Nutrition 0.000 description 1
- 238000004519 manufacturing process Methods 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 238000000926 separation method Methods 0.000 description 1
- 238000010008 shearing Methods 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The invention discloses a terrain modeling method for complex terrain wind field numerical simulation, which comprises the following steps: step 1, determining a microscale calculation domain for mesomicroscale coupling, and acquiring microscale terrains with higher precision; step 2, acquiring a mesoscale topography or determining a flat topography; step 3, determining relevant parameters of terrain fusion, calculating a fusion function, and carrying out mesoscale terrain fusion or flat terrain fusion; and step 4, obtaining a topography file of microscale numerical simulation for microscale calculation domain grid generation. According to the invention, on the basis of the traditional terrain fusion function, the extraction and input process of the mesoscale data is considered, the terrain fusion function is corrected, the obtained mesoscale fusion terrain provides good terrain boundary conditions for microscale calculation, the data information from the mesoscale can be well transmitted, and the precision of the fine resource assessment of the complex terrain wind power plant is improved.
Description
Technical Field
The invention relates to the technical field of wind power generation, in particular to a terrain modeling method for complex terrain wind field numerical simulation.
Background
The wind farm is designed to be in two directions of future wind farm construction, but is popular due to high construction and transmission costs of the deep sea wind farm. The phenomena of high turbulence, high wind shearing, flow separation and the like exist in the complicated mountain wind field, and the method is a great challenge for the assessment of refined wind resources. In complex terrains, due to the influence of terrains, the wind towers are insufficient in representativeness, the local machine position needs to be evaluated by means of mesoscale data, and a mesoscale numerical weather forecast mode (NWP) +a microscale Computational Fluid Dynamics (CFD) mode has become an important means for complex mountain unsteady flow field simulation and unit load and power generation evaluation.
When complex terrain flow simulation is based on mesoscale coupling, terrain modeling is a key step. The terrain modeling needs to consider the following points: firstly, the range of a calculation domain is reasonably evaluated, the accuracy and the efficiency of calculation are guaranteed, secondly, the accuracy of microscale terrains is guaranteed, thirdly, the convenience of microscale calculation domain grid modeling is guaranteed, and thirdly, the accuracy of mesoscale wind speed, temperature and the like information transmission is guaranteed. When the mesoscale simulation is adopted, the common topographic resolution is 30' (about 900 m), and the topographic resolution can reach 10-50m when the microscale calculation is performed, so that the mesoscale topographic fusion is required to be performed so as to ensure the consistency of the mesoscale output and the microscale input. The modeling of the complex terrain wind field grid can adopt commercial software (such as ICEM) or open source codes such as OpenFOAM, and the terrain input formats required by the two software are all stl formats, so that the universality of the fused terrain is ensured, and the fused terrain is stored in the stl format.
At present, most of complex terrain modeling adopts commercial software, format conversion, projection conversion and cutting work are required to be continuously carried out, and the operation flow is quite complicated. And when the grids are drawn, the topography also needs to be cut and repaired for a plurality of times, and a great amount of time is consumed. When the traditional fusion function is adopted for terrain fusion, fusion is directly started from the boundary, and the extraction of mesoscale data is influenced. In order to save computing resources, it is necessary to ensure that the fusion transition area is not too large, and in order to ensure the accuracy of mesoscale data extraction, a part of mesoscale topography needs to be left near the computational domain boundary for extracting the mesoscale data, which is not considered by the conventional fusion function.
Disclosure of Invention
The invention aims to provide a terrain modeling method for complex terrain wind field numerical simulation, which considers the extraction and input process of mesoscale data on the basis of the traditional terrain fusion function, corrects the terrain fusion function, and the obtained mesoscale fusion terrain provides good terrain boundary conditions for microscale calculation, can better transmit data information from the mesoscale, and improves the precision of the fine resource assessment of the complex terrain wind field.
In order to achieve the technical purpose, the invention adopts the following technical scheme:
A terrain modeling method for complex terrain wind field numerical simulation, the terrain modeling method comprising the steps of:
Step 1, obtaining microscale topography:
Determining a microscale calculation domain and a reference point position for microscale calculation according to the actual terrain altitude difference of the researched wind farm area and the position of a wind tower, and obtaining microscale terrain of the microscale calculation domain; the microscale topography is larger than the region of the wind farm under investigation;
Step 2, acquiring a mesoscale terrain or determining a flat terrain: extracting longitude and latitude and altitude data corresponding to nested topography of the innermost layer of the region where the anemometer tower is positioned from mesoscale result data of the researched wind farm region, and extracting mesoscale topography corresponding to a microscale calculation domain; the resolution of the mesoscale topographic data is lower than the resolution of the microscale topographic data; if the method is fused with the flat terrain, determining the altitude of the flat terrain;
step3, merging the mesoscale terrain or merging the mesoscale terrain with the flat terrain:
The reservation ranges in the four directions of the west, the east, the south and the north during the merging of the mesoscale topography are respectively determined according to the mesoscale topography resolution: d w,w、dw,e、dw,s and d w,n,dw,w、dw,e、dw,s and d w,n are both greater than 0.5 x d s and less than d s,ds are mesoscale topographic resolution;
Determining the transition range from low-resolution terrain to high-resolution terrain in the west, east, south and north directions respectively: d c,w、dc,e、dc,s and d c,n,dc,w、dc,e、dc,s and d c,n are both greater than 0.5 x d s less than 1.5d s;
Calculating a terrain fusion function F b:
；
；
；
；
；
Wherein, F b,w、Fb,e、Fb,s and F b,n are respectively fusion functions of the terrain in the west, east, south and north directions, x and y are respectively the x-axis coordinate and the y-axis coordinate of one point of the microscale computational domain, and x min、xmax、ymin and y max are respectively the boundary coordinates of the microscale computational domain on the x-axis and the y-axis; d f is the reserved range of four directions when the plate is fused with the flat terrain;
And (3) fusing the altitudes of the mesoscale terrain to obtain fused altitude distribution:
；
Wherein z b is the altitude of the terrain after fusion; z h is the altitude of the microscale terrain; z l is the altitude of the low-precision terrain, which is the mesoscale terrain or the flat terrain; f b is a terrain fusion function;
and step 4, obtaining a topography file of microscale numerical simulation for microscale calculation domain grid generation.
Further, the calculated domain range is 6h×4h×h, where H is the atmospheric boundary layer thickness, and is determined by the atmospheric stability, 6H corresponds to the direction of flow, 4H corresponds to the direction perpendicular to the direction of flow, and H corresponds to the perpendicular direction.
Further, in the step 1, the micro-scale topographic resolution range is 10m-30m, and after the micro-scale topographic data is obtained, the micro-scale topographic data is converted into a tif format;
In the process of determining the microscale computing domain, the position information of a reference point during mesoscale coupling is determined, the reference point is generally positioned at the center of the microscale computing domain, and microscale topographic data are stored in a geographic coordinate mode.
Further, in step 2, the mesoscale topography is in USGS format.
Further, in the step 3, when the terrains are fused, the microscale terrains are fused with the flat-plate terrains or the mesoscale terrains; d f takes 0 if the microscale topography and the mesoscale topography are fused; if the microscale topography is fused with the slab topography, d w,w、dw,e、dw,s and d w,n take 0.
Further, in step 4, the fused terrain format is stl.
Compared with the prior art, the invention has the following beneficial effects:
Compared with the actual wind farm numerical simulation result and radar measured data, the terrain modeling method for the complex terrain wind farm numerical simulation is simple to operate, convenient to automate, capable of saving a large amount of time cost, high in efficiency and good in precision, the established fusion terrain has good consistency with the mesoscale terrain at the boundary, the mesoscale output and the microscale input have good consistency, and the numerical simulation result has good consistency with the measured data. In addition, the method can be completed by adopting the python language, has high modeling efficiency, and can ensure the consistency of the mesoscale output and the microscale input.
Drawings
FIG. 1 is a comparative schematic diagram of complex terrain obtained by an unfused, old and new fusion method;
FIG. 2 is a schematic view of a microscale topography of a wind farm (UTM projection, altitude cloud);
FIG. 3a is a schematic representation of the topography of a wind farm (altitude cloud) for microscale computational domain modeling;
FIG. 3b is a schematic representation of the topography (three-dimensional topography, top view) of a wind farm for microscale computational domain modeling;
FIG. 4a is a schematic representation of a mesoscale terrain (altitude cloud) for a wind farm;
FIG. 4b is a schematic representation of a mesoscale terrain (three-dimensional terrain, top view) for a wind farm;
FIG. 5 is a graph of a terrain fusion function for a wind farm;
FIG. 6a is a schematic view of a fused terrain (altitude cloud) for a wind farm using the terrain modeling method of the present invention;
FIG. 6b is a schematic view (three-dimensional terrain, top view) of a fused wind farm obtained by the terrain modeling method of the present invention;
FIG. 7a is a schematic representation of the terrain (altitude cloud) of a fused wind farm using prior art methods;
FIG. 7b is a schematic representation of the terrain (three-dimensional terrain, top view) of a fused wind farm using prior art methods;
FIG. 8 is a schematic illustration of the altitude distribution of terrain on the west side obtained by the two fusion methods of the present invention and the prior art;
FIG. 9a is a schematic diagram of a consistency verification result (u-direction, east-west) of an inlet wind speed input and a mesoscale wind speed output at the time of mesoscale coupling;
FIG. 9b is a schematic diagram of the results of a consistency verification of the input of the inlet wind speed and the output of the mesoscale wind speed (v direction, from south to north) at the time of mesoscale coupling;
FIG. 9c is a schematic diagram of the results of a verification of the consistency of the input of the inlet wind speed with the output of the mesoscale wind speed (w direction, bottom to top) at the time of mesoscale coupling;
FIG. 10 is a schematic diagram of a numerical simulation computational grid created using the fused terrain;
FIG. 11a is a schematic diagram of the comparison of the results of a simulation of the radar observations of a wind farm with mesoscale coupling values (wind profile);
FIG. 11b is a graphical representation of the results of a comparison of radar observations of a wind farm with mesoscale coupling numerical simulations (turbulence intensity profile);
FIG. 12 is a flow chart of a terrain modeling method for complex terrain wind field numerical modeling in accordance with the present invention.
Detailed Description
Embodiments of the present invention are described in further detail below with reference to the accompanying drawings.
Referring to fig. 12, an embodiment of the present invention discloses a terrain modeling method for complex terrain wind field numerical simulation, including:
Step 1, obtaining microscale topography: determining a calculation domain for micro-scale calculation according to the actual terrain altitude difference of a specific wind farm and the position of a wind measuring tower, and then acquiring micro-scale terrain with higher precision;
Step 2, acquiring a mesoscale terrain or determining a flat terrain: extracting the mesoscale topography of the area according to the mesoscale result data of the research area; if the microscale terrain is fused with the flat terrain, determining the elevation of the flat terrain;
Step 3, merging the mesoscale terrain or merging the mesoscale terrain with the flat terrain: determining a reserved mesoscale range and a region for mesoscale terrain fusion transition, and calculating a terrain fusion function to further obtain the fused terrain altitude distribution;
；
Wherein z b is the altitude of the terrain after fusion; z h is the altitude of the high-precision microscale terrain, and the resolution of the microscale terrain is generally 10-30m; z l is the altitude of the low-precision terrain, which can be the terrain data in the mesoscale WRF or the flat terrain data; f b is a terrain fusion function.
The formula of F b is as follows:
；
；
；
；
；
Where F b,w、Fb,e、Fb,s and F b,n are fusion functions of the terrain in the west, east, south and north directions, x and y are x-axis coordinates and y-axis coordinates of one point of the micro-scale computation domain, x min、xmax、ymin and y max are coordinate extrema of the micro-scale computation domain on the x-axis and y-axis, d w,w、dw,e、dw,s and d w,n are retention ranges of the terrain in the four directions when fused with the medium-scale terrain, d f is a retention range of the terrain in the flat-plate terrain, and d c,w、dc,e、dc,s and d c,n are transition ranges of the terrain in the four directions from the low-resolution terrain to the high-resolution terrain.
And step 4, obtaining a topography file of microscale numerical simulation for microscale calculation domain grid generation.
1. Obtaining microscale topography
In step 1, the micro-scale topography format is elevation data, which may be SRTM format, or another format, and then converted into tif format.
The micro-scale calculation domain range is comprehensively determined according to the unit position and the terrain complexity of specific attention. The recommended range is 6H×4H×H, wherein the value of H is the thickness of the atmospheric boundary layer, and is determined by the atmospheric stability, 6H corresponds to the direction of flow, 4H corresponds to the direction perpendicular to the direction of flow, and H corresponds to the perpendicular direction. If there is no definite dominant wind direction, 6H may be used in the horizontal direction.
The micro-scale topography can download topography with 30m resolution from the SRTM official network, and can obtain topography with higher precision from measured data. The microscale terrain should be larger (the latitude and longitude of the boundary need to be determined) than the region of the wind farm under study, because a small part needs to be cut off when the grid is generated. After the micro-scale topography is obtained, it needs to be converted into tif format. In the process of determining the microscale calculation domain, the position (including the coordinates and longitude and latitude of the reference point) of the reference point during mesoscale coupling is required to be determined, and the microscale elevation topographic data at the moment only need to be stored in a geographic coordinate mode (EPSG:4326) without projection conversion.
2. Acquiring mesoscale topography or determining slab topography
The mesoscale topographic data may be extracted from the mesoscale computation result file (wrfout x.nc files) to obtain longitude, latitude and altitude data corresponding to the mesoscale topographic data. Because the innermost layer nesting resolution is highest during mesoscale calculation, the topographic data source is also extracted from the innermost layer result file, and the topographic raw data format is USGS. When the method is fused with the flat terrain, the elevation of the flat terrain needs to be determined;
3. mesoscale terrain fusion or slab terrain fusion
Before terrain fusion, the retention range of the mesoscale terrain in four directions, namely four parameters d w,w、dw,e、dw,s and d w,n, are required to be determined according to the mesoscale terrain resolution, wherein the four parameters are larger than 0.5 x d s and smaller than d s（ds and are the mesoscale terrain resolution. The transition range from low to high resolution terrain is then determined, i.e., four parameters d c,w、dc,e、dc,s and d c,n, each of which is greater than 0.5 x d s and less than 1.5d s. The boundary at which the transition starts is preferably located at the center of the mesoscale grid. The research area needs to be far away from the wake area of the fusion area so as to ensure that the fusion terrain does not influence the flow field of the downstream research object. After the parameters are determined, a terrain fusion function can be calculated:
；
；
；
；
；
Where F b,w、Fb,e、Fb,s and F b,n are fusion functions of the terrain in the west, east, south and north directions, respectively, x and y are the x-axis and y-axis coordinates of one of the points of the micro-scale computational domain, and x min、xmax、ymin and y max are the boundary coordinates of the micro-scale computational domain on the x-axis and y-axis, respectively.
After the fusion function is obtained, the fusion of the altitude of the mesoscale terrain can be carried out, and the altitude distribution after fusion is obtained:
；
Where z b is the altitude of the terrain after fusion, z h is the high-precision microscale terrain, z l is the mesoscale WRF terrain data, and F b is the terrain fusion function.
When the topography is fused, the method can only select one mode to be fused with the flat topography or the mesoscale topography. D f takes 0 if the fusion is carried out with the mesoscale topography; if fused with the flat terrain, d w,w、dw,e、dw,s and d w,n take 0.
4. Obtaining a topographic file of microscale numerical simulation
The fused topography is stored in stl format and is consistent with the mesoscale data, and the projection mode of the fused topography is UTM. The computational domain boundary is preferably positioned at the starting position of the terrain fusion transition during grid modeling, and the grid generation can be completed by adopting commercial software or self-programming.
The comparison of the complex topographic schematic diagrams obtained by the unfused and old fusion methods and the new fusion method is shown in fig. 1.
Case (B)
Taking a certain wind farm as an example, mesoscale topography fusion is carried out, and wind farm wind characteristic analysis is carried out.
Step 1, obtaining microscale topography. Firstly, according to the position of the wind measuring tower, high-precision elevation data are obtained, as shown in fig. 1, then a calculation domain (red square frame in fig. 2) for micro-scale calculation is determined according to the characteristics of the researched object, mountain height and the like, and the topography of the area is extracted, as shown in fig. 3a and 3b (EPSG: 32648).
And 2, obtaining the mesoscale topography. Firstly, longitude and latitude and altitude data corresponding to the innermost nested topography of the region where the anemometer tower is located are extracted from the mesoscale result data, and then mesoscale topography corresponding to a microscale calculation domain is obtained, as shown in fig. 4a and 4b (EPSG: 32648).
And 3, merging the mesoscale topography. First, the reserved range of the mesoscale topography in four directions needs to be determined, in this case, d w,w=dw,e=dw,s=dw,n= 750,750, 750 m; then the transition range of the low resolution terrain to the high resolution terrain, in this case, d c,w=dc,e=dc,s=dc,n=1000,df =0; finally, calculating a terrain fusion function F b:
；
；
；
；
；
Where F b,w、Fb,e、Fb,s and F b,n are fusion functions of the terrain in the west, east, south and north directions, x and y are coordinates of the micro-scale computational domain, and x min、xmax、ymin and y max are coordinate extrema of the micro-scale computational domain, in this case, x min=317423.46,xmax=317423.46,ymin=2638514.08,ymax = 2643364.08, which are coordinates under UTM projection (i.e., EPSG: 32648). The calculated fusion function distribution cloud chart is shown in fig. 5.
And further obtaining the fused terrain altitude distribution:
；
Wherein z b is the altitude of the terrain after fusion; z h is the high-precision microscale terrain altitude; z l is the terrain data in low accuracy terrain altitude, i.e., mesoscale WRF. The fused topography is shown in fig. 6a and 6b.
In the old terrain fusion method, d w,w=dw,e=dw,s=dw,n =0 m, the fused terrain is as shown in fig. 7a and 7b.
The biggest difference between the new method and the old method is that the terrain error exists during the extraction of the mesoscale wind ginseng, taking the terrain in fig. 4b, 6b and 7b as an example, the terrain altitude on the west side (near x min) is extracted, as shown in fig. 8.
As can be seen from fig. 8, the unfused terrain and the terrain obtained by the old fusion method have certain errors with the mesoscale terrain at the boundary, and such errors will affect the accuracy of the mesoscale data extraction. The obtained topography adopting the new fusion method has better consistency with the mesoscale topography at the boundary. To further verify the accuracy of the mesoscale data extraction, taking the micro-scale computing domain southwest corner topography as an example, a comparison of wind profiles as shown in fig. 9a to 9c is obtained.
As can be seen from fig. 9a to 9c, the terrain modeling method can ensure the consistency of the transmission from the mesoscale data to the microscale data, thereby ensuring the accuracy of the simulation input of the mesoscale coupling wind field.
And 4, obtaining a topography file of microscale numerical simulation. The fused topography is stored in stl format and is consistent with the mesoscale data, and the projection mode of the fused topography is UTM. The numerical simulation calculation domain and the grid established by the fused terrain are shown in fig. 10, the calculation domain range is slightly smaller than the terrain range (the calculation domain boundary is positioned at the grid center of the mesoscale terrain during numerical simulation, namely the initial position of the terrain fusion transition), and the grid near the ground is encrypted.
As can be seen from an analysis of fig. 11a and 11 b:
The wind speed and turbulence intensity profile obtained by the mesoscale coupling numerical simulation based on the new complex topography fusion method has good consistency with radar observation data, and can simulate the special wind profile and turbulence profile of the complex topography better.
In summary, based on the requirements of terrain modeling in complex terrain numerical simulation, the invention provides a novel complex terrain modeling method for mesoscale coupling. The terrain obtained by the novel fusion method has better consistency with the mesoscale terrain at the boundary, ensures the consistency of mesoscale data to microscale data transmission, can further realize the fine evaluation of wind resources of the wind power plant with complex terrain, improves the generated energy of the wind turbine generator and reduces the economic loss. The invention has the advantages of simple operation, easy coding, simple generation and saving of a great deal of time cost.
It will be appreciated by those skilled in the art that embodiments of the present application may be provided as a method, system, or computer program product. Accordingly, the present application may take the form of an entirely hardware embodiment, an entirely software embodiment or an embodiment combining software and hardware aspects. Furthermore, the present application may take the form of a computer program product embodied on one or more computer-usable storage media (including, but not limited to, disk storage, CD-ROM, optical storage, and the like) having computer-usable program code embodied therein. The scheme in the embodiment of the application can be realized by adopting various computer languages, such as object-oriented programming language Java, an transliteration script language JavaScript and the like.
The present application is described with reference to flowchart illustrations and/or block diagrams of methods, apparatus (systems) and computer program products according to embodiments of the application. It will be understood that each flow and/or block of the flowchart illustrations and/or block diagrams, and combinations of flows and/or blocks in the flowchart illustrations and/or block diagrams, can be implemented by computer program instructions. These computer program instructions may be provided to a processor of a general purpose computer, special purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, create means for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may also be stored in a computer-readable memory that can direct a computer or other programmable data processing apparatus to function in a particular manner, such that the instructions stored in the computer-readable memory produce an article of manufacture including instruction means which implement the function specified in the flowchart flow or flows and/or block diagram block or blocks.
These computer program instructions may 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 such that the instructions which execute on the computer or other programmable apparatus provide steps for implementing the functions specified in the flowchart flow or flows and/or block diagram block or blocks.
Claims (6)
1. A terrain modeling method for complex terrain wind field numerical simulation, characterized in that the terrain modeling method comprises the following steps:
Step 1, obtaining microscale topography:
Determining a microscale calculation domain and a reference point position for microscale calculation according to the actual terrain altitude difference of the researched wind farm area and the position of a wind tower, and obtaining microscale terrain of the microscale calculation domain; the microscale topography is larger than the region of the wind farm under investigation;
Step 2, acquiring a mesoscale terrain or determining a flat terrain: extracting longitude and latitude and altitude data corresponding to nested topography of the innermost layer of the region where the anemometer tower is positioned from mesoscale result data of the researched wind farm region, and extracting mesoscale topography corresponding to a microscale calculation domain; the resolution of the mesoscale topographic data is lower than the resolution of the microscale topographic data;
step3, merging the mesoscale terrain or merging the mesoscale terrain with the flat terrain:
The reservation ranges in the four directions of the west, the east, the south and the north during the merging of the mesoscale topography are respectively determined according to the mesoscale topography resolution: d w,w、dw,e、dw,s and d w,n,dw,w、dw,e、dw,s and d w,n are both greater than 0.5 x d s and less than d s,ds are mesoscale topographic resolution;
Determining the transition range from low-resolution terrain to high-resolution terrain in the west, east, south and north directions respectively: d c,w、dc,e、dc,s and d c,n,dc,w 、dc,e、dc,s and d c,n are both greater than 0.5 x d s less than 1.5d s;
Calculating a terrain fusion function F b:
Wherein, F b,w、Fb,e、Fb,s and F b,n are respectively fusion functions of the terrain in the west, east, south and north directions, x and y are respectively the x-axis coordinate and the y-axis coordinate of one point of the microscale computational domain, and x min、xmax、ymin and y max are respectively the boundary coordinates of the microscale computational domain on the x-axis and the y-axis; d f is the reserved range of four directions when the plate is fused with the flat terrain;
And (3) fusing the altitudes of the mesoscale terrain to obtain fused altitude distribution:
Wherein z b is the altitude of the terrain after fusion; z h is the altitude of the microscale terrain; z l is the altitude of the low-precision terrain, which is the mesoscale terrain or the flat terrain; f b is a terrain fusion function;
and step 4, obtaining a topography file of microscale numerical simulation for microscale calculation domain grid generation.
2. The terrain modeling method for complex terrain wind field numerical simulation according to claim 1, wherein the calculation domain range is 6hχ4hχh, wherein the value of H is the atmospheric boundary layer thickness, and is determined by the atmospheric stability, 6H corresponds to the direction of the flow direction, 4H corresponds to the direction perpendicular to the flow direction, and H corresponds to the perpendicular direction.
3. The terrain modeling method for complex terrain wind field numerical simulation according to claim 1, wherein in step 1, the micro-scale terrain resolution range is 10m-30m, and after the micro-scale terrain data is acquired, the micro-scale terrain data is converted into tif format;
in the process of determining a microscale computing domain, the position information of a reference point during mesoscale coupling is determined, the reference point is positioned at the center of the microscale computing domain, and microscale topographic data are stored in a geographic coordinate mode.
4. The terrain modeling method for complex terrain wind field numerical simulation of claim 1, wherein in step 2, the mesoscale terrain is in USGS format.
5. The terrain modeling method for complex terrain wind field numerical simulation according to claim 1, wherein in step 3, when the terrain is fused, the micro-scale terrain is fused with the flat-plate terrain or the meso-scale terrain; d f takes 0 if the microscale topography and the mesoscale topography are fused; if the microscale topography is fused with the slab topography, d w,w、dw,e、dw,s and d w,n take 0.
6. A terrain modeling method for use in complex terrain wind field numerical simulation according to claim 1, wherein in step 4, the fused terrain format is stl.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410571925.5A CN118135138B (en) | 2024-05-10 | 2024-05-10 | Terrain modeling method for complex terrain wind field numerical simulation |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410571925.5A CN118135138B (en) | 2024-05-10 | 2024-05-10 | Terrain modeling method for complex terrain wind field numerical simulation |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN118135138A true CN118135138A (en) | 2024-06-04 |
| CN118135138B CN118135138B (en) | 2024-07-02 |
Family
ID=91232141
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410571925.5A Active CN118135138B (en) | 2024-05-10 | 2024-05-10 | Terrain modeling method for complex terrain wind field numerical simulation |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118135138B (en) |
Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190338758A1 (en) * | 2017-03-06 | 2019-11-07 | Xinjiang Goldwind Science & Technology Co., Ltd. | Method and device for calculating power generation of wind farm |
| CN114580310A (en) * | 2020-11-28 | 2022-06-03 | 北京零壹橡树科技有限公司 | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) |
| US20230274520A1 (en) * | 2022-02-25 | 2023-08-31 | National Institute of Meteorological Sciences | Method of calculating high-resolution 3d wind numerical information considering detailed topographic information |
| CN117390992A (en) * | 2023-10-25 | 2024-01-12 | 上海能源科技发展有限公司 | Method for predicting refined wind field under complex terrain |
-
2024
- 2024-05-10 CN CN202410571925.5A patent/CN118135138B/en active Active
Patent Citations (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20190338758A1 (en) * | 2017-03-06 | 2019-11-07 | Xinjiang Goldwind Science & Technology Co., Ltd. | Method and device for calculating power generation of wind farm |
| CN114580310A (en) * | 2020-11-28 | 2022-06-03 | 北京零壹橡树科技有限公司 | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) |
| US20230274520A1 (en) * | 2022-02-25 | 2023-08-31 | National Institute of Meteorological Sciences | Method of calculating high-resolution 3d wind numerical information considering detailed topographic information |
| CN117390992A (en) * | 2023-10-25 | 2024-01-12 | 上海能源科技发展有限公司 | Method for predicting refined wind field under complex terrain |
Non-Patent Citations (1)
| Title |
|---|
| 程雪玲;胡非;曾庆存;: "复杂地形风场的精细数值模拟", 气候与环境研究, no. 01, 20 January 2015 (2015-01-20) * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN118135138B (en) | 2024-07-02 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11168667B2 (en) | Method and device for calculating power generation of wind farm | |
| CN110298115B (en) | Wind field power downscaling method based on simplified terrain aerodynamic parameters | |
| CN102663251B (en) | Physical prediction method for wind power station power based on computational fluid mechanics model | |
| CN104331621A (en) | Wind resource computation method | |
| CN110009736B (en) | Method, device, equipment and storage medium for establishing three-dimensional wake model | |
| Malz et al. | Drag-mode airborne wind energy vs. wind turbines: An analysis of power production, variability and geography | |
| CN114580310A (en) | Method for realizing scale reduction processing of WRF (hand-wrenching simulation) wind field based on PALM (PALM fiber laser) | |
| CN116341410B (en) | Typical terrain wind field simulation method and system considering different stabilities | |
| CN111666725A (en) | Anemometer tower planning and site selection method and system suitable for wind power plant with non-complex terrain | |
| CN115329667A (en) | Method, device and equipment for determining point location of fan and storage medium | |
| CN114707437A (en) | Method for establishing three-dimensional full wake flow model of horizontal axis wind turbine | |
| CN105184667B (en) | The method of dual nested simulation wind farm wind velocity distribution | |
| CN114139391A (en) | Wind energy resource assessment method and device based on separation vortex turbulence model and storage medium | |
| Rahmani et al. | Wind power assessment and site matching of wind turbines in Lootak of Zabol | |
| CN118135138B (en) | Terrain modeling method for complex terrain wind field numerical simulation | |
| CN113051845B (en) | Method, system, equipment and storage medium for visually evaluating real-time wind resources of in-service mountain wind power plant | |
| CN115586591B (en) | Automatic forecasting method and device for storm surge | |
| CN111639381A (en) | Urban building solar energy resource evaluation information system and working method thereof | |
| CN116050291A (en) | Wind farm large-scale weather prediction downscaling method and system under typhoon condition | |
| CN113158598B (en) | Foundation type wind power plant flow field CFD partition calculation method, device, equipment and medium | |
| CN110925147B (en) | Method for tracking wake flow of wind driven generator | |
| CN112800697A (en) | Method, system and device for planning and site selection of wind measuring tower in early stage of wind power plant and storage medium | |
| Simões et al. | A first methodology for wind energy resource assessment in urbanised areas in Portugal | |
| KR102385673B1 (en) | Method and apparatus for numerical information production of high-resolusion photovoltaic resources concerning detailed topographical effects | |
| Van Thin et al. | Study to evaluate the effect of terrain surface on performance of a wind farm in Ninh Thuan province, Vietnam |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant |