Pumped storage reservoir leakage parameterized investigation and arrangement method based on reservoir type classification

Abstract

The invention provides a pumped storage reservoir leakage parameterization investigation and arrangement method based on reservoir type classification, and the method comprises the steps: dividing a pumped storage power station reservoir into a plurality of reservoir types according to the position of the reservoir and the landform characteristics of a reservoir region; setting reservoir leakage investigation key parts; outputting a parameterized reservoir leakage exploration and drilling workload scheme by adopting a reservoir leakage exploration and drilling scheme parameterization algorithm according to the reservoir type, topographic features and reservoir leakage exploration key parts of each reservoir type; training and identifying library type classification, and constructing a library type data source; and carrying out reservoir leakage parameterized investigation arrangement based on the reservoir type. According to the invention, for leakage investigation elements emphasized during the investigation and design of the pumped storage power station, pumped storage reservoir type classification and leakage parameterization investigation are combined through intelligent identification and parameterization arrangement rules, the investigation and arrangement efficiency is improved, manual field climbing point selection is greatly reduced, and the working efficiency is improved. And the efficiency and the accuracy of reservoir leakage investigation arrangement work are improved.

Description

Technical Field

[0001] The invention belongs to the technical field of investigation and arrangement of leakage problems in pumped storage power station reservoirs, and in particular relates to a parameterized investigation and arrangement method for leakage of pumped storage reservoirs based on reservoir type classification. Background Art

[0002] Engineering survey is an important link in the early stage of identifying the engineering geological conditions in the construction area of pumped storage power stations and analyzing and evaluating engineering geological problems. Among them, reservoir leakage is the main content of engineering survey.

[0003] The site selection of pumped storage power stations is relatively flexible, and there are many types of reservoirs. In the existing technology, when arranging reservoir leakage surveys for different types of reservoirs, the on-site layout of survey drilling holes mostly relies on manual field climbing in mountainous areas. It is difficult to reasonably optimize the drilling hole positions, which is not only time-consuming and labor-intensive, but also poses a safety hazard.

[0004] Therefore, how to achieve scientific, standardized and efficient survey layout has become a technical problem that needs to be solved urgently. Summary of the invention

[0005] In view of the defects of the prior art, the present invention provides a parameterized leakage survey and layout method for pumped storage reservoirs based on reservoir type classification, which can effectively solve the above problems.

[0006] The technical solution adopted by the present invention is as follows:

[0007] The present invention provides a pumped storage reservoir leakage parameterization survey arrangement method based on reservoir type classification, comprising the following steps:

[0008] Step S1, dividing the pumped storage power station reservoir into multiple reservoir types according to the reservoir location and the topographical features of the reservoir area;

[0009] Step S2, setting the key parts of reservoir leakage investigation for each reservoir type;

[0010] Step S3, according to the reservoir type, topographic features and key parts of reservoir leakage investigation of each reservoir type, a parameterized algorithm for reservoir leakage investigation drilling scheme is adopted to output a parameterized reservoir leakage investigation drilling workload scheme;

[0011] Step S4, training and identifying library type classification, and building library type data source:

[0012] Based on 3D images, orthophotos, oblique images, and satellite images, various reservoir types of pumped storage power stations are marked and the reservoir type data source is trained;

[0013] Step S5, performing parameterized survey and arrangement of reservoir leakage based on the reservoir type:

[0014] Based on the reservoir image data, the reservoir type is intelligently identified using the described reservoir type data source, and with the reservoir dam axis, the reservoir perimeter boundary, the normal storage level elevation value, and the topographic line extracted from the reservoir image data, combined with the key parts for reservoir seepage investigation in step S2, the reservoir seepage investigation boreholes are parametrically arranged in the manner of step S3.

[0015] Preferably, the reservoir types include valley-type reservoirs, gully-type reservoirs, shore-type reservoirs, gully-source-type reservoirs, and mountaintop terrace-type reservoirs;

[0016] The valley-type reservoir refers to a reservoir formed by building a dam across the main river in the project area;

[0017] The gully-type reservoir refers to a reservoir formed by building a dam in a tributary gully of a river, and the reservoir backwater line does not reach the gully source;

[0018] The shore-type reservoir refers to a reservoir formed by a combination of excavation and filling and enclosing with a dam on a river bank terrace or concave plot, without encroaching on the river bed;

[0019] The gully-source-type reservoir refers to a reservoir formed by building a dam within the scope of a gully source;

[0020] The mountaintop terrace-type reservoir refers to a reservoir formed by a combination of excavation and filling and enclosing with a dam on a relatively flat mountaintop terrace in the project area.

[0021] Preferably, the key parts for reservoir seepage investigation of each reservoir type refer to setting the key parts for reservoir seepage investigation of each reservoir type according to the topographic and geomorphic type characteristics, engineering geology and hydrogeological conditions, and reservoir type characteristics of the pumped storage power station reservoir area.

[0022] Preferably, the key parts for reservoir seepage investigation of each reservoir type specifically include:

[0023] For the valley-type reservoir, the key parts for reservoir seepage investigation are the front dam reservoir section, river bends, and interfluve areas;

[0024] For the shore-type reservoir and the mountaintop terrace-type reservoir, the key parts for reservoir seepage investigation are the reservoir bottom and the reservoir perimeter area;

[0025] For the gully-type reservoir and the gully-source-type reservoir, the key parts for reservoir seepage investigation are the front dam reservoir section, the thin watershed, and the low saddle areas.

[0026] Preferably, step S3 is specifically as follows:

[0027] Step S31, parametrically investigate the boundary of the investigation area:

[0028] According to the axis of the reservoir type dam, the boundary of the reservoir perimeter, and the normal storage level elevation value, generate the boundary of the exploration area based on the topographic features of the reservoir type and the distance conditions around the reservoir type; wherein, the distance conditions around the reservoir type refer to taking the dam axis as the fixed boundary, and combining the selected reservoir perimeter boundary or the normal storage level elevation line by the user to independently set the outward expansion distance; the topographic features of the reservoir type refer to the terrain slope, aspect, and terrain curvature within the area of the distance conditions around the reservoir type that meet the distance conditions around the reservoir type.

[0029] Step S32, parameterize the drilling area of the key parts for reservoir leakage exploration:

[0030] According to the key parts for reservoir leakage exploration of each reservoir type in step S2, combined with the topographic feature factors of the exploration area, set the drilling area of the key parts for reservoir leakage exploration in a parameterized manner.

[0031] Step S33, parameterize the drilling layout:

[0032] The overlapping area of the exploration area set in step S31 and the drilling area of the key parts for reservoir leakage exploration set in step S32 is used as the boundary of the drilling layout area; within the boundary of the drilling layout area, parameterize the drilling hole spacing, hole position, and drilling depth, and intelligently set the drilling workload.

[0033] Preferably, step S32 specifically includes:

[0034] Step S32-1, generate the drilling area of the key parts for reservoir leakage exploration according to the reservoir type and the key parts for reservoir leakage exploration of each reservoir type in step S2, including:

[0035] For the front dam reservoir section, combine the reservoir perimeter line or the normal storage level elevation line with the river valley area with an independently set distance condition upstream of the reservoir dam axis as the drilling area of the front dam reservoir section;

[0036] For the river bend plot, according to the terrain curvature of the river valley, when the radius of curvature is less than 5 times the width of the normal storage level river channel, set it in the area with a terrain slope less than 25° and relatively gentle terrain of the river bend as the drilling area of the river bend plot;

[0037] For the interfluve plot, according to the terrain data between adjacent river channels, set it in the area with a terrain slope less than 25° and relatively gentle terrain between the two river channels as the drilling area of the interfluve plot;

[0038] For the reservoir bottom, use the internal range area of the exploration area of the reservoir dam axis and the reservoir perimeter line as the drilling area of the reservoir bottom;

[0039] For the reservoir perimeter, use the reservoir perimeter line area as the drilling area of the reservoir perimeter;

[0040] For a thin watershed, according to the reservoir type topographic features of the surveyed area in step S31, identify the ridge width, ridge peak, and pass value at the normal storage level elevation or the elevation around the reservoir. If the ridge width around the reservoir is lower than the self-set width value, it is a thin watershed, and the ridge peak and pass areas are used as the drilling areas for the thin watershed;

[0041] For the low pass area, use the gently sloping area of the pass on the ridge line of the non-thin watershed around the reservoir within the area that meets the reservoir type perimeter distance condition in the surveyed area of step S31 as the drilling area for the low pass area;

[0042] Step S32-2: Based on the drilling areas for the key parts of the reservoir leakage investigation generated in step S32-1, use the spatial geometric analysis algorithm to generate the boundaries of the drilling areas for the key parts of the reservoir leakage investigation.

[0043] Preferably, step S33 specifically includes:

[0044] Step S33-1: Parameterize the drilling hole spacing:

[0045] Within the boundary of the drilling layout area, self-set the drilling hole spacing and parameterize it;

[0046] Step S33-2: Parameterize the initial drilling hole positions:

[0047] Combined with the boundary of the drilling layout area and the drilling hole spacing, optimize and select the initial drilling hole positions, including:

[0048] (1) For the reservoir section in front of the dam and the reservoir bottom part including the dam axis, that is, the area including the dam axis factor, take the center point of the dam axis as the starting point, generate an arrangement line perpendicular to the dam axis through vector operation, and then, according to the drilling hole spacing condition and the boundary condition of the drilling layout area, gradually parameterize to generate the arrangement lines parallel to the dam axis and the arrangement lines perpendicular to the dam axis; generate grid intersection points for each horizontal and vertical arrangement line, and each of these grid intersection points is an initial drilling hole position; moreover, initial drilling hole positions are also generated on the dam axis section according to the drilling hole spacing condition and the boundary condition of the drilling layout area;

[0049] (2) For the river bend and the interfluve area, and the reservoir bottom part without the dam axis, that is, the area without the dam axis factor, take the center point of the fixed-point area by the spatial geometric algorithm, and generate positive grid lines according to the drilling hole spacing and the boundary condition of the drilling layout area. The intersection points of each grid line are the initial drilling hole positions;

[0050] (3) For the reservoir perimeter part, that is, the multi-segment lines around the reservoir of the shore type and the mountaintop platform type, take the starting point of the multi-segment line as the drilling starting point, and generate each initial drilling hole position according to the drilling hole spacing;

[0051] (4) For the thin watershed and low pass areas, the initial drilling hole positions are arranged at the peak and the ridge pass;

[0052] Meanwhile, according to the self-set parameterized arrangement sequence of the initial drilling hole positions and the condition of the repeated arrangement spacing of the hole positions, the said initial drilling hole positions are filtered to obtain the filtered initial drilling hole positions;

[0053] Step S33-3, optimize the said initial drilling hole positions:

[0054] For the said initial drilling hole positions, adopt adaptive weight optimization by combining the optimization boundary distance threshold of the initial drilling hole positions and the terrain slope, screen out the optimized drilling hole positions that are closest to the initial drilling hole positions, have the most suitable terrain slope and are convenient for drilling arrangement, and identify the elevation values of the optimized drilling hole positions;

[0055] Step S33-4, parameterize the drilling depth value:

[0056] Self-set the drilling depth value under the drilling depth distance condition, where the drilling depth distance condition includes the underlying distance value below the top boundary of the relatively impermeable layer and the underlying distance value below the groundwater level. The distance condition options can be screened and selected, and the drilling depth distance value is self-set for each option; if the drilling depth distance condition factor is not available, self-set the drilling depth value;

[0057] Step S33-5, parameterize the drilling attribute information:

[0058] After parameterizing the drilling hole positions, the elevation values of the drilling hole positions and the drilling depth values, parameterize the drilling number, drilling elevation, and drilling depth attribute information; among them, the drilling number is determined in sequence according to the parameterized drilling arrangement sequence and the sequence of the hole position intersection points. The sequence is defaulted to the terrain polyline node sequence, and the bottom of the reservoir, river bends and interfluvial voids are defaulted to north first, south second, west first, east second, left first, right second, and the front dam section is defaulted to the dam axis first and the upstream dam front area second.

[0059] Preferably, in step S33-3, the said adaptive weight optimization is:

[0060] F(x,y) = w1·D(x,y) + w2·S(x,y)

[0061]

[0062] Where:

[0063] F(x,y) is the objective function, D(x,y) is the normalized distance score between the candidate point and the initial hole position, S(x,y) is the normalized slope score of the candidate point, w1 and w2 are weight coefficients, which are self-set or set as dynamic weight coefficients, satisfying w1 + w2 = 1;

[0064] (x0, y0) is the initial hole position coordinate, and d is the threshold distance for optimizing the boundary of the initial hole position;

[0065] (x, y) is the candidate point coordinate; Slope(x,y) is the terrain slope value of the candidate point, and Slope min is the minimum terrain slope value within the threshold distance range of the boundary, and Slope max is the maximum terrain slope value within the threshold distance range of the boundary.

[0066] Preferably, when parametrically arranging the reservoir leakage exploration boreholes, the elevation and depth attributes of the boreholes are parametrically displayed, and the borehole arrangement can be adjusted independently.

[0067] A parametric exploration layout method for pumped-storage reservoir leakage based on reservoir type classification provided by the present invention has the following advantages:

[0068] (1) The reservoir type classification extracted by the present invention is based on conditions such as terrain, rivers, and valleys. Using this classification method can improve the accuracy of subsequent parametric exploration layout for reservoir leakage;

[0069] (2) Aiming at the leakage exploration elements emphasized during the exploration and design of pumped-storage power stations, the present invention combines the pumped-storage reservoir type classification with parametric leakage exploration through intelligent recognition and parametric layout rules, improving the efficiency of exploration layout, greatly reducing manual field climbing and point selection, and enhancing the efficiency and accuracy of reservoir leakage exploration layout work.

[0070] (3) Regarding the borehole layout, based on multiple factors such as terrain slope, curvature, and borehole spacing, the present invention introduces a multi-objective optimization algorithm to ensure the rationality and accuracy of borehole layout, facilitating the solution of the difficult problem of borehole positioning in steep terrains; the present invention not only improves the efficiency of borehole positioning, but also optimizes the number and position of boreholes, effectively reducing the exploration cost and improving the exploration quality. Description of the Drawings

[0071] Figure 1 is a flowchart of a parametric exploration layout method for pumped-storage reservoir leakage based on reservoir type classification provided by the present invention;

[0072] Figure 2 is a schematic image of a certain pumped-storage upper reservoir in this embodiment;

[0073] Figure 3 is a schematic layout diagram of a certain pumped-storage upper reservoir in this embodiment (including borehole attributes);

[0074] Figure 4 is a schematic terrain layout diagram of a certain pumped-storage upper reservoir in this embodiment (including borehole attributes);

[0075] Figure 5 Schematic diagram of an image of a lower reservoir of a pumped - storage power station according to this embodiment;

[0076] Figure 6 Schematic layout diagram of an image of a lower reservoir of a pumped - storage power station according to this embodiment (including borehole attributes);

[0077] Figure 7 Schematic topographic layout diagram of a lower reservoir of a pumped - storage power station according to this embodiment (including borehole attributes). Detailed implementation manners

[0078] In order to make the technical problems, technical solutions and beneficial effects solved by the present invention clearer and more understandable, the present invention will be further described in detail below with reference to the accompanying drawings and embodiments. It should be understood that the specific embodiments described herein are only used to explain the present invention and are not used to limit the present invention.

[0079] The present invention provides a parametric exploration layout method for leakage of pumped - storage reservoirs based on reservoir type classification. This method can conduct targeted exploration layouts according to the reservoir type, improve the accuracy and reliability of exploration results, and at the same time reduce exploration costs and time costs.

[0080] Refer to Figure 1 , the present invention provides a parametric exploration layout method for leakage of pumped - storage reservoirs based on reservoir type classification, including the following steps:

[0081] Step S1: According to the location of the reservoir and the topographic and geomorphic characteristics of the reservoir area, divide the pumped - storage power station reservoir into multiple reservoir types;

[0082] Specifically, the reservoir types include valley - type reservoirs, gully - type reservoirs, shore - type reservoirs, gully - source - type reservoirs, and mountaintop - terrace - type reservoirs;

[0083] The valley - type reservoir refers to a reservoir formed by building a dam across the main river in the project area;

[0084] The gully - type reservoir refers to a reservoir formed by building a dam in a tributary gully of a river, and the reservoir backwater line does not reach the gully source;

[0085] The shore - type reservoir refers to a reservoir formed by a combination of excavation and filling to build a dam on a river - bank terrace or concave land plot without occupying the river bed;

[0086] The gully - source - type reservoir refers to a reservoir formed by building a dam within the range of a gully source;

[0087] The mountaintop - terrace - type reservoir refers to a reservoir formed by a combination of excavation and filling to build a dam on a relatively flat mountaintop terrace in the project area.

[0088] Step S2: Set the key parts for reservoir leakage exploration of each reservoir type;

[0089] The key areas for reservoir seepage investigation for each reservoir type refer to the key areas for reservoir seepage investigation for each reservoir type set according to the topographic and geomorphic type characteristics, engineering geology and hydrogeological conditions of the reservoir area of the pumped-storage power station, as well as the characteristics of the reservoir type.

[0090] The key areas for reservoir seepage investigation for each reservoir type specifically include:

[0091] For the valley-type reservoir, the key areas for reservoir seepage investigation are the reservoir section in front of the dam, river bends and interfluve areas;

[0092] For the shore-type reservoir and the mountaintop platform-type reservoir, the key areas for reservoir seepage investigation are the reservoir bottom and the surrounding area of the reservoir;

[0093] For the gully-type reservoir and the gully-source-type reservoir, the key areas for reservoir seepage investigation are the reservoir section in front of the dam, the thin watershed and the low saddle area.

[0094] Step S3: According to the reservoir type, topographic characteristics of each reservoir type and the key areas for reservoir seepage investigation, adopt the parametric algorithm for the reservoir seepage investigation drilling plan to output the parametric reservoir seepage investigation drilling workload plan;

[0095] Step S3 is specifically as follows:

[0096] Step S31: Parametrize the boundary of the investigation area:

[0097] According to the dam axis of the reservoir type, the reservoir perimeter boundary and the normal storage level elevation value, generate the boundary of the investigation area based on the topographic characteristics of the reservoir type and the reservoir perimeter distance conditions; among them, the reservoir perimeter distance conditions refer to taking the dam axis as the fixed boundary, combining the reservoir perimeter boundary or the normal storage level elevation line selected by the user, and independently setting the outward expansion distance; the topographic characteristics of the reservoir type refer to the terrain slope, slope direction and terrain curvature within the area of the reservoir perimeter distance conditions that meet the reservoir perimeter distance conditions.

[0098] Step S32: Parametrize the drilling area for the key areas of reservoir seepage investigation:

[0099] According to the key areas for reservoir seepage investigation of each reservoir type in Step S2, combined with the topographic characteristics factors of the investigation area, set the parametric drilling area for the key areas of reservoir seepage investigation;

[0100] Step S32 specifically includes:

[0101] Step S32-1: According to the reservoir type and the key areas for reservoir seepage investigation of each reservoir type in Step S2, generate the drilling area for the key areas of reservoir seepage investigation, including:

[0102] For the reservoir section in front of the dam, in combination with the reservoir perimeter line or the elevation line of the normal storage level, the river valley area with a self-set distance upstream of the reservoir dam axis is used as the drilling area for the reservoir section in front of the dam;

[0103] For the river bend plot, according to the curvature of the river valley topography, where the radius of curvature is less than 5 times the width of the river channel at the normal storage level, it is set in the area with a gentle terrain slope less than 25° of the river bend topography, which is used as the drilling area for the river bend plot;

[0104] For the plot between rivers, according to the topographic data between adjacent river channels, it is set in the area with a gentle terrain slope less than 25° between the two river channels, which is used as the drilling area for the plot between rivers;

[0105] For the reservoir bottom, the internal range area of the exploration area between the reservoir dam axis and the reservoir perimeter line is used as the drilling area for the reservoir bottom;

[0106] For the reservoir perimeter, the reservoir perimeter line area is used as the drilling area for the reservoir perimeter;

[0107] For the thin watershed, according to the reservoir type topographic features of the exploration area in step S31, identify the ridge width, ridge peak, and pass value at the normal storage level elevation or reservoir perimeter elevation. The reservoir perimeter ridge width lower than the self-set width value is the thin watershed, and the ridge peak and pass parts are used as the drilling area for the thin watershed;

[0108] For the low pass part, the gentle terrain area of the pass on the ridge line of the non-thin watershed of the reservoir perimeter within the reservoir type perimeter distance condition area that meets the reservoir type perimeter distance condition in the exploration area of step S31 is used as the drilling area for the low pass part;

[0109] Step S32 - 2, based on the drilling area of the key parts for reservoir leakage exploration generated in step S32 - 1, use the spatial geometric analysis algorithm to generate the boundary of the drilling area of the key parts for reservoir leakage exploration.

[0110] Step S33, parametric layout of drilling:

[0111] The overlapping area of the exploration area set in step S31 and the drilling area of the key parts for reservoir leakage exploration set in step S32 is used as the boundary of the drilling layout area; within the boundary of the drilling layout area, parametrically set the drilling spacing, drilling hole positions, and drilling depths, and intelligently set the drilling workload.

[0112] Step S33 specifically includes:

[0113] Step S33 - 1, parametric drilling spacing:

[0114] Within the boundary of the drilling layout area, independently set the drilling spacing and parameterize it;

[0115] Step S33-2, parameterize the preliminary drilling hole positions:

[0116] Combined with the boundary of the drilling layout area and the drilling hole spacing, optimize and select the preliminary drilling hole positions, including:

[0117] (1) For the pre-dam reservoir section and the reservoir bottom part containing the dam axis, that is, the area containing the dam axis factor, take the center point of the dam axis as the starting point, generate a layout line perpendicular to the dam axis through vector operation, and then, according to the drilling hole spacing condition and the boundary condition of the drilling layout area, gradually parameterize to generate a layout line parallel to the dam axis and a layout line perpendicular to the dam axis; generate grid intersection points for each horizontal and vertical layout line, and each of the grid intersection points is a preliminary drilling hole position; moreover, preliminary drilling hole positions are also generated on the dam axis section according to the drilling hole spacing condition and the boundary condition of the drilling layout area;

[0118] (2) For river bends, interfluve areas, and reservoir bottom parts without the dam axis, that is, areas without the dam axis factor, use spatial geometry to determine the center point of the fixed-point area, and generate positive grid lines according to the drilling hole spacing and the boundary condition of the drilling layout area. The intersection points of each grid line are the preliminary drilling hole positions;

[0119] (3) For the reservoir perimeter part, that is, the multi-segment lines around the shore-type and mountaintop terrace-type reservoirs, take the starting point of the multi-segment line as the drilling starting point, and generate each preliminary drilling hole position according to the drilling hole spacing;

[0120] (4) For thin watersheds and low saddle parts, arrange the preliminary drilling hole positions at the peaks and ridge saddles;

[0121] Meanwhile, according to the independently set parameterized layout sequence of the preliminary drilling hole positions and the hole position repeated layout spacing condition, filter the preliminary drilling hole positions to obtain the filtered preliminary drilling hole positions;

[0122] Step S33-3, optimize the preliminary drilling hole positions:

[0123] For the preliminary drilling hole positions, adopt adaptive weight optimization by combining the optimization boundary distance threshold of the preliminary drilling hole positions and the terrain slope, screen out the optimized drilling hole positions that are closest to the preliminary drilling hole positions, have the most suitable terrain slope, and are convenient for drilling layout, and identify the elevation values of the optimized drilling hole positions;

[0124] The adaptive weight optimization is:

[0125] F(x,y) = w1·D(x,y) + w2·S(x,y)

[0126]

[0127] Where:

[0128] F(x, y) is the objective function, D(x, y) is the normalized distance score between the candidate point and the initial hole position, S(x, y) is the normalized slope score of the candidate point, w1 and w2 are weight coefficients, which are set independently or set as dynamic weight coefficients, satisfying w1 + w2 = 1;

[0129] (x0, y0) is the coordinate of the initial hole position, and d is the threshold of the optimized boundary distance of the initial hole position;

[0130] (x, y) is the coordinate of the candidate point; Slope(x, y) is the terrain slope value of the candidate point, Slope min is the minimum terrain slope value within the range of the boundary distance threshold, Slope max is the maximum terrain slope value within the range of the boundary distance threshold.

[0131] Step S33-4, parameterize the drilling depth value:

[0132] Independently set the drilling depth value according to the drilling depth distance conditions, where the drilling depth distance conditions include the underlay distance value below the top boundary of the relatively impermeable layer and the underlay distance value below the groundwater level. The distance condition options can be screened and selected, and the drilling depth distance value is independently set for each option; if there are no factors of drilling depth distance conditions, independently set the drilling depth value;

[0133] Step S33-5, parameterize the drilling attribute information:

[0134] After parameterizing the drilling hole position, the elevation value of the drilling hole position and the drilling depth value, parameterize the drilling number, drilling elevation, and drilling depth attribute information; among them, the drilling number is determined according to the parameterized drilling layout order and the sequence of the intersection points of the hole positions in turn. The sequence is defaulted to the sequence of the terrain polyline nodes. For the bottom of the reservoir, the river bend and the interfluvial void area, it is defaulted to north first, south second, west first, east second, left first, and right second. For the front dam section, it is defaulted to the dam axis first and then the upstream dam front area.

[0135] Step S4, train to identify the reservoir type classification and construct the reservoir type data source:

[0136] Based on the three-dimensional image, orthophoto image, oblique image, and satellite image data, mark various reservoir types of the pumped-storage power station and train the reservoir type data source;

[0137] Step S5, carry out parametric exploration layout of reservoir leakage based on the reservoir type:

[0138] Based on the reservoir image data, the reservoir type is intelligently identified using the described reservoir type data source, and the reservoir leakage investigation boreholes are parametrically arranged in the following way: using the dam axis line of the reservoir, the boundary line of the reservoir perimeter, the normal storage level elevation value, and the topographic line extracted from the reservoir image data, in combination with the key parts of the reservoir leakage investigation in step S2. When parametrically arranging the reservoir leakage investigation boreholes, the elevation and depth attributes of the boreholes are parametrically displayed, and the borehole arrangement can be adjusted independently.

[0139] The following are examples:

[0140] Example 1:

[0141] S1. The pumped-storage power station reservoir types are carefully classified, that is, according to the location of the reservoir and the topographic and geomorphic characteristics of the reservoir area, the planned reservoir types are classified in detail.

[0142] In the present invention, the reservoir types are divided into five categories, including valley-type reservoirs, gully-type reservoirs, shore-type reservoirs, gully-source-type reservoirs, and mountaintop platform-type reservoirs.

[0143] A valley-type reservoir refers to a reservoir formed by building a dam across the main river in the project area.

[0144] A gully-type reservoir refers to a reservoir formed by building a dam in a tributary gully of a river, and the reservoir backwater line does not reach the gully source.

[0145] A shore-type reservoir refers to a reservoir formed by a combination of excavation and filling and enclosing with a dam on a river bank terrace or concave plot, without encroaching on the riverbed.

[0146] A gully-source-type reservoir refers to a reservoir formed by building a dam within the range of the gully source.

[0147] A mountaintop platform-type reservoir refers to a reservoir formed by a combination of excavation and filling and enclosing with a dam on a relatively flat mountaintop platform in the project area.

[0148] S2. The key points of reservoir leakage investigation for each reservoir type classification are sorted out and set, that is, according to the topographic and geomorphic type characteristics, engineering geology and hydrogeological conditions of the pumped-storage power station reservoir area, the key points of reservoir leakage investigation for each reservoir type are set for subsequent parametric intelligent arrangement of the investigation plan.

[0149] The key points of reservoir leakage investigation for each reservoir type set in step S2 include: for valley-type reservoirs, the key points are the front dam reservoir section, river bends, and interfluve areas; for shore-type and mountaintop platform-type reservoirs, the key points are the reservoir bottom and the reservoir perimeter; for gully-type and gully-source-type reservoirs, the key points are the front dam reservoir section, thin watersheds, low passes, etc.

[0150] S3. Parametric reservoir seepage investigation plan, that is, according to the topographic characteristics of reservoirs of different types and the key points of seepage investigation, parameterize the drilling workload plan for reservoir seepage investigation. The topographic characteristics in step S3 refer to the topographic slope, aspect, and topographic curvature within the area of the distance condition around the reservoir type. Among them, the distance condition around the reservoir type means taking the dam axis as the fixed boundary, combining the selected reservoir perimeter boundary or the normal storage level elevation line by the user, and independently setting the outward expansion distance, with a default value of 500m.

[0151] Step S3 includes the following steps:

[0152] S31. Parametric investigation area boundary, that is, according to the dam axis of the reservoir type, the reservoir perimeter boundary, and the normal storage level elevation value, generate the investigation area boundary line based on the topographic characteristics of the reservoir type and the distance condition around the reservoir type.

[0153] S32. Parametric key points for seepage problem investigation, that is, according to the key points of reservoir seepage investigation for different reservoir types in step S2, combined with the topographic characteristic factors of the investigation area, set the drilling areas for the key points of parametric layout. The steps are as follows:

[0154] S32-1. Generate the key areas for seepage investigation according to the reservoir type and the key points of seepage investigation in step S2, and set the following parametric key points for investigation.

[0155] (1) The drilling area in the front of the dam section of the reservoir is the river valley area combined with the reservoir perimeter line or the normal storage level elevation line and the independently set distance condition upstream of the reservoir dam axis. Among them, the distance condition defaults to the range of 100 - 500m, and the D8 algorithm is used to identify the upstream of the reservoir;

[0156] (2) The drilling area in the river bend is set in the area with a gentle terrain where the terrain slope of the river bend is less than 25° according to the topographic curvature of the river valley, and the radius of curvature is generally less than 5 times the width of the river channel at the normal storage level;

[0157] (3) The drilling area in the interfluve is set in the area with a gentle terrain where the terrain slope between two adjacent river channels is less than 25° according to the topographic data between the adjacent river channels.

[0158] (4) The drilling area at the bottom of the reservoir is the internal range area within the investigation area of the reservoir dam axis and the reservoir perimeter line; the drilling area around the reservoir is the reservoir perimeter line area.

[0159] (5) For the thin watershed, according to the topographic elevation value, topographic slope, and curvature of the investigation area in step S31, identify the ridge width, ridge peak, and pass value at the normal storage level elevation or the reservoir perimeter elevation. The ridge width around the reservoir lower than the independently set width value is the thin watershed, and the investigation area is the ridge peak and pass area, where the independently set width value defaults to 200m.

[0160] (6) The low pass drilling area is the area with gentle terrain on the ridge line of the non-thin watershed on the reservoir perimeter within the self-set distance condition in the survey area of step S31. The distance condition is defaulted to the area within 500m from the center point of the pass in the direction normal to the reservoir perimeter.

[0161] S32-2. Generate the boundary line of the survey area generated according to S32-1 with a spatial geometric analysis algorithm; and the user can independently delimit the boundary of the drilling survey area for subsequent borehole layout.

[0162] S33. Parametrically arrange the drilling, that is, according to the terrain area and the overlapping area boundary of the key drilling area boundary set in the above steps S31 and S32, the overlapping part and the boundary of the drilling layout area, set the parametric borehole spacing, borehole position, and borehole depth conditions, and intelligently set the drilling workload. It includes the following steps:

[0163] S33-1. Parametric borehole spacing, which is the parameter for setting the self-set spacing condition within the above drilling area. The default setting in the pre-feasibility study stage is 100m, and the default setting in the feasibility study stage is 50m. The user can also set the value independently.

[0164] S33-2. Parametric initial borehole position, which is to optimize and select the borehole orifice position by combining the drilling area boundary and borehole spacing factors. According to the above S32 key survey area and reservoir type factors, set the following parametric methods:

[0165] (1) In the front dam reservoir section and the reservoir bottom part including the dam axis, that is, the area including the dam axis factor, take the center point of the dam axis as the starting point, generate the layout line perpendicular to the dam axis through vector operation, and then gradually parametrically generate the layout line parallel to the dam axis and the perpendicular layout line according to the borehole spacing condition and the area boundary condition to generate grid intersection points, and each point is the preliminary borehole position; and the preliminary borehole positions are also generated on the dam axis section according to the borehole spacing condition and the area boundary condition.

[0166] (2) In the river bend and interfluve area, and the reservoir bottom part without the dam axis, that is, the area without the dam axis factor, take the center point of the fixed-point area by spatial geometric algorithm, and generate the positive grid line and intersection points according to the borehole spacing condition and the area boundary condition, and each point is the preliminary borehole position.

[0167] (3) In the reservoir perimeter part, that is, the multi-segment line of the reservoir perimeter of the shore type and mountaintop platform type reservoirs, take the starting point of the multi-segment line as the drilling starting point, and generate each preliminary borehole position according to the borehole spacing condition.

[0168] (4) In parts such as thin watersheds and low passes, as mentioned above, the borehole positions are basically arranged at the peaks and ridge passes.

[0169] Meanwhile, the user sets the parametric layout sequence of the preliminary drilling hole positions and the spacing condition for repeated layout of the hole positions according to their own will. The default value is 50m. That is, when arranging the drilling hole positions in subsequent sequences, if there is a drilling hole arranged in the previous sequence within 50m, the hole position at this intersection point will be automatically cancelled.

[0170] S33-3. Optimize the drilling hole positions. That is, based on the preliminary drilling hole positions determined in S33-2, adopt an adaptive weight optimization by combining the optimization boundary distance threshold from the initial hole positions and the terrain slope to screen out the optimized hole positions that are closest to the initial drilling hole positions, most adaptable to the terrain slope, and convenient for drilling layout, and identify the elevation values of the hole positions.

[0171] Among them: The adaptive weight optimization algorithm is as follows:

[0172] F(x,y) = w1·D(x,y) + w2·S(x,y)

[0173]

[0174] In the formula, F(x,y) is the objective function, D(x,y) is the normalized distance score between the candidate point and the initial hole position, S(x,y) is the normalized slope score of the candidate point, (x0, y0) is the coordinate of the initial hole position, d is the optimization boundary distance threshold of the initial hole position, generally set to 5 - 10m. (x, y) is the coordinate of the candidate point.

[0175] In the formula, w1 and w2 are weight coefficients, which can be set by the user independently or set as dynamic weight coefficients, satisfying w1 + w2 = 1.

[0176] In the formula, Slope(x,y) is the terrain slope value of the candidate point, Slope min is the minimum terrain slope value within the boundary distance threshold, Slope max is the maximum terrain slope value within the boundary distance threshold.

[0177] Furthermore, since there are many factors of steep and unstable terrain within the selected range of pumped - storage power stations, the whale WOA algorithm and the particle swarm PSO algorithm can be used to further improve the selection efficiency.

[0178] S33-4. Parametric drilling depth, which is the drilling depth value set according to the depth distance conditions set independently. The depth distance conditions include the underlying distance value below the top boundary of the relatively impermeable layer and the underlying distance value below the groundwater level. The distance condition options can be screened and selected, and the distance values can be set independently for each option. Among them, the default value of the relative impermeable layer distance condition is 15m, and the default value of the groundwater level distance condition is 20m. If there are no depth distance condition factors, the depth value can be set independently, with a default value of 50m.

[0179] S33-5. Parametric drilling attribute information, that is, after parameterizing the drilling hole positions, elevation values, and depth values according to user parameters, parameterize the drilling number, elevation, and depth attribute information. Among them, the drilling number is determined in sequence according to the above parametric drilling layout order and the sequence of the intersection points of the hole positions. The sequence is defaulted to the node sequence of the topographic polyline. For the bottom of the reservoir, river bends, and interfluvial voids, it is defaulted to north first, south second, west first, east second, left first, and right second. For the front dam section, it is defaulted to the dam axis first and then the upstream front dam area.

[0180] S4. Train and identify the reservoir type classification and construct the reservoir type data source. That is, according to data such as 3D images, orthophotos, oblique images, and satellite images, users mark various reservoir type classifications for pumped-storage, and adopt the YOLOv5 or YOLOv8 algorithm to train the reservoir type data source.

[0181] As is well known, the reservoir image data is limited. In this training and identification, geometric transformation methods such as rotation, flipping, scaling, and cropping are used for the image data, color transformation methods such as adjusting brightness, contrast, and saturation are used, and image data enhancement such as Gaussian blur processing is carried out; as well as the pre-training transfer learning method to effectively achieve high-precision classification of reservoir types.

[0182] In this embodiment, according to data such as satellite images, orthophotos, oblique images, and photo images, geometric transformation methods such as rotation, flipping, scaling, and cropping of the image data are adopted at the same time, and color transformation methods such as adjusting brightness, contrast, and saturation are used to train the reservoir type database based on the YOLOv5 algorithm for the image data.

[0183] S5. Carry out parametric investigation and layout of reservoir seepage based on reservoir type classification, that is, according to the user's reservoir image data, obtain the reservoir classification from the intelligent recognition data source, and use the dam axis of the reservoir, the boundary line of the reservoir perimeter, the normal storage level elevation value, and the topographic line extracted from the user's image to parametrically layout the seepage investigation drilling holes according to the above steps.

[0184] Through the following embodiments, combined with the attached Figures 2 to 4 The technical solution of the present invention will be described in detail.

[0185] During the feasibility study stage of a pumped-storage power station, according to the perimeter line of the upper reservoir and the 3D satellite image of the Ovi software as Figure 2 shown, the upper reservoir is intelligently identified as a gully-source type reservoir, and then the key points of the upper reservoir seepage investigation are sorted out as the front dam section, the thin watershed, the low pass, etc.

[0186] Based on satellite images, terrain elevation lines, reservoir perimeters, dam axes, and normal water level lines are proposed to carry out parameterized reservoir leakage investigation plans. The setting sequence is as follows: the first step is the reservoir section in front of the dam, and the distance condition is set to 100m from the dam axis; the second step is the low pass around the reservoir, and the distance condition is set to be within 500m from the center of the pass to the normal line of the reservoir; the third step is the thin watershed, and the area with a ridge width of 100m at the normal water level elevation is screened out. Figure 3 The two thin watershed areas of the Kuzhou Ridge A and B are shown.

[0187] Then, because leakage investigation had been carried out in the pre-feasibility study stage of the project, the leakage problem was identified in the feasibility study stage, and the borehole spacing was set to 100m, the borehole depth was 100m, and the borehole numbering sequence referred to the above survey plan sequence. After parameterizing the preliminary borehole positions, the adaptive weight optimization of the Whale WOA algorithm was performed with the optimized boundary distance threshold of 50m. 21 leakage investigation boreholes were arranged on the terrain line, with a total drilling workload of 2100m, and the borehole elevation and depth attributes were parameterized, see the attached Figures 3 to 4 , so that subsequent users can adjust the drilling layout independently.

[0188] Through the following examples, combined with the attached Figures 5 to 7 The technical solution of the present invention is described in detail.

[0189] In this embodiment, during the pre-feasibility study phase of a pumped storage power station, Figure 5 The lower reservoir perimeter and Aowei software 3D satellite image shown are identified as a shore-type reservoir or a mountaintop terrace-type reservoir because the image is flat. The user independently selects it as a shore-type reservoir, and then sorts out that the focus of the lower reservoir leakage investigation is the reservoir bottom and the surrounding area.

[0190] According to the satellite images, the terrain elevation line and reservoir perimeter are proposed, and a parameterized reservoir leakage survey plan is carried out. The setting sequence is the first step of the reservoir perimeter, and the second step is the reservoir bottom area, such as Figure 6 .

[0191] Then, because the project is in the pre-feasibility study stage and needs to preliminarily identify the leakage problem, the drilling spacing around the reservoir is set to 150m, the drilling spacing inside the reservoir is set to 200m, and the drilling depth is set to 50m. The drilling numbering sequence refers to the above survey plan sequence and the sequence of polyline nodes around the reservoir. After parameterizing the initial hole positions, the Whale WOA algorithm is used to optimize the hole positions with an optimized boundary distance threshold of 50m. Multiple hole positions are optimized, and 31 leakage exploration holes are arranged on the terrain line, with a total drilling workload of 1550m. The drilling elevation and depth attributes are displayed parameterized, such as Figures 6 - 7 , so that subsequent users can adjust the drilling layout independently.

[0192] The beneficial effects of the present invention are:

[0193] (1) The reservoir type classification extracted by the present invention is based on conditions such as terrain, rivers, and valleys for classifying reservoirs. By adopting this classification method, the accuracy of subsequent parametric exploration layout for reservoir seepage can be improved.

[0194] (2) For the seepage exploration elements emphasized during the exploration and design of pumped-storage power stations, the present invention combines the pumped-storage reservoir type classification with parametric seepage exploration through intelligent recognition and parametric layout rules, improving the efficiency of exploration layout, significantly reducing manual field climbing for point selection, and enhancing the efficiency and accuracy of reservoir seepage exploration layout work.

[0195] (3) Regarding the borehole layout, based on multiple factors such as terrain slope, curvature, and borehole spacing, the present invention introduces a multi-objective optimization algorithm to ensure the rationality and accuracy of borehole layout, facilitating the solution of the difficult problem of borehole positioning in steep terrains. The present invention not only improves the efficiency of borehole positioning but also optimizes the number and location of boreholes, effectively reducing the exploration cost and improving the exploration quality.

[0196] The above are only the preferred embodiments of the present invention. It should be noted that for those of ordinary skill in the art, without departing from the principle of the present invention, several improvements and refinements can be made, and these improvements and refinements should also be regarded as the protection scope of the present invention.

Claims

1. A parametric exploration layout method for leakage of pumped storage reservoirs based on library type classification, characterized in that, It includes the following steps: Step S1: Divide the pumped-storage power station reservoir into multiple reservoir types according to the location of the reservoir and the topographic and geomorphic characteristics of the reservoir area; Step S2: Set the key parts for reservoir seepage investigation for each reservoir type; Step S3: According to the reservoir type, topographic characteristics and key parts for reservoir seepage investigation of each reservoir type, adopt the parametric algorithm for reservoir seepage investigation drilling plan to output the parametric reservoir seepage investigation drilling workload plan; Step S4: Train the recognition of reservoir type classification and construct the reservoir type data source: Mark various reservoir types of the pumped-storage power station according to the three-dimensional image, orthophoto image, oblique image and satellite image data, and train the reservoir type data source; Step S5: Carry out parametric investigation layout of reservoir seepage based on the reservoir type: According to the reservoir image data, intelligently identify the reservoir type by using the reservoir type data source, and combine the dam axis of the reservoir, the boundary line of the reservoir perimeter, the normal storage level elevation value, and the topographic line extracted from the reservoir image data, and combine with the key parts for reservoir seepage investigation in Step S2, and parametrically arrange the reservoir seepage investigation boreholes in the way of Step S3.

2. The parametric exploration layout method for leakage of pumped - storage reservoirs based on library - type classification according to claim 1, characterized in that, The reservoir types include valley-type reservoirs, gully-type reservoirs, shore-type reservoirs, gully-source type reservoirs and mountaintop terrace-type reservoirs; The valley-type reservoir refers to the reservoir formed by building a dam across the main river in the project area; The gully-type reservoir refers to the reservoir formed by building a dam in the tributary gully of the river and the reservoir backwater line does not reach the gully source; The shore-type reservoir refers to the reservoir formed by combining excavation and filling and building a dam on the river bank terrace or concave land plot without occupying the river bed; The gully-source type reservoir refers to the reservoir formed by building a dam within the range of the gully source; The mountaintop terrace-type reservoir refers to the reservoir formed by combining excavation and filling and building a dam on the relatively flat terrace on the mountaintop in the project area.

3. A parametric exploration layout method for leakage of pumped storage reservoirs based on library type classification according to claim 2, characterized in that, The key parts for reservoir seepage investigation of each reservoir type refer to setting the key parts for reservoir seepage investigation of each reservoir type according to the topographic and geomorphic type characteristics, engineering geology and hydrogeological conditions of the pumped-storage power station reservoir area, and the characteristics of the reservoir type.

4. A parametric exploration layout method for leakage of pumped - storage reservoirs based on library - type classification according to claim 2, characterized in that The key parts for reservoir seepage investigation of each reservoir type specifically include: For the valley-type reservoir, the key parts for reservoir seepage investigation are the reservoir section in front of the dam, river bends and interfluve areas; For the shore-type reservoir and the mountaintop terrace-type reservoir, the key parts for reservoir seepage investigation are the reservoir bottom and the reservoir perimeter area; For the gully-type reservoir and the gully-source type reservoir, the key parts for reservoir seepage investigation are the reservoir section in front of the dam, the thin watershed and the low saddle parts.

5. A parametric exploration layout method for leakage of pumped storage reservoirs based on library type classification according to claim 4, characterized in that Step S3 is specifically: Step S31: Parametrize the boundary of the investigation area: According to the dam axis of the reservoir type, the boundary line of the reservoir perimeter and the normal storage level elevation value, generate the boundary of the investigation area based on the topographic characteristics of the reservoir type and the distance conditions around the reservoir type; among them, the distance conditions around the reservoir type refer to taking the dam axis as the fixed boundary, combining the selected boundary line of the reservoir perimeter or the normal storage level elevation line by the user, and independently setting the outward expansion distance; the topographic characteristics of the reservoir type refer to the terrain slope, aspect and terrain curvature within the area of the distance conditions around the reservoir type that meet the distance conditions around the reservoir type; Step S32: Parametrize the drilling area of the key parts for reservoir seepage investigation: According to the key parts of reservoir seepage investigation for each reservoir type in step S2, combined with the topographic characteristics of the investigation area, a parameterized drilling area for the key parts of reservoir seepage investigation is set up; Step S33, parameterized drilling layout: The overlapping area of the investigation area set in step S31 and the drilling area for the key parts of reservoir seepage investigation set in step S32 is used as the boundary of the drilling layout area; within the boundary of the drilling layout area, the drilling workload is intelligently set by parameterizing the drilling hole spacing, hole position, and drilling depth.

6. A parametric exploration layout method for leakage of pumped storage reservoirs based on library type classification according to claim 5, characterized in that Step S32 specifically includes: Step S32-1, generating a drilling area for the key parts of reservoir seepage investigation according to the reservoir type and the key parts of reservoir seepage investigation for each reservoir type in step S2, including: For the pre-dam reservoir section, the valley area with a self-set distance condition upstream of the reservoir dam axis, combined with the reservoir perimeter line or the normal storage level elevation line, is used as the drilling area for the pre-dam reservoir section; For the river bend plot, according to the curvature of the valley terrain, where the radius of curvature is less than 5 times the width of the river channel at the normal storage level, it is set in the area with a gentle terrain slope less than 25° of the river bend terrain as the drilling area for the river bend plot; For the interfluve plot, according to the topographic data between adjacent river channels, it is set in the area with a gentle terrain slope less than 25° between the two river channels as the drilling area for the interfluve plot; For the reservoir bottom, the internal range area of the investigation area between the reservoir dam axis and the reservoir perimeter line is used as the drilling area for the reservoir bottom; For the reservoir perimeter, the reservoir perimeter line area is used as the drilling area for the reservoir perimeter; For the thin watershed, according to the topographic characteristics of the reservoir type in the investigation area of step S31, the ridge width, ridge peak, and pass value at the normal storage level elevation or the reservoir perimeter elevation are identified. The reservoir perimeter ridge width lower than the self-set width value is the thin watershed, and the ridge peak and pass parts are used as the drilling area for the thin watershed; For the low pass part, the gentle terrain area of the pass on the ridge line of the non-thin watershed of the reservoir perimeter within the area that meets the reservoir type perimeter distance condition in the investigation area of step S31 is used as the drilling area for the low pass part; Step S32-2, generating the boundary of the drilling area for the key parts of reservoir seepage investigation by using a spatial geometric analysis algorithm according to the drilling area for the key parts of reservoir seepage investigation generated in step S32-1.

7. A parametric exploration layout method for leakage of pumped storage reservoirs based on library classification according to claim 5, characterized in that Step S33 specifically includes: Step S33-1, parameterized drilling hole spacing: Within the boundary of the drilling layout area, the drilling hole spacing is self-set and parameterized; Step S33-2, parameterized preliminary drilling hole positions: Combined with the boundary of the drilling layout area and the drilling hole spacing, the preliminary drilling hole positions are optimized and selected, including: (1) For the reservoir section in front of the dam and the reservoir bottom part including the dam axis, that is, the area containing the dam axis factor, the center point of the dam axis is taken as the starting point, and the layout line perpendicular to the dam axis is generated through vector operations. Subsequently, according to the borehole spacing conditions and the boundary conditions of the drilling layout area, the layout lines parallel to the dam axis and the layout lines perpendicular to the dam axis are gradually parameterized; grid intersection points are generated for each horizontal and vertical layout line, and each of the grid intersection points is the preliminary borehole position; moreover, preliminary borehole positions are also generated on the dam axis section according to the borehole spacing conditions and the boundary conditions of the drilling layout area. (2) For river bends, interfluve plots, and the reservoir bottom part without the dam axis, that is, the area without the dam axis factor, the center point of the fixed-point area is calculated by spatial geometry, and positive grid lines are generated according to the borehole spacing and the boundary conditions of the drilling layout area. The intersection points of each grid line are the preliminary borehole positions. (3) For the reservoir perimeter part, that is, the multi-segment lines around the shore-type and mountaintop platform-type reservoirs, the starting point of the multi-segment line is taken as the drilling starting point, and each preliminary borehole position is generated according to the borehole spacing. (4) For thin watersheds and low passes, the preliminary borehole positions are arranged at the peak and ridge passes. Meanwhile, according to the independently set parameterized layout sequence of the preliminary borehole positions and the spacing conditions for repeated layout of the hole positions, the preliminary borehole positions are filtered to obtain the filtered preliminary borehole positions. Step S33-3, optimize the preliminary borehole positions: For the preliminary borehole positions, adaptive weight optimization is carried out by combining the optimization boundary distance threshold of the preliminary borehole positions and the terrain slope, screening out the optimized borehole positions that are closest to the preliminary borehole positions, have the most suitable terrain slope, and are convenient for drilling layout, and identifying the elevation values of the optimized borehole positions. Step S33-4, parameterize the borehole depth value: The borehole depth value of the borehole depth distance condition is independently set, where the borehole depth distance condition includes the underlying distance value below the top boundary of the relatively impermeable layer and the underlying distance value below the groundwater level. The distance condition options can be screened and selected, and the borehole depth distance value is independently set for each option; if the borehole depth distance condition factor is not available, the borehole depth value is independently set. Step S33-5, parameterize the borehole attribute information: After parameterizing the borehole positions, the elevation values of the borehole positions, and the borehole depth values, the borehole number, borehole elevation, and borehole depth attribute information are parameterized; among them, the borehole number is determined in sequence according to the parameterized borehole layout order and the sequence of the hole position intersection points. The sequence is defaulted to the node sequence of the terrain multi-segment line. For the reservoir bottom, river bends, and interfluve empty areas, it is defaulted to north first, south second, west first, east second, left first, right second, and for the front dam section, it is defaulted to the dam axis first and then the upstream dam front area.

8. A parametric exploration layout method for leakage of pumped storage reservoirs based on library classification according to claim 7, characterized in that, In step S33-3, the adaptive weight optimization is as follows: F(x,y) = w1·D(x,y) + w2·S(x,y) Where: F(x,y) is the objective function, D(x,y) is the normalized distance score between the candidate point and the initial hole position, S(x,y) is the normalized slope score of the candidate point, w1 and w2 are weight coefficients, which are independently set or set as dynamic weight coefficients, and satisfy w1 + w2 = 1. (x0, y0) is the initial hole position coordinate, and d is the threshold of the optimization boundary distance of the initial hole position; (x, y) is the coordinate of the candidate point; Slope(x, y) is the terrain slope value of the candidate point, Slope min is the minimum terrain slope value within the boundary distance threshold range, Slope max is the maximum terrain slope value within the boundary distance threshold range.

9. A parametric exploration layout method for leakage of pumped-storage reservoirs based on library classification according to claim 1, characterized in that, When parametrically arranging the reservoir leakage investigation boreholes, the elevation and depth attributes of the boreholes are parametrically displayed, and the borehole arrangement can be adjusted independently.

Images