Method and device for generating heliostat field arrangement scheme of mountain terrain

Abstract

The invention provides a method and device for generating a heliostat field arrangement scheme of mountain topography, and relates to the technical field of solar thermal power generation engineering, and the method comprises the steps: constructing an irregular triangular network model of an initial heliostat field arrangement range through employing a Delou triangulation method; establishing a three-dimensional terrain curved surface model according to the irregular triangular net model; regularized grid arrangement is implemented on the three-dimensional terrain curved surface model; performing interpolation analysis on each grid in the regularized grid arrangement by adopting a cubic spline interpolation function, and determining an initial coordinate and a height of the heliostat field; and carrying out multi-objective optimization on the initial coordinates and the height by utilizing a projection area constraint condition and a self-adaptive gravitational search algorithm to obtain a final heliostat field arrangement scheme. According to the method, the technical defects of insufficient terrain adaptability, imperfect shielding effect control, insufficient system energy efficiency optimization and the like in the implementation process of the existing technical scheme can be effectively solved.

Description

Background Art

[0002] In recent years, with the rapid development of renewable energy development and utilization technologies, Concentrated Solar Power (CSP) technology has been promoted and applied globally. In conventional technical solutions, heliostat fields are usually located in flat terrain areas with a slope of no more than 5% to achieve a regular layout of the field. However, in practical engineering applications, there are a large number of mountainous terrain areas with slopes between 5% and 20%. The development and utilization of such terrain face significant technical challenges.

[0003] In existing technical solutions, there are mainly two types of methods for arranging heliostat fields in mountainous terrain: The first type of method directly applies flat terrain field layout algorithms (such as the Campo algorithm, genetic algorithm, etc.) to the mountainous environment, and adjusts the horizontal inclination angle of the heliostat base to maintain the installation angle of the mirror surface relative to the horizontal plane. Although this method inherits the mature technical framework of existing algorithms, it has the defect of insufficient terrain adaptability, specifically manifested as: (a) it cannot effectively compensate for the deviation of the mirror normal vector caused by terrain slope changes; (b) it causes an increase in the mirror surface occlusion effect; (c) the overall light field heat collection efficiency decreases. The second type of method adopts a terrain zoning processing strategy, divides the field area into several sub-areas for independent layout and local optimization. Although this scheme improves the local terrain adaptability, it has the following technical limitations: (i) mutual mirror occlusion is likely to occur in the transition area between sub-areas; (ii) it is difficult to achieve the optimization of global optical efficiency; (iii) the improvement effect of system-level energy efficiency is limited.

[0004] Obviously, existing technical solutions all have technical defects such as insufficient terrain adaptability, imperfect control of occlusion effects, and insufficient optimization of system energy efficiency during implementation. There is an urgent need to develop a new method for arranging heliostat fields in mountainous terrain to solve these problems. Summary of the Invention

[0005] To overcome the problems existing in the related art, the present invention provides a method and device for generating a layout scheme of a heliostat field in mountainous terrain.

[0006] According to the first aspect of the embodiments of the present invention, there is provided a method for generating a layout scheme of a heliostat field in mountainous terrain, the method including: Using the Delaunay triangulation method, constructing an irregular triangular mesh model of the initial layout range of the heliostat field; Establishing a three-dimensional terrain surface model based on the irregular triangular mesh model; Implementing a regular grid layout on the three-dimensional terrain surface model; Using a cubic spline interpolation function to perform interpolation analysis on each grid in the regular grid layout to determine the initial coordinates and height of the heliostat field; The initial coordinates and height are multi-objectively optimized by using the projected area constraint condition and the adaptive gravitational search algorithm to obtain the final layout scheme of the heliostat field.

[0007] In some exemplary embodiments of the present invention, based on the foregoing scheme, interpolation analysis is performed on each grid in the regularized grid layout by using a cubic spline interpolation function, and determining the initial coordinates and height of the heliostat field includes: Calculating the interpolation value of each grid by using a cubic spline interpolation function; Based on the interpolation values of each grid, estimating the initial coordinates and height of the heliostat field corresponding to each grid within the boundary conditions.

[0008] In some exemplary embodiments of the present invention, based on the foregoing scheme, the boundary conditions include: The value of the cubic spline interpolation function defined on the interval at is equal to the value of the cubic spline interpolation function defined on the interval at ; The first derivative value of the cubic spline interpolation function defined on the interval at is equal to the first derivative value of the cubic spline interpolation function defined on the interval at ; and The second derivative value of the cubic spline interpolation function defined on the interval at is equal to the second derivative value of the cubic spline interpolation function defined on the interval at , where represents the coordinate value of a vertex of the th grid, represents the coordinate value of the th grid, represents the coordinate value of the th grid.

[0009] In some exemplary embodiments of the present invention, based on the foregoing scheme, the initial coordinates and height are multi-objectively optimized by using the projected area constraint condition and the adaptive gravitational search algorithm to obtain the final layout scheme of the heliostat field, including: Selecting the solar parameters at a typical moment as reference data, and calculating the projected area of each heliostat field based on the initial coordinates and height; Taking the maximum projected area coverage, the least projection overlap, and no obstruction to the mechanical movement of the heliostat as the objectives, and using the adaptive gravitational search algorithm to perform multi-objective optimization on the initial coordinates and height to obtain the final layout scheme of the heliostat field.

[0010] In some exemplary embodiments of the present invention, based on the foregoing solution, solar parameters at typical moments are selected as reference data, and calculating the projected area of each heliostat field based on the initial coordinates includes: Calculating the vertical height where the center of the current heliostat field is located according to the initial coordinates and height; Based on the vertical height where the center of the heliostat field is located, calculating the distance between the center of the current heliostat field and the receiver; Based on the distance and the actual area of the current heliostat field, calculating the projected area of the current heliostat field.

[0011] In some exemplary embodiments of the present invention, based on the foregoing solution, with the goals of maximizing the projected area coverage, minimizing the projection overlap, and having no obstruction to the mechanical movement of the heliostats, the initial coordinates and height are optimized by using an adaptive gravitational search algorithm, and obtaining the final heliostat field layout plan includes: Setting the number of heliostats in the initial mirror field, encoding the initial mirror field to generate a matrix with N rows and R columns, where N is the number of mirror field layout plans and R is the number of rows of heliostats; Evaluating N mirror field layout plans by using a multi-objective fitness function including the projected area coverage, the projection overlap rate, and the mechanical movement constraint; Calculating the gravitational force according to the Euclidean distance between particles, calculating the acceleration by combining random parameters and an adaptive gravitational constant, and iteratively optimizing the layout plan through a velocity-position update formula; Continuously updating the global optimal solution until the preset convergence accuracy is satisfied, and outputting the heliostat coordinate matrix corresponding to the highest fitness value.

[0012] In some exemplary embodiments of the present invention, based on the foregoing solution, the number of heliostats in the initial mirror field is greater than or equal to 1.1 times the number of heliostats in the expected heliostat field.

[0013] According to the second aspect of the embodiments of the present invention, there is provided a generating device for a heliostat field layout plan for mountainous terrain, including: A first model construction module, configured to construct an irregular triangular mesh model of the initial heliostat field layout range by using the Delaunay triangulation method; A second model construction module, configured to establish a three-dimensional terrain surface model according to the irregular triangular mesh model; A grid layout module, configured to perform regular grid layout on the three-dimensional terrain surface model; An interpolation analysis module, configured to perform interpolation analysis on each grid in the regular grid layout by using a cubic spline interpolation function to determine the initial coordinates and height of the heliostat field; A scheme generation module, configured to perform multi-objective optimization on the initial coordinates and height by using the projected area constraint condition and the adaptive gravitational search algorithm to obtain the final heliostat field layout scheme.

[0014] According to a third aspect of an embodiment of the present invention, there is provided an electronic device, including: a processor; and a memory, where computer-readable instructions are stored on the memory, and when the computer-readable instructions are executed by the processor, the method in the first aspect is implemented.

[0015] According to a fourth aspect of an embodiment of the present invention, there is provided a computer-readable storage medium, on which a computer program is stored, and when the computer program is executed by a processor, the method in the first aspect is implemented.

[0016] The technical solutions provided by the embodiments of the present invention may include the following beneficial effects: In the embodiments of the present invention, through the Triangulated Irregular Network (TIN) model constructed by Delaunay triangulation, the elevation mutation characteristics of the mountain terrain (slope change ≥ 25°) can be effectively captured, and the complex and changeable terrain of the mountain can be accurately fitted.

[0017] Implementing a regular grid layout on the three-dimensional terrain surface model and using the cubic spline interpolation function to analyze and determine the initial coordinates and height can reduce the irregularity of the layout. Compared with random or empirical layouts, it can make more efficient use of the mountain space, reduce the mutual occlusion between heliostats, and improve the overall daylighting efficiency.

[0018] Performing multi-objective optimization by using the projected area constraint condition and the adaptive gravitational search algorithm can, while considering terrain factors, comprehensively balance multiple key performance indicators such as the daylighting area and energy collection efficiency of the heliostats, so that the final heliostat field layout scheme achieves the best performance in a complex mountain environment, improves the collection and conversion efficiency of solar energy, and reduces the construction and operation costs.

[0019] It should be understood that the above general description and subsequent detailed description are only exemplary and explanatory, and cannot limit the present invention. BRIEF DESCRIPTION OF THE DRAWINGS

[0020] The accompanying drawings herein are incorporated into the specification and constitute a part of the present invention, showing embodiments consistent with the present invention, and are used together with the specification to explain the principles of the present invention.

[0021] Figure 1 A schematic diagram of a system architecture showing an exemplary application environment of a method and apparatus to which embodiments of the present invention can be applied; Figure 2 Schematically showing a flowchart of a method according to some embodiments of the present invention; Figure 3 A schematic diagram showing a device according to some embodiments of the present invention; Figure 4 A schematic diagram showing the structure of a computer system of an electronic device according to some embodiments of the present invention; Figure 5 A schematic diagram showing a computer-readable storage medium according to some embodiments of the present invention. Detailed implementation manners

[0022] Here, exemplary embodiments will be described in detail, and examples thereof are shown in the drawings. When the following description refers to the drawings, unless otherwise indicated, the same numbers in different drawings represent the same or similar elements. The implementation manners described in the following exemplary embodiments do not represent all implementation manners consistent with the present invention. On the contrary, they are merely examples of devices and methods consistent with some aspects of the present invention as detailed in the appended claims.

[0023] The terms used in the present invention are for the purpose of describing specific embodiments only and are not intended to limit the present invention. The singular forms "a", "said", and "the" used in the present invention and the appended claims are also intended to include the plural forms unless the context clearly indicates otherwise. It should also be understood that the term "and / or" used herein refers to and includes any or all possible combinations of one or more of the associated listed items.

[0024] It should be understood that although the terms first, second, third, etc. may be used in the present invention to describe various information, such information should not be limited to these terms. These terms are only used to distinguish the same type of information from each other. For example, without departing from the scope of the present invention, the first information may also be referred to as the second information, and similarly, the second information may also be referred to as the first information. Depending on the context, the word "if" as used herein may be interpreted as "when" or "while" or "in response to determining".

[0025] Figure 1 A schematic diagram showing the system architecture of an exemplary application environment of a method and apparatus for generating a heliostat field layout scheme for a mountainous terrain to which embodiments of the present invention can be applied.

[0026] Such as Figure 1As shown in the figure, the system architecture 100 may include one or more terminal devices such as a desktop computer 101, a portable computer 102, a smart phone 103, etc., a network 104, and a server 105. The network 104 is used to provide a medium for communication links between the terminal devices and the server 105. The network 104 may include various connection types, such as wired, wireless communication links, or fiber optic cables, etc. The terminal devices may be various electronic devices with data processing functions, and a display screen is provided on the electronic device, and the display screen is used to display the generation result of the heliostat field layout scheme of the mountain terrain to the user, including but not limited to the above-mentioned desktop computer, portable computer, smart phone, etc. It should be understood that Figure 1 the numbers of the terminal devices, networks, and servers in the figure are only illustrative. According to the implementation requirements, there may be any number of terminal devices, networks, and servers. For example, the server 105 may be a sub-server cluster composed of multiple sub-servers, etc.

[0027] The method for generating the heliostat field layout scheme of the mountain terrain provided by the embodiments of the present invention can generally be executed by a terminal device. Correspondingly, the device for generating the heliostat field layout scheme of the mountain terrain is generally set in the terminal device. However, it is easy for those skilled in the art to understand that the method for generating the heliostat field layout scheme of the mountain terrain provided by the embodiments of the present invention can also be executed by the server 105. Correspondingly, the device for generating the heliostat field layout scheme of the mountain terrain can also be set in the server 105. No special limitation is made in this exemplary embodiment.

[0028] In addition, it should be understood that the method for generating the heliostat field layout scheme of the mountain terrain in the embodiments of the present invention can be configured as a software module. In some implementation scenarios, the method for generating the heliostat field layout scheme of the mountain terrain of the present invention can be deployed separately to adapt to the heliostat field layout of different mountain terrains. In other implementation scenarios, the method for generating the heliostat field layout scheme of the mountain terrain of the present invention can be deployed in other software as a function module of the software, such as being deployed in the analysis software of the mountain terrain. The present invention does not make special restrictions on the application manner of the method for generating the heliostat field layout scheme of the mountain terrain.

[0029] Next, the embodiments of the present invention will be described in detail.

[0030] The heliostat field is the core component of a Concentrated Solar Power (CSP) system, specifically referring to an array system composed of several planar reflectors with a two-axis rotation mechanism. Each mirror body reflects solar radiation to a fixed-position collector (usually a central receiver tower) through real-time angle adjustment. Its core function is to achieve the spatial directional transmission and energy convergence of solar radiation energy.

[0031] like Figure 2 As shown, Figure 2 1 is a flow chart of a method for generating a heliostat field arrangement scheme for mountainous terrain according to an exemplary embodiment of the present invention, comprising the following steps: S210: Using the Delaunay triangulation method, an irregular triangle network model of the initial heliostat field layout is constructed; S220: Establishing a three-dimensional terrain surface model according to the irregular triangle network model; S230: Implementing regularized grid arrangement on the three-dimensional terrain surface model; S240: performing interpolation analysis on each grid in the regularized grid arrangement using a cubic spline interpolation function to determine the initial coordinates and height of the heliostat field; S250: Utilizing the projection area constraint condition and the adaptive gravitational search algorithm to perform multi-objective optimization on the initial coordinates and height to obtain a final heliostat field layout plan.

[0032] In S210, an irregular triangle network model of the initial heliostat field layout range is constructed using the Delaunay triangulation method.

[0033] Delaunay triangulation is an algorithm that constructs a discrete point set into a triangular network. Its characteristic is that the circumscribed circle of any triangle does not contain other points.

[0034] The initial heliostat field layout range can be determined by taking into account multiple factors such as terrain factors, lighting conditions, surrounding environment and construction costs.

[0035] Topographic factors generally need to consider the overall shape, slope, slope direction and terrain undulation of the mountain. For example, areas with steep slopes (generally considered to be unfavorable for heliostat installation and maintenance if the slope exceeds a certain angle, such as 30°), geologically unstable areas such as landslides, and areas with water accumulation and severe obstruction at the bottom of the valley should be avoided. Areas with relatively flat, open and well-lit terrain should be given priority, such as gentle slopes facing south, which can ensure that the heliostat receives sufficient sunlight to the greatest extent.

[0036] The light condition factor is generally determined by analyzing the local solar radiation data, sunshine duration and other light information for many years to clarify the sun's trajectory and light intensity distribution. The layout range should try to cover areas with long sunshine hours and high solar radiation intensity throughout the year to improve the solar energy collection efficiency of the heliostat.

[0037] Peripheral environmental factors involve considering the relationship between the heliostat field and surrounding infrastructure and obstacles. Stay away from large buildings, trees and other objects that may cause shading, and at the same time, it should be convenient to connect to infrastructure such as power transmission lines and roads to reduce construction and operation costs. For example, being close to an existing substation or power transmission line can reduce the cost of laying power transmission lines; being near a traffic road is convenient for equipment transportation and later maintenance.

[0038] Construction cost factors generally include land acquisition or lease costs, site leveling costs, etc. On the premise of meeting the requirements of sunlight and terrain, give priority to areas with lower land costs and less site leveling work to control the overall investment cost of the project.

[0039] When constructing a Triangulated Irregular Network (TIN) model, by connecting the discrete points within the layout range of the heliostat field (these points may be obtained based on topographic surveys, such as topographic feature points, control points, etc., and feature points can be mountaintops, valleys, saddles, etc.) according to the rules of Delaunay triangulation to form triangles, thus forming an irregular triangular network that can reflect the original characteristics of the site.

[0040] For complex and irregular mountain terrains, due to their undulating terrain and diverse slopes, traditional regular grid division methods are difficult to accurately describe terrain features. And the TIN model constructed using the Delaunay triangulation method, since this TIN model is connected based on the actual discrete points of the site and can retain more terrain information, can well adapt to the complexity of mountain terrains, integrate details such as the ups and downs and terrain changes of the mountains into the triangular network, and provide basic data support for subsequent steps such as establishing a three-dimensional terrain surface model, reducing the problem of unreasonable heliostat layout caused by inaccurate terrain modeling.

[0041] In S220, establish a three-dimensional terrain surface model according to the irregular triangular network model.

[0042] In some embodiments, the contour lines of the irregular triangular network model can be extracted, and a three-dimensional terrain curve model can be established using the contour lines.

[0043] For example, for each triangle in the TIN model, according to the coordinates and elevation values of its three vertices, a plane equation can be defined to represent the plane where the triangle is located. In three-dimensional space, the plane equation can generally be expressed in the form of By solving the coefficients of the equation through the triangle vertex coordinates 、 、 、 Combining the plane equations of all triangles can construct a three-dimensional terrain surface model of the entire mirror field area. Each triangle, as a component of the surface, jointly constitutes a continuous surface that approximates the actual terrain.

[0044] The contour lines of the irregular triangular network model can be extracted using algorithms such as contour tracing algorithms, data interpolation, and fitting methods.

[0045] When using the contour tracing algorithm to extract the contour lines of the irregular triangular network model, a starting triangle can be randomly selected from the TIN model. According to the given contour elevation value, judge the intersection situation of the three sides of the triangle with the contour line. If the elevations of the two endpoints of a side are on both sides of the contour elevation, then this side intersects the contour line, and the coordinates of the intersection point are calculated by linear interpolation. Then continue to trace along the triangles adjacent to this triangle and likely to intersect the contour line, repeating the above process until a complete contour line is formed or the tracing termination condition is reached. Repeating this process continuously can extract the contour lines of different elevations, thus obtaining the contour map of the mirror field.

[0046] When using the data interpolation and fitting method to extract the contour lines of the irregular triangular network model, the elevation of any point inside the triangle can also be calculated by the interpolation method. Based on the dense elevation point data obtained from these interpolations, a fitting algorithm such as spline curve fitting is used to connect the points with the same elevation to form smooth contour lines. This method can improve the accuracy and smoothness of the contour lines to a certain extent and better reflect the detailed features of the mirror field terrain.

[0047] In S230, a regular grid layout is implemented on the three-dimensional terrain surface model.

[0048] A regular grid refers to a grid structure evenly arranged in two-dimensional or three-dimensional space at fixed row spacing, column spacing (or layer spacing), such as a rectangular grid (cornfield type) or a hexagonal grid.

[0049] In some embodiments, when implementing a regular grid layout on the three-dimensional terrain surface model, tools such as GIS software (ArcGIS), 3D modeling tools (Blender), or mathematical methods (differential geometric parameterization) can be used to map the three-dimensional terrain surface model to a two-dimensional parameter space (such as UV unwrapping) to facilitate the generation of a regular grid. Then, methods such as surface subdivision (Subdivision) and projection interpolation (such as radial basis functions) are used to generate a regular grid in the parameter space, and then the grid is projected back to the three-dimensional terrain surface model through an inverse mapping.

[0050] In some other embodiments, when implementing a regularized grid layout on the three-dimensional terrain surface model, the grid parameters can be determined first according to the three-dimensional terrain surface model; then, according to the grid parameters, regularly arranged grid nodes can be generated in three-dimensional space; and then the grid nodes are connected.

[0051] Here, the grid parameters include grid resolution, grid range, and grid type. The grid resolution can be determined according to the accuracy requirements of the three-dimensional terrain surface model and the actual application needs. For example, when performing a detailed mirror field analysis, a smaller grid spacing (such as 1 m × 1 m) may be required to accurately capture terrain changes; while when performing a macroscopic terrain display, a larger grid spacing (such as 10 m × 10 m) can be selected. The grid range is used to clarify the coverage range of the regularized grid on the three-dimensional terrain surface model, and generally should cover the entire area to be analyzed. The starting and ending positions of the grid can be determined according to the boundary coordinates of the model. Common regularized grid types include square grids and rectangular grids, etc. Generally, a suitable grid type is selected according to the actual situation and the convenience of subsequent processing. The square grid has the same resolution in all directions, and the calculation and analysis are relatively simple.

[0052] According to the grid parameters, regularly arranged grid nodes are generated in three-dimensional space. Here, taking a square grid as an example, on the horizontal plane, starting from the starting point according to the set grid spacing, nodes are successively generated in the axis and axis directions to form a two-dimensional grid node array.

[0053] Finally, the adjacent grid nodes are connected in a certain order to form a regular grid.

[0054] The regularized grid has a unified structure and a fixed spacing, making the organization and storage of data simpler and more regular. When performing data processing, the computer can quickly access and process the data of each grid point according to a fixed pattern, greatly improving the efficiency of data processing, reducing the calculation time and resource consumption. And due to the regularity of the grid, some efficient data compression algorithms can be used to compress and store the grid data. For example, for elevation data, only the relative elevation change values of each grid point need to be stored, rather than storing the absolute coordinates and elevations of each point, thus saving a large amount of storage space and facilitating data transmission and management.

[0055] In addition, by reasonably setting the grid spacing and selecting the interpolation method, the error propagation in the analysis process can be effectively controlled to ensure the accuracy and reliability of the analysis results. For example, in a mirror field area with complex terrain, the analysis accuracy can be improved and the error can be reduced by reducing the grid spacing.

[0056] In S240, a cubic spline interpolation function is used to perform interpolation analysis on each grid in the regularized grid layout to determine the initial coordinates and height of the heliostat field.

[0057] In some embodiments, for each heliostat in the heliostat field, its initial coordinates are determined according to its position in the regularized grid. Assuming that the heliostat is located within a certain grid, its coordinates can be determined by the node coordinates of the grid and the relative position of the heliostat within the grid. For example, if the heliostat is located at the center of the grid, its coordinates can be taken as the average of the coordinates of the four vertices of the grid.

[0058] Then, the constructed cubic spline interpolation function is used to calculate the elevation value at the heliostat position, that is, the height of the heliostat. Substituting the coordinates of the heliostat into the corresponding cubic spline interpolation function, the value obtained is the height of this point.

[0059] In other embodiments, using a cubic spline interpolation function to perform interpolation analysis on each grid in the regularized grid layout to determine the initial coordinates and height of the heliostat field includes: Using a cubic spline interpolation function to calculate the interpolation value of each grid; Based on the interpolation values of each grid, estimate the initial coordinates and height of the heliostat field corresponding to each grid within the boundary conditions.

[0060] This process can be specifically expressed as: is a series of functions , ,..., 's set, where each function corresponds to the interpolation value of each grid in the regularized grid layout of the TIN model, is the independent variable, representing the interval coordinates of the interpolation. For each small interval , a cubic spline interpolation function is defined, and it is required that this function is three times differentiable, that is, . Here represents that this function has a continuous third-order derivative on the interval .

[0061] The boundary conditions include: Among them, represents 's value at , represents the cubic spline interpolation function defined on the interval , represents the The coordinate value of a vertex of a grid denotes the value at ; denotes a cubic spline interpolation function defined on the interval ; and respectively denote the first derivative and the second derivative of and respectively denote the first derivative and the second derivative of

[0062] That is to say: The value of the cubic spline interpolation function defined on the interval at is equal to the value of the cubic spline interpolation function defined on the interval at ; The first derivative value of the cubic spline interpolation function defined on the interval at is equal to the first derivative value of the cubic spline interpolation function defined on the interval at ; and The second derivative value of the cubic spline interpolation function defined on the interval at is equal to the second derivative value of the cubic spline interpolation function defined on the interval at , where denotes the coordinate value of a vertex of the th grid denotes the coordinate value of the th grid denotes the coordinate value of the th grid

[0063] In S250, the initial coordinates and height are optimized by using the projected area constraint condition and the adaptive gravitational search algorithm to obtain the final heliostat field layout scheme

[0064] The projected area constraint refers to setting certain limitations on the projected area of a heliostat field in a specific direction (such as the horizontal direction or a direction related to the sun's rays). This constraint may be based on various practical requirements, such as land use restrictions and avoiding excessive occlusion between heliostats. By restricting the projected area, the layout of heliostats can be reasonably planned to ensure the efficient operation of the mirror field within a limited space while reducing unnecessary land occupation. As an important constraint, it can exclude layout schemes that do not meet the projected area requirements during the optimization process, narrow the range of feasible solutions, and guide the algorithm to search for the optimal solution in the direction that meets the actual site and usage requirements.

[0065] The Adaptive Gravitational Search Algorithm (AGSA), like the basic gravitational search algorithm, has its core idea originating from Newton's law of universal gravitation and laws of motion. In the algorithm, each solution in the search space is simulated as a particle with mass, and there is a gravitational force between particles. Particles with larger mass have a greater gravitational force on other particles and will attract other particles to approach it, thereby guiding the search process towards a better solution.

[0066] Advantages of AGSA in the optimization of heliostat fields: In the optimization of heliostat field layout, the combination of the coordinates and heights of heliostats forms a complex search space. The adaptive nature of AGSA enables it to better handle this complex situation. It can automatically adjust the search strategy at different search stages, being able to quickly search for areas where optimal solutions may exist in a large range and conduct detailed searches locally to find a more accurate combination of heliostat coordinates and heights. Compared with traditional fixed-parameter algorithms, it can more effectively explore better layout schemes.

[0067] By adaptively adjusting parameters, AGSA can more quickly identify "high-quality regions" in the search space and guide particles to gather in these regions, thereby accelerating the convergence speed of the algorithm. In the optimization of heliostat field layout, this means that a better solution that meets the multi-objective optimization requirements can be found in a shorter time, saving computational resources and time costs and improving the design efficiency.

[0068] Since it can dynamically adjust the search strategy according to the search situation, AGSA is not easily trapped in local optimal solutions. That is to say, in the layout of heliostat fields, there may be multiple locally optimal layout schemes, but only one or a few are globally optimal. AGSA can have a greater chance of jumping out of the local optimal trap by adaptively adjusting the movement and search range of particles, finding the true global optimal or near-global optimal heliostat field layout scheme, thereby improving the overall performance and benefits of the mirror field.

[0069] When working in conjunction with the projected area constraint condition, AGSA continuously checks whether the heliostat field layout plan represented by the particles meets the projected area requirements during the search process. For the plans that do not meet the constraints, an adaptive adjustment mechanism is used to guide the particles to move in the direction that meets the constraints, ensuring that the final obtained plan meets the actual requirements in terms of the projected area.

[0070] Under the framework of multi-objective optimization, AGSA adaptively adjusts the search direction and intensity of the particles according to the importance of each objective and the current optimization situation to balance the relationship between different objectives. For example, when there is a conflict between the two objectives of light energy collection efficiency and cost, AGSA will find a balance point that can optimize both as much as possible during the search process according to the preset weights or dynamically changing priorities, and finally obtain a heliostat field layout plan with the optimal comprehensive performance, meeting the various requirements for the heliostat field in actual engineering.

[0071] In some embodiments, this process may include: Selecting the solar parameters at typical moments as reference data, and calculating the projected area of each heliostat field based on the initial coordinates and heights; Using the adaptive gravitational search algorithm with the goals of maximizing the projected area coverage, minimizing the projection overlap, and having unobstructed mechanical movement of the heliostats, perform multi-objective optimization on the initial coordinates and heights to obtain the final heliostat field layout plan.

[0072] The selection of typical moments usually needs to consider the representativeness of parameters such as the solar altitude angle and azimuth angle and their importance to the operation of the heliostat field. For example, at noon on the winter solstice and summer solstice, the solar altitude angles reach the minimum and relatively large values of the year respectively, which are crucial for evaluating the performance of the heliostat field under extreme seasonal conditions; at sunrise and sunset on the spring equinox and autumn equinox, the solar azimuth angle changes greatly, which can be used to analyze the situation of the heliostat field under different lighting directions. In addition, the moment with the strongest local solar radiation intensity may also be selected as a typical moment.

[0073] Solar parameters mainly include the solar altitude angle , the solar azimuth angle , etc. These parameters can be calculated through astronomical formulas or obtained from professional meteorological databases and solar energy resource assessment software. For example, the solar altitude angle and azimuth angle can be calculated through astronomical algorithms based on the date, time, geographical location (longitude, latitude ), etc.: Among them, is the solar declination, is the hour angle.

[0074] Calculating the projected area of each heliostat field based on the initial coordinates and height includes: Calculating the vertical height where the center of the current heliostat field is located according to the initial coordinates and height ; Based on the vertical height where the center of the heliostat field is located , calculating the distance from the center of the current heliostat field to the receiver ; Based on the distance and the actual area of the current heliostat field , calculating the projected area of the current heliostat field .

[0075] Among them, is the vertical height where the heliostat center is located, is the circumferential radial spacing, is the slope angle, is the distance from the heliostat center to the ground, is the distance from the heliostat center to the receiver, is the height where the receiver is located, is the projected area of the heliostat, is the actual area of the heliostat.

[0076] The main function of the heliostat field is to collect solar energy. A larger projected area coverage means that more light can be reflected and focused on the receiver, thereby improving the solar energy collection efficiency. This goal is directly related to the energy output of the mirror field. During the optimization process, it is necessary to try to arrange the heliostats so that a larger effective illumination area can be formed at the receiver.

[0077] Projection overlap will cause the light reflected by some heliostats to be blocked by other heliostats, resulting in waste of light energy. Reducing projection overlap can improve the utilization efficiency of heliostats, enabling the light reflected by each heliostat to effectively reach the receiver and avoiding unnecessary energy losses.

[0078] The heliostats need to perform mechanical movements during the process of tracking the sun. If the spacing between the heliostats is too small or the layout is unreasonable, it may cause collisions or obstructions during the movement process. Ensuring unobstructed mechanical movement of the heliostats is a basic condition for ensuring the normal operation of the mirror field. When optimizing the layout plan, it is necessary to consider the movement range and space requirements of the heliostats.

[0079] Therefore, aiming at maximizing the projected area coverage, minimizing the projection overlap, and ensuring unobstructed mechanical movement of heliostats, the initial coordinates and height are multi-objectively optimized using the adaptive gravitational search algorithm, and the final heliostat field layout scheme obtained includes: Set the number of heliostats in the initial mirror field, encode the initial mirror field, and generate a matrix with N rows and R columns, where N is the number of mirror field layout schemes and R is the number of rows of heliostats; Evaluate N mirror field layout schemes using a multi-objective fitness function that includes projected area coverage, projection overlap rate, and mechanical movement constraints; Calculate the gravitational force based on the Euclidean distance between particles, calculate the acceleration by combining random parameters and the adaptive gravitational constant, and iteratively optimize the layout scheme through the velocity-position update formula; Continuously update the global optimal solution until the preset convergence accuracy is met, and output the heliostat coordinate matrix corresponding to the highest fitness value.

[0080] When setting the number of heliostats in the initial mirror field, the present invention sets the number of heliostats in the initial mirror field to be greater than or equal to 1.1 times the number of heliostats in the expected heliostat field, so as to reserve sufficient layout space for the final layout scheme. Set the adjacent row spacing to be 1 to 2 times the minimum increment to increase population diversity.

[0081] When evaluating N mirror field layout schemes, mark the highest and lowest fitness values to determine the current optimal layout scheme. The fitness function of the present invention is: Among them, represents at time the relative fitness of the th particle, at time the fitness value of the th particle (corresponding to a heliostat field layout scheme), which reflects the quality of this layout scheme for the multi-objective optimization problem (such as maximizing the projected area coverage and minimizing the projection overlap, etc.); at time represents the worst value among the fitness values of all particles in the particle swarm at time and is used as the reference lower limit for calculating the relative fitness; represents the best value among the fitness values of all particles in the particle swarm at time and is used as the reference upper limit for calculating the relative fitness; at time represents at time, the calculated value related to the mass of the The relative fitness of a particle respectively represent the -th particle and the -th particle, where $N$ represents the total number of particles.

[0082] Here, it is also necessary to perform a boundary check on the population to ensure that all particles are within the search space, avoid overlap between heliostats, and maintain the integrity of the population. Calculate the inertial mass of each particle and adaptively adjust the gravitational constant .

[0083] represents the gravitational constant at time $t$, represents the initial gravitational constant, represents the exponential function with base $e$, $\alpha$ is a coefficient that controls the rate of change of the gravitational constant and is a constant; represents time $t$, where $T$ represents the time period.

[0084] Calculate the gravitational force according to the Euclidean distance between particles, combine the random parameter and the adaptive gravitational constant to calculate the acceleration, and iteratively optimize the layout scheme through the velocity-position update formula, including: represents the gravitational force exerted by particle $i$ on particle $j$ at time $t$, and respectively represent the calculated values related to the masses of the -th particle and the -th particle at time $t$, represents the Euclidean distance between particle $i$ and particle $j$ at time $t$, $\epsilon$ is a very small positive number used to ensure that the denominator is not zero and avoid calculation errors; $\vec{r}_{it}$ and $\vec{r}_{jt}$ are respectively the position vectors of particle $i$ and particle $j$ at time $t$.

[0085] The acceleration calculation formula is: represents the acceleration of particle Acceleration represents a random number between 0 and 1. Introducing randomness can increase the diversity of the algorithm search and avoid falling into local optima.

[0086] Velocity - position update formula: represents the velocity of the th particle at time represents the velocity of the th particle at time represents the acceleration of the th particle at time represents the position of the th particle at time represents the position of the th particle at time represents the velocity of the th particle at time

[0087] Based on the initial mirror field of the three - dimensional mirror field terrain model, the present invention combines the AGSA algorithm to make the projection coverage area of the heliostats as large as possible in the field while satisfying the boundary conditions for the normal operation of the mirror field, with as little projection overlap area as possible, and ensuring unobstructed mechanical movement of the heliostats. Through parameter adjustment, better arrangement effects can be obtained within a determined area, and finally an optimal mirror field layout scheme can be achieved.

[0088] According to the second aspect of the embodiments of the present invention, there is also provided a generating device for the layout scheme of the heliostat field on mountainous terrain. Referring to Figure 3 shown, the generating device for the layout scheme of the heliostat field on mountainous terrain includes: A first model construction module 310, configured to construct an irregular triangular mesh model of the initial layout range of the heliostat field by using the Delaunay triangulation method; A second model construction module 320, configured to establish a three - dimensional terrain surface model according to the irregular triangular mesh model; A grid layout module 330, configured to perform regular grid layout on the three - dimensional terrain surface model; An interpolation analysis module 340, configured to perform interpolation analysis on each grid in the regular grid layout by using a cubic spline interpolation function to determine the initial coordinates and height of the heliostat field; A solution generation module 350 is configured to perform multi-objective optimization on the initial coordinates and height by using a projected area constraint condition and an adaptive gravitational search algorithm to obtain a final heliostat field layout solution.

[0089] It should be noted that although several modules and sub-modules of the apparatus for generating a heliostat field layout solution for mountainous terrain are mentioned in the above detailed description, such a division is not mandatory. In fact, according to an embodiment of the present invention, the features and functions of two or more of the above-described modules or sub-modules can be embodied in one module or unit. Conversely, the features and functions of one module or sub-module described above can be further divided and embodied by multiple modules or sub-modules.

[0090] In addition, in an exemplary embodiment of the present invention, an electronic device capable of implementing a method for generating a heliostat field layout solution for mountainous terrain is also provided.

[0091] Those skilled in the art can understand that various aspects of the present invention can be implemented as a system, a method, or a program product. Therefore, various aspects of the present invention can be specifically implemented in the following forms, namely: a complete hardware embodiment, a complete software embodiment (including firmware, microcode, etc.), or an embodiment combining hardware and software aspects, which can be collectively referred to as "circuit", "module", or "system" here.

[0092] Next, reference is made to Figure 4 to describe the electronic device 400 according to such an embodiment of the present invention. Figure 4 The illustrated electronic device 400 is merely an example and should not impose any limitation on the functions and usage scope of the embodiments of the present invention.

[0093] As Figure 4 shown, the electronic device 400 is presented in the form of a general-purpose computing device. The components of the electronic device 400 may include, but are not limited to: at least one of the above-mentioned processing units 410, at least one of the above-mentioned storage units 420, a bus 430 connecting different system components (including the storage unit 420 and the processing unit 410), and a display unit 440.

[0094] Among them, the storage unit stores program codes, and the program codes can be executed by the processing unit 410, so that the processing unit 410 executes the steps according to various exemplary embodiments of the present invention described in the above "Exemplary Method" section of the present invention. For example, the processing unit 410 can execute as Figure 2S210 shown in FIG. [ID]: Using the Delaunay triangulation method, construct an irregular triangular network model of the initial heliostat field layout range; S220: Establish a three-dimensional terrain surface model according to the irregular triangular network model; S230: Implement a regular grid layout on the three-dimensional terrain surface model; S240: Use a cubic spline interpolation function to perform interpolation analysis on each grid in the regular grid layout to determine the initial coordinates and height of the heliostat field; S250: Use the projection area constraint condition and the adaptive gravitational search algorithm to perform multi-objective optimization on the initial coordinates and height to obtain the final heliostat field layout plan.

[0095] The storage unit 420 may include a readable medium in the form of a volatile storage unit, such as a random access storage unit (RAM) 421 and / or a cache storage unit 422, and may further include a read-only storage unit (ROM) 423.

[0096] The storage unit 420 may also include a program / utility 424 having a set (at least one) of program modules 425. Such program modules 425 include, but are not limited to: an operating system, one or more application programs, other program modules, and program data. Each or some combination of these examples may include an implementation of a network environment.

[0097] The bus 430 may represent one or more of several types of bus structures, including a memory bus or memory controller, a peripheral bus, a graphics acceleration port, a processing unit, or a local bus using any of a variety of bus structures.

[0098] The electronic device 400 may also communicate with one or more external devices 470 (such as a keyboard, a pointing device, a Bluetooth device, etc.), may also communicate with one or more devices that enable a user to interact with the electronic device 400, and / or communicate with any device that enables the electronic device 400 to communicate with one or more other computing devices (such as a router, a modem, etc.). Such communication may be through the input / output (I / O) interface 450. And, the electronic device 400 may also communicate with one or more networks (such as a local area network (LAN), a wide area network (WAN), and / or a public network, such as the Internet) through the network adapter 460. As shown in the figure, the network adapter 460 communicates with other modules of the electronic device 400 through the bus 430. It should be understood that, although not shown in the figure, other hardware and / or software modules may be used in conjunction with the electronic device 400, including but not limited to: microcode, device drivers, redundant processing units, external disk drive arrays, RAID systems, tape drives, and data backup storage systems, etc.

[0099] Through the description of the above embodiments, those skilled in the art can easily understand that the exemplary embodiments described herein can be implemented by software or by a combination of software and necessary hardware. Therefore, the technical solutions according to the embodiments of the present invention can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (such as a CD-ROM, a USB flash drive, a mobile hard disk, etc.) or on a network, including several instructions to enable a computing device (such as a personal computer, a server, a terminal device, or a network device, etc.) to execute the method according to the embodiments of the present invention.

[0100] In an exemplary embodiment of the present invention, there is also provided a computer-readable storage medium having stored thereon a program product capable of implementing the above method of the present invention. In some possible embodiments, various aspects of the present invention can also be implemented in the form of a program product, which includes program code. When the program product runs on a terminal device, the program code is used to cause the terminal device to execute the steps according to various exemplary embodiments of the present invention described in the above "Exemplary Method" section of the present invention.

[0101] Reference Figure 5 As shown, a program product 500 for implementing the method for generating the heliostat field layout scheme for the above mountain terrain according to the embodiments of the present invention is described. It can be a portable compact disc read-only memory (CD-ROM) and includes program code, and can run on a terminal device, such as a personal computer. However, the program product of the present invention is not limited thereto. In the present invention, the readable storage medium can be any tangible medium that contains or stores a program, and the program can be used by or in combination with an instruction execution system, apparatus, or device.

[0102] The program product can adopt any combination of one or more readable storage media. The readable storage medium can be, for example, but not limited to, an electrical, magnetic, optical, electromagnetic, infrared, or semiconductor system, apparatus, or device, or any combination of the above. More specific examples (non-exhaustive list) of the readable storage medium include: an electrical connection having one or more wires, a portable disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), an optical fiber, a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the above.

[0103] The program code for performing the operations of the present invention can be written in any combination of one or more programming languages, including object-oriented programming languages such as Java, C++, etc., and also including conventional procedural programming languages such as the "C" language or similar programming languages. The program code can be executed entirely on the user's computing device, partially on the user's device, executed as a stand-alone software package, partially on the user's computing device and partially on a remote computing device, or entirely on a remote computing device or server. In the case of a remote computing device, the remote computing device can be connected to the user's computing device through any type of network, including a local area network (LAN) or a wide area network (WAN), or, alternatively, can be connected to an external computing device (e.g., by using an Internet service provider to connect through the Internet).

[0104] In addition, the above-mentioned drawings are only schematic illustrations of the processes included in the method according to the exemplary embodiments of the present invention, and are not for limiting purposes. It is easy to understand that the processes shown in the above-mentioned drawings do not indicate or limit the chronological order of these processes. Additionally, it is also easy to understand that these processes can be executed, for example, synchronously or asynchronously in multiple modules.

[0105] From the description of the above embodiments, those skilled in the art can easily understand that the example embodiments described herein can be implemented by software, or can be implemented by a combination of software and necessary hardware. Therefore, the technical solutions according to the embodiments of the present invention can be embodied in the form of a software product, which can be stored in a non-volatile storage medium (which can be a CD-ROM, a USB flash drive, a portable hard disk, etc.) or on a network, and includes several instructions to enable a computing device (which can be a personal computer, a server, a touch terminal, or a network device, etc.) to execute the method according to the embodiments of the present invention.

[0106] After considering the specification and practicing the invention disclosed herein, those skilled in the art will readily conceive of other embodiments of the present invention. The present invention is intended to cover any variations, uses, or adaptations of the present invention, which follow the general principles of the present invention and include known common general knowledge or conventional technical means in the technical field not disclosed by the present invention. The specification and embodiments are only regarded as exemplary, and the true scope and spirit of the present invention are pointed out by the claims.

[0107] It should be understood that the present invention is not limited to the exact structures described above and shown in the drawings, and various modifications and changes can be made without departing from its scope. The scope of the present invention is only limited by the appended claims.

Claims

1. A method for generating a heliostat field layout plan for mountainous terrain, characterized in that: include: The irregular triangle network model of the initial heliostat field layout is constructed using the Delaunay triangulation method. Establishing a three-dimensional terrain surface model based on the irregular triangle network model; Implementing regularized grid arrangement on the three-dimensional terrain surface model; Using a cubic spline interpolation function to perform interpolation analysis on each grid in the regularized grid arrangement to determine the initial coordinates and height of the heliostat field; The initial coordinates and height are optimized by multi-objective using projection area constraints and an adaptive gravitational search algorithm to obtain a final heliostat field layout plan.

2. The method for generating a heliostat field arrangement plan for mountainous terrain according to claim 1, characterized in that: Using a cubic spline interpolation function to perform interpolation analysis on each grid in the regularized grid arrangement to determine the initial coordinates and height of the heliostat field includes: The interpolation value of each grid is calculated using the cubic spline interpolation function; Based on the interpolation value of each grid, the initial coordinates and height of the heliostat field corresponding to each grid within the boundary condition are estimated.

3. The method for generating a heliostat field arrangement plan for mountainous terrain according to claim 2, characterized in that: The boundary conditions include: Defined in the interval The cubic spline interpolation function on The value at is defined in the interval The cubic spline interpolation function on The values ​​at are equal; Defined in the interval The cubic spline interpolation function on The first-order derivative value at is the same as that defined on the interval The cubic spline interpolation function on The values ​​of the first-order derivatives at are equal; and Defined in the interval The cubic spline interpolation function on The second-order derivative value at is the same as that defined on the interval The cubic spline interpolation function on The second-order derivative values ​​at are equal, where Indicates The coordinate value of a vertex of a mesh, Indicates The coordinate values ​​of the grid, Indicates The coordinate values ​​of the grid.

4. The method for generating a heliostat field arrangement plan for mountainous terrain according to claim 1, characterized in that: The initial coordinates and height are optimized by multi-objective using the projection area constraint and the adaptive gravitational search algorithm, and the final heliostat field layout scheme is obtained, including: Selecting solar parameters at typical times as reference data, and calculating the projection area of ​​each heliostat field based on the initial coordinates and height; With the goal of maximizing the projection area coverage, minimizing the projection overlap, and unhindered mechanical movement of the heliostats, the initial coordinates and heights are optimized by using an adaptive gravitational search algorithm to obtain a final heliostat field layout plan.

5. The method for generating a heliostat field arrangement plan for mountainous terrain according to claim 4, characterized in that: Selecting the solar parameters at a typical time as reference data, and calculating the projection area of ​​each heliostat field based on the initial coordinates includes: Calculate the vertical height of the current heliostat field center according to the initial coordinates and height; Calculate the distance between the current heliostat field center and the receiver based on the vertical height of the heliostat field center; Based on the distance and the actual area of ​​the current heliostat field, a projected area of ​​the current heliostat field is calculated.

6. The method for generating a heliostat field arrangement plan for mountainous terrain according to claim 4, characterized in that: With the goal of maximizing the projection area coverage, minimizing the projection overlap, and ensuring unimpeded mechanical movement of the heliostats, the initial coordinates and heights are optimized using an adaptive gravitational search algorithm to obtain the final heliostat field layout plan, including: Setting the number of heliostats in an initial mirror field, encoding the initial mirror field, and generating a matrix with N rows and R columns, where N is the number of mirror field layout schemes, and R is the number of heliostat rows; A multi-objective fitness function including projection area coverage, projection overlap and mechanical motion constraints is used to evaluate N mirror field layout schemes. The gravitational force is calculated based on the Euclidean distance between particles, the acceleration is calculated by combining random parameters and adaptive gravitational constants, and the layout scheme is iteratively optimized through the velocity-position update formula; The global optimal solution is continuously updated until the preset convergence accuracy is met, and the heliostat coordinate matrix corresponding to the highest fitness value is output.

7. The method for generating a heliostat field arrangement plan for mountainous terrain according to claim 6, characterized in that: The number of heliostats in the initial mirror field is greater than or equal to 1.1 times the number of heliostats in the expected heliostat field.

8. A device for generating a heliostat field arrangement plan for mountainous terrain, characterized in that: include: The first model building module is used to build an irregular triangle network model of the initial heliostat field layout range by using the Delaunay triangulation method; A second model building module is used to build a three-dimensional terrain surface model based on the irregular triangle network model; A grid arrangement module, used for implementing regularized grid arrangement on the three-dimensional terrain surface model; An interpolation analysis module, used to perform interpolation analysis on each grid in the regularized grid arrangement using a cubic spline interpolation function to determine the initial coordinates and height of the heliostat field; The scheme generation module is used to perform multi-objective optimization on the initial coordinates and height by using the projection area constraint condition and the adaptive gravitational search algorithm to obtain the final heliostat field layout scheme.

9. An electronic device, characterized in that: include: processor; as well as A memory having computer-readable instructions stored thereon, wherein the computer-readable instructions, when executed by the processor, implement the method for generating a heliostat field arrangement plan for a mountainous terrain according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and when the computer program is executed by a processor, the method for generating a heliostat field arrangement plan for a mountainous terrain as claimed in any one of claims 1 to 7 is implemented.

Images