CN113742637A - Method and device for calculating annual average silt loss rate of reservoir, electronic equipment and storage medium - Google Patents
Method and device for calculating annual average silt loss rate of reservoir, electronic equipment and storage medium Download PDFInfo
- Publication number
- CN113742637A CN113742637A CN202110937105.XA CN202110937105A CN113742637A CN 113742637 A CN113742637 A CN 113742637A CN 202110937105 A CN202110937105 A CN 202110937105A CN 113742637 A CN113742637 A CN 113742637A
- Authority
- CN
- China
- Prior art keywords
- reservoir
- target
- annual average
- loss rate
- basin
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 76
- 238000004364 calculation method Methods 0.000 claims abstract description 35
- 238000004422 calculation algorithm Methods 0.000 claims abstract description 23
- 238000004062 sedimentation Methods 0.000 claims description 16
- 230000008021 deposition Effects 0.000 claims description 13
- 238000012544 monitoring process Methods 0.000 description 13
- 238000012806 monitoring device Methods 0.000 description 11
- 238000010586 diagram Methods 0.000 description 9
- 238000012545 processing Methods 0.000 description 7
- 239000013049 sediment Substances 0.000 description 6
- 230000008569 process Effects 0.000 description 5
- 238000004590 computer program Methods 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 238000004891 communication Methods 0.000 description 3
- 238000011156 evaluation Methods 0.000 description 3
- 230000006870 function Effects 0.000 description 3
- 238000005516 engineering process Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 238000009825 accumulation Methods 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000001413 cellular effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000012790 confirmation Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000013500 data storage Methods 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000003993 interaction Effects 0.000 description 1
- 239000011159 matrix material Substances 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 238000000611 regression analysis Methods 0.000 description 1
- 238000009877 rendering Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Mathematical Physics (AREA)
- Data Mining & Analysis (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Computational Mathematics (AREA)
- Pure & Applied Mathematics (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Algebra (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
The application discloses a method and a device for calculating the annual average silt loss rate of a reservoir, electronic equipment and a storage medium. The method comprises the following steps: the method comprises the steps of obtaining siltation data of a plurality of reservoirs in N target drainage basins, wherein the N target drainage basins are located in a target area; calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir; calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir; determining the actual total storage capacity of the reservoirs in each target drainage basin based on the storage capacity data of each reservoir in each target drainage basin; calculating the annual average silt loss rate of the reservoir in the target area by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total reservoir capacity of the reservoir in each target basin, wherein N is a positive integer; the problem that the calculation error of the direct method is large due to the fact that the sample is missing and the difference of the drainage basin is large at present is solved.
Description
Technical Field
The invention relates to the technical field of hydraulic engineering, in particular to a method and a device for calculating the annual average silt loss rate of a reservoir, electronic equipment and a storage medium.
Background
Because the river sand content in China is high, the storage capacity lost due to sedimentation in China is large every year, a quantitative method for the sedimentation characteristic of the reservoir is provided according to the sedimentation characteristic of the reservoir in China, great theoretical guidance significance is provided for the sedimentation evaluation and treatment work of the reservoir, and the dynamic sedimentation rate of the reservoir is generally reflected by the annual average sedimentation loss rate of the reservoir at present.
In the actual process of calculating the annual average silt loss rate of the reservoir, the annual average silt loss rate of the reservoir is calculated by adopting a direct method at present, the annual average silt loss rate of a single reservoir in a determined time period is calculated according to the collected reservoir sediment deposition and the total reservoir capacity of the reservoir, and the annual average silt loss rate of the national reservoir is calculated by combining reservoir capacity data and annual average silt loss rate data. However, the related technology lacks comprehensive reservoir silt loss data, and the difference of the number of different basin samples is large, so that the annual average silt loss rate of the national reservoirs is not high in accuracy.
Disclosure of Invention
In view of this, the present application provides a method and an apparatus for calculating an annual average silt loss rate of a reservoir, an electronic device, and a storage medium, so as to improve the calculation accuracy of the annual average silt loss rate of the reservoir.
In a first aspect, an embodiment of the present application provides a method for calculating an annual average silt loss rate of a reservoir, including: the method comprises the steps of obtaining siltation data of a plurality of reservoirs in N target drainage basins, wherein the N target drainage basins are located in a target area;
calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir;
calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the annual average silt loss rate of each reservoir and the storage capacity data corresponding to each reservoir;
determining the actual total storage capacity of the reservoirs in each target drainage basin based on the storage capacity data of each reservoir in each target drainage basin;
and calculating the annual average silt loss rate of the reservoir in the target region within a preset time period by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in the preset time period of each target basin and the actual total reservoir capacity of the reservoir in each target basin, wherein N is a positive integer.
In a second aspect, an embodiment of the present application provides a device for calculating an annual average silt loss rate of a reservoir, including: the device comprises an acquisition module, a storage module and a storage module, wherein the acquisition module is used for acquiring the sedimentation data of a plurality of reservoirs in N target drainage basins, and the N target drainage basins are located in a target area. And the first calculation module is used for calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir. And the second calculation module is used for calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir. And the third calculation module is used for determining the actual total storage capacity of the reservoirs in each target drainage basin based on the storage capacity data of each reservoir in each target drainage basin. And the fourth calculation module is used for calculating the annual average silt loss rate of the reservoir in the target area by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total reservoir capacity of the reservoir in each target basin, wherein N is a positive integer.
In a third aspect, an embodiment of the present application provides an electronic device, including: at least one processor and storage;
the processor is used for executing the computer program stored in the storage to realize the method for calculating the annual average silt loss rate of the reservoir as described in any one of the embodiments of the first aspect.
In a fourth aspect, embodiments of the present application provide a computer storage medium, where one or more programs are stored, and the one or more programs may be executed by the electronic device described in the third aspect, so as to implement the method for calculating an annual average silt loss rate of a reservoir described in any one of the embodiments of the first aspect.
The application has the following advantages and beneficial effects:
according to the method, the device, the electronic equipment and the storage medium for calculating the annual average silt loss rate of the reservoir, the annual average silt loss rate of each basin in each target basin in a target area (such as the whole country) is calculated by using a stock weighted average method, the integral characteristic of reservoir siltation distribution in each target basin is obtained, the calculation result has basin representativeness, and the annual average silt loss rate of the reservoir in the target basin is weighted and averaged by combining the total reservoir capacity of the reservoir in each target basin, so that the annual average silt loss rate of the reservoir in the target area is calculated, and the calculation accuracy of the annual average silt loss rate of the reservoir in the target area is improved.
Drawings
The accompanying drawings, which are included to provide a further understanding of the embodiments of the invention and are incorporated in and constitute a part of this application, illustrate embodiment(s) of the invention and together with the description serve to explain the principles of the invention. In the drawings:
fig. 1 is a schematic view illustrating an application scenario of a method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present application;
fig. 2 is a schematic flow chart illustrating a method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present application;
FIG. 3 is a schematic flow chart illustrating a method for calculating the annual average silt loss rate of a reservoir according to another embodiment of the present application;
FIG. 4 is a schematic flow chart illustrating a method for calculating the annual average silt loss rate of a reservoir according to another embodiment of the present application;
FIG. 5 is a flow chart illustrating a method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present application;
FIG. 6 is a schematic diagram showing calculation results of annual average silt loss rates of reservoirs in various domestic watersheds, which is provided in the embodiment of the application;
fig. 7 is a block diagram illustrating a computing device for an annual average silt loss rate of a reservoir according to an embodiment of the present application;
fig. 8 is a structural block diagram of a second computing module in the computing apparatus for an annual average silt loss rate of a reservoir according to an embodiment of the present application;
fig. 9 is a block diagram illustrating a fourth computing module in a computing apparatus for an annual average silt loss rate of a reservoir according to an embodiment of the present application;
fig. 10 is a block diagram showing a configuration of an electronic device for executing a method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present application;
fig. 11 shows a storage unit for storing or carrying program codes for implementing the method for calculating the annual average silt loss rate of a reservoir according to the embodiment of the present application.
Detailed Description
In order to make the objects, technical solutions and advantages of the present invention more apparent, the present invention is further described in detail below with reference to examples and accompanying drawings, and the exemplary embodiments and descriptions thereof are only used for explaining the present invention and are not meant to limit the present invention.
The technical solutions in the embodiments of the present application will be described below in a clear and complete manner with reference to the drawings in the embodiments of the present application, and it is obvious that the described embodiments are only a part of the embodiments of the present application, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments given herein without making any creative effort, shall fall within the protection scope of the present application.
The reservoir siltation loss refers to the fact that a large amount of river sand silts in a reservoir area, the reservoir capacity is damaged, the reservoir siltation loss seriously affects the operation efficiency and the stability of an ecological system of a reservoir, at present, China is the country with the largest quantity of reservoirs in the world, the reservoir siltation situation of China is very severe at the same time, in view of the characteristic that the basin of China has high sand content, a quantification method is provided for reservoir siltation characteristics, important theoretical guidance significance is provided for reservoir siltation evaluation and treatment work, and at present, the annual average siltation rate of the reservoir is usually used as an index of the annual average reservoir capacity loss degree of the reservoir in a certain period.
In the research on the quantitative calculation of the current related reservoir silt loss, the applicant finds that the calculation of the annual average silt loss rate of the reservoir needs to be concerned with the following steps: the quantity and the operation stage of the reservoirs are in a dynamic change process along with the increase of time, and the difference of reservoir samples in different watersheds is large; a large amount of reservoir sedimentation data lack large-scale statistics, so the annual average sedimentation loss rate of reservoirs in the watershed and the whole country can be analyzed only by combining part of representative reservoirs. In view of the above reasons, the calculation accuracy of the annual average silt loss rate of the reservoir in China is greatly influenced by statistical samples. The calculation result is lack of demonstration in a macroscopic level, and the problem of low calculation precision exists.
Specifically, the annual average silt loss rate is used as an evaluation index, and the annual average silt loss rate of a reservoir in China is calculated by a direct method in the related technology, wherein the method specifically adopts the following method for calculation:
calculating the annual average silt loss rate of the single reservoir in the determined time period according to the collected reservoir sediment deposition and the total reservoir capacity, wherein the specific calculation formula is as follows:
wherein k represents the number of the whole reservoir in the statistical sample, RkThe annual average silt loss rate, Δ V, of the kth reservoirkThe silting and loss storage capacity V of the kth reservoir in the statistical timekIs the storage capacity of the kth reservoir, tkThe year was counted for the deposition in the kth reservoir.
And calculating the annual average silt loss rate of the national reservoirs by counting the number of all the silted reservoirs in the sample silted data, wherein the calculation formula is as follows:
wherein c is the number of all the reservoir deposits in the statistical sample.
Therefore, the embodiment of the application provides a method for calculating the annual average silt loss rate of a reservoir, the method includes the steps of firstly obtaining siltation data of a plurality of reservoirs in N target basins, wherein the N target basins are located in the target areas, then calculating the silt loss rate of each reservoir according to the siltation data of each reservoir, calculating the annual average silt loss rate of each reservoir in each target basin by adopting a weighted average algorithm according to the silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir, then determining the actual total reservoir capacity of each reservoir in each target basin according to the reservoir annual average silt loss rate of each target basin and the actual total reservoir capacity of each reservoir in each target basin, and finally calculating the reservoir silt loss rate of a preset time period in each target area by adopting a weighted average algorithm according to the annual average silt loss rate of each reservoir in each target basin and the actual total reservoir capacity of each reservoir in each target basin, wherein N is a positive integer. The annual average silt loss rate of the reservoir in the target area can be calculated, and the calculation precision of the annual average silt loss rate of the reservoir in the target area is improved.
The following introduces an application environment of the method for calculating the annual average silt loss rate of the reservoir provided by the embodiment of the application:
referring to fig. 1, fig. 1 is a schematic view of an application scenario of a method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present disclosure, where the method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present disclosure may be applied to a system 100 for monitoring an annual average silt loss rate of a reservoir, and taking a server as an example of an electronic device in the present disclosure, the system 100 for monitoring an annual average silt loss rate of a reservoir includes the server 110, a router 120, and a reservoir intelligent monitoring device 130. The server 110, the router 120 and the intelligent reservoir monitoring device 130 may be connected to each other via the internet to form the internet of things.
The intelligent reservoir monitoring devices 130 may include a plurality of devices, and the intelligent reservoir monitoring devices 130 may be acquisition devices for acquiring data related to the annual average silt loss rate of the reservoir, such as terminal devices including a reservoir silt amount data acquisition device, a reservoir sediment discharge amount data acquisition device, and a total reservoir capacity data acquisition device. The server 110 may be a server for managing the reservoir intelligent monitoring device 130, for example, a server for statistically collecting data, and specifically, the server 110 may be a cloud server. The server 110 and the intelligent reservoir monitoring device 130 form an internet of things, and the information acquired by the intelligent reservoir monitoring device 130 is acquired to perform online control monitoring and the like on the intelligent reservoir monitoring device 130, for example, the online time length and other related information of the data acquired by the intelligent reservoir monitoring device 130 controls the working state of the intelligent reservoir monitoring device 130, and pushes the related information of the intelligent reservoir monitoring device 130 to corresponding staff of a reservoir workstation and the like.
It should be noted that the configuration shown in fig. 1 does not constitute a limitation on the system 100 for monitoring the annual average loss rate of a reservoir, and in other embodiments, the system 100 may include fewer or more devices or components than those shown, or may constitute some devices, or a different arrangement of components.
Embodiments of the present application will be described in detail below with reference to the accompanying drawings.
Referring to fig. 2, fig. 2 is a schematic flow chart of a method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present application, where the method for calculating an annual average silt loss rate of a reservoir according to an embodiment of the present application is applied to a system 100 for monitoring an annual average silt loss rate of a reservoir, and the method includes:
step S10, the sedimentation data of a plurality of reservoirs in N target basins are obtained, wherein the N target basins are located in the target area.
In the embodiment of the application, the N target watersheds are adaptively adjusted according to the selected target area, specifically, the N target watersheds can be classified according to the ground water collecting area and the underground water collecting area, and can be selected and adjusted according to the actual data use requirement of the target area. The types of the deposition data of a plurality of reservoirs are more, each reservoir comprises different reservoir parameters, construction age, sediment discharge amount, sediment deposition amount and other parameters, and the annual average deposition and loss rate correlation parameters of the reservoirs in different target areas are different, for example, the reservoir deposition of China is mainly based on the sediment content, and the deposition data of each reservoir can comprise water area ecological parameters related to the deposition and loss of each country and the like.
For example, in China, the number of reservoirs is the largest in China, and the storage capacity is 10 ten thousand meters3The reservoir is provided with 9.8 thousands seats. Meanwhile, reservoir silting conditions in China are very severe, and particularly river sand is used as the reservoir silting characteristic in China to select reservoir data of each basin.
As an implementation manner, the deposition data of a plurality of reservoirs in the N target watersheds are stored in advance in servers corresponding to the reservoir workstations of each watershed, the intelligent reservoir monitoring equipment continuously collects real-time data, and the servers acquire the real-time data and update the deposition data.
And step S20, calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir.
In the embodiment of the application, the calculation time of the annual average silt loss rate can be automatically acquired through the time of updating the connection information of the router or the intelligent reservoir monitoring equipment. The embodiments of the present application are not further limited thereto, which may be in accordance with certain practical sections. In addition, the preset time period can also be preset by a hydropower station worker through a server, for example, the worker can index the siltation data of each reservoir in a certain time period in the server and calculate the siltation rate.
Whether the router or the intelligent reservoir monitoring equipment is automatically set and acquired or preset by a worker, the preset time period can be an absolute time point interval. For example, the time information and the accumulation data information are updated and transmitted to the server for a fixed period of time (e.g., 2010, 05 months) consisting of a year, a week, a month, etc.
In addition, when the siltation data and the time information are transmitted between the server and the server, the server sends a server identifier or a password or simultaneously contains the server identifier and the password, and specifically, the server identifier corresponds to a certain reservoir of different watersheds.
And step S30, calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the silt loss rate of each reservoir and the storage capacity data corresponding to each reservoir.
In the embodiment of the application, the silt loss rate of each reservoir and the storage capacity data corresponding to each reservoir are counted by the server and corresponding associated information is carried out. The information association can be manual association of workers through a server, or automatic association can be performed according to server identification and corresponding storage data, the server comprises a monitoring node in the information association, and the monitoring node can be set to be intelligent algorithm judgment or artificial judgment. In an example, regression analysis can be performed according to historical data to judge the annual average silt loss rate value range of the reservoir.
Step S40, determining the actual total storage capacity of the reservoirs in each target basin based on the storage capacity data of each reservoir in each target basin.
In the embodiment of the application, the server performs accounting on the storage capacity data of different reservoirs of each target drainage basin, and summarizes the actual total storage capacity of the reservoirs.
And step S50, calculating the annual average silt loss rate of the reservoir in the target area within a preset time period by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total storage capacity of the reservoir in each target basin, wherein N is a positive integer.
According to the method for calculating the annual average silt loss rate of the reservoir, the annual average silt loss rate of each basin in each target basin in a target area (such as the whole country) is calculated by using a inventory weighted average method, the integral characteristics of the silting distribution of the reservoir in each target basin are obtained, the calculation result has basin representativeness, and the annual average silt loss rate of the reservoir in the target basin is weighted and averaged by combining the total reservoir capacity of the reservoir in each target basin, so that the annual average silt loss rate of the reservoir in the target area is calculated, and the calculation precision of the annual average silt loss rate of the reservoir in the target area is improved.
In one embodiment, the fouling profile comprises: each target basin identification, after step S10, the method further comprising: and classifying the sedimentation data of the plurality of reservoirs based on the identification of each target basin to obtain the sedimentation data of each reservoir in each target basin.
In this embodiment, the server performs identification processing on a plurality of reservoirs in different watersheds, and sums up the reservoirs identified by the servers in step S20 of the embodiment to corresponding watersheds for secondary identification processing, specifically, the identification processing includes setting a first character set and a second character set, the server processes each watershed according to a first character set label, and marks different reservoirs according to serial numbers on the basis of the first character set, specifically, the first character set and the second character set include symbols such as numbers and letters.
Preferably, the fouling data includes: the siltation loss storage capacity and the total storage capacity of internal statistics, the siltation loss rate of each reservoir is calculated according to the siltation data of each reservoir, and the siltation loss rate comprises the following steps:
calculating the silt loss rate of each reservoir by adopting the following formula:
wherein R isijThe annual average silt loss rate of the jth reservoir in the ith target watershed, i is the target watershed to which the reservoir belongs, j is the serial number of the reservoir in the target watershed, and delta VijIs the silt loss storage capacity V of the jth reservoir in the ith target flow fieldijIs the total storage capacity, t, of the jth reservoir in the ith target flow fieldijAnd (3) calculating the deposition time of the jth reservoir in the ith basin, wherein i is less than or equal to N.
Referring to fig. 3, fig. 3 is a schematic flow chart of another method for calculating an annual average silt loss rate of a reservoir according to the embodiment of the present application, which is applied to a system 100 for monitoring an annual average silt loss rate of a reservoir, in step S30: calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir, and the method comprises the following steps:
s302: determining the total storage capacity of the reservoirs of each target basin based on the total storage capacity of each reservoir;
s304: determining the silt loss storage capacity of each target basin based on the annual average silt loss rate of each reservoir and the total storage capacity of each reservoir;
s306: and dividing the silt loss storage capacity of each target basin by the total storage capacity of the corresponding reservoir of the target basin to obtain the annual average silt loss rate of the reservoir of the target basin.
As a preference of the above embodiment, the step S306: dividing the silt loss storage capacity of each target basin by the total storage capacity of the corresponding reservoir of the target basin to obtain the annual average silt loss rate of the reservoir of the target basin, and the method comprises the following steps:
calculating the annual average silt loss rate of the reservoir of each target watershed by adopting the following formula:
wherein R isiIs the annual average silt loss rate of the reservoir of the ith target watershed, m is the total number of investigated reservoirs in the ith target watershed, j is the serial number of the reservoir in the target watershed, wherein j is a positive integer,
is the total storage capacity of the reservoir of the ith target basin,
the silt storage capacity of the ith target basin.
Referring to fig. 4, fig. 4 is a schematic flow chart of another method for calculating an annual average silt loss rate of a reservoir according to the embodiment of the present application, which is applied to a system 100 for monitoring an annual average silt loss rate of a reservoir, in the step S50: calculating the annual average silt loss rate of the reservoir in the target region in a preset time period by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total storage capacity of the reservoir in each target basin, wherein the method comprises the following steps:
s502: determining the silt loss storage capacity in the target area based on the annual average silt loss rate of the reservoir of each target basin and the actual total storage capacity of the reservoir of each target basin;
s504: determining the actual total storage capacity of the reservoir in the target area based on the actual total storage capacity of the reservoir in each target drainage basin;
s506: and dividing the actual total storage capacity of the reservoir in the target area by the storage capacity of the silt in the target area to determine the annual average silt rate of the reservoir in the target area within a preset time period.
As a preference of the above embodiment, the step S506: dividing the silt loss storage capacity in the target area by the actual total storage capacity of the reservoir in the target area to determine the annual average silt loss rate of the reservoir in the target area within a preset time period, wherein the method comprises the following steps:
calculating the annual average silt loss rate of the reservoir in a preset time period in the target area by adopting the following formula:
wherein R represents the annual average silt loss rate of the reservoir in a preset time period in the target area, Vi' indicates the actual reservoir total storage capacity contained in the target basin i,
is the actual total storage capacity of the reservoir in the target area,
the silt storage capacity in the target area is obtained.
For example, the time period is set to year, and domestic is taken as a target area, please refer to fig. 5, fig. 5 is a schematic flow chart of a method for calculating an annual average silt loss rate of a reservoir in an embodiment of the present application, and the attached drawings and the above-mentioned examples are combinedThe embodiment step is that S1, the storage capacity V of the data reservoir is obtained according to the obtaining moduleijAnd the silt loss storage capacity delta V of the reservoirijAnd the statistical time delta t of reservoir sedimentation, S2, calculating the annual average silt loss rate of a single reservoir in the basin through a first calculation module
The annual average silt loss rate of the national reservoir can be calculated, S3, the total annual average silt loss rate of the basin i is calculated through a second calculation module
S4, calculating the actual total storage capacity V of the basin i through a third calculation modulei', S5 calculating the annual average silt loss rate of the national reservoir through a fourth calculation module
The application provides a method, a device, electronic equipment and a storage medium for calculating the annual average silt loss rate of a reservoir, which are exemplified by taking the annual average silt loss rate as a preset time period and taking a domestic reservoir as a calculation target of the annual average silt loss rate of the reservoir, firstly obtaining the siltation data of a plurality of reservoirs in a domestic target basin, then calculating the annual average silt loss rate of each reservoir according to the reservoir capacity data of each reservoir, calculating the annual average silt loss rate of each target basin by adopting a weighted average algorithm according to the annual average silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir, then determining the actual total reservoir capacity of each target basin according to the reservoir capacity data of each reservoir in each target basin, and finally calculating the annual average silt loss rate of the domestic reservoir by adopting a weighted average reservoir algorithm according to the annual average silt loss rate of each target basin and the actual total reservoir capacity of each target basin, the method and the device reduce the large deviation in calculating the annual average silt loss rate of the reservoir caused by basin difference and the quantity proportion of the collected samples of each basin, and provide a reliable processing method for calculating and applying the annual average silt loss rate of the reservoir in the target area.
Please refer to fig. 6, fig. 6 is a schematic diagram of an annual average silt loss rate calculation result of each basin in China according to an embodiment of the present application, where a label 1 is an annual average silt loss rate value of a reservoir calculated by a direct method mentioned in the related art, and a label 2 is an annual average silt loss rate value of a reservoir obtained by a method for calculating an annual average silt loss rate of a reservoir provided by the present application, and specifically, the annual average silt loss rate of each basin is calculated according to siltation data of 6702 seats of seven basins in China, where the annual average silt loss rate of a reservoir calculated by a direct method is 0.49%, and the result calculated by a weighted average method is 0.41%, so that the present application can improve the calculation accuracy of the annual average silt loss rate of a reservoir, and has a wide application prospect.
Referring to fig. 7, an apparatus 700 for calculating an annual average silt loss rate of a reservoir, provided in the embodiment of the present application, is applied to a system 100 for monitoring an annual average silt loss rate of a reservoir, the apparatus includes:
an obtaining module 710, configured to obtain siltation data of a plurality of reservoirs in N target basins, where the N target basins are located in a target area;
the first calculating module 720 is used for calculating the silt loss rate of each reservoir according to the storage capacity data of each reservoir;
the second calculating module 730 is used for calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir;
the third calculation module 740 is configured to determine the actual total storage capacity of the reservoirs in each target drainage basin based on the storage capacity data of each reservoir in each target drainage basin;
and a fourth calculating module 750, configured to calculate the annual average silt loss rate of the reservoir in the target area in a preset time period by using a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target drainage basin and the actual total reservoir capacity of the reservoir in each target drainage basin, where N is a positive integer.
Referring to fig. 8, fig. 8 is a structural block diagram of a second calculating module in the calculating apparatus for an annual average silt loss rate of a reservoir according to the embodiment of the present application, where the second calculating module 730 includes:
a first confirmation module 731, configured to determine a total storage capacity of the reservoirs of each target watershed based on the total storage capacity of each reservoir;
a second confirming module 732, configured to determine the silt storage capacity of each target basin based on the silt rate of each reservoir and the total storage capacity of each reservoir;
the second sub-calculation module 733, configured to divide the silt storage capacity of each target drainage basin by the total storage capacity of the reservoir of the corresponding target drainage basin to obtain an annual average silt rate of the reservoir of the target drainage basin.
Referring to fig. 9, fig. 9 is a structural block diagram of a fourth calculating module in the calculating apparatus for an annual average silt loss rate of a reservoir according to the embodiment of the present application, where the fourth calculating module 750 includes:
the first acquisition module 751 is used for determining the silt loss storage capacity in the target area based on the annual average silt loss rate of the reservoirs in each target basin and the actual total storage capacity of the reservoirs in each target basin;
the second obtaining module 752 is configured to determine an actual total storage capacity of the reservoir in the target area based on the actual total storage capacity of the reservoir in each target drainage basin;
and a third sub-calculation module 753, configured to divide the silt storage capacity in the target area by the actual total storage capacity of the reservoir in the target area to determine an annual average silt rate of the reservoir in the target area within a preset time period.
It should be noted that the device embodiment in the present application corresponds to the foregoing method embodiment, and specific principles in the device embodiment may refer to the contents in the foregoing method embodiment, which is not described herein again.
Referring to fig. 10, based on the above, an electronic device 200 including the method for calculating the annual average silt loss rate of the reservoir may be further provided, where the electronic device 200 may be a smart phone, a tablet computer, a portable computer, or the like.
The electronic device 200 also includes a processor 202 and memory 204. The storage 204 stores programs that can execute the content of the foregoing embodiments, and the processor 202 can execute the programs stored in the storage 204.
The Memory 204 may include a Random Access Memory (RAM) or a Read-Only Memory (Read-Only Memory). The storage 204 may be used to store instructions, programs, code sets, or instruction sets. Memory 204 may include a program storage area and a data storage area, where the program storage area may store instructions for implementing an operating system, instructions for implementing at least one function (e.g., instructions for a user to obtain a random number), instructions for implementing the various method embodiments described below, and the like. The stored data area may also store data (e.g., random numbers) created by the terminal in use, and the like.
The electronic device 200 may further include a network module for receiving and transmitting electromagnetic waves, and implementing interconversion between the electromagnetic waves and the electrical signals, so as to communicate with a communication network or other devices, for example, an audio playing device. The network module may include various existing circuit elements for performing these functions, such as an antenna, a radio frequency transceiver, a digital signal processor, an encryption/decryption chip, a Subscriber Identity Module (SIM) card, memory, and so forth. The network module may communicate with various networks such as the internet, an intranet, a wireless network, or with other devices via a wireless network. The wireless network may comprise a cellular telephone network, a wireless local area network, or a metropolitan area network. The screen can display the interface content and perform data interaction.
Referring to fig. 11, a block diagram of a computer-readable storage medium according to an embodiment of the present application is shown. The computer readable medium 900 has stored therein a program code that can be called by a processor to execute the method described in the above method embodiments.
The computer-readable storage medium may be an electronic memory such as a flash memory, an EEPROM (electrically erasable programmable read only memory), an EPROM, a hard disk, or a ROM. Alternatively, the computer-readable storage medium includes a non-volatile computer-readable storage medium. The computer readable storage medium has a storage space for program code for performing any of the method steps of the above-described method. The program code may be read from or written to one or more computer program products. The program code may be compressed, for example, in a suitable form.
Embodiments of the present application also provide a computer program product or computer program comprising computer instructions stored in a computer-readable storage medium. The processor of the computer device reads the computer instructions from the computer readable storage medium, and the processor executes the computer instructions to enable the computer device to execute the method for calculating the annual average silt loss rate of the reservoir described in the various optional implementation modes.
The above embodiments are only used to illustrate the technical solutions of the present application, and not to limit the same; although the present application has been described in detail with reference to the foregoing embodiments, it will be understood by those of ordinary skill in the art that: the technical solutions described in the foregoing embodiments may still be modified, or some technical features may be equivalently replaced; such modifications and substitutions do not necessarily depart from the spirit and scope of the corresponding technical solutions in the embodiments of the present application.
Claims (10)
1. A method for calculating the annual average silt loss rate of a reservoir is characterized by comprising the following steps:
the method comprises the steps of obtaining siltation data of a plurality of reservoirs in N target drainage basins, wherein the N target drainage basins are located in a target area;
calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir;
calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the annual average silt loss rate of each reservoir and the storage capacity data corresponding to each reservoir;
determining the actual total storage capacity of the reservoirs in each target drainage basin based on the storage capacity data of each reservoir in each target drainage basin;
and calculating the annual average silt loss rate of the reservoir in the target area by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total reservoir capacity of the reservoir in each target basin, wherein N is a positive integer.
2. The method of claim 1, wherein the siltation data includes: each target basin is identified, and the method further comprises the following steps:
and classifying the sedimentation data of the plurality of reservoirs based on the identification of each target basin to obtain the sedimentation data of each reservoir in each target basin.
3. The method of claim 2, wherein the siltation data includes: the silt loss storage capacity and the total storage capacity which are counted in the preset time period are used for calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir, and the method comprises the following steps:
the annual average silt loss rate of each reservoir is calculated by the following formula:
wherein R isijThe annual average silt loss rate of the jth reservoir in the ith target watershed, i is the target watershed to which the reservoir belongs, j is the serial number of the reservoir in the target watershed, and delta VijIs the silt loss storage capacity V of the jth reservoir in the ith target flow fieldijIs the total storage capacity, t, of the jth reservoir in the ith target flow fieldijAnd (3) the deposition statistical time of the jth reservoir in the ith drainage basin is year, and i is less than or equal to N.
4. The method for calculating the annual average silt loss rate of the reservoir according to claim 3, wherein the calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the annual average silt loss rate of each reservoir and the reservoir capacity data corresponding to each reservoir comprises:
determining the total storage capacity of the reservoirs of each target basin based on the total storage capacity of each reservoir;
determining the annual average silt loss storage capacity of each target basin based on the annual average silt loss rate of each reservoir and the total storage capacity of each reservoir;
and dividing the annual average silt loss storage capacity of each target basin by the total storage capacity of the corresponding target basin to obtain the annual average silt loss rate of the reservoir of the target basin.
5. The method according to claim 4, wherein the step of dividing the silt storage capacity of each target basin by the total storage capacity of the corresponding target basin to obtain the annual average silt loss rate of the reservoir of the target basin comprises:
calculating the annual average silt loss rate of the reservoir of each target watershed by adopting the following formula:
wherein R isiIs the annual average silt loss rate of the reservoir of the ith target watershed, m is the total number of investigated reservoirs in the ith target watershed, j is the serial number of the reservoir in the target watershed, wherein j is a positive integer,
is the total storage capacity of the reservoir of the ith target basin,
the storage capacity of the silt in the ith target basin is the annual average silt loss.
6. The method for calculating the annual average silt loss rate of the reservoir according to claim 5, wherein the calculating the annual average silt loss rate of the reservoir in the target area by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total reservoir capacity of the reservoir in each target basin comprises:
determining the annual average silt loss storage capacity in the target area based on the annual average silt loss rate of the reservoir in each target basin and the actual total storage capacity of the reservoir in each target basin;
determining the actual total storage capacity of the reservoir in the target area based on the actual total storage capacity of the reservoir in each target drainage basin;
and dividing the annual average silt loss storage capacity in the target area by the actual total storage capacity of the reservoir in the target area to determine the annual average silt loss rate of the reservoir in the target area.
7. The method for calculating the annual average silt rate of the reservoir according to claim 6, wherein the step of dividing the annual average silt rate storage capacity in the target area by the actual total storage capacity of the reservoir in the target area to determine the annual average silt rate of the reservoir in the target area comprises the following steps:
calculating the annual average silt loss rate of the reservoir in the target area by adopting the following formula:
wherein R represents the annual average silt loss rate of the reservoir in the target area, Vi' indicates the actual reservoir total storage capacity contained in the target basin i,
is the actual total storage capacity of the reservoir in the target area,
the annual average silt loss storage capacity in the target area is shown.
8. A calculating device for the annual average silt loss rate of a reservoir is characterized by comprising:
the system comprises an acquisition module, a storage module and a storage module, wherein the acquisition module is used for acquiring the sedimentation data of a plurality of reservoirs in N target drainage basins, and the N target drainage basins are positioned in a target area;
the first calculation module is used for calculating the annual average silt loss rate of each reservoir according to the siltation data of each reservoir;
the second calculation module is used for calculating the annual average silt loss rate of the reservoir of each target basin by adopting a weighted average algorithm based on the annual average silt loss rate of each reservoir and the storage capacity data corresponding to each reservoir;
the third calculation module is used for determining the actual total storage capacity of the reservoirs in each target drainage basin based on the storage capacity data of each reservoir in each target drainage basin;
and the fourth calculation module is used for calculating the annual average silt loss rate of the reservoir in the target area by adopting a weighted average algorithm according to the annual average silt loss rate of the reservoir in each target basin and the actual total reservoir capacity of the reservoir in each target basin, wherein N is a positive integer.
9. An electronic device, comprising:
one or more processors;
a reservoir;
one or more programs, wherein the one or more programs are stored in the storage and configured to be executed by the one or more processors, the one or more programs configured to perform the method of calculating an annual average silt rate of a reservoir of any of claims 1-7.
10. A computer-readable storage medium storing program code that can be invoked by a processor to perform the method of calculating an annual average reservoir silt rate according to any one of claims 1 to 7.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110937105.XA CN113742637B (en) | 2021-08-16 | 2021-08-16 | Calculation method and device for annual average silt loss rate of reservoir, electronic equipment and storage medium |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110937105.XA CN113742637B (en) | 2021-08-16 | 2021-08-16 | Calculation method and device for annual average silt loss rate of reservoir, electronic equipment and storage medium |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN113742637A true CN113742637A (en) | 2021-12-03 |
| CN113742637B CN113742637B (en) | 2024-05-07 |
Family
ID=78731211
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202110937105.XA Active CN113742637B (en) | 2021-08-16 | 2021-08-16 | Calculation method and device for annual average silt loss rate of reservoir, electronic equipment and storage medium |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN113742637B (en) |
Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100563757B1 (en) * | 2006-01-23 | 2006-03-24 | 김은주 | Methods to measure pondage |
| CN108320095A (en) * | 2018-01-31 | 2018-07-24 | 黄河水利委员会黄河水利科学研究院 | A kind of reservoir sedimentation methods of risk assessment |
| CN111027226A (en) * | 2019-12-25 | 2020-04-17 | 中国电建集团中南勘测设计研究院有限公司 | Pumped storage power station water-sand numerical simulation calculation method |
| CN111723995A (en) * | 2020-06-24 | 2020-09-29 | 中国长江三峡集团有限公司 | Optimized scheduling method for sediment in flood season of reservoir under combined scheduling of cascade reservoirs |
| CN111967157A (en) * | 2020-08-14 | 2020-11-20 | 中国水利水电科学研究院 | Reservoir group water and sand optimization scheduling method and system considering staged multi-objective function |
| CN112287528A (en) * | 2020-10-16 | 2021-01-29 | 黄河勘测规划设计研究院有限公司 | Flood control high water level determination method for sandy river reservoir based on high beach trough |
| CN112287436A (en) * | 2020-10-16 | 2021-01-29 | 黄河勘测规划设计研究院有限公司 | Method and system for designing sediment section and effective storage capacity of reservoir of sandy river |
| CN112418491A (en) * | 2020-11-06 | 2021-02-26 | 黄河勘测规划设计研究院有限公司 | Dynamic configuration method for residual sediment storage capacity of reservoir |
| CN112561213A (en) * | 2021-02-23 | 2021-03-26 | 浙江贵仁信息科技股份有限公司 | Reservoir flood control scheduling method, system, equipment and medium based on improved POA |
-
2021
- 2021-08-16 CN CN202110937105.XA patent/CN113742637B/en active Active
Patent Citations (9)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| KR100563757B1 (en) * | 2006-01-23 | 2006-03-24 | 김은주 | Methods to measure pondage |
| CN108320095A (en) * | 2018-01-31 | 2018-07-24 | 黄河水利委员会黄河水利科学研究院 | A kind of reservoir sedimentation methods of risk assessment |
| CN111027226A (en) * | 2019-12-25 | 2020-04-17 | 中国电建集团中南勘测设计研究院有限公司 | Pumped storage power station water-sand numerical simulation calculation method |
| CN111723995A (en) * | 2020-06-24 | 2020-09-29 | 中国长江三峡集团有限公司 | Optimized scheduling method for sediment in flood season of reservoir under combined scheduling of cascade reservoirs |
| CN111967157A (en) * | 2020-08-14 | 2020-11-20 | 中国水利水电科学研究院 | Reservoir group water and sand optimization scheduling method and system considering staged multi-objective function |
| CN112287528A (en) * | 2020-10-16 | 2021-01-29 | 黄河勘测规划设计研究院有限公司 | Flood control high water level determination method for sandy river reservoir based on high beach trough |
| CN112287436A (en) * | 2020-10-16 | 2021-01-29 | 黄河勘测规划设计研究院有限公司 | Method and system for designing sediment section and effective storage capacity of reservoir of sandy river |
| CN112418491A (en) * | 2020-11-06 | 2021-02-26 | 黄河勘测规划设计研究院有限公司 | Dynamic configuration method for residual sediment storage capacity of reservoir |
| CN112561213A (en) * | 2021-02-23 | 2021-03-26 | 浙江贵仁信息科技股份有限公司 | Reservoir flood control scheduling method, system, equipment and medium based on improved POA |
Non-Patent Citations (9)
| Title |
|---|
| JIANGGUO CHEN 等: "An analysis of the factors affecting the annual average sedimentation loss rates of reservoirs", 《E-PROCEEDINGS OF THE 38TH IAHR WORLD CONGRESS》, 6 December 2019 (2019-12-06), pages 3098 - 3102 * |
| S.K.马宗德 , 刘敏: "喜马拉雅山区一座水库的使用期限", 水利水电快报, no. 23, 14 December 2001 (2001-12-14), pages 5 - 6 * |
| ZHAOYIN WANG 等: "Strategied for managing reservoir sedimentation", 《INTERNATIONAL JOURNAL OF SEDIMENT RESEARCH》, no. 4, 31 December 2009 (2009-12-31), pages 369 - 284 * |
| 何录合: "碧口水库泥沙淤积规律及控制", 西北电力技术, no. 04, 30 August 2000 (2000-08-30), pages 49 - 51 * |
| 吴虹娟, 杨国录, 余明辉: "清江隔河岩水库泥沙计算分析", 人民长江, no. 08, 20 August 2000 (2000-08-20), pages 27 - 28 * |
| 归力佳;顾圣平;林乐曼;咸京;李江林;: "基于组合赋权-理想点法的梯级水库蓄水研究", 人民黄河, no. 05, 10 May 2018 (2018-05-10), pages 50 - 54 * |
| 王育杰 等: "水库淤积物密实变化对冲淤计算与库容影响的研究", 《泥沙研究》, vol. 45, no. 3, 29 May 2020 (2020-05-29), pages 46 - 52 * |
| 赵定平;汤小云;何海涛;: "水库泥沙淤积量计算方法探讨", 科技资讯, no. 35, 13 December 2008 (2008-12-13), pages 241 - 242 * |
| 陈建国 等: "水库年平均淤损率影响因素分析", 《泥沙研究》, vol. 44, no. 4, 6 August 2019 (2019-08-06), pages 1 - 5 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN113742637B (en) | 2024-05-07 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN109948669A (en) | A kind of abnormal deviation data examination method and device | |
| CN110570133A (en) | building engineering dynamic cost analysis and control system based on big data | |
| CN103971289A (en) | Mine drawing data processing method and device | |
| CN112330065A (en) | Runoff forecasting method based on basic flow segmentation and artificial neural network model | |
| CN112668773A (en) | Method and device for predicting warehousing traffic and electronic equipment | |
| CN110836683A (en) | River channel management method, river channel management device and terminal | |
| CN107256231A (en) | A kind of Team Member's identification equipment, method and system | |
| CN106845831A (en) | A kind of Irrigation Forecast method, apparatus and system | |
| CN113688542A (en) | Intelligent optimization water resource configuration method and device, computer equipment and storage medium | |
| CN111428197B (en) | Data processing method, device and equipment | |
| CN112040507B (en) | Azimuth angle prediction method based on terminal sampling data AOA | |
| CN113742637A (en) | Method and device for calculating annual average silt loss rate of reservoir, electronic equipment and storage medium | |
| CN115842847B (en) | Intelligent control method, system and medium for water meter based on Internet of things | |
| CN113222245A (en) | Method and system for checking monthly electric quantity and electricity charge abnormity of residential user and storage medium | |
| CN109828253B (en) | Multi-station radar quantitative fusion target detection method | |
| CN114358548A (en) | Method and device for determining evaluation index and electronic equipment | |
| CN106550387B (en) | A kind of wireless sensor network routing layer QoS evaluating method | |
| CN114971053A (en) | Training method and device for online prediction model of network line loss rate of low-voltage transformer area | |
| CN112860996A (en) | Interest point processing method and device, electronic equipment and medium | |
| CN108093423A (en) | A kind of discovery method of base station location exception in user bill big data based on Ransac algorithms | |
| CN117391644B (en) | Parameter adjustment method, device, equipment and medium in contract management process | |
| CN114511225A (en) | Section identification method, device, equipment and storage medium | |
| CN116911466A (en) | Method, device, equipment and medium for determining park construction degree | |
| CN116011689A (en) | Path planning method and device, electronic equipment and storage medium | |
| CN116777674A (en) | Power distribution network data processing method and device, electronic equipment and storage medium |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |