A pumped storage power station site selection method based on digital elevation model
The pumped storage power station site selection method based on digital elevation model solves the problem of low computational efficiency of traditional site selection methods in non-mountainous terrain areas, realizes efficient site selection over a large area, finds more candidate dam sites, and improves computational efficiency and applicability.
Patent Information
- Authority / Receiving Office
- CN · China
- Patent Type
- Applications(China)
- Current Assignee / Owner
- CHINA WATER NORTHEASTERN INVESTIGATION DESIGN & RES
- Filing Date
- 2026-06-25
- Publication Date
- 2026-07-24
AI Technical Summary
In Northeast China, North China, and Eastern Inner Mongolia, the terrain is mainly plains and hills with few mountains. Traditional pumped storage power station site selection methods are inefficient in small-area site selection, have low computational efficiency, and lack effective screening mechanisms, resulting in a waste of computational resources.
A site selection method based on digital elevation model (DEM) is adopted. By preprocessing the DEM file, unsuitable terrain areas are eliminated, contour lines and river network lines are generated, and intersection points are calculated to generate candidate dam points and reservoir basins, thereby reducing invalid calculations and improving calculation efficiency.
It can quickly find more candidate dam sites in a large area, improve calculation efficiency, reduce human intervention, enhance the objectivity and applicability of site selection, and is suitable for non-mountainous terrain areas.
Smart Images

Figure CN122454094A_ABST
Abstract
Description
Technical Field
[0001] This invention belongs to the field of power plant site selection technology, specifically relating to the design of a pumped storage power plant site selection method based on a digital elevation model. Background Technology
[0002] Pumped storage power stations are important energy storage facilities that utilize surplus electricity during off-peak hours to pump water from a lower reservoir to an upper reservoir, releasing it to generate electricity during peak hours, thus storing and releasing electrical energy. Site selection is a crucial aspect of pumped storage power station construction; the rationality and efficiency of the site selection directly affect the economic viability and feasibility of the power station.
[0003] Traditional methods for selecting pumped-storage power station sites primarily target areas with significant topographic relief, such as mountains and hills, determining candidate dam sites through manual reconnaissance and terrain analysis. However, in regions like Northeast China, North China, and Eastern Inner Mongolia, the terrain is predominantly plains and hills, with fewer mountains and relatively flat terrain. Selecting pumped-storage power station sites under these terrain conditions presents the following technical challenges: (1) Low site selection efficiency in small areas: Due to the small topographic relief, when selecting a site in a small area, there are few candidate hydropower station sites that can meet the construction conditions of pumped storage power stations (such as the need for sufficient height difference and reservoir capacity), resulting in low site selection efficiency.
[0004] (2) Low computational efficiency of large-scale site selection: In order to find enough candidate dam sites, it is necessary to expand the site selection area, for example, within the area of a city (about 50,000 square kilometers). However, traditional site selection methods require detailed topographic analysis and calculation of the entire area, which involves a huge amount of calculation and has extremely low computational efficiency, making it difficult to complete the site selection work within a reasonable time.
[0005] (3) Lack of effective screening mechanism: Existing technology lacks an effective method to quickly screen out terrain areas that meet the basic conditions in a large area, resulting in a large amount of invalid calculations and wasting computing resources.
[0006] Therefore, there is an urgent need for a technical solution that can efficiently select sites for pumped storage power stations in large areas of non-mountainous terrain, which can both expand the site selection range to find enough candidate dam sites and improve computational efficiency and reduce invalid calculations. Summary of the Invention
[0007] The purpose of this invention is to propose a method for selecting pumped storage power stations based on digital elevation models (DEMs). This method enables large-area pumped storage power station site selection in non-mountainous terrain areas using DEMs, which can both expand the site selection range to find sufficient candidate dam sites and improve computational efficiency, avoiding inefficient detailed calculations for the entire large area.
[0008] The technical solution of this invention is: a method for site selection of pumped storage power stations based on digital elevation models, comprising the following steps: S1. Process the DEM file according to preset constraints, remove terrain areas that do not meet the conditions, and obtain the filtered terrain areas.
[0009] S2. Generate contour lines based on the selected terrain regions and preset contour intervals.
[0010] S3. Generate river network lines based on the selected terrain regions.
[0011] S4. Generate candidate dam points based on the intersections of contour lines and river network lines.
[0012] S5. Generate candidate reservoir basins based on candidate dam points.
[0013] S6. Evaluate and select sites for candidate reservoir basins to complete the site selection for pumped storage power stations.
[0014] Further, step S1 includes the following sub-steps: S11, using each raster cell in the DEM file Define a search window centered on a given location, and define the radius of the search window. for: in Indicates the search distance. Indicates the spatial resolution of DEM data. This represents the function for rounding up.
[0015] S12. Calculate the maximum elevation value of each grid cell. and minimum elevation value : in Represents grid cells Elevation value, Represented by grid cells Centered on, with radius as The collection of all grid cells within the search window.
[0016] S13, Based on the maximum elevation value and minimum elevation value Computational raster cells Maximum elevation difference within the search window : S14. Based on the maximum elevation difference and the preset elevation difference threshold For each grid cell Make a judgment: in Indicates the judgment value. Represents grid cells Meets the requirements. Represents grid cells It does not meet the requirements.
[0017] S15, Reservation The grid cells will The elevation values of the raster cells are set to invalid values or removed from the dataset to obtain the filtered terrain regions.
[0018] Furthermore, step S2 includes the following sub-steps: S21. Based on the minimum elevation value in the selected terrain regions. Maximum elevation value and preset contour interval Determine the sequence of contour line elevation values : in Indicates the first Each contour line elevation value, Indicates the number of contour line elevation values. This represents the floor function.
[0019] S22, For each contour line elevation value Track all elevation values equal to in DEM data The raster cells are used to generate contour lines, for elevation values that are not equal to... The grid cells are used to determine the positions through which contour lines pass using linear interpolation.
[0020] S23. Smooth the generated contour lines to eliminate jagged edges caused by the resolution limitations of the DEM data, and obtain smooth contour vector data.
[0021] Furthermore, step S4 includes the following sub-steps: S41. Calculate each contour line. and each river network line geometric intersection : in This represents the geometric intersection operation.
[0022] S42. For any two intersection points and , in response to the distance between them Then only one intersection point is retained to achieve the deduplication of intersection points. Indicates the minimum spacing between dam points. Indicates the intersection point coordinates Indicates the intersection point The coordinates.
[0023] S43. Generate a set of candidate dam points based on the intersections after deduplication and filtering. ,in This represents the total number of candidate dam sites.
[0024] Furthermore, step S5 includes the following sub-steps: S51. Determining candidate dam sites The river network line where it is located Calculate the river network line At candidate dam sites Tangent direction vector at point ,in .
[0025] S52, based on the tangent direction vector At candidate dam sites A perpendicular line is generated at the location that is perpendicular to the river network segment. : in This represents the direction vector of the vertical line. Indicates the calculated parameters. Indicates the length of the vertical line.
[0026] S53, Calculate the vertical line and contour lines intersection Intersection Vertical line Points that satisfy the following conditions: in Represents vertical line Points on The corresponding DEM elevation.
[0027] S54, vertical line and contour lines All intersections Connecting regions in ascending order of elevation, forming closed polygonal areas, serves as candidate reservoir basins: in Indicates the first A closed region surface, This represents the set of intersection points sorted by elevation from low to high. This represents the basin surface generated from a point.
[0028] S55. Repeat steps S51 to S54 until all candidate reservoirs corresponding to candidate dam points are generated.
[0029] Furthermore, the site selection assessment in step S6 includes reservoir basin screening, geological condition assessment, and environmental impact assessment.
[0030] The beneficial effects of this invention are: (1) The present invention eliminates a large number of terrain areas that do not meet the basic conditions in the early stage of calculation through the preliminary terrain screening in step S1, which greatly reduces the amount of subsequent calculation and significantly improves the calculation efficiency of large-scale pumped storage power station site selection.
[0031] (2) Due to the improved computational efficiency, the present invention can be used to select sites in a larger area (such as 50,000 square kilometers), thereby finding more candidate dam sites and solving the problem of insufficient number of candidate dam sites in small areas.
[0032] (3) In this invention, key parameters such as search distance range, elevation difference threshold, and contour interval can be configured and can be flexibly adjusted according to the terrain characteristics and site selection requirements of different regions, which improves the applicability and practicality of the method.
[0033] (4) The entire site selection process in this invention is automatically completed based on DEM data, which reduces manual intervention and improves the objectivity and consistency of site selection.
[0034] (5) This invention is applicable to non-mountainous terrain areas with more plains and hills and fewer mountains, filling the technical gap in the site selection of large-scale pumped storage power stations in such areas. Attached Figure Description
[0035] Figure 1 The diagram shown is a flowchart of a pumped storage power station site selection method based on a digital elevation model provided by an embodiment of the present invention. Detailed Implementation
[0036] Exemplary embodiments of the present invention will now be described in detail with reference to the accompanying drawings. It should be understood that the embodiments shown and described in the drawings are merely exemplary and are intended to illustrate the principles and spirit of the invention, and are not intended to limit the scope of the invention.
[0037] This invention provides a method for site selection of pumped storage power stations based on digital elevation models, such as... Figure 1As shown, the process includes the following steps S1 to S6: S1. Process the DEM file according to preset constraints, remove terrain areas that do not meet the conditions, and obtain the filtered terrain areas to reduce the amount of subsequent calculations.
[0038] In this embodiment of the invention, the preset limiting condition includes the existence of a required height difference (e.g., 50 meters) within a specified distance range (e.g., 300 meters). The specified distance range (e.g., 300 meters) and the required height difference (e.g., 50 meters) are configurable parameters and can be adjusted according to actual site selection requirements.
[0039] Step S1 includes the following sub-steps S11 to S15: S11, using each raster cell in the DEM file Define a search window centered on a given location, and define the radius of the search window. for: in This indicates the search distance (set to 300 meters in this embodiment of the invention). Indicates the spatial resolution of DEM data. This represents the function for rounding up.
[0040] S12. Calculate the maximum elevation value of each grid cell. and minimum elevation value : in Represents grid cells Elevation value, Represented by grid cells Centered on, with radius as The collection of all grid cells within the search window.
[0041] S13, Based on the maximum elevation value and minimum elevation value Computational raster cells Maximum elevation difference within the search window : S14. Based on the maximum elevation difference and the preset elevation difference threshold (In this embodiment of the invention, it is set to 50m) for each grid cell Make a judgment: in Indicates the judgment value. Represents grid cells Meets the requirements. Represents grid cells It does not meet the requirements.
[0042] S15. Generate a new DEM data file based on the screening results, and retain it. The grid cells will The elevation values of the raster cells are set to invalid values (such as NoData) or removed from the dataset to obtain the filtered terrain regions.
[0043] Through the above-mentioned preliminary terrain screening process, a large number of areas that do not meet the basic terrain conditions (i.e., lack sufficient height difference in a local area) can be quickly eliminated in the early stage of calculation, which greatly reduces the amount of calculation in subsequent steps such as contour line generation, river network extraction, and dam point generation, and significantly improves the calculation efficiency of large-scale site selection.
[0044] S2. Generate contour lines based on the selected terrain regions and preset contour intervals.
[0045] Step S2 includes the following sub-steps S21 to S23: S21. Based on the minimum elevation value in the selected terrain regions. Maximum elevation value and preset contour interval Determine the sequence of contour line elevation values : in Indicates the first Each contour line elevation value, Indicates the number of contour line elevation values. This represents the floor function. In this embodiment of the invention, the contour interval... It can be set, for example, to 5 meters, 10 meters, 20 meters, etc., and can be selected according to the terrain characteristics and accuracy requirements.
[0046] S22, For each contour line elevation value Track all elevation values equal to in DEM data The raster cells are used to generate contour lines, for elevation values that are not equal to... The grid cells are used to determine the positions through which contour lines pass using linear interpolation.
[0047] S23. Smooth the generated contour lines to eliminate jagged edges caused by the resolution limitations of the DEM data, and obtain smooth contour vector data.
[0048] In this embodiment of the invention, the generated contour data is stored in vector format, and each contour line contains its corresponding elevation value attribute.
[0049] S3. Generate river network lines based on the selected terrain regions.
[0050] In this embodiment of the invention, the open-source GIS software package Whitebox is used to generate river network lines. The generation efficiency for large-scale terrain is within 10 minutes. The generated river network line data is stored in vector format and includes attribute information such as river network level and flow direction.
[0051] S4. Generate candidate dam points based on the intersections of contour lines and river network lines.
[0052] The intersection of contour lines and river network lines usually presents good conditions for dam construction because there is both a waterway passing through and a suitable terrain elevation.
[0053] Step S4 includes the following sub-steps S41 to S43: S41. Calculate each contour line. (corresponding elevation value) ) and each river network line geometric intersection : in This represents geometric intersection operations. Since contour lines and river network lines are both vector segments, intersection point calculation can be achieved using a line segment intersection algorithm. Let the contour lines... Consisting of a series of ordered points Composition, river network line From ordered points Composition, for each intersection point Its elevation value is equal to the elevation value of the corresponding contour line. .
[0054] S42. For any two intersection points and , in response to the distance between them Then only one intersection point is retained to achieve the deduplication of intersection points. This indicates the minimum distance between dam points (100 meters in this embodiment of the invention). Indicates the intersection point coordinates Indicates the intersection point The coordinates.
[0055] S43. Generate a set of candidate dam points based on the intersections after deduplication and filtering. ,in This represents the total number of candidate dam sites.
[0056] S5. Generate candidate reservoir basins based on candidate dam points.
[0057] Step S5 includes the following sub-steps S51 to S55: S51. Determining candidate dam sites The river network line where it is located Calculate the river network line At candidate dam sites Tangent direction vector at point ,in .
[0058] S52, based on the tangent direction vector At candidate dam sites A perpendicular line is generated at the location that is perpendicular to the river network segment. : in The direction vector of the vertical line (and) vertical), Indicates the calculated parameters. Indicates the length of the vertical line.
[0059] S53, Calculate the vertical line and contour lines intersection Intersection Vertical line Points that satisfy the following conditions: in Represents vertical line Points on The corresponding DEM elevation.
[0060] S54, vertical line and contour lines All intersections Connecting regions in ascending order of elevation, forming closed polygonal areas, serves as candidate reservoir basins: in Indicates the first A closed region surface, This represents the set of intersection points sorted by elevation from low to high. This represents the basin surface generated from a point.
[0061] S55. Repeat steps S51 to S54 until all candidate reservoirs corresponding to candidate dam points are generated.
[0062] S6. Evaluate and select sites for candidate reservoir basins to complete the site selection for pumped storage power stations.
[0063] In this embodiment of the invention, the site selection assessment includes reservoir basin screening, geological condition assessment, and environmental impact assessment.
[0064] Those skilled in the art will recognize that the embodiments described herein are intended to help the reader understand the principles of the invention, and should be understood that the scope of protection of the invention is not limited to such specific statements and embodiments. Those skilled in the art can make various other specific modifications and combinations based on the technical teachings disclosed in this invention without departing from the spirit of the invention, and these modifications and combinations are still within the scope of protection of this invention.
Claims
1. A method for site selection of pumped storage power stations based on digital elevation models, characterized in that, Includes the following steps: S1. Process the DEM file according to preset constraints, remove terrain areas that do not meet the conditions, and obtain the filtered terrain areas. S2. Generate contour lines based on the selected terrain regions and preset contour intervals; S3. Generate river network lines based on the filtered terrain areas; S4. Generate candidate dam points based on the intersections of contour lines and river network lines; S5. Generate candidate reservoir basins based on candidate dam sites; S6. Evaluate and select sites for candidate reservoir basins to complete the site selection for pumped storage power stations.
2. The method for selecting a pumped storage power station based on a digital elevation model according to claim 1, characterized in that, Step S1 includes the following sub-steps: S11, using each raster cell in the DEM file A search window is defined centered on a point, the radius of which is... for: in Indicates the search distance. Indicates the spatial resolution of DEM data. This represents the floor function; S12. Calculate the maximum elevation value of each grid cell. and minimum elevation value : in Represents grid cells Elevation value, Represented by grid cells Centered on, with radius as The collection of all grid cells within the search window; S13, Based on the maximum elevation value and minimum elevation value Computational raster cells Maximum elevation difference within the search window : S14. Based on the maximum elevation difference and the preset elevation difference threshold For each grid cell Make a judgment: in Indicates the judgment value. Represents grid cells Meets the requirements. Represents grid cells Does not meet the requirements; S15, Reservation The grid cells will The elevation values of the raster cells are set to invalid values or removed from the dataset to obtain the filtered terrain regions.
3. The method for selecting a pumped storage power station based on a digital elevation model according to claim 1, characterized in that, Step S2 includes the following sub-steps: S21. Based on the minimum elevation value in the selected terrain regions. Maximum elevation value and preset contour interval Determine the sequence of contour line elevation values : in Indicates the first Each contour line elevation value, Indicates the number of contour line elevation values. This represents the floor function; S22, For each contour line elevation value Track all elevation values equal to in DEM data The raster cells are used to generate contour lines, for elevation values that are not equal to... The grid cells are used to determine the positions through which contour lines pass using linear interpolation. S23. Smooth the generated contour lines to eliminate jagged edges caused by the resolution limitations of the DEM data, and obtain smooth contour vector data.
4. The method for selecting a pumped storage power station based on a digital elevation model according to claim 3, characterized in that, Step S4 includes the following sub-steps: S41. Calculate each contour line. and each river network line geometric intersection : in This represents the geometric intersection operation; S42. For any two intersection points and In response to the distance between them Then only one intersection point is retained to achieve the deduplication of intersection points. Indicates the minimum spacing between dam points. Indicates the intersection point coordinates Indicates the intersection point The coordinates; S43. Generate a set of candidate dam points based on the intersections after deduplication and filtering. ,in This represents the total number of candidate dam sites.
5. The method for selecting a pumped storage power station based on a digital elevation model according to claim 4, characterized in that, Step S5 includes the following sub-steps: S51. Determining candidate dam sites The river network line Calculate the river network line At candidate dam sites Tangent direction vector at point ,in ; S52, based on the tangent direction vector At candidate dam sites A perpendicular line is generated at the location that is perpendicular to the river network segment. : in This represents the direction vector of the vertical line. Indicates the calculated parameters. Indicates the length of the vertical line; S53, Calculate the vertical line and contour lines intersection The intersection point Vertical line Points that satisfy the following conditions: in Represents vertical line Points on Corresponding DEM elevation; S54, vertical line and contour lines All intersections Connect the regions in ascending order of elevation to form closed polygonal areas, which will serve as candidate reservoir basins: in Indicates the first A closed region surface, This represents the set of intersection points sorted by elevation from low to high. This represents the basin surface generated from a point; S55. Repeat steps S51 to S54 until all candidate reservoirs corresponding to candidate dam points are generated.
6. The method for selecting a pumped storage power station based on a digital elevation model according to claim 1, characterized in that, The site selection assessment in step S6 includes reservoir / basin screening, geological condition assessment, and environmental impact assessment.