A pumped storage power station dam engineering simulation method and system considering terrain

Abstract

The present disclosure provides a method and system for simulating pumped storage power station dam engineering considering terrain, comprising: obtaining DEM data of a reservoir area where a pumped storage power station dam is located; based on the DEM data, generating a reservoir surface and a horizontal transverse center line of a dam top plane by using a spatial data analysis method; based on the obtained horizontal transverse center line of the dam top plane and a preset width of the dam top plane, generating a dam top plane; and based on the constructed slope equations of both sides of the dam and the obtained boundary of the slope of both sides of the dam, constructing a dam body geometric shape and obtaining the volume of the dam body; according to the obtained DEM data, generating a thickness planing surface boundary of the mountain bodies on both ends of the dam; based on the obtained reservoir surface, dam body geometric shape, dam body volume and thickness planing surface boundary of the mountain bodies on both ends of the dam, constructing a dam body model, and based on the constructed dam body model, simulating pumped storage power station dam engineering.

Description

Technical Field

[0001] This disclosure belongs to the field of dam engineering simulation technology, and in particular relates to a simulation method and system for pumped storage power station dam engineering that takes into account the terrain. Background Technology

[0002] The statements in this section are merely background information relating to this disclosure and do not necessarily constitute prior art.

[0003] Dams are fundamental infrastructure for pumped storage power stations, and dam design and construction are crucial aspects of their development. However, the inventors discovered a lack of research on dam engineering simulation for pumped storage power stations in existing technologies. Furthermore, current dam modeling methods for general scenarios typically only consider the dam itself, neglecting the impact of the surrounding environment on the simulation and consequently overlooking the dam's safety. Moreover, existing dam modeling methods are significantly inadequate in terms of modeling accuracy and efficiency, making them unsuitable for direct application in the simulation of pumped storage power station dam projects. Summary of the Invention

[0004] To address the aforementioned problems, this disclosure provides a simulation method and system for pumped storage power station dam projects that takes into account topography. The scheme is based on DEM data of the reservoir area where the pumped storage power station dam is located. It generates horizontal centerlines of the reservoir surface and the dam crest plane using spatial data analysis methods, and constructs the dam's geometric shape using these horizontal centerlines as a reference. This effectively improves the accuracy and efficiency of dam modeling. Furthermore, in the simulation of the pumped storage power station dam project, the analysis of the thickness profiles of the mountains on which the dam rests is located fully considers the dam's safety, achieving accurate and comprehensive simulation analysis of the pumped storage power station dam project.

[0005] According to a first aspect of the present disclosure, a simulation method for pumped storage power station dam engineering that takes into account terrain is provided, including:

[0006] Obtain DEM data of the reservoir area where the pumped storage power station dam is located;

[0007] Based on the DEM data, a spatial data analysis method is used to generate the horizontal centerline of the reservoir surface and the dam crest plane. The spatial data analysis method includes: extracting contour lines and runoff lines from the DEM data of the reservoir area; generating perpendicular segments of the runoff lines based on the obtained runoff lines; calculating the intersection points of the perpendicular segments and contour lines, and using the line connecting the two sides of the contour lines closest to the center point of the contour line as the interception line segment; and determining the horizontal centerline of the reservoir surface and the dam crest plane based on the interception line segment and the contour lines.

[0008] Based on the obtained horizontal centerline of the dam crest plane and the preset width of the dam crest plane, the dam crest plane is generated; and based on the constructed equations of the dam's two sides and the obtained boundaries of the dam's two sides, the geometry of the dam body is constructed, and the volume of the dam body is obtained.

[0009] Based on the obtained DEM data, the thickness profile boundary of the mountains on which the dam is supported is generated;

[0010] Based on the obtained reservoir surface, dam geometry, dam volume, and the thickness profile boundaries of the mountains on both ends of the dam, a dam model is constructed, and a pumped storage power station dam engineering simulation is performed based on the constructed dam model.

[0011] Furthermore, the horizontal centerline of the reservoir surface and the dam crest plane is determined based on the interception line segment and contour line. Specifically, the closed area enclosed by the contour line with the highest elevation and its corresponding interception line segment is selected as the reservoir surface, and the interception line segment is used as the horizontal centerline of the dam crest plane.

[0012] Furthermore, the generation of the dam top plane specifically involves: determining a first plane line and a second plane line parallel to the horizontal center line of the dam top plane, wherein the first plane line and the second plane line are located on both sides of the center line; determining the edge lines at both ends of the dam based on the intersection points of the first plane line and the second plane line with the contour line at the uppermost edge of the reservoir; and determining the top plane of the dam based on the edge lines at both ends of the dam, the first plane line, and the second plane line.

[0013] Furthermore, the boundaries of the sloping sides of the dam are obtained by using the horizontal scan line method to trace the boundaries based on the constructed equations of the sloping sides of the dam, thereby obtaining the boundaries of the inner and outer sloping sides of the dam and the volume covered by the inner and outer sloping sides of the dam.

[0014] Furthermore, the volume of the dam body includes the volume covered by the inner and outer slopes of the dam and the volume covered by the top plane of the dam, wherein the volume covered by the top plane of the dam is obtained based on the horizontal scan line method.

[0015] Furthermore, the step of generating the thickness profile boundary of the mountain body on which the two ends of the dam rely based on the obtained DEM data specifically involves: determining the elevation value at each layer based on the preset vertical layering parameters, the minimum elevation of the dam crest surface, and the elevation of the inner and outer side boundary; drawing a line parallel to the center line of the dam crest plane at the elevation of each layer; extending both ends of the parallel line; and generating the mountain thickness boundary line based on the obtained extension line.

[0016] Furthermore, the DEM data is point cloud data obtained by scanning with an airborne LiDAR. After filtering, smoothing, and removal of ground features, DEM data with a sub-decimeter-level regular grid is generated using a spatial interpolation method.

[0017] According to a second aspect of the present disclosure, a simulation system for pumped storage power station dam engineering that takes into account terrain is provided, comprising:

[0018] The data acquisition unit is used to acquire DEM data of the reservoir area where the pumped storage power station dam is located.

[0019] The reservoir surface and centerline generation unit is used to generate the horizontal centerline of the reservoir surface and dam crest plane based on the DEM data using spatial data analysis methods. The spatial data analysis method includes: extracting contour lines and runoff lines from the DEM data of the reservoir area; generating perpendicular segments of the runoff lines based on the obtained runoff lines; calculating the intersection points of the perpendicular segments and contour lines, and using the line connecting the two sides of the contour lines closest to the center point of the contour line as the interception line segment; and determining the horizontal centerline of the reservoir surface and dam crest plane based on the interception line segment and the contour lines.

[0020] The geometric shape construction unit is used to generate the dam top plane based on the obtained horizontal centerline of the dam top plane and the preset width of the dam top plane; and to construct the geometric shape of the dam body based on the constructed equations of the dam's two sides and the obtained boundaries of the dam's two sides, and to obtain the volume of the dam body.

[0021] The terrain generation unit is used to generate the thickness profile boundary of the mountains on which the dam is supported based on the obtained DEM data.

[0022] The simulation unit is used to construct a dam model based on the obtained reservoir surface, dam geometry, dam volume, and the thickness profile boundary of the mountains on both ends of the dam, and to perform pumped storage power station dam engineering simulation based on the constructed dam model.

[0023] According to a third aspect of the present disclosure, an electronic device is provided, including a memory, a processor, and a computer program stored in the memory and running on the memory. When the processor executes the program, it implements the aforementioned simulation method for pumped storage power station dam engineering that takes into account terrain.

[0024] According to a fourth aspect of the present disclosure, a non-transitory computer-readable storage medium is provided, on which a computer program is stored, which, when executed by a processor, implements the aforementioned simulation method for pumped storage power station dam engineering that takes into account terrain.

[0025] Compared with the prior art, the beneficial effects of this disclosure are:

[0026] (1) This disclosure provides a simulation method and system for pumped storage power station dam projects that takes into account the terrain. The scheme is based on the DEM data of the reservoir area where the pumped storage power station dam is located. It generates the horizontal center line of the reservoir surface and the dam top plane through spatial data analysis methods, and constructs the geometric shape of the dam body based on the horizontal center line of the dam top plane, which effectively improves the accuracy and efficiency of dam body modeling. At the same time, in the simulation of pumped storage power station dam projects, by introducing the analysis of the thickness profile of the mountains on both ends of the dam, the safety of the dam body is fully considered, and accurate and comprehensive simulation analysis of pumped storage power station dam projects is realized.

[0027] (2) The proposed scheme adopts an intelligent extraction method for the horizontal transverse centerline of the dam crest plane. Based on the DEM data of the reservoir area, contour lines and drainage lines are extracted; based on the obtained drainage lines, perpendicular segments of the drainage lines are generated; the intersection points of the perpendicular segments and contour lines are calculated, and the line connecting the intersection points on both sides of the contour lines closest to the center point of the contour lines is used as the interception line segment; based on the interception line segment and the contour lines, the horizontal transverse centerline of the reservoir surface and the dam crest plane is determined; at the same time, in the process of extracting drainage lines, a depression filling processing algorithm based on block matrix parallel operation is adopted. By using the method of CPU+GPU parallel data processing, the computation efficiency is significantly improved compared with the traditional algorithm.

[0028] (3) The proposed scheme adopts a horizontal scanning line method to trace the boundary. The scheme uses the horizontal scanning line method to trace the boundary based on the constructed equation of the slopes on both sides of the dam. While obtaining the boundaries of the inner and outer slopes of the dam, it can quickly and accurately determine the volume covered by the inner and outer slopes of the dam.

[0029] Advantages of this disclosure in additional aspects will be set forth in part in the description which follows, and in part will be obvious from the description, or may be learned by practice of this disclosure. Attached Figure Description

[0030] The accompanying drawings, which form part of this disclosure, are used to provide a further understanding of this disclosure. The illustrative embodiments of this disclosure and their descriptions are used to explain this disclosure and do not constitute an undue limitation of this disclosure.

[0031] Figure 1 This is a flowchart of a simulation method for a pumped storage power station dam project that takes into account the terrain, as described in an embodiment of this disclosure.

[0032] Figure 2 This is a schematic diagram of the one-click reservoir parameter extraction setting interface described in this embodiment of the disclosure;

[0033] Figure 3 This is a schematic diagram of the generated contour data described in the embodiments of this disclosure;

[0034] Figure 4 This is a schematic diagram of the extracted waterline data described in the embodiments of this disclosure;

[0035] Figure 5 This is a schematic diagram of the vertical line of the water catchment line described in the embodiments of this disclosure;

[0036] Figure 6 This is a schematic diagram of the extracted centerline of the dam crest plane as described in the embodiments of this disclosure;

[0037] Figure 7 This is a schematic diagram of the generated dam top plane as described in the embodiments of this disclosure;

[0038] Figure 8 This is a schematic diagram of the outer boundary of the dam body as described in the embodiments of this disclosure;

[0039] Figure 9(a) is a schematic diagram of the thickness of the left mountain body as described in the embodiment of this disclosure;

[0040] Figure 9(b) is a schematic diagram of the thickness of the two mountain bodies described in the embodiment of this disclosure;

[0041] Figure 10 This is a schematic diagram simulating half of the dam structure as described in the embodiments of this disclosure;

[0042] Figure 11(a) is a frontal view of the visualization of the 3D model of the dam described in the embodiment of this disclosure;

[0043] Figure 11(b) is a schematic diagram of the tilted view of the 3D model simulation visualization of the dam described in the embodiments of this disclosure. Detailed Implementation

[0044] The present disclosure will be further described below with reference to the accompanying drawings and embodiments.

[0045] It should be noted that the following detailed descriptions are illustrative and intended to provide further explanation of this disclosure. Unless otherwise specified, all technical and scientific terms used herein have the same meaning as commonly understood by one of ordinary skill in the art to which this disclosure pertains.

[0046] It should be noted that the terminology used herein is for the purpose of describing particular embodiments only and is not intended to limit the exemplary embodiments according to this disclosure. As used herein, the singular form is intended to include the plural form as well, unless the context clearly indicates otherwise. Furthermore, it should be understood that when the terms “comprising” and / or “including” are used in this specification, they indicate the presence of features, steps, operations, devices, components, and / or combinations thereof.

[0047] Where there is no conflict, the embodiments and features described herein can be combined with each other.

[0048] Example 1:

[0049] The purpose of this embodiment is to provide a simulation method for pumped storage power station dam projects that takes into account the terrain.

[0050] A simulation method for pumped storage power station dam engineering that takes into account terrain features includes:

[0051] Obtain DEM data of the reservoir area where the pumped storage power station dam is located;

[0052] Based on the DEM data, a spatial data analysis method is used to generate the horizontal centerline of the reservoir surface and the dam crest plane. The spatial data analysis method includes: extracting contour lines and runoff lines from the DEM data of the reservoir area; generating perpendicular segments of the runoff lines based on the obtained runoff lines; calculating the intersection points of the perpendicular segments and contour lines, and using the line connecting the two sides of the contour lines closest to the center point of the contour line as the interception line segment; and determining the horizontal centerline of the reservoir surface and the dam crest plane based on the interception line segment and the contour lines.

[0053] Based on the obtained horizontal centerline of the dam crest plane and the preset width of the dam crest plane, the dam crest plane is generated; and based on the constructed equations of the dam's two sides and the obtained boundaries of the dam's two sides, the geometry of the dam body is constructed, and the volume of the dam body is obtained.

[0054] Based on the obtained DEM data, the thickness profile boundary of the mountains on which the dam is supported is generated;

[0055] Based on the obtained reservoir surface, dam geometry, dam volume, and the thickness profile boundaries of the mountains on both ends of the dam, a dam model is constructed, and a pumped storage power station dam engineering simulation is performed based on the constructed dam model.

[0056] In specific implementation, the horizontal centerline of the reservoir surface and the dam crest plane is determined based on the interception line segment and contour line. Specifically, the closed area enclosed by the contour line with the highest elevation and its corresponding interception line segment is selected as the reservoir surface, and the interception line segment is used as the horizontal centerline of the dam crest plane.

[0057] In specific implementation, the generation of the dam top plane is as follows: based on the horizontal centerline of the dam top plane, a first plane line and a second plane line parallel to the centerline are determined, wherein the first plane line and the second plane line are located on both sides of the centerline; based on the intersection of the first plane line and the second plane line with the contour line of the uppermost edge of the reservoir, the edge lines at both ends of the dam are determined; based on the edge lines at both ends of the dam, the first plane line, and the second plane line, the top plane of the dam is determined.

[0058] In specific implementation, the boundary of the sloping sides of the dam is obtained by using the horizontal scan line method to trace the boundary based on the constructed equation of the sloping sides of the dam, thereby obtaining the boundary of the inner and outer sloping sides of the dam and the volume covered by the inner and outer sloping sides of the dam.

[0059] In specific implementation, the volume of the dam body includes the volume covered by the inner and outer slopes of the dam and the volume covered by the top plane of the dam. The volume covered by the top plane of the dam is obtained based on the horizontal scan line method.

[0060] In specific implementation, the step of generating the thickness profile boundary of the mountain on both ends of the dam based on the obtained DEM data is as follows: the elevation value of each layer is determined according to the preset vertical layering parameters, the minimum value of the dam crest elevation and the inner and outer side boundary elevations; parallel lines are drawn at the elevations of each layer to the center line of the dam crest plane; and the two ends of the parallel lines are extended; the thickness boundary line of the mountain is generated based on the obtained extension lines.

[0061] In practice, the DEM data is point cloud data obtained by scanning with an airborne LiDAR. After filtering, smoothing and removing ground features, DEM data with a sub-decimeter-level regular grid is generated by spatial interpolation.

[0062] Specifically, for ease of understanding, the following detailed description of the solution in this embodiment is provided in conjunction with the accompanying drawings:

[0063] A simulation method for pumped storage power station dam engineering that takes into account terrain includes the following:

[0064] S1: Based on high-precision DEM (Digital Elevation Model) data, spatial data analysis methods are applied to achieve automated location of pumped storage power station dam sites;

[0065] Specifically, in step S1, this embodiment provides a method for one-click intelligent extraction of the centerline of the dam top surface, which includes the following process:

[0066] S1.1: Parameter Preset

[0067] A parameter preset interface was designed and developed. Users can set parameters through this interface based on data characteristics and power plant requirements. These parameters mainly include DEM data source, baseline elevation, contour interval, output contour line name, water catchment threshold, output water catchment line name, maximum dam length (km), and expected reservoir capacity (10,000 m³). 3 The function outputs the reservoir surface name, etc., and then starts the automatic reservoir surface extraction process to achieve one-click reservoir surface extraction. Figure 2 As shown.

[0068] S1.2: Extract contour lines using sub-decimeter level DEM data

[0069] (1) The current method for producing high-precision DEMs in reservoir areas mainly involves using airborne LiDAR to scan and obtain point cloud data. After filtering, smoothing, and removing ground features, spatial interpolation methods are applied to generate DEM data with a sub-decimeter-level regular grid. The process of extracting contour lines from the DEM data is as follows:

[0070] Depending on the size of the reservoir area, the characteristics of the Lidar data, and the accuracy of the reservoir capacity calculation, it is necessary to specify the base elevation and contour interval. Generally, the base elevation is set to the rounded value of the lowest point elevation in the Lidar scan data, and the contour interval is set to 1 decimeter to meet the requirements.

[0071] Using an adaptive search algorithm, contour lines can be automatically generated based on the set base elevation and contour interval parameters, such as... Figure 3 As shown, since this algorithm is not the inventive point of this patent, it will not be described in detail.

[0072] S1.3: Extract the water catchment line using sub-decimeter level DEM data.

[0073] (1) Filling depressions to generate flow direction data

[0074] This embodiment proposes a depression-filling algorithm based on block matrix parallel computation for high-resolution, large-area DEM data. It employs a CPU+GPU parallel data processing method, significantly improving computation speed compared to traditional algorithms. The main process is as follows:

[0075] 1) DEM data partitioning

[0076] To address the challenges of processing large datasets, DEM data is divided into blocks. The number of blocks is determined by the computer's CPU, CPU memory, and number of cores, such as 1024*1024 blocks. However, the number of blocks should not be too large, as this would affect computational efficiency.

[0077] Because the depression detection is performed using 8-directional calculation, two rows or two columns of data need to be overlapped at the adjacent edges of the block to save the block boundary processing with pixels in all 8 directions.

[0078] 2) Depression detection

[0079] In general, depressions account for a very small proportion of the entire DEM data. Therefore, algorithms that quickly detect depressions and then fill them in a concentrated manner are significantly more efficient than traditional pixel-by-pixel iterative depression filling algorithms.

[0080] Then, 8-directional detection is performed on each DEM block data, as follows:

[0081] Leftward direction: DEM block data is denoted as M0. Copy M0 and remove the last column to get M. L0 Copy M0 from column 2 to column 3 to obtain matrix M. L1 M L0 Subtract M L1 The resulting matrix is ​​denoted as M. L .

[0082] Rightward direction: Copy M0 and remove the first column to get M R0 Copy M0 from column 1 to column 2 to obtain matrix M. R1 M R0 Subtract M R1 The resulting matrix is ​​denoted as M. R .

[0083] Top direction: Copy M0 and remove the last line to get M U0 Copy M0 from row 2 to row 3 to obtain matrix M. U1 M U0 Subtract M U1 The resulting matrix is ​​denoted as M. U .

[0084] Downwards: Copy M0 and remove the first line to get M D0 Copy M0 from row 1 to row 2-3 to obtain matrix M. D1 M D0 Subtract M D1 The resulting matrix is ​​denoted as M. D .

[0085] Top left: Copy M0 and remove the last column and last row to get M. LU0 Copy M0 from column 2 to the last column and row 2 to the last row to obtain matrix M. LU1 M L0 Subtract M LU1 The resulting matrix is ​​denoted as M. LU .

[0086] Top right direction: Copy M0 and remove the first column and last row to get M. RU0 Copy M0 from column 1 to column 2 to the last, and row 2 to row 3 to obtain matrix M.RU1 M RU0 Subtract M RU1 The resulting matrix is ​​denoted as M. RU .

[0087] Bottom left: Copy M0 and remove the last column and the first row to get M. LD0 From column 2 to the last column, and row 1 to the second-to-last row, copy M0 to obtain matrix M. LD1 M LD0 Subtract M LD1 The resulting matrix is ​​denoted as M. LD .

[0088] Bottom right direction: Copy M0 and remove the first column and first row to get M RD0 Copy M0 from column 1 to column 2 to the second to last, and from row 1 to row 2 to the second to last, to obtain matrix M. RD1 M RD0 Subtract M RD1 The resulting matrix is ​​denoted as M. RD .

[0089] M L M R M U M D M LU M RU M LD M RD Eight matrices are aligned to form a three-dimensional matrix. Local statistical methods are applied to calculate the maximum value, yielding the statistical matrix M. Max It also records which matrix the flow originates from, indicating the direction of the flow, denoted as matrix M. FD Then for M Max Threshold segmentation is performed, and pixels with values ​​less than or equal to 0 are set to NoData, while other pixel values ​​remain unchanged. This represents the lowest point of the depression, yielding the depression matrix M. De Clear M from memory L M R M U M D M LU M RU M LD M RD Eight matrices free up storage space.

[0090] The data of each DEM block is processed in parallel to obtain the depression matrix and flow direction matrix data of each block.

[0091] 3) Modify the flow direction

[0092] Search for the lowest cell in the depression matrix. If it is an independent cell, check if all 8 of its neighboring cells flow to that depression. If not all of them flow to that depression, then in M... FD The value of the depression is changed to flow out of the neighborhood into 8 neighboring cells and in M Max The cell with the largest median value; if all flows are directed to the depression, it means that its 8 neighboring cells are also depressions. Continue searching the adjacent cells in the outer ring until a cell that does not flow to the depression area is found. Also, search for the cell with the smallest elevation in M0, which is the outlet of the depression. Modify the flow direction of all depressions to the direction with the smallest distance to the outlet.

[0093] If the lowest pixel in the depression is a continuous sheet-like area, continue searching the adjacent pixels on the outer ring until a pixel that does not flow into the depression area is found. Also, search for the pixel with the smallest elevation in M0, which will be the outlet of the depression, and modify the flow direction of all depressions to the direction with the smallest distance to the outlet.

[0094] If a depression is located at the edge of a block, depressions in adjacent blocks will be treated together, so as not to affect the flow direction of the local watershed.

[0095] The depression matrix data of each block is processed in parallel to obtain the final flow direction matrix data of each block.

[0096] (3) Calculation of water catchment

[0097] Based on the final flow matrix M FD The data is used to calculate the cumulative number of pixels that flow into each pixel from its eight neighboring pixels, generating a drainage matrix M. CQ .

[0098] However, since the cumulative amount cannot be calculated in parallel, the cumulative amount must be calculated sequentially (e.g., starting from the top left corner). When calculating the cumulative amount of the current block, the cumulative amount of the overlapping pixels at the edges of adjacent blocks is copied over.

[0099] (4) Drainage line extraction

[0100] Application of the water catchment matrix M CQ Water catchment lines can be easily extracted using a threshold segmentation method. A larger threshold yields coarser-grained catchment lines, while a smaller threshold results in denser lines. Based on reservoir extraction requirements, extensive experiments have demonstrated that setting a uniform threshold of 1 / 10 of the maximum catchment volume for the entire watershed is sufficient for the application. Values ​​greater than or equal to 1 / 10 are assigned a value of 1, and values ​​less than 1 / 10 are assigned a value of 0. This yields a binary matrix M of the water catchment lines. CL .

[0101] A raster-to-vector conversion method is used to generate a surface graphic, denoted as L. C_Lyr, using the method of extracting the centerline of a surface to extract the watercourse, such as Figure 4 As shown, this method is significantly more efficient than the traditional image-by-image element search algorithm.

[0102] S1.4: The perpendicular segment that generates the water catchment line

[0103] The upstream and downstream endpoints of each straight segment of the catchment line can be determined by the amount of water collected at the two endpoints, denoted as P. up and P down At the downstream endpoint P down At that point, draw a perpendicular line segment along the normal direction of the straight line segment, with a length equal to the maximum dam length, such as... Figure 5 As shown, the midpoint of the line segment lies on the watercourse, and the formulas for calculating the plane coordinates of the two endpoints of the perpendicular line segment are as follows: (1) and (2):

[0104] when When using this formula, the following formula is used:

[0105]

[0106] when When using this formula, the following formula is used:

[0107]

[0108] In the formula, x1 and y1 are the plane coordinates of one endpoint of the perpendicular segment, x2 and y2 are the plane coordinates of the other endpoint of the perpendicular segment, θ is the tangent of the slope of the water catchment line, x and y are the plane coordinates of the midpoint of the perpendicular segment (which are also points on the water catchment line), and L is the maximum dam length.

[0109] S1.5: Generate reservoir surface

[0110] To calculate the intersection of a vertical line segment or its extension from the center to both ends with a contour line, assuming the plane coordinates of the two ends of a line segment on the contour line are (x3, y3) and (x4, y4) respectively, the formula for calculating the intersection point is as shown in equation (3):

[0111]

[0112] In the formula, x1 and y1 are the plane coordinates of one endpoint of the perpendicular line segment, x2 and y2 are the plane coordinates of the other endpoint of the perpendicular line segment, and x and y are the coordinates of the intersection point. If the intersection point is located between (x3, y3) and (x4, y4), it is retained; otherwise, it is discarded. The calculation is repeated for all line segments on the contour line. If there are intersection points on both sides of the center point of the perpendicular line segment, then the intersection point (x1, y1) closest to the center point is selected. dc1 y dc1 ) and (x dc2 y dc2The line segment between two points (x, y) is called the intercepting line segment. Calculate the line segment between the two points (x, y). dc1 y dc1 ) and (x dc2 y dc2 The distance D between the two contour lines is retained if it is less than the maximum dam length L; otherwise, it is discarded. The process continues with the next contour line until all contour lines intersecting the vertical line segment or its extension have been processed. The closed area enclosed by the contour line with the highest elevation and its corresponding interception line segment is the reservoir surface. The interception line segment is the horizontal transverse centerline of the dam crest plane, denoted as L. DC0 This contour line is denoted as the uppermost contour line L0 of the reservoir, and its elevation is the elevation of the reservoir surface, which is also the elevation of the top plane of the dam, denoted as H. For example... Figure 6 As shown.

[0113] Repeat the process for the next perpendicular line segment until all perpendicular line segments or their extensions form the reservoir surface and interception line segments.

[0114] S2: Based on parameters such as dam length and inclination angles of inner and outer slopes, the coordinates of the dam boundary line are calculated using a data model to construct the dam's geometric shape, taking into account topographical factors to meet the requirement of building according to the terrain.

[0115] Specifically, in step S2, according to the "Design Code for Pumped Storage Power Stations (NB / T 10072-2018)," based on high-precision DEM data, and according to parameters such as dam length and inner and outer slope inclination angles, the boundary of the pumped storage power station dam is automatically extracted to construct the dam's geometric shape. This includes the following steps:

[0116] S2.1: Generate the dam crest plane

[0117] According to dam design specifications, the width W of the dam crest plane is set, based on the plane coordinates (x, y) of the two endpoints of the horizontal centerline of the dam crest plane. dc1 y dc1 ) and (x dc2 y dc2 The equation for the centerline can be written as shown in equation (4):

[0118]

[0119] Pass (x) dc1 y dc1 ) and (x dc2 y dc2 Draw the center line L DC0 perpendicular line L P1 and L P2 , respectively in L P1 and L P2 Take two points on each side of the center line, located on the perpendicular line L. P1The plane coordinates of two points on the x-axis are (x-y) LP1_i y LP1_1 ) and (x LP1_2 y LP1_2 ), the midpoint is (x d1 y d1 ), located on the perpendicular line L P2 The plane coordinates of two points on the x-axis are (x-y) LP2_1 y LP2_1 ) and (x LP2_2 y LP2_2 ), the midpoint is (x d2 y d2 ), making L P1 The distance between two points on the same surface is W, L P2 The distance between two points on the x-axis is also W. LP1_1 y L1_1 ) and (x LP2_1 y L2_1 The straight line in a plane between two points is denoted as L. DS1 Its equation is as shown in equation (5):

[0120]

[0121] (x LP1_2 y L1_2 ) and (x LP2_2 y L2_2 The straight line in a plane between two points is denoted as L. DS2 Its equation is as shown in equation (6):

[0122]

[0123] L DC0 L DS1 and L Ds2 The three plane lines are parallel lines, L DS1 and L DS2 The elevation of the side line of the dam crest plane is H.

[0124] Calculate L using formula (3) DS1 and L DS2 Contour line L at the uppermost edge of the reservoir C The intersection point can be used to obtain the distance from the midpoint (x). d1 y d1 ) and (x d2 y d2 The four nearest intersections are located at L. DS1 The two points on the graph are denoted as (x) ds1_1 y ds1_1 ) and (x ds1_2 y ds1_2 ), L CThe portion between these two points is the boundary line at one end of the dam, denoted as DEL1; the portion at L... DS2 The two points on the graph are denoted as (x) ds2_1 y ds2_1 ) and (x ds2_2 y ds2_2 ), L C The portion between these two points is the boundary line at the other end of the dam, denoted as DEL2.

[0125] By DEL1, L DS1 DEL2 and L DS2 The top plane of the dam is formed by four parts, as follows: Figure 7 As shown.

[0126] S2.2: Establish the equations for the two inclined planes

[0127] According to dam design specifications, the inclination angle of the inner side of the dam is set as θ. in The outer tilt angle θ ex In general, θ in =θ ex ,

[0128] By L DS1 and θ in The equation for the inner slope of the dam can be obtained from L. DS2 and θ ex The equations for the outer slope of the dam can be obtained. The detailed process for establishing the equations for the inner slope is as follows:

[0129] Through line segment L DS1 midpoint As L DS1 By finding the vertical line, we can determine the intersection point of the vertical line and the horizontal plane. Then pass P ds1_0 With line segment L DS1 Parallel three-dimensional straight lines are denoted as L. ds1_0 Its equation is as shown in equation (7):

[0130]

[0131] Pass P ds1_0 With line segment L ds1_0 A perpendicular three-dimensional straight line is denoted as L. dsp1_0 Its equation is as shown in equation (8):

[0132]

[0133] Find line L dsp1_0 A little bit P dsin_1 Make point P ds1_0 The distance is d = H × ctanθ in And it is located inside the dam. To find P...ds1_1 The coordinates of the point are used to establish the system of equations as shown in equation (9):

[0134]

[0135] Solving the system of equations (9) yields the coordinates of two points. Then, the distance from these two points to point P is calculated. up The distance from P up The point with the smaller distance is the inner point, denoted as P. dsin_1 (x dsin_1 y dsin_1 ,0).

[0136] According to L DS1 Two points (x) on ds1_1 y ds1_1 , z H ) and (x ds1_2 y ds1_2 , z H ) and P dsin_1 (x dsin_1 y dsin_1 The equation of the inner inclined plane can be determined as shown in equation (10): , 0),

[0137] A×(xx ds1_1 )+B×(yy ds1_1 )+C×(zz H )+D=0 (10)

[0138] In the formula B = (x ds1_2 -x ds1_1 )×z H , D = -(A × x ds1 _1+B×y ds1_1 +C×z H )

[0139] In the formula, || represents the determinant operation.

[0140] The detailed process of establishing the equation for the outer inclined plane is as follows:

[0141] Through line segment L DS2 midpoint As L DS2 By finding the vertical line, we can determine the intersection point of the vertical line and the horizontal plane. Then pass P ds2_0 With line segment L DS2 Parallel three-dimensional straight lines are denoted as L. ds2_0 Its equation is as shown in equation (11):

[0142]

[0143] Pass P ds2_0 With line segment L ds2_0 A perpendicular three-dimensional straight line is denoted as L. dsp2_0 Its equation is as shown in equation (12):

[0144]

[0145] Find line L dsp2_0 A little bit P dsex_1 Make point P ds2_0 The distance is d = H × ctanθ ex And it is located inside the dam. To find P... ds2_1 The coordinates of the point are used to establish the system of equations as shown in equation (13):

[0146]

[0147] Solving the system of equations (9) yields the coordinates of two points. Then, the distance from these two points to point P is calculated. down The distance from P down The point with the larger distance is the outermost point, denoted as P. dsex_1 (x dsex_1 y dsex_1 ,0).

[0148] According to L DS2 Two points (x) on ds2_1 y ds2_1 , z H ) and (x ds2_2 y ds2_2 , z H ) and P dsex_1 (x dsex_1 y dsex_1 The equation of the inner inclined plane can be determined as shown in equation (14): , 0),

[0149] A×(xx ds2_1 )+B×(yy ds2_1 )+C×(zz H )+D=0 (14)

[0150] In the formula, A = (y dsex_1 -y ds2_1 )×(-z H ), B = (x ds2_2 -x ds2_1 )×z H , D = -(A × x ds2_1 +B×y ds2_1 +C×z H )

[0151] S2.3: Tracing the boundaries of the two inclined planes

[0152] To achieve the function of building dams according to the terrain, it is necessary to trace the boundary lines of the two slopes based on the DEM terrain features. This embodiment proposes an algorithm for tracing the boundary using the horizontal scan line method. Taking the tracing of the inner side boundary as an example, the specific process is as follows:

[0153] Through line segment L DS1 midpoint Draw an L on the inner side DS1 perpendicular line L DS1_P Then L DS1_P The equation is as shown in equation (15):

[0154]

[0155] In the formula, A, B, C, and D are the same as in equation (10). G = 0,

[0156] From the midpoint P ds1 Along L DS1_P North Korea P up each direction Take a point at intervals (Δx is the size of the DEM cell), and calculate the coordinates of each point according to equation (9). DS1_P1 y DS1_P1 , z DS1_P1 ), ...,(x DS1_Pn y DS1_Pn , z DS1_Pn ). Perform L through each point DS1 Parallel lines are used as scan lines, starting from each point. Points are taken at intervals from both ends of the scan line. The coordinates of the two points can be calculated according to formulas (9) and (16). The row and column numbers of the DEM pixels at the coordinates of the two points can be calculated using formula (16).

[0157]

[0158] In the formula, x and y are the coordinates of the current point, i and j represent the column number and row number of the DEM cell corresponding to the current point, x0 and y0 are the starting coordinates of the DEM data, and Δx and Δy are the X-direction length and Y-direction length of the DEM cell, respectively, which are generally equal.

[0159] The elevation of the terrain at point i and j can be obtained, denoted as z. DS1_P1_DEM If z DS1_P1 >z DS1_P1_DEM This indicates that the side of this point has not yet touched the terrain surface. The volume of this point can be estimated using formula (17):

[0160]

[0161] If z DS1_P1 ≤z DS1_P1_DEM This indicates that the side of the point has already touched the terrain surface. The two intersections of the scan line and the boundary of the pixel are calculated, and the line drawn is the partial boundary line on the side. The volume block at the boundary can be ignored.

[0162] Iterate through all points on the scan line, capturing all boundary lines located on the scan line; process each scan line in the same way. Finally, from L... DS1 Starting from one endpoint, connect all the boundary lines in sequence to form a complete inner boundary line, such as... Figure 8 As shown; the estimated volume blocks are summed up to obtain the volume covered by the inner surface.

[0163] The same method can be used to generate the outer boundary line, such as Figure 8 As shown, the volume under the outer surface cover is calculated, which will not be elaborated further. During the scanning calculation, the minimum elevation of the inner and outer surface boundaries is recorded and denoted as z. DS_Pmin .

[0164] S3: Apply the cumulative calculation algorithm to estimate the volume of the dam body, and then estimate the earthwork required for the construction of the dam body;

[0165] Specifically, in step S3, in addition to estimating the volume under the inner and outer surfaces of the dam body as in S2.3, the volume under the top plane is also estimated. This part of the volume can be estimated using the same scanning line method, which will not be elaborated further.

[0166] S4: To assess the safety of the dam, it is crucial to calculate the thickness of the mountains on both ends of the dam, as this is an important parameter for evaluating the dam's bearing capacity.

[0167] Dam safety is a crucial aspect that must be considered in dam design. Calculating the thickness of the mountainsides on both ends of the dam is an important parameter for assessing the dam's bearing capacity. This process includes the following steps:

[0168] S4.1: Set layer parameters

[0169] To draw the thickness profile of the mountains on both ends of the dam, this embodiment sets vertical layering parameters, such as n layers, based on the layering parameters and the minimum value z of the dam crest elevation H and the inner and outer side boundary elevations. DS_Pmin The elevation value at each layer can be determined, denoted as z. PS_1 , ..., z PS_n The elevation calculation formula is as shown in equation (18):

[0170]

[0171] In the formula, i represents the layer number from top to bottom.

[0172] S4.2: Generate parallel lines

[0173] Draw the centerline L of the dam crest plane at the elevation of each layer. DC0 Parallel lines, including L DC0 There are a total of n+1 parallel lines, and their plane equations are all the same as formula (4).

[0174] Based on the elevation of each parallel line, select the contour line that is closest to that elevation, and extend the lines from the two endpoints of the parallel lines respectively. The first intersection point with the selected contour line is a boundary point of the profile line. The coordinates of the intersection point can be obtained by applying formula (3).

[0175] Calculate the intersection points of the extended line with other contour lines in a loop until the elevation of the intersection point is lower than the elevation of the first intersection point.

[0176] S4.3: Generate mountain thickness boundary lines

[0177] By connecting the intersection points in sequence, the boundary lines of the mountain thickness profile are formed. A mountain thickness boundary line is generated at each end of the dam. At the same time, the thickness of each layer of the mountain can be calculated, as shown in Figure 9(a) and Figure 9(b).

[0178] S5: Using the constructed geometric coordinates of the dam body, apply the automatic spatial 3D modeling method to generate a 3D model of the dam body, and display the simulation effect in a 3D scene.

[0179] Step S5 specifically involves: using the constructed geometric coordinates of the dam body, and according to the specifications and methods of 3D model building, automatically constructing a 3D model of the dam body, generating a 3D model file in the .obj universal format, and displaying the simulation effect in a 3D scene. Wherein:

[0180] S5.1: Create a set of points

[0181] (1) Generate the set of coordinates of the points

[0182] Based on the top plane coordinates, two side boundary coordinates, and bottom point coordinates obtained from processes S2.1, S2.3, and S3 above, a coordinate set of points for the dam's 3D model is created, with the following structure:

[0183] v...x1, y1, z1

[0184] v...x², y², z²

[0185] v...x3, y3, z3

[0186] v...x4, y4, z4

[0187] v...x5, y5, z5

[0188] v...x6, y6, z6

[0189] v...x7, y7, z7

[0190] v...x8, y8, z8

[0191] Where V refers to the point coordinates of the 3D model, and each coordinate includes three values: x, y, and z.

[0192] (2) Normal set of the generated points

[0193] 1) Algorithm for determining the normal direction

[0194] Let the average point be c and the normal vector be n. The problem can be transformed into finding a direction n such that the distribution of the projection points of all neighborhood points on direction n is most concentrated. Such an optimization problem is shown in equation (19):

[0195]

[0196] Clearly, we seek a direction n such that the distribution of all neighborhood points projected onto direction n is most concentrated, which means that the variance of the projection of a point onto that direction is minimized.

[0197] Choose c as the mean of this group of points, then x i -c represents the decentralized point. Therefore, the problem becomes the following problem, as shown in equation (20):

[0198]

[0199] Expanding the parentheses and simplifying, we get the result, as shown in equation (21):

[0200]

[0201] The n we are looking for is The eigenvector corresponding to the smallest eigenvalue.

[0202] 2) The set of normals of the generated points

[0203] For each point in the coordinate set generated in step (1), the algorithm for determining the normal is applied to obtain the following set of normals:

[0204] vn...x1,y1,z1

[0205] vn...x2,y2,z2

[0206] vn...x3,y3,z3

[0207] vn...x4,y4,z4

[0208] vn...x5,y5,z5

[0209] vn...x6,y6,z6

[0210] vn...x7,y7,z7

[0211] vn...x8,y8,z8

[0212] Here, vn represents the normal direction.

[0213] S5.2: Face Set

[0214] In the half-edge data structure, mesh edges are explicitly stored by representing each edge with a pair of directed half-edge twins, each pointing in opposite directions. A half-edge stores a reference to its twin, as well as references to the preceding and following half-edges along the same face or hole. Vertices store their position and references to any half-edges originating from that vertex, while faces store any half-edges belonging to that face. The half-edge data structure stores an array of vertex, face, and half-edge records.

[0215] The top plane, front face, rear face, and bottom face are sorted clockwise according to the half-side structure, as follows: Figure 10 The resulting dam model conforming to the half-side structure is as follows: Figure 10 As shown.

[0216] #Top Plane

[0217] f 1 2 3 4

[0218] #Front-end

[0219] f 4 3 5 6

[0220] #Backend

[0221] f 2 1 7 8

[0222] #bottom

[0223] f 2 8 7 1 4 6 5 3

[0224] S5.3: 3D Simulation Visualization

[0225] The generated Obj format model was converted into model data in a format supported by SuperMap software. A 3D simulation visualization system was developed using iSuperObject, which integrates basic data such as DEM, contour lines, and reservoir surface into the 3D scene, and can display the 3D model of the dam in 360 degrees, as shown in Figure 11(a) and Figure 11(b).

[0226] Example 2:

[0227] The purpose of this embodiment is to provide a simulation system for pumped storage power station dam projects that takes into account the terrain.

[0228] A simulation system for pumped storage power station dam engineering that takes into account terrain features includes:

[0229] The data acquisition unit is used to acquire DEM data of the reservoir area where the pumped storage power station dam is located.

[0230] The reservoir surface and centerline generation unit is used to generate the horizontal centerline of the reservoir surface and dam crest plane based on the DEM data using spatial data analysis methods. The spatial data analysis method includes: extracting contour lines and runoff lines from the DEM data of the reservoir area; generating perpendicular segments of the runoff lines based on the obtained runoff lines; calculating the intersection points of the perpendicular segments and contour lines, and using the line connecting the two sides of the contour lines closest to the center point of the contour line as the interception line segment; and determining the horizontal centerline of the reservoir surface and dam crest plane based on the interception line segment and the contour lines.

[0231] The geometric shape construction unit is used to generate the dam top plane based on the obtained horizontal centerline of the dam top plane and the preset width of the dam top plane; and to construct the geometric shape of the dam body based on the constructed equations of the dam's two sides and the obtained boundaries of the dam's two sides, and to obtain the volume of the dam body.

[0232] The terrain generation unit is used to generate the thickness profile boundary of the mountains on which the dam is supported based on the obtained DEM data.

[0233] The simulation unit is used to construct a dam model based on the obtained reservoir surface, dam geometry, dam volume, and the thickness profile boundary of the mountains on both ends of the dam, and to perform pumped storage power station dam engineering simulation based on the constructed dam model.

[0234] Furthermore, the system described in this embodiment corresponds to the method described in Embodiment 1, and its technical details have been described in detail in Embodiment 1, so they will not be repeated here.

[0235] In further embodiments, the following is also provided:

[0236] An electronic device includes a memory and a processor, as well as computer instructions stored in the memory and running on the processor. When executed by the processor, the computer instructions perform the method described in Embodiment 1. For brevity, further details are omitted here.

[0237] It should be understood that in this embodiment, the processor can be a central processing unit (CPU), or it can be other general-purpose processors, digital signal processors (DSPs), application-specific integrated circuits (ASICs), programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. The general-purpose processor can be a microprocessor or any conventional processor, etc.

[0238] Memory may include read-only memory and random access memory, and provides instructions and data to the processor. A portion of memory may also include non-volatile random access memory. For example, memory may also store information about the device type.

[0239] A computer-readable storage medium for storing computer instructions, which, when executed by a processor, perform the method described in Embodiment 1.

[0240] The method in Embodiment 1 can be directly implemented by a hardware processor, or implemented by a combination of hardware and software modules within the processor. The software modules can reside in readily available storage media in the art, such as random access memory, flash memory, read-only memory, programmable read-only memory, electrically erasable programmable memory, or registers. This storage medium is located in memory; the processor reads information from the memory and, in conjunction with its hardware, completes the steps of the above method. To avoid repetition, a detailed description is not provided here.

[0241] Those skilled in the art will recognize that the units, i.e., algorithm steps, of the various examples described in connection with this embodiment can be implemented in electronic hardware or a combination of computer software and electronic hardware. Whether these functions are implemented in hardware or software depends on the specific application and design constraints of the technical solution. Those skilled in the art can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of this disclosure.

[0242] The above-described embodiment provides a simulation method and system for pumped storage power station dam engineering that takes into account terrain, which is feasible and has broad application prospects.

[0243] The above description is merely a preferred embodiment of this disclosure and is not intended to limit this disclosure. Various modifications and variations can be made to this disclosure by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this disclosure should be included within the scope of protection of this disclosure.

Claims

1. A simulation method for pumped storage power station dam engineering that takes into account terrain, characterized in that, include: Obtain DEM data of the reservoir area where the pumped storage power station dam is located; Based on the DEM data, a spatial data analysis method is used to generate the horizontal centerline of the reservoir surface and the dam crest plane; wherein, the spatial data analysis method is as follows: extracting contour lines and water catchment lines based on the DEM data of the reservoir area; generating vertical segments of the water catchment lines based on the obtained water catchment lines; Calculate the intersection of the vertical line segment and the contour line, and use the line connecting the two sides of the contour line closest to the center point of the contour line as the interception line segment; determine the horizontal centerline of the reservoir surface and the dam crest plane based on the interception line segment and the contour line. Based on the obtained horizontal centerline of the dam crest plane and the preset width of the dam crest plane, the dam crest plane is generated; and based on the constructed equations of the dam's two sides and the obtained boundaries of the dam's two sides, the geometry of the dam body is constructed, and the volume of the dam body is obtained. Based on the obtained DEM data, the thickness profile boundary of the mountains on which the dam is supported is generated; Based on the obtained reservoir surface, dam geometry, dam volume, and the thickness profile boundaries of the mountains on both ends of the dam, a dam model is constructed, and a pumped storage power station dam engineering simulation is performed based on the constructed dam model.

2. The simulation method for pumped storage power station dam engineering considering terrain as described in claim 1, characterized in that, The horizontal centerline of the reservoir surface and the dam crest plane is determined based on the interception line segment and contour line. Specifically, the closed area enclosed by the contour line with the highest elevation and its corresponding interception line segment is selected as the reservoir surface, and the interception line segment is used as the horizontal centerline of the dam crest plane.

3. The simulation method for pumped storage power station dam engineering considering terrain as described in claim 1, characterized in that, The generation of the dam top plane specifically involves: determining a first plane line and a second plane line parallel to the horizontal center line of the dam top plane, wherein the first plane line and the second plane line are located on both sides of the center line; determining the edge lines at both ends of the dam based on the intersection points of the first plane line and the second plane line with the contour line at the uppermost edge of the reservoir; and determining the top plane of the dam based on the edge lines at both ends of the dam, the first plane line, and the second plane line.

4. The simulation method for pumped storage power station dam engineering considering terrain as described in claim 1, characterized in that, The boundaries of the sloping sides of the dam are obtained by using the horizontal scan line method to trace the boundaries based on the constructed equations of the sloping sides of the dam, thereby obtaining the boundaries of the inner and outer sloping sides of the dam and the volume covered by the inner and outer sloping sides of the dam.

5. The simulation method for pumped storage power station dam engineering considering terrain as described in claim 1, characterized in that, The volume of the dam body includes the volume covered by the inner and outer slopes of the dam and the volume covered by the top plane of the dam, wherein the volume covered by the top plane of the dam is obtained based on the horizontal scan line method.

6. The simulation method for pumped storage power station dam engineering considering terrain as described in claim 1, characterized in that, The process of generating the thickness profile boundary of the mountain body on which the two ends of the dam rely based on the obtained DEM data is as follows: the elevation value at each layer is determined based on the preset vertical layering parameters, the elevation of the dam top surface, and the minimum value of the elevation of the inner and outer side boundary. Draw a line parallel to the centerline of the dam crest plane at the elevation of each layer; and extend both ends of the parallel line. The mountain thickness boundary line is generated based on the obtained extension line.

7. The simulation method for pumped storage power station dam engineering considering terrain as described in claim 1, characterized in that, The DEM data is point cloud data obtained by scanning with an airborne LiDAR. After filtering, smoothing and removing ground features, DEM data with a sub-decimeter-level regular grid is generated by spatial interpolation.

8. A simulation system for pumped storage power station dam engineering that takes into account terrain, characterized in that, include: The data acquisition unit is used to acquire DEM data of the reservoir area where the pumped storage power station dam is located. The reservoir surface and centerline generation unit is used to generate the horizontal centerline of the reservoir surface and dam crest plane based on the DEM data using spatial data analysis methods. The spatial data analysis method includes: extracting contour lines and runoff lines from the DEM data of the reservoir area; generating perpendicular segments of the runoff lines based on the obtained runoff lines; calculating the intersection points of the perpendicular segments and contour lines, and using the line connecting the two sides of the contour lines closest to the center point of the contour line as the interception line segment; and determining the horizontal centerline of the reservoir surface and dam crest plane based on the interception line segment and the contour lines. The geometric shape construction unit is used to generate the dam top plane based on the obtained horizontal centerline of the dam top plane and the preset width of the dam top plane; and to construct the geometric shape of the dam body based on the constructed equations of the dam's two sides and the obtained boundaries of the dam's two sides, and to obtain the volume of the dam body. The terrain generation unit is used to generate the thickness profile boundary of the mountains on which the dam is supported based on the obtained DEM data. The simulation unit is used to construct a dam model based on the obtained reservoir surface, dam geometry, dam volume, and the thickness profile boundary of the mountains on both ends of the dam, and to perform pumped storage power station dam engineering simulation based on the constructed dam model.

9. An electronic device, comprising a memory, a processor, and a computer program stored in the memory and running thereon, characterized in that, When the processor executes the program, it implements a simulation method for pumped storage power station dam engineering that takes into account the terrain, as described in any one of claims 1-7.

10. A non-transitory computer-readable storage medium having a computer program stored thereon, characterized in that, When executed by the processor, the program implements a simulation method for pumped storage power station dam engineering that takes into account the terrain, as described in any one of claims 1-7.

Images