Contour line method storage capacity determination method and system considering river center island

Abstract

The invention discloses a contour method reservoir capacity determination method and system considering a river center island. The method comprises the following steps: acquiring elevation information in a reservoir area and a closed contour area corresponding to the elevation information; performing self-intersection processing on the closed contour line area under each piece of elevation information to obtain a completely overlapped area in each closed contour line area; determining a river center island contour line region corresponding to each piece of elevation information according to each completely overlapped region; and determining the reservoir capacity according to each elevation information, the corresponding closed contour line area and the corresponding river center island contour line area. According to the method, under the condition of the same elevation, the contour line area of the river center island is quickly identified and positioned by utilizing the spatial relationship between the contour line area of the river center island and the closed contour line area, and accurate subtraction is carried out on the contour line area of the river center island through a spatial analysis technology, so that the area of the effective contour line area is obtained; and the accuracy of reservoir capacity calculation is ensured.

Description

Technical Field

[0001] The present invention relates to the fields of reservoir management and water resources assessment, and in particular to a method and system for determining reservoir capacity using a contour method taking river islands into consideration. Background Art

[0002] In the field of reservoir management and water resources assessment, reservoir capacity is a core indicator for measuring the water storage capacity of a reservoir. Specifically, it refers to the total water storage capacity or storage volume of the reservoir below a certain water level or between two water levels.

[0003] The contour method is widely used in the calculation or verification of reservoir capacity. The basic principle of this method is to calculate the reservoir's storage capacity at different water levels by measuring the horizontal area enclosed by adjacent contour lines and combining the elevation difference between them (i.e., the contour interval).

[0004] However, in practical applications, the contour method often treats the entire reservoir area as a continuous and uniform water storage space. This assumption is generally feasible when the reservoir terrain is relatively simple and there are no significant non-water storage areas. However, when river islands exist within the reservoir area, as land within the reservoir area, the river islands do not have water storage capacity, resulting in biased calculation results and failing to accurately reflect the actual water storage capacity of the reservoir. Summary of the Invention

[0005] In order to solve the above problems, the present invention discloses a method and system for determining reservoir capacity using the contour method taking river islands into consideration.

[0006] The present invention discloses a method for determining reservoir capacity by using a contour method taking river islands into consideration, comprising the following steps:

[0007] Obtain each elevation information within the reservoir area and its corresponding closed contour area;

[0008] Perform self-intersection processing on the closed contour area under each elevation information to obtain the completely overlapping area within each closed contour area;

[0009] Determine the river island contour area corresponding to each elevation information based on each completely overlapping area;

[0010] The reservoir capacity is determined based on each elevation information and the corresponding closed contour area and river island contour area.

[0011] Preferably, the reservoir capacity is determined based on each elevation information and the corresponding closed contour area and river island contour area, specifically:

[0012] Determine the effective contour size information corresponding to each elevation information according to the closed contour area and river island contour area corresponding to each elevation information;

[0013] The reservoir capacity is determined based on each elevation information and its corresponding effective contour size information.

[0014] Preferably, the effective contour size information corresponding to each elevation information is determined based on the closed contour area and the river island contour area corresponding to each elevation information, specifically:

[0015] Determine the valid contour area corresponding to each elevation information according to the closed contour area and river island contour area corresponding to each elevation information;

[0016] The effective contour line size information corresponding to each elevation information is determined according to each effective contour line area.

[0017] Preferably, the effective contour size information corresponding to each elevation information is determined based on the closed contour area and the river island contour area corresponding to each elevation information, specifically:

[0018] Determine the closed contour size information corresponding to each elevation information according to the closed contour area corresponding to each elevation information;

[0019] Determine the river island contour size information corresponding to each elevation information according to the river island contour area corresponding to each elevation information;

[0020] The effective contour line size information corresponding to each elevation information is determined according to the closed contour line size information and the river island contour line size information corresponding to each elevation information.

[0021] Preferably, each elevation information in the reservoir area and its corresponding closed contour area are obtained, specifically:

[0022] Obtain the overall closed contour area within the reservoir area;

[0023] The overall closed contour area is segmented according to the elevation field to obtain each elevation information in the reservoir area and its corresponding closed contour area.

[0024] Preferably, the reservoir capacity is determined based on each elevation information and its corresponding effective contour size information, including:

[0025] Integrate each elevation information and its corresponding effective contour size information to obtain the actual overall closed contour information;

[0026] Determine the volume between each group of adjacent contour lines based on the actual overall closed contour line information;

[0027] The storage capacity is determined based on the volume between multiple adjacent contour lines.

[0028] Preferably, the overall closed contour area within the reservoir area is obtained, specifically:

[0029] Acquire the reservoir area contour area, wherein the reservoir area contour area includes closed contour area and unclosed contour area under all elevation information;

[0030] Performing closing processing on the unclosed contour area;

[0031] The closed contour area and the unclosed contour area after the closure process are regarded as the overall closed contour area in the reservoir area.

[0032] The present invention also discloses a system for determining reservoir capacity using a contour method taking river islands into consideration, comprising:

[0033] The acquisition module is used to obtain each elevation information in the reservoir area and its corresponding closed contour area;

[0034] A processing module is used to perform self-intersection processing on the closed contour area under each elevation information to obtain a completely overlapping area within each closed contour area;

[0035] A preset area determination module is used to determine the river island contour area corresponding to each elevation information based on each completely overlapping area;

[0036] The reservoir capacity determination module is used to determine the reservoir capacity based on each elevation information and the corresponding closed contour area and river island contour area.

[0037] Compared with the prior art, the present invention has the following beneficial effects:

[0038] (1) Under the condition of the same elevation, the present invention uses the spatial relationship between the river island contour area and the closed contour area to quickly identify and locate the river island contour area, and uses spatial analysis technology to accurately reduce the river island contour area, thereby obtaining the area of ​​the effective contour area and ensuring the accuracy of the reservoir capacity calculation;

[0039] (2) The method of the present invention does not require manual measurement. Through automated data processing and spatial analysis technology, it realizes the rapid processing of reservoir contour data and the accurate calculation of reservoir capacity, improves the efficiency of calculating area and reservoir capacity, and reduces labor and time costs. BRIEF DESCRIPTION OF THE DRAWINGS

[0040] Figure 1 is a flow chart of the method of the present invention;

[0041] Figure 2 It is a schematic diagram of identifying the river island contour area from the closed contour area under a single elevation of the present invention;

[0042] Figure 3 It is a schematic diagram of subtracting the river island contour area from the closed contour area according to the present invention;

[0043] Figure 4 (a) is the display effect of the unclosed contour area and the dam axis after physical edge connection in Example 1 in the CAD environment;

[0044] Figure 4 The surface filling area in (b) is the unclosed contour area after the closing process in Example 1;

[0045] Figure 5 (a) in the figure is the unclosed contour area after logical edge joining in Example 2;

[0046] Figure 5 (b) in the figure is the unclosed contour area after physical edge joining in Example 2;

[0047] Figure 6 (a) is the GIS format graphic display of the unclosed contour lines and dam axis output in Example 4;

[0048] Figure 6 (b) is the result of surface filling display of the unclosed contour lines of the inner circle and the dam axis surface in Example 4;

[0049] Figure 6 (c) in the figure is the surface filling display of the unclosed contour lines of the outer circle and the dam axis surface in Example 4.

[0050] In the figure, 1 is the closed contour area under each elevation information; 2 is the river island contour area under each elevation information; 3 is the effective contour area. DETAILED DESCRIPTION

[0051] In the following description, specific details such as particular system structures and techniques are provided for purposes of illustration, not limitation, to facilitate a thorough understanding of the embodiments of the present invention. However, it will be apparent to those skilled in the art that the present invention may be practiced in other embodiments without these specific details. In other cases, detailed descriptions of well-known systems, devices, circuits, and methods are omitted so as not to obscure the description of the present invention with unnecessary detail.

[0052] like Figure 1 As shown, the present invention discloses a method for determining reservoir capacity by contour line method taking river island into consideration, comprising the following steps:

[0053] S1. Obtain each elevation information in the reservoir area and the closed contour area 1 under each elevation information; in the present invention, the elevation information includes the elevation value.

[0054] Preferably, S1 is specifically:

[0055] S11. Obtain the overall closed contour area within the reservoir area; the overall closed contour area includes the range enclosed by closed contour lines at all elevations within the reservoir area.

[0056] In the present invention, the region specifically includes the line type and coordinates corresponding to the closed contour line. The concept of region appearing in the subsequent description is the same as that mentioned this time and will not be repeated.

[0057] The overall closed contour area is a large and comprehensive concept, which is usually represented on a map as a continuous area spanning multiple elevation values.

[0058] Preferably, S11 is specifically:

[0059] S111. Obtaining the reservoir area contour area, where the reservoir area contour area includes closed contour areas and unclosed contour areas under all elevation information;

[0060] S112, closing the unclosed contour area;

[0061] According to relevant data from the National Geographic System, when drawing contour lines within the reservoir area, the lines are usually broken at the dam site. This requires using the dam axis as an auxiliary line to close all unclosed contour lines and ensure the continuity of the contour lines.

[0062] In the GIS system, the corresponding dimension information can be directly generated based on the closed contour area 1 under each elevation information. In other words, the basis for obtaining the dimension information of each contour area is to ensure that the contour is closed, that is, the unclosed contour areas must be logically connected or closed. The specific closing process is: using the dam axis as an auxiliary line, physical and logical connection are performed. The closing process can be implemented in the CAD environment or other software environments. For example, if the data file is relatively large, it can be processed in the GIS system.

[0063] Specifically, the size information in the present invention includes the area or perimeter of the closed contour area 1 under each elevation information.

[0064] Several embodiments of the sealing process of the present invention are as follows:

[0065] In Example 1, physical edge connection is first performed to ensure that the unclosed contour lines at the non-dam axis are complete. The physical edge connection process can be directly realized by automatic connection in a CAD environment, or it can be processed in other software environments. For example, in the case of a large data file, the DISSOLVE tool can be used in a GIS system. The next step is to connect the contour lines at the dam axis [i.e. Figure 4 (a) and Figure 4 The short vertical straight line in (b) will close the unclosed contour line. Figure 4As shown in (b), you can use the boundary command in CAD to generate a closed area, and then reassign attributes to the closed area.

[0066] In Example 2, auxiliary lines can also be drawn manually. Figure 5 In (a)], first ensure the logical connection; Figure 5 As shown in (b), the edges are then physically connected to form a whole.

[0067] In Example 3, the dam axis can also be copied multiple times, interrupted at the intersections with the unclosed contour lines, and then matched with the original unclosed contour lines, and then physically connected into a whole.

[0068] like Figure 6 As shown in (a) in the figure, Example 4 selects all unclosed contour lines and the dam axis at the dam axis under GIS, constructs faces with the dam axis in turn through iterative feature selection, and finally integrates all closed faces and outputs them in CAD format to replace the original unclosed contour lines, as shown in Figure 4. Figure 6 (b) and Figure 6 As shown in (c) in the figure.

[0069] S113. The closed contour areas under all elevation information and the unclosed contour areas after closing are regarded as the overall closed contour area in the reservoir area.

[0070] After the closing process is completed, the relevant information of the overall closed contour area in the reservoir area is finally output as the original CAD format contour file. Specifically, the CAD format is or .

[0071] Load the original CAD format contour files into the GIS system. These files usually contain the overall closed contour area within the reservoir area (such as line type, coordinates, etc.) and the corresponding attribute information (elevation value).

[0072] S12. Segment the entire closed contour area according to the elevation field to obtain each elevation information in the reservoir area and the closed contour area 1 under each elevation information.

[0073] S12 specifically includes: in the GIS system, directly outputting the Polygon feature layer as a GIS format file, and storing the result in a database (GDB / MDB).

[0074] Output the polygonal area as a database-formatted area vector file A. This file is used to store the area and attribute information of the Polygon feature layer. Each polygon in area vector file A corresponds to a specific elevation surface and contains the geometric characteristics (area, perimeter) and corresponding attribute information (elevation value) of that elevation surface.

[0075] Use the elevation field to split the surface vector A in the database format. In the GIS system, use the GIS tool "Split ByAttributes" to split the surface vector A according to the elevation field to obtain the elevation feature layer Bi corresponding to each elevation value.

[0076] The present invention utilizes the conversion relationship between CAD and GIS data to realize the rapid reading and processing of reservoir area contour data; and by segmenting the overall closed contour area, the closed contour area 1 corresponding to each elevation information is accurately obtained, providing a basis for subsequent area calculation and reservoir capacity determination.

[0077] S2, performing self-intersection processing on the closed contour area 1 under each elevation information to obtain a completely overlapping area within the closed contour area 1 under each elevation information;

[0078] S3, determining the river island contour area 2 under each elevation information according to each completely overlapping area;

[0079] like Figure 2 As shown, at each individual elevation value, the river island is completely spatially contained within the closed contour area 1 corresponding to each elevation information. Using the Intersect tool to perform self-intersection processing on each elevation feature layer Bi, internal contour patches can be identified, resulting in the corresponding completely overlapping area, namely the river island contour area. In the GIS system, each river island contour area is displayed as the river island location layer Ci. Ci is batch output using the Intersect tool combined with the "Iterate Feature Class" iterator, and naming is achieved by adding the wildcard character %Name% to the specified string.

[0080] like Figure 3 As shown, the present invention uses the closed contour area 1 under each elevation value to accurately subtract the river island contour area 2 under each elevation information, avoiding the problem of identifying the closed contour areas under other elevation values ​​as river island contour areas, ensuring that the contour areas under each elevation value are correctly processed and calculated, thereby greatly improving the accuracy and reliability of reservoir capacity calculation.

[0081] Furthermore, the treatment method of the present invention is not affected by the complexity of the reservoir terrain and can be applied to various types of reservoirs, including reservoirs with complex terrain and multiple river islands.

[0082] S4. Determine the reservoir capacity according to each elevation information, the closed contour area 1 under each elevation information, and the river island contour area 2 under each elevation information.

[0083] Preferably, S4 is specifically:

[0084] S41, determining effective contour size information corresponding to each elevation information according to the closed contour area 1 under each elevation information and the river island contour area 2 under each elevation information;

[0085] In one embodiment, S41 specifically includes:

[0086] S411, determining a valid contour area 3 corresponding to each elevation information according to the closed contour area 1 under each elevation information and the river island contour area 2 under each elevation information;

[0087] S412. Determine the effective contour size information corresponding to each elevation information according to each effective contour area 3.

[0088] S41 is specifically as follows: in the GIS system, the erase tool is used, with each elevation feature layer Bi as the input feature and the corresponding river island location layer Ci as the erase feature, to obtain the effective contour area 3 after deducting the river island, named layer Di. At the same time, the GIS system will generate the size information of the effective contour area 3. In this embodiment, the size information of the effective contour area 3 refers to the area of ​​the effective contour area 3.

[0089] Use the merge tool in the GIS system to merge all Di layers into a new layer E.

[0090] Perform summary statistics on the attribute table of the new layer E, with the statistical field being area and the grouping field being elevation, and obtain the summed area record corresponding to each elevation value.

[0091] In another embodiment, S41 is specifically:

[0092] S411, determining the closed contour size information corresponding to each elevation information according to the closed contour area 1 under each elevation information;

[0093] S412, determining the river island contour size information corresponding to each elevation information according to the river island contour area 2 under each elevation information;

[0094] S413. Determine the effective contour line size information corresponding to each elevation information according to the closed contour line size information and the river island contour line size information corresponding to each elevation information.

[0095] In this embodiment, the size information refers to the area; combined with S12, it can be seen that in the GIS system, the corresponding closed contour line size information can be directly obtained based on the closed contour line area 1 under each elevation information. Similarly, the river island contour line is closed to obtain the river island contour line size information at each elevation value.

[0096] The difference between the closed contour area and the river island contour area at the same elevation is the effective contour area at that elevation. In the present invention, the number of river island contour areas 2 corresponding to each elevation is 0 to n, and correspondingly, the number of river island contour size information corresponding to each elevation is 0 to n. That is to say, in addition to river islands, the present invention is also applicable to terrain with local protrusions underwater but not exposed to the water surface. The self-intersection process can determine whether there are contour areas of raised parts of underwater terrain or river islands; if there are, they will be subtracted; if not, the subtraction operation can still be performed, but an empty layer is generated, and the original contour area remains unchanged after subtraction.

[0097] S42. Determine the reservoir capacity based on each elevation information and its corresponding effective contour dimension information.

[0098] Preferably, S42 includes:

[0099] S421. Integrate each elevation information and its corresponding valid contour size information to obtain the actual overall closed contour information; the actual overall closed contour information refers to all elevation information and its corresponding valid contour size information. This step can be regarded as a data sorting and merging process.

[0100] S422, determining the volume between each group of adjacent contour lines based on the actual overall closed contour line information;

[0101] S423. Determine the storage capacity based on the volume between multiple adjacent contour lines.

[0102] Calculate the segmented storage capacity (the storage capacity between two adjacent contour lines) according to the following contour line formula:

[0103]

[0104] in, and are the effective contour areas of two adjacent closed contour areas, is the contour interval, that is, the elevation difference between two adjacent contour lines.

[0105] Then sum up the segmented storage capacity to calculate the total storage capacity (the storage capacity below a certain elevation value).

[0106] The following is an example of using an editor in Excel to implement steps S422 and S423 of the present invention according to the formula in S423:

[0107]

[0108] In the table above, column C represents the storage capacity between adjacent volume contours, and the values ​​in column D represent the total storage capacity below each level elevation;

[0109] The calculation of segmented storage capacity is edited as follows:

[0110] Considering the lowest point elevation, the formula corresponding to cell C2 is: ; The formula corresponding to cell C3: ;

[0111] Accordingly, the storage capacity summary calculation is edited as follows:

[0112] The formula corresponding to cell D2 is: ; The formula corresponding to cell D3 is: .

[0113] The method of the present invention does not require manual measurement. Through automated data processing and spatial analysis technology, it realizes rapid processing of reservoir area contour data and accurate calculation of reservoir capacity; improves the efficiency of calculating area and reservoir capacity, and reduces labor and time costs.

[0114] The present invention also discloses a system for determining reservoir capacity using a contour method taking river islands into consideration, comprising:

[0115] An acquisition module is used to acquire each elevation information in the reservoir area and a closed contour area 1 under each elevation information;

[0116] A processing module is used to perform self-intersection processing on the closed contour area 1 under each elevation information to obtain a completely overlapping area within the closed contour area 1 under each elevation information;

[0117] A preset area determination module is used to determine the river island contour area 2 under each elevation information based on each completely overlapping area;

[0118] The reservoir capacity determination module is used to determine the reservoir capacity based on each elevation information and the closed contour area 1 under each elevation information and the river island contour area 2 under each elevation information.

[0119] Compared with the prior art, the present invention has the following beneficial effects:

[0120] (1) Under the condition of the same elevation, the present invention uses the spatial relationship between the river island contour area and the closed contour area to quickly identify and locate the river island contour area, and uses spatial analysis technology to accurately reduce the river island contour area, thereby obtaining the area of ​​the effective contour area and ensuring the accuracy of the reservoir capacity calculation;

[0121] (2) The method of the present invention does not require manual measurement. Through automated data processing and spatial analysis technology, it realizes the rapid processing of reservoir contour data and the accurate calculation of reservoir capacity, improves the efficiency of calculating area and reservoir capacity, and reduces labor and time costs.

[0122] The above descriptions are merely a few embodiments of the present application and do not constitute any form of limitation to the present application. Although the present application discloses the preferred embodiments as above, they are not intended to limit the present application. Any technical personnel familiar with the present profession, without departing from the scope of the technical solution of the present application, using the technical content disclosed above to make slight changes or modifications are equivalent to equivalent implementation cases and fall within the scope of the technical solution.

Claims

1. A method for determining reservoir capacity by contour method taking river island into consideration, characterized in that: The following steps are involved: Obtain each elevation information in the reservoir area and its corresponding closed contour area; Perform self-intersection processing on the closed contour area under each elevation information to obtain the completely overlapping area in each closed contour area; Determine the river island contour area corresponding to each elevation information according to each completely overlapping area; The reservoir capacity is determined based on each elevation information and the corresponding closed contour area and river island contour area.

2. The method for determining reservoir capacity by contour method taking river island into consideration according to claim 1, characterized in that: The reservoir capacity is determined based on each elevation information and the corresponding closed contour area and river island contour area, specifically: Determine the effective contour size information corresponding to each elevation information according to the closed contour area and the river island contour area corresponding to each elevation information; The reservoir capacity is determined based on each elevation information and its corresponding effective contour size information.

3. The method for determining reservoir capacity by contour method taking river island into consideration according to claim 2, characterized in that: According to the closed contour area and river island contour area corresponding to each elevation information, the effective contour size information corresponding to each elevation information is determined, specifically: Determine the valid contour area corresponding to each elevation information according to the closed contour area and the river island contour area corresponding to each elevation information; The effective contour line size information corresponding to each elevation information is determined according to each effective contour line area.

4. The method for determining reservoir capacity by contour method taking river island into consideration according to claim 2, characterized in that: According to the closed contour area and river island contour area corresponding to each elevation information, the effective contour size information corresponding to each elevation information is determined, specifically: Determine the closed contour size information corresponding to each elevation information according to the closed contour area corresponding to each elevation information; Determine the river island contour size information corresponding to each elevation information according to the river island contour area corresponding to each elevation information; The effective contour line size information corresponding to each elevation information is determined according to the closed contour line size information and the river island contour line size information corresponding to each elevation information.

5. The method for determining reservoir capacity by contour method taking river island into consideration according to claim 1, characterized in that: Get each elevation information in the reservoir area and its corresponding closed contour area, specifically: Obtain the overall closed contour area within the reservoir area; The overall closed contour area is segmented according to the elevation field to obtain each elevation information in the reservoir area and its corresponding closed contour area.

6. The method for determining reservoir capacity by contour method taking river island into consideration according to claim 2, characterized in that: The reservoir capacity is determined based on each elevation information and its corresponding effective contour size information, including: Integrate each elevation information and its corresponding effective contour size information to obtain the actual overall closed contour information; Determine the volume between each group of adjacent contour lines based on the actual overall closed contour line information; The reservoir capacity is determined based on the volume between multiple adjacent contour lines.

7. The method for determining reservoir capacity by contour method taking river island into consideration according to claim 5, characterized in that: Get the overall closed contour area within the reservoir area, specifically: Acquire the reservoir area contour area, wherein the reservoir area contour area includes closed contour area and unclosed contour area under all elevation information; Closing the unclosed contour area; The closed contour area and the unclosed contour area after the closure process are taken as the overall closed contour area in the reservoir area.

8. A system for determining reservoir capacity by contour method taking river island into consideration, characterized in that: include: The acquisition module is used to obtain each elevation information in the reservoir area and its corresponding closed contour area; A processing module is used to perform self-intersection processing on the closed contour area under each elevation information to obtain a completely overlapping area in each closed contour area; A preset area determination module is used to determine the river island contour area corresponding to each elevation information according to each completely overlapping area; The reservoir capacity determination module is used to determine the reservoir capacity based on each elevation information and the corresponding closed contour area and river island contour area.

Images