CN109583096A - A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling - Google Patents

A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling Download PDF

Info

Publication number
CN109583096A
CN109583096A CN201811467690.6A CN201811467690A CN109583096A CN 109583096 A CN109583096 A CN 109583096A CN 201811467690 A CN201811467690 A CN 201811467690A CN 109583096 A CN109583096 A CN 109583096A
Authority
CN
China
Prior art keywords
mesoscale
model
minute yardstick
wind
wrf
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
CN201811467690.6A
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.)
Huarun Power Technology Research Institute Co Ltd
Original Assignee
Huarun Power Technology Research Institute 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 Huarun Power Technology Research Institute Co Ltd filed Critical Huarun Power Technology Research Institute Co Ltd
Priority to CN201811467690.6A priority Critical patent/CN109583096A/en
Publication of CN109583096A publication Critical patent/CN109583096A/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
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate 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

The wind-resources calculation method based on mesoscale model and minute yardstick models coupling that the present invention provides a kind of, comprising: obtain the anemometer tower measured data of target area, and obtain data of analyzing again corresponding with the target area as ambient field;Simulation calculating is carried out to ambient field by mesoscale model, obtains WRF Mesoscale Meteorology data;Simulation calculating is carried out to ambient field by minute yardstick model, obtains OpenFOAM minute yardstick calculated result;Statistical relationship is established according to anemometer tower measured data and WRF Mesoscale Meteorology data;OpenFOAM minute yardstick calculated result is corrected using statistical relationship, obtains wind-resources amendment meter result.Present invention can be implemented in mesoclimate elements and the area of regional climate element correlation difference to carry out operation and prediction, improves computational efficiency, reduces calculation amount and substantially increases the accuracy to climatic simulation.

Description

A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling
Technical field
The present invention relates to wind-resources assessment technical fields, are based on mesoscale model and micro- ruler more specifically to one kind Spend the wind-resources calculation method of models coupling.
Background technique
WRF (research of Weather Research and Forecast weather and prediction) is the gas being widely used As model.WRF mode has numerical method the most advanced and data assimilation now, using improved physical process Scheme, while there is multinest and be easily positioned in the ability of diverse geographic location, emphasis considers several kilometers to tens kilometers The horizontal grid of resolution ratio such as improves from cloud scale to synoptic scale at the simulation and forecast essence of different scales significant weather feature Degree is well adapted for the needs from idealization research to different applications such as operational forecasts, and oneself is widely used in wind-resources mould It fits in the business of forecasting wind speed.Since WRF belongs to mesoscale model, emphasis simulation is from several kilometers to hundreds of kilometer scale Weather phenomenon, the general horizontal grid for only considering 1-10 kilometers of resolution ratio, and be unable to reach the small scale of tens meters of resolution ratio Fining simulation.Especially under complicated landform, landform is difficult to be embodied in WRF for the forced change effect of wind speed. Therefore the WRF low resolution air speed data exported is converted high-resolution fining number by the method for generally using NO emissions reduction According to.
NO emissions reduction method is broadly divided into statistics NO emissions reduction and power NO emissions reduction at present.Wherein, Statistical downscaling is to utilize The observational data of many years establishes the statistical relationship between mesoclimate situation and regional climate element, and is provided with independent observation Material examines this relationship, this relationship is finally applied to the mesoclimate information of global climate model output again, to predict The climatic change trend of regional feature;And power NO emissions reduction is exactly that mesoscale model drives small scale model more to be refined Wind field, as Mesoscale Numerical Simulation result provides whole region really entry condition of the meteorological background as CFD model;And CFD mould Type refines Mesoscale Numerical Simulation result.
Currently, the shortcomings that Statistical downscaling is that enough observational datas is needed to establish statistical model, it is impossible to be used in middle ruler Spend the area of climatic elements and regional climate element correlation difference;And the shortcomings that traditional power NO emissions reduction method is calculation amount Greatly, Fei Jishi, by the Boundary Condition Effect that global climate model provides, dynamic mode compares the systematic error of climatic simulation Greatly.
Summary of the invention
In view of this, the present invention provides a kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling So as to solve the deficiencies in the prior art.
To solve the above problems, the present invention provide it is a kind of based on the wind-resources of mesoscale model and minute yardstick models coupling Calculation method, the mesoscale model are WRF mesoscale model, and the minute yardstick model is OpenFOAM minute yardstick model, comprising:
The anemometer tower measured data of target area is obtained, and obtains data of analyzing again corresponding with the target area and makees For ambient field;
Simulation calculating is carried out to the ambient field by the mesoscale model, obtains WRF Mesoscale Meteorology data; Also, simulation calculating is carried out to the ambient field by the minute yardstick model, obtains OpenFOAM minute yardstick calculated result;
Statistical relationship is established according to the anemometer tower measured data and the WRF Mesoscale Meteorology data;
The OpenFOAM minute yardstick calculated result is corrected using the statistical relationship, obtains wind-resources corrected Calculation knot Fruit.
Preferably, described " simulation calculating to be carried out to the ambient field by the mesoscale model, obtains WRF mesoscale Meteorological Models data " include:
According to the mesoscale model, the kilometer grade grid in the target area of WRF mesoscale model is constructed;
Based on the ambient field, the wind speed of the kilometer grade grid is calculated for the target area, obtains WRF mesoscale Meteorological Models data.
Preferably, described " according to the mesoscale model, to construct the public affairs in the target area of WRF mesoscale model In grade grid " include:
The target area is split as unit of pre-determined distance, according to mesoscale model building based on segmentation Three layers of nested grid of the corresponding WRF mesoscale model in the target area afterwards, as kilometer grade grid.
Preferably, the grid distance of three layers of nested grid is respectively 27km, 9km and 3km.
Preferably, described " it is based on the ambient field, the wind speed of the kilometer grade grid is calculated for the target area, Obtain WRF Mesoscale Meteorology data " include:
Acquisition is combined with the matched Parameterization Scheme in the target area;
It is combined according to the Parameterization Scheme, determines simulated time;
It is directed to the numerical model simulation of the kilometer grade grid of the target area, is obtained in the simulated time The corresponding wind speed of each kilometer grade grid, as WRF Mesoscale Meteorology data.
Preferably, described " simulation calculating to be carried out to the ambient field by the minute yardstick model, it is micro- to obtain OpenFOAM Dimension calculation result " includes:
According to the minute yardstick model, the fining net in the target area of OpenFOAM minute yardstick mode is constructed Lattice;
Obtain the boundary condition corresponding with the target area and initial value under OpenFOAM minute yardstick mode;
According to the boundary condition and the initial value, solved by the OpenFOAM of the OpenFOAM minute yardstick mode Device calculates the flow field of the fining grid of the target area, obtains the OpenFOAM minute yardstick calculated result.
Preferably, described " statistics to be established according to the anemometer tower measured data and the WRF Mesoscale Meteorology data Relationship " includes:
The WRF Mesoscale Meteorology data are corrected by the anemometer tower measured data, obtain mesoscale correction phase Relationship number;Also, related coefficient is corrected according to the mesoscale and establishes the statistical relationship.
In addition, to solve the above problems, the present invention also provides a kind of based on mesoscale model and minute yardstick models coupling Wind-resources computing device, the mesoscale model are WRF mesoscale model, and the minute yardstick model is OpenFOAM minute yardstick mould Type, comprising:
Module is obtained, for obtaining the anemometer tower measured data of target area, and is obtained corresponding with the target area Analyze data again as ambient field;
Computing module obtains WRF mesoscale for carrying out simulation calculating to the ambient field by the mesoscale model Meteorological Models data;Also, simulation calculating is carried out to the ambient field by the minute yardstick model, obtains the micro- ruler of OpenFOAM Spend calculated result;
Module is established, for establishing system according to the anemometer tower measured data and the WRF Mesoscale Meteorology data Meter relationship;
Correction module obtains wind money for correcting the OpenFOAM minute yardstick calculated result using the statistical relationship Source corrected Calculation result.
In addition, to solve the above problems, the computer equipment includes storage the present invention also provides a kind of computer equipment Device and processor, the memory calculate journey based on the wind-resources of mesoscale model and minute yardstick models coupling for storing Sequence, processor operation it is described by the wind-resources calculation procedure of mesoscale model and minute yardstick models coupling so that based on described It calculates machine equipment and executes the wind-resources calculation method as described above based on mesoscale model and minute yardstick models coupling.
In addition, to solve the above problems, the present invention also provides a kind of computer readable storage medium, it is described computer-readable The wind-resources calculation procedure based on mesoscale model and minute yardstick models coupling is stored on storage medium, it is described to be based on mesoscale The wind-resources calculation procedure of model and minute yardstick models coupling is realized as described above based on mesoscale mould when being executed by processor The wind-resources calculation method of type and minute yardstick models coupling.
A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling provided by the invention, by dividing Simulation calculating is not carried out by mesoscale model and minute yardstick model, respectively obtain WRF Mesoscale Meteorology data and OpenFOAM minute yardstick calculated result, and then statistics is established according to anemometer tower measured data and WRF Mesoscale Meteorology data Relationship recycles statistical relationship to be modified OpenFOAM minute yardstick calculated result, obtains wind-resources corrected Calculation result.This Invention effectively combines statistics NO emissions reduction and power drop by utilizing mesoscale model and minute yardstick models coupling calculation method The advantages of two kinds of NO emissions reduction methods of scale, can either obtain the wind-resources of higher resolution in the case where surveying the insufficient situation of wind data amount Map, to quickly carry out corresponding previous work;Also the region for having a small amount of anemometer tower can be carried out beneficial to complementary and ginseng It examines, missing data situation is made up and carries out long-term correlation, it can be achieved that in mesoclimate element and region gas The area for waiting element correlation difference carries out operation and prediction, improves computational efficiency, reduces calculation amount and substantially increase Accuracy to climatic simulation.
Detailed description of the invention
Fig. 1 is that the present invention is based on the wind-resources calculation method example schemes of mesoscale model and minute yardstick models coupling to relate to And hardware running environment structural schematic diagram;
Fig. 2 is that the present invention is based on mesoscale models and the wind-resources calculation method first embodiment of minute yardstick models coupling Flow diagram;
Fig. 3 is that the present invention is based on mesoscale models and the wind-resources calculation method second embodiment of minute yardstick models coupling Flow diagram;
Fig. 4 is that the present invention is based on mesoscale models and the wind-resources calculation method second embodiment of minute yardstick models coupling The flow diagram of step S210 and step S220 refinement;
Fig. 5 is that the present invention is based on mesoscale models and the wind-resources calculation method 3rd embodiment of minute yardstick models coupling Flow diagram;
Fig. 6 is that the present invention is based on mesoscale models and the wind-resources calculation method fourth embodiment of minute yardstick models coupling Flow diagram;
Fig. 7 is that the present invention is based on mesoscale models and the 5th embodiment of wind-resources calculation method of minute yardstick models coupling Journey schematic diagram;
Fig. 8 is that the present invention is based on the functional modules of mesoscale model and the wind-resources computing device of minute yardstick models coupling to show It is intended to.
The embodiments will be further described with reference to the accompanying drawings for the realization, the function and the advantages of the object of the present invention.
Specific embodiment
The embodiment of the present invention is described below in detail, in which the same or similar labels are throughly indicated same or like Element or element with the same or similar functions.
In addition, term " first ", " second " are used for descriptive purposes only and cannot be understood as indicating or suggesting relative importance Or implicitly indicate the quantity of indicated technical characteristic.Define " first " as a result, the feature of " second " can be expressed or Implicitly include one or more of the features.In the description of the present invention, the meaning of " plurality " is two or more, Unless otherwise specifically defined.
In the present invention unless specifically defined or limited otherwise, term " installation ", " connected ", " connection ", " fixation " etc. Term shall be understood in a broad sense, for example, it may be being fixedly connected, may be a detachable connection, or integral;It can be mechanical connect It connects, is also possible to be electrically connected;It can be directly connected, can also can be in two elements indirectly connected through an intermediary The interaction relationship of the connection in portion or two elements.It for the ordinary skill in the art, can be according to specific feelings Condition understands the concrete meaning of above-mentioned term in the present invention.
It should be appreciated that the specific embodiments described herein are merely illustrative of the present invention, it is not intended to limit the present invention.
As shown in Figure 1, being the structural schematic diagram of the hardware running environment for the terminal that the embodiment of the present invention is related to.
The PC that computer equipment of the embodiment of the present invention can be is also possible to smart phone, tablet computer or with one Determine computing capability and includes the E-book reader of biomedical information acquisition equipment (such as image collecting device), MP3 broadcasting The packaged types terminal devices such as device, MP4 player, portable computer.As shown in Figure 1, the computer equipment may include: processing Device 1001, such as CPU, network interface 1004, user interface 1003, memory 1005, communication bus 1002.Wherein, communication is total Line 1002 is for realizing the connection communication between these components.User interface 1003 may include display screen, input unit such as Keyboard, remote controler, optional user interface 1003 can also include standard wireline interface and wireless interface.Network interface 1004 can Choosing may include standard wireline interface and wireless interface (such as WI-FI interface).Memory 1005 can be high-speed RAM storage Device is also possible to stable memory, such as magnetic disk storage.Memory 1005 optionally can also be independently of aforementioned processing The storage device of device 1001.Optionally, computer equipment can also include RF (Radio Frequency, radio frequency) circuit, audio Circuit, WiFi module etc..It is passed in addition, computer equipment can also configure gyroscope, barometer, hygrometer, thermometer, infrared ray The other sensors such as sensor, details are not described herein.
It will be understood by those skilled in the art that computer equipment shown in Fig. 1 does not constitute the limit to computer equipment It is fixed, it may include perhaps combining certain components or different component layouts than illustrating more or fewer components.Such as Fig. 1 institute Show, as may include in a kind of memory 1005 of computer readable storage medium operating system, data-interface control program, Network attachment procedure and wind-resources calculation procedure based on mesoscale model and minute yardstick models coupling.
The present invention provides a kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling, passes through utilization Mesoscale model and minute yardstick models coupling calculation method effectively combine statistics two kinds of NO emissions reductions of NO emissions reduction and power NO emissions reduction The advantages of method, can either obtain the wind-resources map of higher resolution in the case where surveying the insufficient situation of wind data amount, thus quickly Carry out corresponding previous work;Also the region for having a small amount of anemometer tower can be carried out beneficial to complementary and reference, to missing data Situation is made up and carries out long-term correlation, it can be achieved that poor in mesoclimate element and regional climate element correlation Area carry out operation and prediction, improve computational efficiency, reduce calculation amount and substantially increase the standard to climatic simulation Exactness.
Embodiment 1:
Referring to Fig. 2, first embodiment of the invention provides a kind of based on the wind of mesoscale model and minute yardstick models coupling money Source calculation method, wherein the mesoscale model is WRF mesoscale model, and the minute yardstick model is OpenFOAM minute yardstick Model, which comprises
Step S100 obtains the anemometer tower measured data of target area, and acquisition is corresponding with the target area again Data are analyzed as ambient field;
It should be noted that small scale model is generally divided into two kinds, one is the linear models based on the conservation of mass, a kind of It is CFD (Computational Fluid Dynamics) mode based on Fluid Mechanics Computation.Linear model calculating speed is fast, But it is appropriate only for level terrain, CFD mode is preferred for complicated landform.CFD mode be by originally in time-domain and The field discretization of continuous physical quantity in spatial domain, then calculating solution fluid motion fundamental equation by numerical value, (such as quality is moved Amount, energy equation etc.), the relationship for obtaining the fundamental physical quantity in flow field on each position and its changing over time.Since CFD is asked Solution is original fluid Fundamental Equation of Motion, contains the nonlinear effect in equation, is based on so calculated result will be much better than The linear model of the conservation of mass.
CFD mode of the OpenFOAM as open source, not only has the computing capability of commercial CFD software, can also keep away for user Exempt from expensive business software expense.And OpenFOAM is there are also opening and flexibility that commercial CFD software does not have, inside Code full disclosure, user can arbitrarily modify solver as needed;The closure, numeric format, parallel etc. of flowing automatically schemes and items Parameter carries out depth customization, facilitates user and carries out calculating fail-safe analysis and be further improved numerical procedure.
Therefore, coupling OpenFOAM mode and WRF mode are conducive to promote wind-resources assessment and wind power prediction business Horizontal and raising prediction accuracy.
Target area is the relative Repeat relative to actual measurement region, and the general range for surveying region is less than target area, For example, actual measurement regional scope be A, target area be then include survey region A range B.For including actual measurement region The target area of range is further calculated and is predicted, the accuracy of the climatic prediction for the region is helped to improve.
It is above-mentioned, then data are analyzed, it is the global analysis of data of ERA-Interim0.5 × 0.5, wherein data is to each The observational data of kind source (ground, ship, radio sounding, pilot balloon, aircraft, satellite etc.) carries out quality control and assimilation Processing, obtains the analysis of data collection again of complete set, and the element that it not only includes is more, and range is wide, and the when segment length extended, It is a comprehensive data set.It is primary that 6h is divided between the update of data.
Above-mentioned, anemometer tower is a kind of for measuring the High-rising Tower Structures of wind energy parameter, i.e., a kind of for gas near the ground The tower that stream motion conditions are observed, record.It was mostly built, was used by wind-power electricity generation enterprise, meteorology, environmental protection administration in the past In meteorological observation and atmosphere environment supervision.
Above-mentioned, anemometer tower measured data is the measured data of all anemometer towers within a certain period of time in target area.
Step S200 carries out simulation calculating to the ambient field by the mesoscale model, obtains WRF Meso-scale meteorology Mode data;Also, simulation calculating is carried out to the ambient field by the minute yardstick model, obtains OpenFOAM minute yardstick meter Calculate result;
Step S300 establishes statistics according to the anemometer tower measured data and the WRF Mesoscale Meteorology data and closes System;
Above-mentioned, WRF Mesoscale Meteorology data are prediction result, wherein the wind speed in a large amount of target areas can be obtained The prediction test result of test point, can up to tens of thousands of;And anemometer tower measured data is the anemometer tower actual measurement in target area Data, can achieve several hundred.By the anemometer tower measured data surveyed using anemometer tower to WRF Meso-scale meteorology mould Formula number is corrected, and is established relevant statistical relationship, is obtained related coefficient, and the standard of mesoscale model data can be improved by amendment Exactness.
Step S400 corrects the OpenFOAM minute yardstick calculated result using the statistical relationship, obtains wind-resources and repair Positive calculated result.
It is above-mentioned, after the statistical relationship for establishing anemometer tower measured data and WRF Mesoscale Meteorology data, that is, pass through reality After the anemometer tower measured data of survey is to WRF Mesoscale Meteorology data correction, using the obtained related coefficient of statistical relationship, The OpenFOAM minute yardstick calculated result of minute yardstick is modified to get wind-resources corrected Calculation result is arrived.
A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling provided in this embodiment, passes through Simulation calculating is carried out by mesoscale model and minute yardstick model respectively, respectively obtain WRF Mesoscale Meteorology data and OpenFOAM minute yardstick calculated result, and then statistics is established according to anemometer tower measured data and WRF Mesoscale Meteorology data Relationship recycles statistical relationship to be modified OpenFOAM minute yardstick calculated result, obtains wind-resources corrected Calculation result.This Embodiment effectively combines statistics NO emissions reduction and power drop by utilizing mesoscale model and minute yardstick models coupling calculation method The advantages of two kinds of NO emissions reduction methods of scale, can either obtain the wind-resources of higher resolution in the case where surveying the insufficient situation of wind data amount Map, to quickly carry out corresponding previous work;Also the region for having a small amount of anemometer tower can be carried out beneficial to complementary and ginseng It examines, missing data situation is made up and carries out long-term correlation, it can be achieved that in mesoclimate element and region gas The area for waiting element correlation difference carries out operation and prediction, improves computational efficiency, reduces calculation amount and substantially increase Accuracy to climatic simulation.
Embodiment 2:
Referring to Fig. 3-4, second embodiment of the invention provides a kind of wind based on mesoscale model and minute yardstick models coupling Resource Calculation method is based on above-mentioned first embodiment shown in Fig. 2, the step S200, " by the mesoscale model to institute State ambient field and carry out simulation calculating, obtain WRF Mesoscale Meteorology data " include:
Step S210 constructs the kilometer grade in the target area of WRF mesoscale model according to the mesoscale model Grid;
It should be noted that mesoscale WRF model or minute yardstick OpenFOAM model are all to belong to numerical value calculating.And numerical value The basic conception of calculating is to solve for the math equation of a physics or Chemical Problem, and math equation here is usually one group inclined The differential equation (including boundary condition, the partial differential equations of boundary condition are not incomplete).Because of partial differential equations It is difficult directly to calculate as a result, so being assumed an initial value with numerical method, then being started to iterate to calculate.Iteration needs Problem is discrete (disassembling).For example, calculate the weather in certain city, then need the region segmentation in the city is only at very Vertical and continuous subelement goes iteration, this subelement is exactly grid.The grid of resolution ratio based on kilometer grade, as kilometer Grade grid.
Step S220 is based on the ambient field, the wind speed of the kilometer grade grid is calculated for the target area, is obtained WRF Mesoscale Meteorology data.
The step S210, " according to the mesoscale model, in the target area for constructing WRF mesoscale model Kilometer grade grid " includes:
Step S211 is split the target area as unit of pre-determined distance, according to the mesoscale model structure Three layers of nested grid of the corresponding WRF mesoscale model in the target area after building based on segmentation, as kilometer grade grid.
The grid distance of three layers of nested grid is respectively 27km, 9km and 3km.
It is above-mentioned, three layers of nested grid of WRF mode are arranged according to target area, grid distance is respectively 27km, 9km, 3km.
Above-mentioned, grid is the basis of iterative calculation.By target area by 1 kilometer × 1 kilometer segmentation, just at kilometer grade net Lattice carry out the calculating of mesoscale WRF on this grid, because of the relationship of WRF model, it cannot calculate more acurrate.So If necessary to predict and know the weather on certain building doorway in target area, it is necessary to target area be pressed 1 meter × 1 meter of grid Segmentation, is further calculated with minute yardstick OpenFOAM model.
The step S220 " is based on the ambient field, the wind of the kilometer grade grid is calculated for the target area Speed obtains WRF Mesoscale Meteorology data " include:
Step S221, acquisition are combined with the matched Parameterization Scheme in the target area;
Step S222 combines according to the Parameterization Scheme, determines simulated time;
Above-mentioned, Parameterization Scheme may include microphysical processes, radiative process, road surface process, boundary layer and cumulus convection Equal parameterized procedures.In the present embodiment, the Parameterization Scheme combination for being suitble to target area, i.e., matched Parameterization Scheme are chosen Combination is set mode beginning and ending time (simulated time), carries out numerical model simulation, thus.Carry out the WRF mode of target area It calculates.It is above-mentioned, be suitble to target area Parameterization Scheme combination, be according to measuring and calculating personnel the experience in this field there are also industries Documents and materials obtain.
Step S223 is directed to the numerical model simulation of the kilometer grade grid of the target area, obtains the mould The corresponding wind speed of the kilometer grade grid each of in pseudotime, as WRF Mesoscale Meteorology data.
Above-mentioned, method provided in the present embodiment is to carry out simulation and forecast, i.e. calculation of wind speed for wind-resources.This In embodiment, by the wind speed for calculating kilometer grade grid (the usually position of lattice point) with WRF mesoscale.Calculating process essence On be mathematical iterations process.Simulation calculating is carried out to the ambient field by the mesoscale model of WRF, obtains WRF mesoscale gas As mode data, wherein constructing the kilometer grade grid of three layers of nested grid of target area, and obtained within the beginning and ending time The wind speed of each kilometer of grade grid then acquires WRF Mesoscale Meteorology data, further improves the accurate of wind speed test Property.
Embodiment 3:
Referring to Fig. 5, third embodiment of the invention provides a kind of based on the wind of mesoscale model and minute yardstick models coupling money Source calculation method is based on above-mentioned first embodiment shown in Fig. 2, the step S200, " by the minute yardstick model to described Ambient field carries out simulation calculating, obtains OpenFOAM minute yardstick calculated result " include:
Step S230, according to the minute yardstick model, in the target area for constructing OpenFOAM minute yardstick mode Refine grid;
It is above-mentioned, the corresponding fining grid in target area is constructed by trellis algorithm.The step can be with construction step S200 In two processes of kilometer grade grid carry out simultaneously, can also have sequencing.
Step S240, obtain OpenFOAM minute yardstick mode under boundary condition corresponding with the target area and initially Value;
Above-mentioned, boundary condition and initial value are according to practical experience and to calculate the ginseng inputted by user in calculating process Number.Wherein, boundary condition is a function, and the inside needs to set several parameters, and parameter has certain range.It is calculated in wind-resources The inside, this boundary condition are the Wind outline of zoning inlet, and Wind outline changes constantly, so being change in location 's.In the present embodiment to solve this problem, it goes to replace using the function of a standard.Initial value also has range.Perimeter strip The effect of part is to define the border condition of Solve problems, and the effect of initial value is to provide an initial value to iterative solution.This two A parameter that is to solve for device and must input, if without boundary condition and initial value, it is subsequent to calculate solution.
Wherein, the function of the standard, for the boundary condition for the partial differential equation that solver solves, so-called canonical function It is an approximate equation obtained in the atmosphere boundary theory, it can rough description Wind outline.
Step S250 passes through the OpenFOAM minute yardstick mode according to the boundary condition and the initial value OpenFOAM solver calculates the flow field of the fining grid of the target area, obtains the OpenFOAM minute yardstick and calculates knot Fruit.
OpenFOAM solver is the function tool for calculating and solving, and calculating to solve itself can be abstracted into one group of partial differential Equation, in addition the boundary condition and initial value of front, this group of partial differential equation can be asked by the method for mathematical iterations Solution, solver are just the function for being iterated solution.In the present embodiment, by constructing the fining grid of target area, into And according to boundary condition and initial value, it is iterated by OpenFOAM solver and calculates its flow field, it is final available OpenFOAM minute yardstick calculated result, thus in the resolution levels of small scale, or calculated value is obtained, in order to further It is coupled, improves the accuracy of operation.
Embodiment 4:
Referring to Fig. 6, fourth embodiment of the invention provides a kind of based on the wind of mesoscale model and minute yardstick models coupling money Source calculation method, is based on above-mentioned first embodiment shown in Fig. 2, the step S300, " according to the anemometer tower measured data and The WRF Mesoscale Meteorology data establish statistical relationship " include:
Step S310 corrects the WRF Mesoscale Meteorology data by the anemometer tower measured data, obtains middle ruler Degree correction related coefficient;Also, related coefficient is corrected according to the mesoscale and establishes the statistical relationship.
It is above-mentioned, because the anemometer tower limited amount in target area, such asTarget area is that a city is (above county level City) then anemometer tower quantity is probably several hundred, but the wind speed of the prediction for the target area being calculated with WRF mesoscale can be with Obtain tens of thousands of a points as a result, so in the present embodiment, first in tens of thousands of of the measured data correction of above-mentioned several hundred a anemometer towers Scale as a result, further go to correct that several hundred million minute yardsticks are calculated as a result, can be obtained with the result of mesoscale again More accurate correction value.In the present embodiment, correction is gone to predict in obtained WRF by the anemometer tower measured data using actual measurement Scale Meteorological Models data, so that the accuracy of WRF Mesoscale Meteorology data is improved, in order to further pass through in WRF Scale Meteorological Models data go to correct finally obtained OpenFOAM minute yardstick calculated result, to improve computational efficiency, drop Low calculation amount and substantially increase accuracy to climatic simulation.Also, in the present embodiment, it can be achieved that in mesoscale gas The area for waiting element and regional climate element correlation difference carries out operation and prediction, avoids provided method in the prior art The drawbacks of operation and prediction can not being carried out for above-mentioned area.
Embodiment 5:
It is provided a kind of based on mesoscale model and minute yardstick model knot in order to better illustrate the present invention referring to Fig. 6 The wind-resources calculation method of conjunction, provides following steps:
1, downloading analyzes data as ambient field again;
2, " kilometer " grade grid of WRF mesoscale model is made;
3, it is calculated with WRF mesoscale model;
4, it exports and handles Mesoscale Meteorology WRF data;
5, the fining grid of production open source CFD mode OpenFOAM;
6, OpenFOAM boundary condition and initial value are set;
7, operation OpenFOAM solver calculates flow field;
8, it exports and handles minute yardstick OpenFOAM calculated result;
9, statistical relationship is established using WRF mesoscale result and anemometer tower data;
10, minute yardstick calculated result is corrected using the statistical relationship of foundation.
This method effectively combines Mesoscale Numerical Simulation and calculates the advantages of calculating with minute yardstick simulation.
The resolution ratio that Mesoscale Numerical Simulation calculates in " kilometer " rank, using with macroscopical scope.For example, to the whole nation or the whole province Whether the wind-resources overall distribution trend investigation in region meets the requirements, but is directed to the specific Microplanner of Wind Power Project, especially It is when requiring to carry out the just calculation of Capacity Assessment and generated energy, this requires the wind-resources maps of higher resolution to expire The requirement of sufficient wind-resources engineer.
On the other hand, the spatial resolution of minute yardstick model is generally 20~100m, can accurately reflect local landform, The influence that the variation of looks generates microcosmic flow field, to the wind power plant production prediction and unit under complicated landform, complicated environmental condition Microcosmic structure provides foundation.Also, two kinds of sides are compensated for by the combination of two methods of mesoscale and minute yardstick to a certain extent The deficiency of method.
In addition, the present invention also provides a kind of based on the wind-resources of mesoscale model and minute yardstick models coupling with reference to Fig. 7 Device is calculated, the mesoscale model is WRF mesoscale model, and the minute yardstick model is OpenFOAM minute yardstick model, comprising:
Module 10 is obtained, for obtaining the anemometer tower measured data of target area, and is obtained and the target area pair That answers analyzes data as ambient field again;
Computing module 20 obtains ruler in WRF for carrying out simulation calculating to the ambient field by the mesoscale model Spend Meteorological Models data;Also, simulation calculating is carried out to the ambient field by the minute yardstick model, it is micro- to obtain OpenFOAM Dimension calculation result;
Module 30 is established, for establishing according to the anemometer tower measured data and the WRF Mesoscale Meteorology data Statistical relationship;
Correction module 40 obtains wind for correcting the OpenFOAM minute yardstick calculated result using the statistical relationship Resource corrected Calculation result.
In addition, the computer equipment includes memory and processor, institute the present invention also provides a kind of computer equipment Memory is stated for storing the wind-resources calculation procedure based on mesoscale model and minute yardstick models coupling, the processor operation The wind-resources calculation procedure based on mesoscale model and minute yardstick models coupling is so that the computer equipment executes as above State the wind-resources calculation method based on mesoscale model and minute yardstick models coupling.
In addition, being stored on the computer readable storage medium the present invention also provides a kind of computer readable storage medium There is the wind-resources calculation procedure based on mesoscale model and minute yardstick models coupling, it is described to be based on mesoscale model and minute yardstick mould The wind-resources calculation procedure that type combines is realized as described above based on mesoscale model and minute yardstick model when being executed by processor In conjunction with wind-resources calculation method.
The serial number of the above embodiments of the invention is only for description, does not represent the advantages or disadvantages of the embodiments.
Through the above description of the embodiments, those skilled in the art can be understood that above-described embodiment side Method can be realized by means of software and necessary general hardware platform, naturally it is also possible to by hardware, but in many cases The former is more preferably embodiment.Based on this understanding, technical solution of the present invention substantially in other words does the prior art The part contributed out can be embodied in the form of software products, which is stored in one as described above In storage medium (such as ROM/RAM, magnetic disk, CD), including some instructions are used so that terminal device (it can be mobile phone, Computer, server or network equipment etc.) execute method described in each embodiment of the present invention.The above is only of the invention Preferred embodiment is not intended to limit the scope of the invention, all using made by description of the invention and accompanying drawing content Equivalent structure or equivalent flow shift is applied directly or indirectly in other relevant technical fields, and is similarly included in this hair In bright scope of patent protection.

Claims (10)

1. a kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling, the mesoscale model is WRF Mesoscale model, the minute yardstick model are OpenFOAM minute yardstick model, which is characterized in that the described method includes:
The anemometer tower measured data of target area is obtained, and obtains data of analyzing again corresponding with the target area as back Jing Chang;
Simulation calculating is carried out to the ambient field by the mesoscale model, obtains WRF Mesoscale Meteorology data;And And simulation calculating is carried out to the ambient field by the minute yardstick model, obtain OpenFOAM minute yardstick calculated result;
Statistical relationship is established according to the anemometer tower measured data and the WRF Mesoscale Meteorology data;
The OpenFOAM minute yardstick calculated result is corrected using the statistical relationship, obtains wind-resources corrected Calculation result.
2. the wind-resources calculation method based on mesoscale model and minute yardstick models coupling, feature exist as described in claim 1 In described " carrying out simulation calculating to the ambient field by the mesoscale model, obtain WRF Mesoscale Meteorology data " Include:
According to the mesoscale model, the kilometer grade grid in the target area of WRF mesoscale model is constructed;
Based on the ambient field, the wind speed of the kilometer grade grid is calculated for the target area, obtains WRF Meso-scale meteorology Mode data.
3. the wind-resources calculation method based on mesoscale model and minute yardstick models coupling, feature exist as claimed in claim 2 In " according to the mesoscale model, constructing the kilometer grade grid in the target area of the WRF mesoscale model " packet It includes:
The target area is split as unit of pre-determined distance, after being constructed according to the mesoscale model based on segmentation Three layers of nested grid of the corresponding WRF mesoscale model in the target area, as kilometer grade grid.
4. the wind-resources calculation method based on mesoscale model and minute yardstick models coupling, feature exist as claimed in claim 3 In the grid distance of three layers of nested grid is respectively 27km, 9km and 3km.
5. the wind-resources calculation method based on mesoscale model and minute yardstick models coupling, feature exist as claimed in claim 2 In described " to be based on the ambient field, the wind speed of the kilometer grade grid is calculated for the target area, obtains WRF mesoscale Meteorological Models data " include:
Acquisition is combined with the matched Parameterization Scheme in the target area;
It is combined according to the Parameterization Scheme, determines simulated time;
Each of it is directed to the numerical model simulation of the kilometer grade grid of the target area, obtain in the simulated time The corresponding wind speed of the kilometer grade grid, as WRF Mesoscale Meteorology data.
6. the wind-resources calculation method based on mesoscale model and minute yardstick models coupling, feature exist as described in claim 1 In described " to carry out simulation calculating to the ambient field by the minute yardstick model, obtain OpenFOAM minute yardstick and calculate knot Fruit " includes:
According to the minute yardstick model, the fining grid in the target area of OpenFOAM minute yardstick mode is constructed;
Obtain the boundary condition corresponding with the target area and initial value under OpenFOAM minute yardstick mode;
According to the boundary condition and the initial value, pass through the OpenFOAM solver meter of the OpenFOAM minute yardstick mode The flow field for calculating the fining grid of the target area, obtains the OpenFOAM minute yardstick calculated result.
7. the wind-resources calculation method based on mesoscale model and minute yardstick models coupling, feature exist as described in claim 1 In " the establishing statistical relationship according to the anemometer tower measured data and the WRF Mesoscale Meteorology data " includes:
The WRF Mesoscale Meteorology data are corrected by the anemometer tower measured data, obtain mesoscale correction phase relation Number;Also, related coefficient is corrected according to the mesoscale and establishes the statistical relationship.
8. a kind of wind-resources computing device based on mesoscale model and minute yardstick models coupling, the mesoscale model is WRF Mesoscale model, the minute yardstick model are OpenFOAM minute yardstick model characterized by comprising
Module is obtained, for obtaining the anemometer tower measured data of target area, and acquisition is corresponding with the target area again Data are analyzed as ambient field;
Computing module obtains WRF Meso-scale meteorology for carrying out simulation calculating to the ambient field by the mesoscale model Mode data;Also, simulation calculating is carried out to the ambient field by the minute yardstick model, obtains OpenFOAM minute yardstick meter Calculate result;
Module is established, is closed for establishing statistics according to the anemometer tower measured data and the WRF Mesoscale Meteorology data System;
Correction module obtains wind-resources and repairs for correcting the OpenFOAM minute yardstick calculated result using the statistical relationship Positive calculated result.
9. a kind of computer equipment, which is characterized in that the computer equipment includes memory and processor, the memory For storing the wind-resources calculation procedure based on mesoscale model and minute yardstick models coupling, it is based on described in the processor operation The wind-resources calculation procedure of mesoscale model and minute yardstick models coupling is so that the computer equipment executes such as claim 1-7 Any one of described in the wind-resources calculation method based on mesoscale model and minute yardstick models coupling.
10. a kind of computer readable storage medium, which is characterized in that be stored on the computer readable storage medium in being based on The wind-resources calculation procedure of Scale Model and minute yardstick models coupling, it is described based on mesoscale model and minute yardstick models coupling It is realized when wind-resources calculation procedure is executed by processor and is based on mesoscale model and micro- ruler as described in any one of claim 1-7 Spend the wind-resources calculation method of models coupling.
CN201811467690.6A 2018-12-03 2018-12-03 A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling Pending CN109583096A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201811467690.6A CN109583096A (en) 2018-12-03 2018-12-03 A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201811467690.6A CN109583096A (en) 2018-12-03 2018-12-03 A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling

Publications (1)

Publication Number Publication Date
CN109583096A true CN109583096A (en) 2019-04-05

Family

ID=65925930

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201811467690.6A Pending CN109583096A (en) 2018-12-03 2018-12-03 A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling

Country Status (1)

Country Link
CN (1) CN109583096A (en)

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110533347A (en) * 2019-09-10 2019-12-03 浙江运达风电股份有限公司 A kind of wind energy turbine set wind-resources calculation method, device, equipment and readable medium
CN111339714A (en) * 2020-02-17 2020-06-26 珠江水利委员会珠江水利科学研究院 Multi-scale hydrodynamic coupling method based on FVCOM and OpenFOAM models
CN111680408A (en) * 2020-05-26 2020-09-18 中国能源建设集团广东省电力设计研究院有限公司 Wind resource map drawing method and device for offshore wind power
CN111814345A (en) * 2020-07-17 2020-10-23 中国船舶重工集团海装风电股份有限公司 Method for interpolating and observing missing wind speed data by utilizing various mesoscale wind speed data
CN112632792A (en) * 2020-12-30 2021-04-09 中交第一航务工程局有限公司 Shore-approaching wave reanalysis simulation system based on measured data and numerical simulation

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103778572A (en) * 2014-02-24 2014-05-07 南方电网科学研究院有限责任公司 Sea wind resource evaluation method based on WRF mode
CN104331621A (en) * 2014-11-05 2015-02-04 中国大唐集团新能源股份有限公司 Wind resource computation method
CN106326625A (en) * 2015-07-09 2017-01-11 中国电力科学研究院 Method for simulating wind field by coupling WRF (weather research and forecasting) and OpenFOAM modes
CN108536881A (en) * 2017-03-06 2018-09-14 新疆金风科技股份有限公司 Method and apparatus for calculating wind power plant generated energy
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 (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103778572A (en) * 2014-02-24 2014-05-07 南方电网科学研究院有限责任公司 Sea wind resource evaluation method based on WRF mode
CN104331621A (en) * 2014-11-05 2015-02-04 中国大唐集团新能源股份有限公司 Wind resource computation method
CN106326625A (en) * 2015-07-09 2017-01-11 中国电力科学研究院 Method for simulating wind field by coupling WRF (weather research and forecasting) and OpenFOAM modes
CN108536881A (en) * 2017-03-06 2018-09-14 新疆金风科技股份有限公司 Method and apparatus for calculating wind power plant generated energy
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 (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110533347A (en) * 2019-09-10 2019-12-03 浙江运达风电股份有限公司 A kind of wind energy turbine set wind-resources calculation method, device, equipment and readable medium
CN111339714A (en) * 2020-02-17 2020-06-26 珠江水利委员会珠江水利科学研究院 Multi-scale hydrodynamic coupling method based on FVCOM and OpenFOAM models
CN111680408A (en) * 2020-05-26 2020-09-18 中国能源建设集团广东省电力设计研究院有限公司 Wind resource map drawing method and device for offshore wind power
CN111814345A (en) * 2020-07-17 2020-10-23 中国船舶重工集团海装风电股份有限公司 Method for interpolating and observing missing wind speed data by utilizing various mesoscale wind speed data
CN111814345B (en) * 2020-07-17 2022-08-02 中国船舶重工集团海装风电股份有限公司 Method for interpolating and observing missing wind speed data by utilizing various mesoscale wind speed data
CN112632792A (en) * 2020-12-30 2021-04-09 中交第一航务工程局有限公司 Shore-approaching wave reanalysis simulation system based on measured data and numerical simulation
CN112632792B (en) * 2020-12-30 2022-03-04 中交第一航务工程局有限公司 Shore-approaching wave reanalysis simulation system based on measured data and numerical simulation

Similar Documents

Publication Publication Date Title
CN109583096A (en) A kind of wind-resources calculation method based on mesoscale model and minute yardstick models coupling
Sailor et al. A neural network approach to local downscaling of GCM output for assessing wind power implications of climate change
Lundquist et al. An immersed boundary method enabling large-eddy simulations of flow over complex terrain in the WRF model
Gall et al. The hurricane forecast improvement project
Zadra et al. The next Canadian regional climate model
Quan et al. Urban data and building energy modeling: A GIS-based urban building energy modeling system using the urban-EPC engine
US9310517B2 (en) Method and system for determining accuracy of a weather prediction model
Barker et al. A three-dimensional variational data assimilation system for MM5: Implementation and initial results
Mohan et al. WRF model performance analysis for a suite of simulation design
Rauscher et al. Exploring a global multiresolution modeling approach using aquaplanet simulations
Zhang et al. Assimilation of precipitation-affected radiances in a cloud-resolving WRF ensemble data assimilation system
CN108563867A (en) A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM
CN102945508B (en) Model correction based wind power forecasting method
Zupanski et al. A prototype WRF-based ensemble data assimilation system for dynamically downscaling satellite precipitation observations
Veronesi et al. Statistical learning approach for wind resource assessment
CN106443833B (en) A kind of numerical weather forecast method based on cloud computing
CN106326625A (en) Method for simulating wind field by coupling WRF (weather research and forecasting) and OpenFOAM modes
CN105760970A (en) Method for predicting AQI
CN105243435A (en) Deep learning cellular automaton model-based soil moisture content prediction method
JP4227889B2 (en) Concentration prediction method and apparatus for air pollutants, program, and apparatus for providing concentration prediction information for air pollutants
CN109543907B (en) Complex terrain wind resource assessment method and device
CN107784165A (en) Surface temperature field multi-scale data assimilation method based on photovoltaic plant
Teixeira et al. Surface to boundary layer coupling in the urban area of Lisbon comparing different urban canopy models in WRF
CN103413036A (en) Continuous forest fire weather level forecasting model and application thereof
Zhang et al. Large-eddy simulation of traffic-related air pollution at a very high resolution in a mega-city: evaluation against mobile sensors and insights for influencing factors

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

Application publication date: 20190405