CN114169576A - Wind resource calculation method and device and electronic equipment - Google Patents

Wind resource calculation method and device and electronic equipment Download PDF

Info

Publication number
CN114169576A
CN114169576A CN202111344542.7A CN202111344542A CN114169576A CN 114169576 A CN114169576 A CN 114169576A CN 202111344542 A CN202111344542 A CN 202111344542A CN 114169576 A CN114169576 A CN 114169576A
Authority
CN
China
Prior art keywords
calculation
power plant
wind
wind power
sectors
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
CN202111344542.7A
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.)
Guodian United Power Technology Co Ltd
Original Assignee
Guodian United Power Technology 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 Guodian United Power Technology Co Ltd filed Critical Guodian United Power Technology Co Ltd
Priority to CN202111344542.7A priority Critical patent/CN114169576A/en
Publication of CN114169576A publication Critical patent/CN114169576A/en
Pending legal-status Critical Current

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q10/00Administration; Management
    • G06Q10/04Forecasting or optimisation specially adapted for administrative or management purposes, e.g. linear programming or "cutting stock problem"
    • G06Q10/043Optimisation of two dimensional placement, e.g. cutting of clothes or wood
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
    • G06Q50/00Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
    • G06Q50/06Energy or water supply

Landscapes

  • Business, Economics & Management (AREA)
  • Engineering & Computer Science (AREA)
  • Human Resources & Organizations (AREA)
  • Economics (AREA)
  • Strategic Management (AREA)
  • Physics & Mathematics (AREA)
  • Health & Medical Sciences (AREA)
  • Theoretical Computer Science (AREA)
  • Marketing (AREA)
  • General Physics & Mathematics (AREA)
  • General Business, Economics & Management (AREA)
  • Tourism & Hospitality (AREA)
  • Game Theory and Decision Science (AREA)
  • Quality & Reliability (AREA)
  • Operations Research (AREA)
  • Entrepreneurship & Innovation (AREA)
  • Development Economics (AREA)
  • Public Health (AREA)
  • Water Supply & Treatment (AREA)
  • General Health & Medical Sciences (AREA)
  • Primary Health Care (AREA)
  • Wind Motors (AREA)

Abstract

The invention provides a wind resource calculation method, a wind resource calculation device and electronic equipment, wherein the method comprises the following steps: acquiring a basic example file and a parameter set of a wind power plant; the wind power plant comprises a plurality of sectors, and the parameter set comprises parameters corresponding to each sector; automatically calculating a plurality of sectors through OpenFOAM based on the basic example file and the parameters corresponding to each sector to obtain the calculation result of each sector; and obtaining wind resource distribution information of the wind power plant according to the calculation results of the plurality of sectors. In the method, the multiple sectors of the wind power plant are automatically calculated through the OpenFOAM, and compared with the method that parameters are manually set before simulation calculation of each sector in the existing method, the calculation efficiency of the multiple sectors and the calculation efficiency of wind resource distribution information of the wind power plant are improved, so that the evaluation efficiency of wind resources is improved, the micro-site selection of the wind power plant in the later period is facilitated, and the method has important practical value.

Description

Wind resource calculation method and device and electronic equipment
Technical Field
The invention relates to the technical field of new energy wind power generation, in particular to a wind resource calculation method, a wind resource calculation device and electronic equipment.
Background
With the continuous development of the wind power generation industry, the micro site selection of the wind power plant is gradually developed from a simple flat terrain to a complex multi-hill gully terrain, and the success or failure of the micro site selection directly influences the economy of each wind power plant, so the micro site selection of the wind power plant has important significance. The wind resource assessment is an important support for micro-site selection of the wind power plant, the wind resource assessment is quantitative assessment of wind energy distribution and turbulence intensity characteristics in a wind power plant flow field, and CFD (Computational Fluid Dynamics) simulation is an important and mature means for wind resource assessment, so that not only can higher computing precision be realized, but also higher computing speed can be realized, a large amount of construction cost can be saved, and the wind resource assessment is widely applied.
In practical application, the range of a wind power plant is generally from dozens of square kilometers to hundreds of square kilometers, and for CFD simulation in such a large range, the conventional method mainly divides the wind power plant into a plurality of sectors, and sets simulation parameters for each sector manually before simulation, so that a large amount of time cost is consumed, the speed of wind resource evaluation is reduced, and the requirements of practical application cannot be met.
Disclosure of Invention
In view of the above, the present invention provides a wind resource calculation method, a wind resource calculation device, and an electronic device, so as to alleviate the above problems, a plurality of sectors of a wind farm are automatically calculated through OpenFOAM, and compared with the existing method in which parameters are manually set before simulation calculation of each sector, the calculation efficiency of the plurality of sectors and the calculation efficiency of wind resource distribution information of the wind farm are improved, so that the evaluation efficiency of wind resources is improved, further, micro-site selection of the wind farm in a later period is facilitated, and the method has an important practical value.
In a first aspect, an embodiment of the present invention provides a wind resource calculation method, which is applied to an electronic device configured with an OpenFOAM, and includes: acquiring a basic example file and a parameter set of a wind power plant; the wind power plant comprises a plurality of sectors, the basic sample files are used for representing basic description information of the wind power plant, the parameter sets comprise parameters corresponding to each sector, and the parameters comprise at least one of the following parameters: the incoming flow direction, boundary parameters and turbulent kinetic energy dissipation rate; automatically calculating a plurality of sectors through OpenFOAM based on the basic example file and the parameters corresponding to each sector to obtain the calculation result of each sector; wherein the calculation results comprise a wind speed result and a turbulence energy result; and obtaining wind resource distribution information of the wind power plant according to the calculation results of the plurality of sectors.
Preferably, an embodiment of the present invention provides a first possible implementation manner of the first aspect, where the parameters further include: calculating the number of threads in parallel; the step of automatically calculating the plurality of sectors through OpenFOAM based on the basic example file and the parameter corresponding to each sector includes: copying a basic arithmetic example file for each sector, and dividing the copied basic arithmetic example file into a plurality of blocks; and based on OpenFOAM, performing parallel computation on a plurality of blocks according to the number of parallel computation threads corresponding to the sector, and obtaining the computation result of the sector according to the computation sub-result of each block.
Preferably, an embodiment of the present invention provides a second possible implementation manner of the first aspect, where before the step of obtaining the basic example file of the wind farm, the method further includes: acquiring topographic data of a wind power plant; determining a corresponding calculation domain of the wind power plant according to the topographic data; carrying out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant; acquiring basic setting information of a wind power plant; the basic setting information comprises boundary surface information, boundary conditions and monitoring positions of the wind power plant; and generating a basic operation example file of the wind power plant according to the calculation domain grid and the basic setting information.
Preferably, an embodiment of the present invention provides a third possible implementation manner of the first aspect, where the step of performing grid division on the computation domain to generate a computation domain grid corresponding to the wind farm includes: a blockMesh module based on OpenFOAM carries out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant; or carrying out grid division on the calculation domain according to a preset mode to generate an initial calculation domain grid, and converting the initial calculation domain grid into a grid which can be identified by OpenFOAM to obtain a calculation domain grid corresponding to the wind power plant; wherein, the preset mode comprises one of the following modes: gridggen, pointwise, gambit, icemcfd, tetgen, gmesh, and ansys.
Preferably, an embodiment of the present invention provides a fourth possible implementation manner of the first aspect, wherein the step of acquiring the topographic data of the wind farm includes: acquiring original topographic data of a wind power plant; carrying out operation processing on the original topographic data to obtain topographic data; wherein the operation process comprises at least one of: cleaning operation and repairing operation.
Preferably, an embodiment of the present invention provides a fifth possible implementation manner of the first aspect, where each sector carries priority information, and the step of automatically calculating the plurality of sectors through OpenFOAM includes: and according to the priority information of the plurality of sectors, automatically calculating the plurality of sectors in sequence through OpenFOAM.
Preferably, an embodiment of the present invention provides a sixth possible implementation manner of the first aspect, where the step of automatically calculating the plurality of sectors through OpenFOAM further includes: acquiring a calculation sequence of a plurality of sectors set by a user; and according to the calculation sequence, automatically calculating the plurality of sectors in sequence through OpenFOAM.
In a second aspect, an embodiment of the present invention further provides a wind resource calculation apparatus, which is applied to an electronic device configured with OpenFOAM, and the apparatus includes: the acquisition module is used for acquiring basic example files and parameter sets of the wind power plant; the wind power plant comprises a plurality of sectors, the basic sample files are used for representing basic description information of the wind power plant, the parameter sets comprise parameters corresponding to each sector, and the parameters comprise at least one of the following parameters: the incoming flow direction, boundary parameters and turbulent kinetic energy dissipation rate; the first calculation module is used for automatically calculating a plurality of sectors through OpenFOAM based on the basic example file and the parameters corresponding to each sector to obtain the calculation result of each sector; wherein the calculation results comprise a wind speed result and a turbulence energy result; and the second calculation module is used for obtaining wind resource distribution information of the wind power plant according to the calculation results of the plurality of sectors.
In a third aspect, an embodiment of the present invention further provides an electronic device, which includes a memory, a processor, and a computer program stored in the memory and executable on the processor, where the processor implements the steps of the method in the first aspect when executing the computer program.
In a fourth aspect, the present invention further provides a computer-readable storage medium, on which a computer program is stored, where the computer program is executed by a processor to perform the steps of the method in the first aspect.
The embodiment of the invention has the following beneficial effects:
the embodiment of the invention provides a wind resource calculation method, a device and electronic equipment, wherein a basic calculation example file and parameters of each sector are preset, and a plurality of sectors of a wind power plant are automatically calculated through OpenFOAM.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and drawings.
In order to make the aforementioned and other objects, features and advantages of the present invention comprehensible, preferred embodiments accompanied with figures are described in detail below.
Drawings
In order to more clearly illustrate the embodiments of the present invention or the technical solutions in the prior art, the drawings used in the description of the embodiments or the prior art will be briefly described below, and it is obvious that the drawings in the following description are some embodiments of the present invention, and other drawings can be obtained by those skilled in the art without creative efforts.
Fig. 1 is a flowchart of a wind resource calculation method according to an embodiment of the present invention;
FIG. 2 is a flow chart of another wind resource calculation method according to an embodiment of the present invention;
FIG. 3 is a schematic diagram of a wind resource calculating device according to an embodiment of the present invention;
fig. 4 is a schematic structural diagram of an electronic device according to an embodiment of the present invention.
Detailed Description
To make the objects, technical solutions and advantages of the embodiments of the present invention clearer, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings, and it is apparent that the described embodiments are some, but not all embodiments of the present invention. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present invention.
Aiming at the problem that the evaluation efficiency of wind resources is low because parameters need to be set manually before simulation calculation of each sector in the existing method, the embodiment of the invention provides a wind resource calculation method, a wind resource calculation device and electronic equipment, which improves the calculation efficiency of a plurality of sectors, improves the calculation efficiency of the plurality of sectors and the calculation efficiency of wind resource distribution information of a wind power plant, improves the evaluation efficiency of the wind resources, is convenient for micro site selection of the wind power plant in the later period, and has important practical value.
To facilitate understanding of the present embodiment, first, a wind resource calculation method provided by the present embodiment is described in detail below. The execution main body is the electronic device configured with the OpenFOAM, and in practical application, the OpenFOAM is used as an open-source CFD mode, so that the method has better computing capacity, avoids expensive software cost, reduces application cost and is convenient to popularize and apply in practical application. In addition, the electronic device may be a PC (Personal Computer), or may also be a mobile intelligent device, such as a smart phone and a palm Computer, which may be specifically set according to actual situations, and the embodiment of the present invention does not limit the description.
Based on the above electronic device, an embodiment of the present invention provides a wind resource calculation method, as shown in fig. 1, the method includes the following steps:
step S102, acquiring a basic example file and a parameter set of the wind power plant; the wind power plant comprises a plurality of sectors, the basic sample files are used for representing basic description information of the wind power plant, and the parameter sets comprise parameters corresponding to each sector; the parameters include at least one of: the incoming flow direction, boundary parameters and turbulent kinetic energy dissipation rate;
specifically, before simulation calculation, a basic example file and a parameter set of the wind power plant are obtained, wherein the basic example file is compiled according to the setting operation of a user or an operator and is used for representing basic description information of the wind power plant; in addition, the wind power plant is further divided into a plurality of sectors according to wind directions, and an operator sets parameters of each sector in advance so as to automatically calculate the plurality of sectors by using a simpleFoam solver and taking a wind power plant flow field as a simulation object based on open source computational fluid dynamics software OpenFOAM.
Before the step of obtaining the basic example file of the wind power plant, the method also needs to compile the basic example file, and the specific process comprises the following steps:
(I) acquiring topographic data of a wind power plant; specifically, the method comprises the steps of firstly, obtaining original topographic data of a wind power plant; if a wind power plant file in stl format is output by adopting a GlobalMapper, and the wind power plant file is extended to the lowest point of a calculation domain to obtain original topographic data of the wind power plant; then, carrying out operation processing on the original topographic data to obtain topographic data; wherein the operation process comprises at least one of: cleaning operation and repairing operation, such as checking and repairing a missing surface, a broken surface, a repeated surface and the like of original topographic data to obtain topographic data;
(II) determining a corresponding calculation domain of the wind power plant according to the topographic data; specifically, an initial calculation domain of the wind power plant is determined according to topographic data of the wind power plant, the initial calculation domain is subjected to external development according to a wind power plant position coordinate range in the topographic data, the calculation domain range which is reasonably designed is used as a calculation domain corresponding to the wind power plant, the calculation domain of the wind power plant is enriched, the evaluation precision of wind resources is guaranteed, and the precision of micro site selection of the wind power plant in the later period is improved;
(III) carrying out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant; specifically, one possible way of meshing is: a blockMesh module based on OpenFOAM carries out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant; for example, a blockMesh module is adopted to generate an initial grid, and an area where a wind power plant is located is encrypted by adjusting a grid scale factor; and generating a mixed grid (comprising a terrain grid and a boundary grid) by adopting a snappyHexMesh module, and generating a boundary layer grid near the ground, wherein the total height of the boundary layer grid needs to cover the height of the hub. Another possible way of meshing is: carrying out grid division on the calculation domain according to a preset mode to generate an initial calculation domain grid, and converting the initial calculation domain grid into a grid which can be identified by OpenFOAM to obtain a calculation domain grid corresponding to the wind power plant; wherein, the preset mode comprises one of the following modes: gridgen, pointwise, gambit, icemcfd, tetgen, gmesh and ansys, which can be specifically set according to actual conditions;
(IV) acquiring basic setting information of the wind power plant; specifically, the basic setting information includes, but is not limited to: boundary surface information, boundary conditions and monitoring positions;
the boundary surface information comprises boundary surface names and boundary surface types, wherein the boundary surface names are respectively east, west, south, normal, top and terrain according to the azimuth; in the boundary surface type setting, the side1 and the side2 are added as spare boundary surfaces, the type is a slip boundary, in practical application, if the incoming flow direction is west wind, the attribute of the spare boundary surface needs to be set, and if the incoming flow direction is north-south wind, the attribute of the spare boundary surface does not need to be set. Boundary conditions include, but are not limited to: a reference height, a reference wind speed, a reference flow direction, a reference height, and a turbulence kinetic energy dissipation ratio; if the reference height is set to be 100 meters, the reference wind speed is 10m/s, a reference flow direction (namely an incoming flow direction) is set to be southwest wind, namely a wind direction vector is set to be (110), the reference height is set to be a coordinate at the lowest part of the terrain, the turbulent kinetic energy dissipation rate is set to be a function of the reference wind speed, and the like; and, the monitoring position can adopt a sample tool to monitor the speed, the turbulence energy value and the like at the hub height of the machine position and the anemometer tower.
In practical application, the parameters of each sector can be set through the parameter dictionary; wherein the parameter dictionary comprises k, epsilon, p and U; where k denotes boundary surface information, epsilon denotes turbulence energy dissipation ratio, p denotes wind pressure, and U denotes wind speed, as in the parameter dictionary, the boundary conditions corresponding to east inlet, east outlet, west inlet, west outlet, south inlet, south outlet, north inlet, and north outlet may be set and named easting, eastoutlet, westinglet, westoutlet, souuthinelet, souuthoutlet, northitlet, northoutlet, etc., respectively, and the incoming wind direction and turbulence energy dissipation ratio of each sector may also be set, and the like.
(V) generating a basic example file of the wind power plant according to the calculation domain grid and the basic setting information; furthermore, a shell script file is compiled by combining the calculation domain grid and the basic setting information, so that the basic example file of the wind power plant is driven to copy, modify parameters, calculate and the like.
Further, in the shell script file, the wind farm is further divided into a plurality of sectors according to the wind direction, for example, 12 wind directions are set, the wind farm can be divided into 12 sectors, and the number of the specific sector divisions can be set according to the actual situation, which is not limited to be described in the embodiment of the present invention.
Step S104, automatically calculating a plurality of sectors through OpenFOAM based on the basic example file and the parameters corresponding to each sector to obtain the calculation result of each sector;
specifically, during calculation for each sector, OpenFOAM acquires a basic algorithm file and parameters corresponding to the sector, modifies the basic algorithm file based on the parameters of the sector, and then performs simulation calculation on the modified basic algorithm file, for example, the modified basic algorithm file corresponding to each sector is input into OpenFOAM, and simulation calculation is performed on each sector based on a simpleFoam solver in combination with a wall function, a turbulence model, a profile calculation model, and the like, so as to obtain a calculation result of the sector, where the calculation result includes but is not limited to a wind speed result and a turbulence result, for example, wind pressure p and wind speed U of the calculated sector are displayed in a parameter dictionary, and the calculation result may be a result corresponding to each of a plurality of machine sites in the sector. Therefore, when the multiple sectors are automatically calculated through OpenFOAM, each sector modifies the basic calculation example file according to the corresponding parameters, so that the calculation accuracy and independence of each sector are guaranteed, the wind resource assessment accuracy is further guaranteed, and the micro site selection accuracy of the wind power plant in the later period is improved.
In one possible automatic calculation mode, each sector carries priority information, and the process of automatically calculating the plurality of sectors through OpenFOAM is as follows: and according to the priority information of the plurality of sectors, automatically calculating the plurality of sectors in sequence through OpenFOAM. For example, priority information of the sector is determined according to the number information of the sector, and OpenFOAM automatically calculates a plurality of sectors one by one according to the number information until all sectors are calculated. In another possible automatic calculation manner, the process of automatically calculating the plurality of sectors through OpenFOAM is as follows: acquiring a calculation sequence of a plurality of sectors set by a user; and according to the calculation sequence, automatically calculating the plurality of sectors in sequence through OpenFOAM. For example, before the calculation is started, a user sets a calculation sequence of a plurality of sectors in advance, and in the calculation process, OpenFOAM automatically calculates the plurality of sectors in sequence according to the calculation sequence until the calculation of all the sectors is finished. It should be noted that, a user may also set only a partial sector of the plurality of sectors for calculation, so that OpenFOAM automatically calculates the set partial sector, and may specifically perform setting according to actual situations.
In addition, in the actual calculation process, because the range of the wind farm is generally dozens of square kilometers to hundreds of square kilometers, for the CFD simulation in such a large range, if the calculation accuracy needs to be ensured, a grid with sufficient compactness needs to be drawn, and the number of the grids is generally in the level of ten million, so that parallel calculation needs to be adopted in order to improve the calculation efficiency. If serial calculation is adopted, the calculation time is increased by geometric multiples, the wind power plant simulation generally needs to perform multi-sector calculation, and if manual setting is performed on each sector, a large amount of time cost is consumed.
Further, in order to improve the calculation efficiency of each sector, the basic algorithm file copied by each sector is also decomposed into a plurality of blocks, and the number of the divided blocks can be specifically set according to the actual situation. At this time, the parameter corresponding to the sector further includes the number of parallel computing threads, where the number of parallel computing threads can be obtained in advance based on the basic algorithm file test. For example, when parallel testing is performed on the basic example file, the parallel computing number with the highest computing efficiency is used as the parallel computing thread number, and it should be noted that the parallel computing thread number of each sector may be the same or different.
Specifically, for each sector, copying a basic algorithm file, and dividing the copied basic algorithm file into a plurality of blocks; and on the basis of OpenFOAM, performing parallel computation on a plurality of blocks according to the number of parallel computation threads corresponding to the sector, and obtaining the computation result of the sector according to the computation sub-result of each block, so that the computation efficiency of the sector can be improved by performing parallel computation on a plurality of blocks divided by basic computation example files corresponding to the same sector, and the evaluation efficiency of wind resources is improved.
And S106, obtaining wind resource distribution information of the wind power plant according to the calculation results of the plurality of sectors.
Specifically, after the calculation result of each sector is obtained through the above calculation, analysis is performed according to the calculation results of the plurality of sectors to obtain wind resource distribution information of the wind farm, where the wind resource distribution information includes, but is not limited to, wind turbine related data in the wind farm, such as average wind speed, wind direction, and turbulence energy throughout the year, and the turbulence intensity may also be calculated according to the turbulence energy, for example, the turbulence intensity is calculated according to the following formula:
Figure BDA0003352329580000101
wherein, IpExpressing turbulence intensity, i.e. the intensity of the directional ambient turbulence calculated at the aircraft position, k expressing turbulence pulsation kinetic energy, i.e. the turbulence kinetic energy calculated at the height of the hub of the aircraft position, Ux、UyAnd UzIndicating the position of the aircraftThree velocity components of the hub height. It should be noted that, for the specific calculation process of each parameter in the wind resource distribution information, reference may be made to the prior art, and details of the embodiment of the present invention are not described herein again.
In addition, in the analysis process, the wind resource distribution information of the wind farm can be automatically generated according to the calculation results of the plurality of sectors, or the wind resource distribution information of the wind farm can be generated by combining the analysis operation of an operator, for example, combining the experience information input by the operator, such as inputting the experience information through a command line, so that the micro-site selection of the wind farm can be better performed according to the wind resource distribution information of the wind farm.
The embodiment of the invention provides a wind resource calculation method, which is characterized in that a basic calculation example file and parameters of each sector are preset, and a plurality of sectors of a wind power plant are automatically calculated through OpenFOAM.
For ease of understanding, 12 sectors are illustrated here. As shown in fig. 2, the wind resource calculation method provided in the embodiment of the present invention specifically includes the following processes:
(1) acquiring topographic data of the wind power plant, and processing the topographic data; if the Globalmapper is adopted to output as an stl format file, extending to the lowest point of a calculation domain, and checking and repairing a missing face, a broken face, a repeated face and the like to enable the data to be complete topographic data;
(2) naming boundaries and generating a computational domain grid; specifically, the method comprises the steps of carrying out external expansion according to a coordinate range of a wind power plant position, designing a reasonable calculation domain range, generating an initial grid by adopting blockMesh, and encrypting an area where the wind power plant is located by adjusting a grid scale factor; a snappyHexmesh is adopted to generate a mixed grid, a boundary layer grid is generated near the ground, and the total height of the boundary layer grid needs to cover the height of a hub; setting the names of boundary surfaces as east, west, south, normal, top and terrain according to the azimuth; it should be noted that the computational domain grid is a set of grids for representing the computational domain range, and the specific grid size is determined according to various factors such as the computational domain range;
(3) setting a boundary condition; specifically, the corresponding boundary conditions are set according to the wind direction, including but not limited to: reference height, reference wind speed, reference flow direction, reference height, turbulent kinetic energy dissipation ratio, and the like may be specifically referred to the above embodiments, and the embodiments of the present invention are not described in detail herein.
(4) Setting machine position parameter monitoring; namely, setting a monitoring position, such as adopting a sample tool to monitor the speed, the turbulence energy value and the like at the machine position and the hub height of the anemometer tower; therefore, a basic example file of the wind power plant is generated by setting boundary conditions and setting machine position parameter monitoring and combining the computational domain grid;
(5) setting parallel computation; if the basic calculation example file is subjected to parallel calculation test, the number of parallel calculation threads with the highest calculation efficiency is obtained, and a parallel calculation dictionary file is set according to the number of the parallel calculation threads; here, the number of parallel computing threads per sector may be set by a parallel computing dictionary file; in addition, in practical application, the number of parallel computing threads of each sector may also be consistent with the number of blocks divided by the basic example file corresponding to the sector, and may be specifically set according to practical situations;
(6) setting sectors in the basic example; if 12 sectors are set: 0 °, 30 °, 60 °, 90 °, 120 °, 150 °, 180 °, 210 °, 240 °, 270 °, 300 °, and 330 °;
(7) compiling an execution script; compiling a script file by adopting Shell language, compiling parameters (such as wind direction angle) of 12 sectors into the script file, and the like;
(8) automatically calculating a plurality of sectors through OpenFOAM; specifically, firstly, sector naming is carried out, and incoming flow directions and boundary names of different sectors are modified; then copying a basic example file and executing dictionary file modification; finally, carrying out partition calculation; taking a 90-degree Sector as an example, firstly naming the 90-degree Sector as Sector90, copying a basic example file as Sector90, setting side1 as south and side2 as normal in Sector90, and changing the incoming flow direction of the Sector from (110) in the basic example file to (-100) because the incoming flow direction is normal; and changing eastrinet in the U, p, k and epsilon parameter dictionaries to east, and changing westoutlet to west. Partitioning the basic algorithm file corresponding to the 90-degree sector, for example, re-partitioning the computational domain grid in the basic algorithm file to obtain a plurality of blocks, performing parallel computation according to the number of corresponding parallel computing threads until a convergence condition set by the basic algorithm file is reached, merging the computation sub-results of each block, deleting the partition file, and only retaining the computation sub-results of each block so as to obtain the computation result of the 90-degree sector according to the computation sub-results of each block; it should be noted that the parameters of each sector further include iteration step number or convergence condition, and in the calculation process, each block corresponding to the sector performs iteration calculation until the iteration step number or the convergence condition is reached, and the calculation is stopped;
(9) repeating the step (8), changing the parameters of different sectors into corresponding incoming flow direction vectors, changing corresponding boundary names, and sequentially executing the calculation of all the sectors;
(10) extracting machine position parameters; specifically, the last line of data in the turbulence kinetic energy monitoring file and the last line of data in the wind speed vector monitoring file are respectively extracted in a circulating mode, namely the turbulence kinetic energy result and the wind speed result in the calculation result of each sector are extracted, so that the parameters (including the wind speed and the turbulence kinetic energy data) of the wind power plant site are obtained according to the turbulence kinetic energy result and the wind speed result in the calculation result of each sector, and the micro-site selection of the wind power plant at the later stage is facilitated. It should be noted that the machine location point includes each wind turbine, a anemometer tower, a boundary point of a wind farm, and the like, and may be specifically set according to an actual situation.
To sum up, the wind resource calculation method provided by the embodiment of the invention adopts Shell language to compile a script file based on open source software OpenFOAM, and compiles wind direction angles of a plurality of sectors into the script file, sets a parallel calculation and result merging command, improves the calculation efficiency of the sectors, thereby improving the evaluation efficiency of the wind resource; meanwhile, in the sector calculation process, redundant data is also deleted, so that data redundancy is reduced, and data processing time is saved; and the wind speed and turbulence kinetic energy data of each sector machine site can be extracted, so that the micro site selection of the wind power plant in the later period is facilitated, and the method has important practical value.
Corresponding to the above method embodiment, an embodiment of the present invention further provides a wind resource calculation apparatus, which is applied to an electronic device configured with an OpenFOAM, as shown in fig. 3, the apparatus includes an obtaining module 31, a first calculation module 32, and a second calculation module 33, which are connected in sequence; the functions of each module are as follows:
the obtaining module 31 is configured to obtain a basic example file and a parameter set of the wind farm; the wind power plant comprises a plurality of sectors, basic example files are used for representing basic description information of the wind power plant, parameter sets comprise parameters corresponding to each sector, and the parameters comprise at least one of the following parameters: the incoming flow direction, boundary parameters and turbulent kinetic energy dissipation rate;
a first calculation module 32, configured to automatically calculate, based on the basic example file and the parameter corresponding to each sector, a plurality of sectors through OpenFOAM to obtain a calculation result of each sector; wherein the calculation results comprise a wind speed result and a turbulence energy result;
and a second calculating module 33, configured to obtain wind resource distribution information of the wind farm according to the calculation results of the multiple sectors.
The embodiment of the invention provides a wind resource calculation device, which is characterized in that a basic calculation example file and parameters of each sector are preset, and a plurality of sectors of a wind power plant are automatically calculated through OpenFOAM.
In one possible embodiment, the parameters further include: calculating the number of threads in parallel; the first calculating module 32 is further configured to: copying a basic arithmetic example file for each sector, and dividing the copied basic arithmetic example file into a plurality of blocks; and based on OpenFOAM, performing parallel computation on a plurality of blocks according to the number of parallel computation threads corresponding to the sector, and obtaining the computation result of the sector according to the computation sub-result of each block.
In another possible embodiment, before the obtaining of the basic example file of the wind farm, the apparatus further includes: acquiring topographic data of a wind power plant; determining a corresponding calculation domain of the wind power plant according to the topographic data; carrying out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant; acquiring basic setting information of a wind power plant; the basic setting information comprises boundary surface information, boundary conditions and monitoring positions of the wind power plant; and generating a basic operation example file of the wind power plant according to the calculation domain grid and the basic setting information.
In another possible embodiment, the above meshing the computation domains to generate the computation domain mesh corresponding to the wind farm includes: a blockMesh module based on OpenFOAM carries out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant; or carrying out grid division on the calculation domain according to a preset mode to generate an initial calculation domain grid, and converting the initial calculation domain grid into a grid which can be identified by OpenFOAM to obtain a calculation domain grid corresponding to the wind power plant; wherein, the preset mode comprises one of the following modes: gridggen, pointwise, gambit, icemcfd, tetgen, gmesh, and ansys.
In another possible embodiment, the acquiring the topographic data of the wind farm includes: acquiring original topographic data of a wind power plant; carrying out operation processing on the original topographic data to obtain topographic data; wherein the operation process comprises at least one of: cleaning operation and repairing operation.
In another possible embodiment, each sector carries priority information, and the automatically calculating the plurality of sectors by OpenFOAM includes: and according to the priority information of the plurality of sectors, automatically calculating the plurality of sectors in sequence through OpenFOAM.
In another possible embodiment, the automatically calculating the plurality of sectors through OpenFOAM further includes: acquiring a calculation sequence of a plurality of sectors set by a user; and according to the calculation sequence, automatically calculating the plurality of sectors in sequence through OpenFOAM.
The wind resource calculation device provided by the embodiment of the invention has the same technical characteristics as the wind resource calculation method provided by the embodiment, so that the same technical problems can be solved, and the same technical effects can be achieved.
The embodiment of the invention also provides electronic equipment which comprises a processor and a memory, wherein the memory stores machine executable instructions capable of being executed by the processor, and the processor executes the machine executable instructions to realize the wind resource calculation method.
Referring to fig. 4, the electronic device comprises a processor 40 and a memory 41, the memory 41 storing machine executable instructions capable of being executed by the processor 40, the processor 40 executing the machine executable instructions to implement the wind resource calculation method described above.
Further, the electronic device shown in fig. 4 further includes a bus 42 and a communication interface 43, and the processor 40, the communication interface 43 and the memory 41 are connected through the bus 42.
The Memory 41 may include a high-speed Random Access Memory (RAM) and may also include a non-volatile Memory (non-volatile Memory), such as at least one disk Memory. The communication connection between the network element of the system and at least one other network element is realized through at least one communication interface 43 (which may be wired or wireless), and the internet, a wide area network, a local network, a metropolitan area network, etc. may be used. The bus 42 may be an ISA (Industry Standard Architecture) bus, a PCI (Peripheral Component Interconnect) bus, an EISA (Enhanced Industry Standard Architecture) bus, or the like. The above-mentioned bus may be classified into an address bus, a data bus, a control bus, and the like. For ease of illustration, only one double-headed arrow is shown in FIG. 4, but that does not indicate only one bus or one type of bus.
The processor 40 may be an integrated circuit chip having signal processing capabilities. In implementation, the steps of the above method may be performed by integrated logic circuits of hardware or instructions in the form of software in the processor 40. The Processor 40 may be a general-purpose Processor, and includes a Central Processing Unit (CPU), a Network Processor (NP), and the like; the device can also be a Digital Signal Processor (DSP), an Application Specific Integrated Circuit (ASIC), a Field Programmable Gate Array (FPGA) or other Programmable logic device, a discrete Gate or transistor logic device, or a discrete hardware component. The various methods, steps and logic blocks disclosed in the embodiments of the present invention may be implemented or performed. A general purpose processor may be a microprocessor or the processor may be any conventional processor or the like. The steps of the method disclosed in connection with the embodiments of the present invention may be directly implemented by a hardware decoding processor, or implemented by a combination of hardware and software modules in the decoding processor. The software module may be located in ram, flash memory, rom, prom, or eprom, registers, etc. storage media as is well known in the art. The storage medium is located in a memory 41, and the processor 40 reads the information in the memory 41 and completes the steps of the method of the foregoing embodiment in combination with the hardware thereof.
The present embodiments also provide a machine-readable storage medium having stored thereon machine-executable instructions that, when invoked and executed by a processor, cause the processor to implement the wind resource calculation method described above.
The wind resource calculation method, the wind resource calculation device, and the computer program product of the electronic device provided in the embodiments of the present invention include a computer-readable storage medium storing a program code, where instructions included in the program code may be used to execute the method described in the foregoing method embodiments, and specific implementations may refer to the method embodiments and are not described herein again.
It is clear to those skilled in the art that, for convenience and brevity of description, the specific working processes of the system and the apparatus described above may refer to the corresponding processes in the foregoing method embodiments, and are not described herein again.
In addition, in the description of the embodiments of the present invention, unless otherwise explicitly specified or limited, the terms "mounted," "connected," and "connected" are to be construed broadly, e.g., as meaning either a fixed connection, a removable connection, or an integral connection; can be mechanically or electrically connected; they may be connected directly or indirectly through intervening media, or they may be interconnected between two elements. The specific meanings of the above terms in the present invention can be understood in specific cases to those skilled in the art.
The functions, if implemented in the form of software functional units and sold or used as a stand-alone product, may be stored in a non-volatile computer-readable storage medium executable by a processor. Based on such understanding, the technical solution of the present invention may be embodied in the form of a software product, which is stored in a storage medium and includes instructions for causing a computer device (which may be a personal computer, a server, or a network device) to execute all or part of the steps of the method according to the embodiments of the present invention. And the aforementioned storage medium includes: a U-disk, a removable hard disk, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk, and other various media capable of storing program codes.
In the description of the present invention, it should be noted that the terms "center", "upper", "lower", "left", "right", "vertical", "horizontal", "inner", "outer", etc., indicate orientations or positional relationships based on the orientations or positional relationships shown in the drawings, and are only for convenience of description and simplicity of description, but do not indicate or imply that the device or element being referred to must have a particular orientation, be constructed and operated in a particular orientation, and thus, should not be construed as limiting the present invention. Furthermore, the terms "first," "second," and "third" are used for descriptive purposes only and are not to be construed as indicating or implying relative importance.
Finally, it should be noted that: the above-mentioned embodiments are only specific embodiments of the present invention, which are used for illustrating the technical solutions of the present invention and not for limiting the same, and the protection scope of the present invention is not limited thereto, although the present invention is described in detail with reference to the foregoing embodiments, those skilled in the art should understand that: any person skilled in the art can modify or easily conceive the technical solutions described in the foregoing embodiments or equivalent substitutes for some technical features within the technical scope of the present disclosure; such modifications, changes or substitutions do not depart from the spirit and scope of the embodiments of the present invention, and they should be construed as being included therein. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.

Claims (10)

1. A wind resource calculation method is applied to an electronic device configured with OpenFOAM, and comprises the following steps:
acquiring a basic example file and a parameter set of a wind power plant; the wind power plant comprises a plurality of sectors, the basic operation example file is used for characterizing basic description information of the wind power plant, the parameter set comprises parameters corresponding to each sector, and the parameters comprise at least one of the following parameters: the incoming flow direction, boundary parameters and turbulent kinetic energy dissipation rate;
based on the basic example file and the parameters corresponding to each sector, automatically calculating a plurality of sectors through the OpenFOAM to obtain a calculation result of each sector; wherein the calculation results include a wind speed result and a turbulence energy result;
and obtaining wind resource distribution information of the wind power plant according to the calculation results of the plurality of sectors.
2. The wind resource calculation method of claim 1, wherein the parameters further comprise: calculating the number of threads in parallel; the step of automatically calculating a plurality of the sectors through the OpenFOAM based on the basic example file and the parameter corresponding to each of the sectors includes:
for each sector, copying the basic algorithm file, and dividing the copied basic algorithm file into a plurality of blocks; and based on the OpenFOAM, performing parallel computation on the blocks according to the number of parallel computation threads corresponding to the sector, and obtaining the computation result of the sector according to the computation sub-result of each block.
3. The wind resource calculation method of claim 1, wherein the step of obtaining a base algorithm file for a wind farm is preceded by the method further comprising:
acquiring topographic data of the wind power plant;
determining a calculation domain corresponding to the wind power plant according to the topographic data;
carrying out grid division on the calculation domain to generate a calculation domain grid corresponding to the wind power plant;
acquiring basic setting information of the wind power plant; the basic setting information comprises boundary surface information, boundary conditions and monitoring positions of the wind power plant;
and generating a basic operation example file of the wind power plant according to the calculation domain grid and the basic setting information.
4. The wind resource calculation method according to claim 3, wherein the step of performing grid division on the calculation domain to generate the calculation domain grid corresponding to the wind farm includes:
carrying out grid division on the calculation domain based on a blockMesh module of the OpenFOAM to generate a calculation domain grid corresponding to the wind power plant; alternatively, the first and second electrodes may be,
carrying out grid division on the calculation domain according to a preset mode to generate an initial calculation domain grid, and converting the initial calculation domain grid into a grid which can be identified by the OpenFOAM to obtain a calculation domain grid corresponding to the wind power plant; wherein the preset mode comprises one of the following modes: gridggen, pointwise, gambit, icemcfd, tetgen, gmesh, and ansys.
5. The wind resource calculation method of claim 3, wherein said step of obtaining terrain data for said wind farm comprises:
acquiring original topographic data of the wind power plant;
carrying out operation processing on the original topographic data to obtain the topographic data; wherein the operation process comprises at least one of: cleaning operation and repairing operation.
6. The wind resource calculation method according to claim 1, wherein each of the sectors carries priority information, and the step of automatically calculating the plurality of sectors by the OpenFOAM comprises:
and according to the priority information of the plurality of sectors, automatically calculating the plurality of sectors in sequence through the OpenFOAM.
7. The wind resource calculation method of claim 1, wherein the step of automatically calculating the plurality of sectors by the OpenFOAM further comprises:
acquiring a calculation sequence of a plurality of sectors set by a user;
and according to the calculation sequence, automatically calculating a plurality of sectors in sequence through the OpenFOAM.
8. A wind resource computing apparatus applied to an OpenFOAM-configured electronic device, the apparatus comprising:
the acquisition module is used for acquiring basic example files and parameter sets of the wind power plant; the wind power plant comprises a plurality of sectors, the basic operation example file is used for characterizing basic description information of the wind power plant, the parameter set comprises parameters corresponding to each sector, and the parameters comprise at least one of the following parameters: the incoming flow direction, boundary parameters and turbulent kinetic energy dissipation rate;
a first calculation module, configured to automatically calculate, based on the basic example file and a parameter corresponding to each sector, a plurality of sectors through the OpenFOAM to obtain a calculation result of each sector; wherein the calculation results include a wind speed result and a turbulence energy result;
and the second calculation module is used for obtaining wind resource distribution information of the wind power plant according to the calculation results of the plurality of sectors.
9. An electronic device comprising a memory, a processor and a computer program stored on the memory and executable on the processor, wherein the steps of the method of any of the preceding claims 1-7 are implemented when the computer program is executed by the processor.
10. A computer-readable storage medium, characterized in that a computer program is stored on the computer-readable storage medium, which computer program, when being executed by a processor, performs the steps of the method of any of the preceding claims 1-7.
CN202111344542.7A 2021-11-12 2021-11-12 Wind resource calculation method and device and electronic equipment Pending CN114169576A (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202111344542.7A CN114169576A (en) 2021-11-12 2021-11-12 Wind resource calculation method and device and electronic equipment

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202111344542.7A CN114169576A (en) 2021-11-12 2021-11-12 Wind resource calculation method and device and electronic equipment

Publications (1)

Publication Number Publication Date
CN114169576A true CN114169576A (en) 2022-03-11

Family

ID=80478941

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202111344542.7A Pending CN114169576A (en) 2021-11-12 2021-11-12 Wind resource calculation method and device and electronic equipment

Country Status (1)

Country Link
CN (1) CN114169576A (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117151352A (en) * 2023-11-01 2023-12-01 北京大学长沙计算与数字经济研究院 Wind resource evaluation method, system, computer storage medium and terminal equipment

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104331621A (en) * 2014-11-05 2015-02-04 中国大唐集团新能源股份有限公司 Wind resource computation method
CN108563867A (en) * 2018-04-13 2018-09-21 华中科技大学 A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM
CN109598064A (en) * 2018-12-04 2019-04-09 华润电力技术研究院有限公司 A kind of wind-resources zoning optimization method based on OpenFOAM

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN104331621A (en) * 2014-11-05 2015-02-04 中国大唐集团新能源股份有限公司 Wind resource computation method
CN108563867A (en) * 2018-04-13 2018-09-21 华中科技大学 A method of WRF and CFD coupled simulation wind fields are realized based on OpenFOAM
CN109598064A (en) * 2018-12-04 2019-04-09 华润电力技术研究院有限公司 A kind of wind-resources zoning optimization method based on OpenFOAM

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
叶昭良;闫姝;史绍平;曾崇济;陈新明;张波;马俊祥;: "基于Openfoam的怀来复杂地形风电场流动特性研究", 可再生能源, no. 04, 16 April 2020 (2020-04-16), pages 471 - 476 *
张德胜: "基于OpenFOAM的风电场流场的数值计算", 中国优秀硕士学位论文全文数据库 工程科技Ⅱ辑, no. 2, 15 December 2013 (2013-12-15), pages 1 - 10 *

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN117151352A (en) * 2023-11-01 2023-12-01 北京大学长沙计算与数字经济研究院 Wind resource evaluation method, system, computer storage medium and terminal equipment
CN117151352B (en) * 2023-11-01 2024-01-30 北京大学长沙计算与数字经济研究院 Wind resource evaluation method, system, computer storage medium and terminal equipment

Similar Documents

Publication Publication Date Title
CN107346357B (en) Offshore wind turbine fatigue analysis system based on integral coupling model
Connolly et al. Development of a computer program to locate potential sites for pumped hydroelectric energy storage
CN106650086A (en) Indirect air cooling system numerical simulation platform based on FLUENT software
CN114169576A (en) Wind resource calculation method and device and electronic equipment
CN114091146A (en) Method and device for block form design based on wind environment simulation and storage medium
CN113937793A (en) Stability analysis method based on impedance segmentation reduced model zero point identification
CN109635317B (en) CFD (computational fluid dynamics) simulation terrain topology method for high-altitude area
CN109684723A (en) A kind of two-dimensional structure internal acoustic method for analyzing performance
CN112949227A (en) Method, system, equipment and storage medium suitable for determining turbulence intensity of wind field in complex terrain
CN113158598B (en) Foundation type wind power plant flow field CFD partition calculation method, device, equipment and medium
CN115375867A (en) Method, system, equipment and medium for calculating geothermal resource quantity by using grid model
CN112967391A (en) Topographic map determining method and device and electronic equipment
CN117195016B (en) Sewage treatment mode determining method and device, computer equipment and storage medium
CN117151352B (en) Wind resource evaluation method, system, computer storage medium and terminal equipment
CN111291464A (en) Dynamic equivalence method and device for power system
CN118211521B (en) Wave eliminating effect simulation method and relevant equipment for wave eliminating hole wharf
CN118135138B (en) Terrain modeling method for complex terrain wind field numerical simulation
CN117725807B (en) Method, system and equipment for calculating heterogeneous reservoir fracture stress based on boundary element
CN111506987B (en) Method and system for lean parameter management and stable intelligent calculation of generator set
CN117706045B (en) Combined control method and system for realizing atmospheric ozone monitoring equipment based on Internet of things
CN116738641B (en) Running water circulation method and system for sewage river network
CN110648399A (en) Wind power plant terrain entity modeling method and device based on Surfer software
CN117454801A (en) Method for selecting dynamic capacity-increasing monitoring line grade of overhead transmission line
CN107247835B (en) Library capacity calculation precision verification method
CN117787007A (en) Method and device for determining radiator parameters, computer equipment and storage medium

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