CN109165476A - A kind of modeling method and simulation of wind method of modularization wind-field model - Google Patents

A kind of modeling method and simulation of wind method of modularization wind-field model Download PDF

Info

Publication number
CN109165476A
CN109165476A CN201811203630.3A CN201811203630A CN109165476A CN 109165476 A CN109165476 A CN 109165476A CN 201811203630 A CN201811203630 A CN 201811203630A CN 109165476 A CN109165476 A CN 109165476A
Authority
CN
China
Prior art keywords
wind
model
data
submodel
information
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.)
Pending
Application number
CN201811203630.3A
Other languages
Chinese (zh)
Inventor
肖凯
谢文平
罗啸宇
聂铭
雷旭
姚博
曾志鹏
李朝
肖仪清
赵子涵
钟继
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Guangdong Power Grid Co Ltd
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Original Assignee
Guangdong Power Grid Co Ltd
Electric Power Research Institute of Guangdong Power Grid Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Guangdong Power Grid Co Ltd, Electric Power Research Institute of Guangdong Power Grid Co Ltd filed Critical Guangdong Power Grid Co Ltd
Priority to CN201811203630.3A priority Critical patent/CN109165476A/en
Publication of CN109165476A publication Critical patent/CN109165476A/en
Pending legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • Computer Hardware Design (AREA)
  • Evolutionary Computation (AREA)
  • Geometry (AREA)
  • General Engineering & Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

This application discloses a kind of modeling method of modularization wind-field model and simulation of wind methods, the application is solved by batch and post-processing, establish the standard wind field database model of the air speed data comprising target area difference wind angle, any position, the standard wind field database model established based on the application, real-time short-term or long-term general mood marquis's characterisitic parameter in combining target area, can quickly obtain the wind speed field distribution simulation result of target area.Cumbersome Modeling Calculation step is greatly simplified by the modularization of wind speed field model, solves the analog simulation complex steps of existing complex topography wind field analogy method, and low efficiency promotes difficult technical problem.

Description

A kind of modeling method and simulation of wind method of modularization wind-field model
Technical field
This application involves the modeling methods and simulation of wind side of Wind Engineering field more particularly to a kind of modularization wind-field model Method.
Background technique
China's Coastal Areas terrain environment is complicated, mostly mountain and hill landforms.The gradient, the height of Complex Mountain landform are planted Quilt and occlusion effect can all influence to flow through the wind field distribution situation in the region.The influence of this landform, geomorphological environment to wind field can To be indicated with the wind speed of nondimensionalization than size.Computational fluid dynamics (computation fluid dynamics, CFD) Can the air-flow in simulation of atmospheric boundary layer generated in complicated landform environment acceleration, separation, around phenomena such as, at present by Wind energy assesses unit, Electric Design unit, and city planning design unit uses.However, predicting certain using CFD analogy method at present When the wind field in one area, generally requires since terrain modeling, then to carry out grid drawing, boundary condition definition, solves parameter Setting and post-processing etc. a series of step, the analog simulation step for resulting in existing complex topography wind field analogy method are numerous Trivial, low efficiency promotes difficult technical problem.
Summary of the invention
It is existing for solving this application provides a kind of modeling method of modularization wind-field model and simulation of wind method The analog simulation complex steps of complex topography wind field analogy method, low efficiency promote difficult technical problem.
The application first aspect provides a kind of modeling method of modularization wind-field model, comprising:
S1: being based on global digital elevation model, and target area is divided into several continuous square subregions, is established The corresponding three-dimensional terrain model of the subregion;
S2: it is based on the three-dimensional terrain model, foundation includes the spatial data of the three-dimensional terrain model, wind direction The standard database submodel of data and wind speed than data;
S3: circulation executes step S2, establishes the corresponding standard database of each three-dimensional terrain model respectively Submodel;
S4: each standard database submodel is integrated, and obtains the standard wind field data of the target area Library model.
Preferably, the step S1 is specifically included:
S11: being based on global digital elevation model, and target area is divided into several continuous square subregions, and The computational domain of the subregion is vertically stretched to sustained height threshold value according to preset blocking rate, obtains the subregion pair The three-dimensional terrain model answered.
Preferably, the step S11 is specifically included:
S111: being based on global digital elevation model, and target area is divided into several continuous square subregions, In, the center of the subregion is divided into core space, expansion area and transition region outward;
S112: according to the global digital elevation model, the core space, the expansion area and the transition are determined respectively Latitude and longitude coordinates information, absolute elevation information, ground mulching type information and the leaf area density information in area;
S113: merge latitude and longitude coordinates information, the absolute elevation letter of the core space, the expansion area and the transition region Breath, ground mulching type information and leaf area density information, generate dimensional topography curved surface, and according to the dimensional topography curved surface and Preset blocking rate determines the height threshold of the computational domain of the subregion, and by the computational domain be vertically stretched to it is described The equal height of height threshold obtains the corresponding three-dimensional terrain model of the subregion.
Preferably, the step S2 is specifically included:
S21: it is based on the three-dimensional terrain model, obtains the roughness of ground surface information of the computational domain of the subregion, and root According to fluid dynamics governing equation, the corresponding relationship of the roughness of ground surface information and aerodynamic resistance factor is established;
S22: the calculating parameter of Aerodynamics Model is obtained, and simultaneously output model is calculated according to the calculating parameter batch Wind speed compares data, wherein the calculating parameter includes: turbulence model type, boundary condition, convergence threshold value;
S23: according to the model wind speed than data, foundation include the three-dimensional terrain model spatial data, The standard database submodel of wind direction data and wind speed than data.
The application second aspect provides a kind of wind of standard wind field database model established based on foregoing invention content Field analogy method, comprising:
S6: according to the latitude and longitude coordinates to simulated domain, called from the standard wind field database model with it is described to The corresponding standard database submodel of simulated domain;
S7: obtaining the practical wind parameter to simulated domain, and the practical wind parameter is substituted into the standard database Submodel carries out operation, obtains the simulation of wind result data to simulated domain, wherein the practical wind parameter includes: The shot and long term wind climate parameter to simulated domain refers to wind speed parameter to the practical of simulated domain with described.
The application third aspect provides a kind of model building device of modularization wind-field model, comprising:
It is continuous to be divided into several for being based on global digital elevation model by initial landform modeling unit for target area Square subregion, establish the corresponding three-dimensional terrain model of the subregion;
Submodel generation unit, for being based on the three-dimensional terrain model, foundation includes the three-dimensional terrain model The standard database submodel of spatial data, wind direction data and wind speed than data;
Submodel integral unit obtains the target area for integrating each standard database submodel The standard wind field database model in domain.
Preferably, the initial landform modeling unit is specifically used for: being based on global digital elevation model, target area is drawn It is divided into several continuous square subregions, and is vertically stretched the computational domain of the subregion according to preset blocking rate To sustained height threshold value, the corresponding three-dimensional terrain model of the subregion is obtained.
Preferably, the initial landform modeling unit specifically includes:
It is continuous to be divided into several for being based on global digital elevation model by sub-zone dividing subelement for target area Square subregion, wherein the center of the subregion is divided into core space, expansion area and transition region outward;
Terrain information obtains subelement, for determining the core space, institute respectively according to the global digital elevation model Latitude and longitude coordinates information, absolute elevation information, ground mulching type information and the leaf area for stating expansion area and the transition region are close Spend information;
Terrain modeling subelement, for merging the latitude and longitude coordinates of the core space, the expansion area and the transition region Information, absolute elevation information, ground mulching type information and leaf area density information generate dimensional topography curved surface, and according to institute State dimensional topography curved surface and preset blocking rate determine the subregion computational domain height threshold, and by the computational domain edge It vertically is stretched to the height equal with the height threshold, obtains the corresponding three-dimensional terrain model of the subregion.
Preferably, the submodel generation unit specifically includes:
Hydrodynamic model models subelement, for being based on the three-dimensional terrain model, obtains the calculating of the subregion The roughness of ground surface information in domain, and according to fluid dynamics governing equation, establish the roughness of ground surface information and air force Learn the corresponding relationship of resistance coefficient;
Hydrokinetic parameter obtains subelement, for obtaining the calculating parameter of Aerodynamics Model, and according to the meter It calculates the calculating of parameter batch and output model wind speed compares data, wherein the calculating parameter includes: turbulence model type, perimeter strip Part, convergence threshold value;
Submodel generates subelement, and according to the model wind speed than data, foundation includes the three-dimensional terrain model The standard database submodel of spatial data, wind direction data and wind speed than data.
The application fourth aspect provides a kind of wind of standard wind field database model established based on foregoing invention content Field simulator, comprising:
Submodel call unit, the latitude and longitude coordinates for basis to simulated domain, from the standard wind field database mould It is called in type with described to the corresponding standard database submodel of simulated domain;
Simulation of wind arithmetic element, for obtaining the practical wind parameter to simulated domain, by the practical wind parameter It substitutes into the standard database submodel and carries out operation, obtain the wind speed field distribution simulation result data to simulated domain, Wherein, the practical wind parameter includes: the shot and long term wind climate parameter and the reality to simulated domain to simulated domain Border refers to wind speed parameter.
As can be seen from the above technical solutions, the application has the following advantages:
The application is solved by batch and post-processing, establishes the wind speed comprising target area difference wind angle, any position The standard wind field database model of data, based on the standard wind field database model that the application establishes, the reality in combining target area When short-term or long-term general mood marquis's characterisitic parameter, can quickly obtain the wind speed field distribution simulation result of target area.Pass through wind speed The modularization of field model greatly simplifies cumbersome Modeling Calculation step, solves existing complex topography wind field analogy method Analog simulation complex steps, low efficiency promotes difficult technical problem.
Detailed description of the invention
In order to illustrate the technical solutions in the embodiments of the present application or in the prior art more clearly, to embodiment or will show below There is attached drawing needed in technical description to be briefly described, it should be apparent that, the accompanying drawings in the following description is only this Some embodiments of application without any creative labor, may be used also for those of ordinary skill in the art To obtain other attached drawings according to these attached drawings.
Fig. 1 is a kind of process signal of one embodiment of the modeling method of modularization wind-field model provided by the present application Figure;
Fig. 2 is a kind of process signal of second embodiment of the modeling method of modularization wind-field model provided by the present application Figure;
Fig. 3 is a kind of one embodiment of the simulation of wind method based on standard wind field database model provided by the present application Flow diagram;
Fig. 4 is a kind of structural representation of one embodiment of the model building device of modularization wind-field model provided by the present application Figure;
Fig. 5 is a kind of one embodiment of the simulation of wind device based on standard wind field database model provided by the present application Structural schematic diagram;
Fig. 6 is a kind of sub-zone dividing top view of standard wind field database model provided by the present application;
Fig. 7 is a kind of sub-zone dividing side view of standard wind field database model provided by the present application;
Fig. 8 is that the three-dimensional terrain model of the application applies the computational domain schematic diagram obtained after boundary condition.
Specific embodiment
Traditional wind speed field analysis method needs to undergo from initial landform surface data to final wind speed profile result Cumbersome step, and the cross-platform operation between different software can take a substantial amount of time, it is extensive that these factors limit it Application.
The embodiment of the present application provides the modeling method and simulation of wind method of a kind of modularization wind-field model, for solving The analog simulation complex steps of existing complex topography wind field analogy method, low efficiency promote difficult technical problem.
To enable present invention purpose, feature, advantage more obvious and understandable, below in conjunction with the application Attached drawing in embodiment, the technical scheme in the embodiment of the application is clearly and completely described, it is clear that disclosed below Embodiment be only some embodiments of the present application, and not all embodiment.Based on the embodiment in the application, this field Those of ordinary skill's all other embodiment obtained without making creative work belongs to the application protection Range.
Referring to Fig. 1, the embodiment of the present application provides a kind of modeling method of modularization wind-field model, comprising:
101: being based on global digital elevation model, target area is divided into several continuous square subregions, is built Found the corresponding three-dimensional terrain model of the subregion;
It should be noted that needing first to obtain the landform of target area before the simulation of wind method for implementing the application Information establishes the three-dimensional terrain model of each sub-regions.
102: be based on three-dimensional terrain model, foundation include the spatial data of three-dimensional terrain model, wind direction data and Standard database submodel of the wind speed than data;
It should be noted that establish include subregion terrain information three-dimensional terrain model after, in dimensional topography On the basis of model, by the spatial data, wind direction data and wind speed of each sub-regions than data be inserted in it is corresponding dimensionally Shape model obtains the standard database submodel of the subregion.
103: judging whether the quantity of standard database submodel is equal with the quantity of three-dimensional terrain model, if then Step 104 is executed, if it is not, then return step 102.
It should be noted that circulation executes step 102, until all subregions of target area establish corresponding mark Step 104 is executed again after quasi- database subsystem module.
104: each standard database submodel being integrated, the standard wind field data of the target area are obtained Library model.
The embodiment of the present application is solved by batch and post-processing, and establishing includes target area difference wind angle, any position Air speed data blocking criteria wind field database model, based on the application establish standard wind field database model, in conjunction with Real-time short-term or long-term general mood marquis's characterisitic parameter of objective area, can quickly obtain the wind speed field distribution simulation knot of target area Fruit.Cumbersome Modeling Calculation step is greatly simplified by the modularization of wind speed field model, solves existing complicated landform The analog simulation complex steps of simulation of wind method, low efficiency promote difficult technical problem.
A kind of above detailed description of one embodiment of the modeling method of modularization wind-field model provided by the present application, It is below a kind of detailed description of second embodiment of the modeling method of modularization wind-field model provided by the present application.
Referring to Fig. 2, the embodiment of the present application provides a kind of modeling method of modularization wind-field model, comprising:
201: it is based on global digital elevation model, target area is divided into several continuous square subregions, In, the center of the subregion is divided into core space, expansion area and transition region outward;
It should be noted that please referring to Fig. 6 and Fig. 7, subregion is carried out according to the square of preset side length, which is to extract Wind field information and the core space for establishing database.The outer layer of core space is provided that expansion area, and setting expansion area is to reduce Influence of the boundary to core space wind field result, wherein the terrain information in expansion area is real terrain elevation, also from Digital Radio It is obtained in elevation model.The outer layer of expansion area is transition region, and setting transition region is in order to avoid the elevation of target area boundaries is prominent Become, guarantee arrives stream wind speed section simultaneously avoid actual complex boundary generate reflux, wherein transition region elevation by expansion area boundary gradually It is reduced to the lowest elevation of core space.
In addition, expansion area and transition region are the square loop with core space common point, and core space, expansion area and mistake The size for crossing area can be adjusted by computing capability, and the optimal case for Practical Project considers, the core space of the present embodiment Side length is preferably 10 kilometers, and the width of expansion area is preferably 1 kilometer, and the width of transition region is preferably 5 kilometers.
202: according to the global digital elevation model, determining the core space, the expansion area and the transition respectively Latitude and longitude coordinates information, absolute elevation information, ground mulching type information and the leaf area density information in area;
It should be noted that then obtaining mesh from global digital elevation model after target area has been divided subregion Each sub-regions marked under region include the corresponding latitude and longitude coordinates information of core space, expansion area and transition region of subregion, sea Pull out elevation information, ground mulching type information and leaf area density information;
203: merging latitude and longitude coordinates information, the absolute elevation letter of the core space, the expansion area and the transition region Breath, ground mulching type information and leaf area density information, generate dimensional topography curved surface, and according to the dimensional topography curved surface and Preset blocking rate determines the height threshold of the computational domain of the subregion, and by the computational domain be vertically stretched to it is described The equal height of height threshold obtains the corresponding three-dimensional terrain model of the subregion.
It should be noted that including the latitude and longitude coordinates and its correspondence of target area discrete point according to the landform altitude Absolute elevation, the roughness of ground surface data include ground mulching type, vegetation pattern and leaf area density information.Merge Core, the latitude and longitude coordinates of extension and transitional region, absolute elevation and roughness of ground surface data generate dimensional topography curved surface.
Then according to the dimensional topography curved surface and preset blocking rate, topographic projection's area under each wind direction is calculated, with And meeting the maximum height of blocking rate requirement is computational domain height threshold, and the computational domain is vertically stretched to and the height The equal height of threshold value is spent, the corresponding three-dimensional terrain model of the subregion is obtained.
Wherein, the height of computational domain takes the computational domain height of the maximum wind direction of the downward blocking rate of full blast.Blocking rate=stream Object maximum front face area/basin cross-sectional area, the maximum front face area for streaming object take target area and extend out area windward Projected area on face, the cross-sectional area in basin are projection of the computational domain in windward side.In addition, to avoid the obstruction in computational domain Effect, according to wind tunnel test experience, the preset blocking rate of the present embodiment can be preferably 3~5%.
204: being based on the three-dimensional terrain model, obtain the roughness of ground surface information of the computational domain of the subregion, and root According to fluid dynamics governing equation, the corresponding relationship of the roughness of ground surface information and aerodynamic resistance factor is established;
It should be noted that in the corresponding relationship for establishing the roughness of ground surface information and aerodynamic resistance factor When, the horizontal grid spacing parameter, vertical grid first floor height parameter and grid for obtaining the three-dimensional terrain model of calculating first increase Long scale parameter.
Meanwhile premised on the independence for guaranteeing numerical solution.The present embodiment preferred levels size of mesh opening parameter is 50 meters, is erected It is 1 meter to first floor height parameter, grid growth ratio parameter is preferably smaller than 1.1.
Then by applying source item in the dynamic governing equations of wind field near the ground, the air force that coarse earth's surface is provided It learns resistance and introduces flow field calculation model;Resistance source item indicates in computational fluid dynamics governing equation are as follows:
su=-CdA|U|u
sk=CdA(βp|U|3d|U|k)
Wherein, su、sk、sεIt is illustrated respectively in the sources of resistance applied in the equation of momentum, tubulence energy equation and dissipative shock wave equation , CdIt indicates that the corresponding resistance coefficient of ground surface type, A are leaf area density coefficient, indicates leaf area contained in unit volume Size, U is mean wind speed, and u is instantaneous wind speed, and k is tubulence energy, and ε is dissipative shock wave, αp、βp、αd、βdFor the warp in turbulence model Test constant.
And the influence of surface vegetation stream field is related with the type of vegetation, closeness.In the present embodiment, resistance source item passes through The leaf area density coefficient of surface vegetation and its corresponding resistance coefficient determine.However the vegetation pattern of the same area and growth shape Condition does not ensure that uniformity, therefore the corresponding roughness of ground surface data information of the present embodiment combination discrete point, by analog domain with Kilometer is unit subregion, and the average leaf area density coefficient in each area is respectively adopted and resistance coefficient calculates.
205: obtaining the calculating parameter of Aerodynamics Model, and simultaneously output model is calculated according to the calculating parameter batch Wind speed compares data, wherein the calculating parameter includes: turbulence model type, boundary condition, convergence threshold value;
It should be noted that obtaining the calculating parameter of Aerodynamics Model, and according to the calculating parameter using parallel Calculation method batch solves wind field movement and output model wind speed compares data.
Specifically, turbulence model type preferably uses the k- ω model of sst, pressure is coupled with speed to be calculated using SIMPLIC Method, convective term use QUICK format, and diffusion term uses central difference schemes.Convergence is that the residual error of all variables restrains, when The residual error of x, y, z momentum sum reaches 10-5, the residual error of ω and continuity equation reaches 10-3When, all residual errors all reach convergence.
The boundary condition is speed entrance, exit condition.Referring to Fig. 8, Fig. 8 is the present embodiment in specified registering weather vane When calculation, used boundary condition form.The form is not needed to regenerate computational domain, be improved when carrying out batch calculating Computational efficiency.
Wherein, wind speed calculating than data of the model wind speed than the single discrete point in data specifically includes: with wind speed ratio Characterize the influence to wind speed profile of landform altitude and roughness of ground surface, under θ wind angle, the wind speed of target area midpoint (x, y, z) It may be expressed as: than S
S (x, y, z, θ)=U (x, y, z, θ)/Uin(z,θ)
Wherein, U (x, y, z, θ) is the wind speed solving result of target area point (x, y, z), Uin(z, θ) is border entries z high Arrives stream wind speed size at degree.
Finally, target area wind speed is converted into latitude coordinate than exporting the cartesian coordinate of result, the sub-district is obtained The wind speed of discrete point in domain compares data.
206: according to the model wind speed than data, foundation include the three-dimensional terrain model spatial data, The standard database submodel of wind direction data and wind speed than data.
It should be noted that obtaining and counting the wind speed of the discrete point being calculated by step 205 than data, and will system The wind speed for the various discrete point counted in conjunction with the three-dimensional terrain model that step 203 is established, obtains the dimensional topography mould than data The corresponding standard database submodel of type.
207: judging whether the quantity of standard database submodel is equal with the quantity of three-dimensional terrain model, if then Step 208 is executed, if it is not, then return step 204
It should be noted that circulation executes step 204 to step 206, until all subregions of target area are established Step 208 is executed again after corresponding standard database submodule.
208: each standard database submodel being integrated, the standard wind field data of the target area are obtained Library model.
In addition, the present embodiment is also based on global digital elevation model frame, the sub- entry in expanding data library, including it is multiple Spatial altitude, multiple wind directions and corresponding wind speed ratio.
The embodiment of the present application is solved by batch and post-processing, and establishing includes target area difference wind angle, any position Air speed data blocking criteria wind field database model, based on the application establish standard wind field database model, in conjunction with Real-time short-term or long-term general mood marquis's characterisitic parameter of objective area, can quickly obtain the wind speed field distribution simulation knot of target area Fruit.Cumbersome Modeling Calculation step is greatly simplified by the modularization of wind speed field model, solves existing complicated landform The analog simulation complex steps of simulation of wind method, low efficiency promote difficult technical problem.
A kind of above detailed description of second embodiment of the modeling method of modularization wind-field model provided by the present application, It is below a kind of the detailed of one embodiment of the simulation of wind method based on standard wind field database model provided by the present application Explanation.
Referring to Fig. 3, the embodiment of the present application provides a kind of standard wind field database mould established based on the above embodiment The simulation of wind method of type, comprising:
301: according to the latitude and longitude coordinates to simulated domain, called from the standard wind field database model with it is described to The corresponding standard database submodel of simulated domain;
302: obtaining the practical wind parameter to simulated domain, the practical wind parameter is substituted into the standard database Submodel carries out operation, obtains the simulation of wind result data to simulated domain.
Wherein, the practical wind parameter include: the shot and long term wind climate parameter to simulated domain with described wait simulate The practical of region refers to wind speed parameter.
It should be noted that present embodiments providing a kind of standard wind field database model established based on previous embodiment Simulation of wind method, first according to coordinate position to simulated domain, called from the standard wind field database model with It is described to the corresponding standard database submodel of simulated domain, then using the standard database submodel selected as template, will The practical wind parameter is substituted into the standard database submodel and carries out operation by the practical wind parameter to simulated domain, Obtain the simulation of wind result data to simulated domain.
In addition, the quantity for the standard database submodel that the present embodiment selects includes but is not limited to 1.
The above are a kind of implementations of the simulation of wind method based on standard wind field database model provided by the present application The detailed description of example, here is a kind of the detailed of one embodiment of the model building device of modularization wind-field model provided by the present application Explanation.
Referring to Fig. 4, the embodiment of the present application provides a kind of model building device of modularization wind-field model, comprising:
Initial landform modeling unit 401, for being based on global digital elevation model, is divided into several companies for target area Continuous square subregion establishes the corresponding three-dimensional terrain model of the subregion;
Submodel generation unit 402, for being based on the three-dimensional terrain model, foundation includes the three-dimensional terrain model Spatial data, wind direction data and wind speed than data standard database submodel;
Submodel integral unit 403 obtains the target for integrating each standard database submodel The standard wind field database model in region.
Preferably, the initial landform modeling unit 401 is specifically used for: global digital elevation model is based on, by target area Domain is divided into several continuous square subregions, and according to preset blocking rate by the computational domain of the subregion vertically It is stretched to sustained height threshold value, obtains the corresponding three-dimensional terrain model of the subregion.
Preferably, the initial landform modeling unit 401 specifically includes:
Target area is divided into several for being based on global digital elevation model by sub-zone dividing subelement 4011 Continuous square subregion, wherein the center of the subregion is divided into core space, expansion area and transition region outward;
Terrain information obtains subelement 4012, for determining the core respectively according to the global digital elevation model Latitude and longitude coordinates information, absolute elevation information, ground mulching type information and the blade face in area, the expansion area and the transition region Product density information;
Terrain modeling subelement 4013, for merging the longitude and latitude of the core space, the expansion area and the transition region Coordinate information, absolute elevation information, ground mulching type information and leaf area density information generate dimensional topography curved surface, and root Determine the height threshold of the computational domain of the subregion according to the dimensional topography curved surface and preset blocking rate, and by the calculating Domain is vertically stretched to the height equal with the height threshold, obtains the corresponding three-dimensional terrain model of the subregion.
Preferably, the submodel generation unit 402 specifically includes:
Hydrodynamic model models subelement 4021, for being based on the three-dimensional terrain model, obtains the subregion The roughness of ground surface information of computational domain, and according to fluid dynamics governing equation, establish the roughness of ground surface information and air The corresponding relationship of dynamical resistance coefficient;
Hydrokinetic parameter obtains subelement 4022, for obtaining the calculating parameter of Aerodynamics Model, and according to institute It states the calculating of calculating parameter batch and output model wind speed compares data, wherein the calculating parameter includes: turbulence model type, side Boundary's condition, convergence threshold value;
Submodel generates subelement 4023, and according to the model wind speed than data, foundation includes the dimensional topography mould Standard database submodel of spatial data, wind direction data and the wind speed of type than data.
The above are a kind of detailed description of one embodiment of the model building device of modularization wind-field model provided by the present application, It is below a kind of the detailed of one embodiment of the simulation of wind method based on standard wind field database model provided by the present application Explanation.
Referring to Fig. 5, the embodiment of the present application provides a kind of standard wind field database established based on foregoing invention content The simulation of wind device of model, comprising:
Submodel call unit 501, the latitude and longitude coordinates for basis to simulated domain, from the standard wind field database It is called in model with described to the corresponding standard database submodel of simulated domain;
Simulation of wind arithmetic element 502 joins the practical wind for obtaining the practical wind parameter to simulated domain Number substitutes into the standard database submodel and carries out operation, obtains the wind speed field distribution simulation number of results to simulated domain It according to, wherein the practical wind parameter include: the shot and long term wind climate parameter to simulated domain and described to simulated domain It is practical to refer to wind speed parameter.
It is apparent to those skilled in the art that for convenience and simplicity of description, the system of foregoing description, The specific work process of device and unit, can refer to corresponding processes in the foregoing method embodiment, and details are not described herein.
In several embodiments provided herein, it should be understood that disclosed system, device and method can be with It realizes by another way.For example, the apparatus embodiments described above are merely exemplary, for example, the unit It divides, only a kind of logical function partition, there may be another division manner in actual implementation, such as multiple units or components It can be combined or can be integrated into another system, or some features can be ignored or not executed.Another point, it is shown or The mutual coupling, direct-coupling or communication connection discussed can be through some interfaces, the indirect coupling of device or unit It closes or communicates to connect, can be electrical property, mechanical or other forms.
The description of the present application and term " first " in above-mentioned attached drawing, " second ", " third ", " the 4th " etc. are (if deposited ) it is to be used to distinguish similar objects, without being used to describe a particular order or precedence order.It should be understood that use in this way Data are interchangeable under appropriate circumstances, so that embodiments herein described herein for example can be in addition to illustrating herein Or the sequence other than those of description is implemented.In addition, term " includes " and " having " and their any deformation, it is intended that Cover it is non-exclusive include, for example, containing the process, method, system, product or equipment of a series of steps or units need not limit In step or unit those of is clearly listed, but may include be not clearly listed or for these process, methods, produce The other step or units of product or equipment inherently.
The unit as illustrated by the separation member may or may not be physically separated, aobvious as unit The component shown may or may not be physical unit, it can and it is in one place, or may be distributed over multiple In network unit.It can select some or all of unit therein according to the actual needs to realize the mesh of this embodiment scheme 's.
It, can also be in addition, the functional units in various embodiments of the present invention may be integrated into one processing unit It is that each unit physically exists alone, can also be integrated in one unit with two or more units.Above-mentioned integrated list Member both can take the form of hardware realization, can also realize in the form of software functional units.
The above, above embodiments are only to illustrate the technical solution of the application, rather than its limitations;Although referring to before Embodiment is stated the application is described in detail, those skilled in the art should understand that: it still can be to preceding Technical solution documented by each embodiment is stated to modify or equivalent replacement of some of the technical features;And these It modifies or replaces, the spirit and scope of each embodiment technical solution of the application that it does not separate the essence of the corresponding technical solution.

Claims (10)

1. a kind of modeling method of modularization wind-field model characterized by comprising
S1: being based on global digital elevation model, target area is divided into several continuous square subregions, described in foundation The corresponding three-dimensional terrain model of subregion;
S2: it is based on the three-dimensional terrain model, foundation includes the spatial data of the three-dimensional terrain model, wind direction data And wind speed is than the standard database submodel of data;
S3: circulation executes step S2, establishes the corresponding standard database submodule of each three-dimensional terrain model respectively Type;
S4: each standard database submodel is integrated, and obtains the standard wind field database mould of the target area Type.
2. modeling method according to claim 1, which is characterized in that the step S1 is specifically included:
S11: being based on global digital elevation model, and target area is divided into several continuous square subregions, and according to The computational domain of the subregion is vertically stretched to sustained height threshold value by preset blocking rate, and it is corresponding to obtain the subregion Three-dimensional terrain model.
3. modeling method according to claim 2, which is characterized in that the step S11 is specifically included:
S111: being based on global digital elevation model, and target area is divided into several continuous square subregions, wherein The center of the subregion is divided into core space, expansion area and transition region outward;
S112: according to the global digital elevation model, the core space, the expansion area and the transition region are determined respectively Latitude and longitude coordinates information, absolute elevation information, ground mulching type information and leaf area density information;
S113: merge the latitude and longitude coordinates information of the core space, the expansion area and the transition region, absolute elevation information, Ground mulching type information and leaf area density information generate dimensional topography curved surface, and according to the dimensional topography curved surface and in advance The blocking rate set determines the height threshold of the computational domain of the subregion, and the computational domain is vertically stretched to and the height Spend the equal height of threshold value.
4. modeling method according to claim 1, which is characterized in that the step S2 is specifically included:
S21: it is based on the three-dimensional terrain model, obtains the roughness of ground surface information of the computational domain of the subregion, and according to stream Body dynamics governing equation establishes the corresponding relationship of the roughness of ground surface information and aerodynamic resistance factor;
S22: the calculating parameter of Aerodynamics Model is obtained, and simultaneously output model wind speed is calculated according to the calculating parameter batch Compare data, wherein the calculating parameter includes: turbulence model type, boundary condition, convergence threshold value;
S23: according to the model wind speed than data, foundation includes the spatial data of the three-dimensional terrain model, wind direction The standard database submodel of data and wind speed than data.
5. a kind of simulation of wind method for the standard wind field database model established based on such as Claims 1-4 any one, It is characterized in that, comprising:
S6: it according to the latitude and longitude coordinates to simulated domain, calls with described from the standard wind field database model wait simulate The corresponding standard database submodel in region;
S7: obtaining the practical wind parameter to simulated domain, and the practical wind parameter is substituted into the standard database submodule Type carries out operation, obtains the simulation of wind result data to simulated domain, wherein the practical wind parameter includes: described Shot and long term wind climate parameter to simulated domain refers to wind speed parameter to the practical of simulated domain with described.
6. a kind of model building device of modularization wind-field model characterized by comprising
Target area is being divided into several continuously just for being based on global digital elevation model by initial landform modeling unit The corresponding three-dimensional terrain model of the subregion is established in square subregions domain;
Submodel generation unit, for being based on the three-dimensional terrain model, foundation includes the space of the three-dimensional terrain model The standard database submodel of coordinate data, wind direction data and wind speed than data;
Submodel integral unit obtains the target area for integrating each standard database submodel Standard wind field database model.
7. model building device according to claim 6, which is characterized in that the initial landform modeling unit is specifically used for: base In global digital elevation model, target area is divided into several continuous square subregions, and according to preset obstruction The computational domain of the subregion is vertically stretched to sustained height threshold value by rate, obtains the corresponding dimensional topography mould of the subregion Type.
8. model building device according to claim 7, which is characterized in that the initial landform modeling unit specifically includes:
Target area is being divided into several continuously just for being based on global digital elevation model by sub-zone dividing subelement Square subregions domain, wherein the center of the subregion is divided into core space, expansion area and transition region outward;
Terrain information obtains subelement, for determining the core space, the expansion respectively according to the global digital elevation model The latitude and longitude coordinates information of exhibition section and the transition region, absolute elevation information, ground mulching type information and leaf area density letter Breath;
Terrain modeling subelement, for merge the core space, the expansion area and the transition region latitude and longitude coordinates information, Absolute elevation information, ground mulching type information and leaf area density information generate dimensional topography curved surface, and according to the three-dimensional Terrain surface and preset blocking rate determine the height threshold of the computational domain of the subregion, and the computational domain is vertically drawn The height equal with the height threshold is extended to, the corresponding three-dimensional terrain model of the subregion is obtained.
9. model building device according to claim 6, which is characterized in that the submodel generation unit specifically includes:
Hydrodynamic model models subelement, for being based on the three-dimensional terrain model, obtains the computational domain of the subregion Roughness of ground surface information, and according to fluid dynamics governing equation, it establishes the roughness of ground surface information and aerodynamics hinders The corresponding relationship of force coefficient;
Hydrokinetic parameter obtains subelement, joins for obtaining the calculating parameter of Aerodynamics Model, and according to the calculating Number batch calculates and output model wind speed compares data, wherein the calculating parameter includes: turbulence model type, boundary condition, receipts Hold back level threshold value;
Submodel generates subelement, and according to the model wind speed than data, foundation includes the space of the three-dimensional terrain model The standard database submodel of coordinate data, wind direction data and wind speed than data.
10. a kind of simulation of wind device for the standard wind field database model established based on such as claim 6 to 9 any one, It is characterised by comprising:
Submodel call unit, the latitude and longitude coordinates for basis to simulated domain, from the standard wind field database model It calls with described to the corresponding standard database submodel of simulated domain;
Simulation of wind arithmetic element substitutes into the practical wind parameter for obtaining the practical wind parameter to simulated domain The standard database submodel carries out operation, obtains the simulation of wind result data to simulated domain, wherein the reality Border wind parameter includes: that the shot and long term wind climate parameter to simulated domain and the practical reference wind speed to simulated domain are joined Number.
CN201811203630.3A 2018-10-16 2018-10-16 A kind of modeling method and simulation of wind method of modularization wind-field model Pending CN109165476A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811203630.3A CN109165476A (en) 2018-10-16 2018-10-16 A kind of modeling method and simulation of wind method of modularization wind-field model

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811203630.3A CN109165476A (en) 2018-10-16 2018-10-16 A kind of modeling method and simulation of wind method of modularization wind-field model

Publications (1)

Publication Number Publication Date
CN109165476A true CN109165476A (en) 2019-01-08

Family

ID=64878296

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811203630.3A Pending CN109165476A (en) 2018-10-16 2018-10-16 A kind of modeling method and simulation of wind method of modularization wind-field model

Country Status (1)

Country Link
CN (1) CN109165476A (en)

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110111420A (en) * 2019-04-12 2019-08-09 东北林业大学 Skifield wind field modeling method
CN111143999A (en) * 2019-12-26 2020-05-12 广东电网有限责任公司电力科学研究院 Method, device and equipment for calculating regional surface roughness
CN112364502A (en) * 2020-11-09 2021-02-12 广东海洋大学寸金学院 Method and device for generating area planning scheme for relieving heat island effect and storage medium
CN113065184A (en) * 2021-03-09 2021-07-02 中国三峡建设管理有限公司 Method and device for predicting wind field of arch dam construction area
CN115577297A (en) * 2022-10-21 2023-01-06 武汉理工大学 Visual sea area wind field simulation and prediction system and method
CN115828684A (en) * 2022-11-28 2023-03-21 中交建筑集团有限公司 Wind resistance analysis method for main girder large cantilever of asymmetric cable-stayed bridge
CN116452739A (en) * 2023-04-18 2023-07-18 中广核风电有限公司 Method and device for realizing three-dimensional visualization of wind field model

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104699936A (en) * 2014-08-18 2015-06-10 沈阳工业大学 Sector management method based on CFD short-term wind speed forecasting wind power plant
US20150302117A1 (en) * 2012-11-09 2015-10-22 Korea Institute Of Energy Research Numerical simulation system and numerical simulation method for atmospheric flow by computational fluid dynamics
CN107390298A (en) * 2017-07-19 2017-11-24 云南电网有限责任公司电力科学研究院 A kind of analogy method and device of Complex Mountain underlying surface strong wind
CN108152867A (en) * 2017-12-29 2018-06-12 南方电网科学研究院有限责任公司 Typhoon early warning system and method for power transmission line
CN108197376A (en) * 2017-12-28 2018-06-22 中铁二局集团有限公司 A kind of temporary building design and modeling method and system based on BIM
CN108563867A (en) * 2018-04-13 2018-09-21 华中科技大学 A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20150302117A1 (en) * 2012-11-09 2015-10-22 Korea Institute Of Energy Research Numerical simulation system and numerical simulation method for atmospheric flow by computational fluid dynamics
CN104699936A (en) * 2014-08-18 2015-06-10 沈阳工业大学 Sector management method based on CFD short-term wind speed forecasting wind power plant
CN107390298A (en) * 2017-07-19 2017-11-24 云南电网有限责任公司电力科学研究院 A kind of analogy method and device of Complex Mountain underlying surface strong wind
CN108197376A (en) * 2017-12-28 2018-06-22 中铁二局集团有限公司 A kind of temporary building design and modeling method and system based on BIM
CN108152867A (en) * 2017-12-29 2018-06-12 南方电网科学研究院有限责任公司 Typhoon early warning system and method for power transmission line
CN108563867A (en) * 2018-04-13 2018-09-21 华中科技大学 A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM

Cited By (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110111420A (en) * 2019-04-12 2019-08-09 东北林业大学 Skifield wind field modeling method
CN110111420B (en) * 2019-04-12 2023-08-01 东北林业大学 Ski field wind field modeling method
CN111143999A (en) * 2019-12-26 2020-05-12 广东电网有限责任公司电力科学研究院 Method, device and equipment for calculating regional surface roughness
CN112364502A (en) * 2020-11-09 2021-02-12 广东海洋大学寸金学院 Method and device for generating area planning scheme for relieving heat island effect and storage medium
CN112364502B (en) * 2020-11-09 2024-04-16 广东海洋大学寸金学院 Region planning scheme generation method, device and storage medium for relieving heat island effect
CN113065184A (en) * 2021-03-09 2021-07-02 中国三峡建设管理有限公司 Method and device for predicting wind field of arch dam construction area
CN113065184B (en) * 2021-03-09 2024-03-12 中国三峡建设管理有限公司 Method and device for predicting wind field of arch dam construction area
CN115577297A (en) * 2022-10-21 2023-01-06 武汉理工大学 Visual sea area wind field simulation and prediction system and method
CN115828684A (en) * 2022-11-28 2023-03-21 中交建筑集团有限公司 Wind resistance analysis method for main girder large cantilever of asymmetric cable-stayed bridge
CN116452739A (en) * 2023-04-18 2023-07-18 中广核风电有限公司 Method and device for realizing three-dimensional visualization of wind field model
CN116452739B (en) * 2023-04-18 2023-12-15 中广核风电有限公司 Method and device for realizing three-dimensional visualization of wind field model

Similar Documents

Publication Publication Date Title
CN109165476A (en) A kind of modeling method and simulation of wind method of modularization wind-field model
CN109446691B (en) Living standing tree wind resistance performance analysis method based on laser point cloud and aerodynamics
CN108536881B (en) Method and device for calculating the power generation of a wind farm
CN108563867A (en) A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM
CN105513133B (en) A kind of city wind environment numerical map makes and display methods
CN111767597B (en) City model verification method, device, equipment and storage medium
CN114997541B (en) Urban waterlogging prediction method and early warning platform based on digital twin technology
CN102663827A (en) Three-dimensional dynamic whole-process simulation method for storm surge and flood routing in complex flooding areas
CN108664705B (en) OpenFOAM-based method for simulating surface roughness of complex terrain
CN112163381B (en) Lateral boundary condition setting method suitable for complex terrain wind field flow numerical simulation
CN109215112A (en) A kind of mask method of unilateral side point cloud model
CN109948214A (en) The multiple dimensioned wind environment method for numerical simulation in city
CN106874602A (en) Process meteorological data method and apparatus
CN108072897A (en) It is a kind of to mix computational methods when two-dimensionally seismic wave is walked
CN109635317A (en) A kind of high altitude localities CFD artificially generated terrain topological method
KR102537488B1 (en) Wind speed analyzing system in urban area based on inflow wind speed and morphological parameters, and wind speed analyzing method using the same
CN107274448B (en) Variable weight cost aggregation stereo matching algorithm based on horizontal tree structure
CN109190187A (en) A kind of wind power plant wake flow calculation method and system based on more physical models
CN112862969A (en) Broken line contributor in civil engineering
CN114937137B (en) BIM and GIS-based building environment intelligent analysis method
CN104517299A (en) Method for restoring and resimulating physical video fluid driving model
CN108534982A (en) A kind of easement curve for mountain area simulation of wind
CN111898819B (en) Space grid dividing method and device
CN112100839B (en) Method and system for constructing high-fidelity single storm analytic wind speed field in cylindrical coordinate system
CN104299265A (en) Group behavior control technology under fluid environment influence

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
RJ01 Rejection of invention patent application after publication
RJ01 Rejection of invention patent application after publication

Application publication date: 20190108