Method, device and equipment for correcting water level-storage capacity relationship based on reservoir accumulation

Abstract

The application discloses a method, device and equipment for correcting water level-storage capacity relationship based on reservoir siltation, which comprises the following steps: obtaining initial storage capacity values corresponding to each water level layer of a reservoir; obtaining running state change data of each flow level of each water level layer of the reservoir within a period of time; calculating a relative siltation parameter of the reservoir according to a proportional relationship between the running state change data and initial running state data; calculating a siltation storage capacity according to the initial storage capacity values and the relative siltation parameter of the reservoir; and correcting actual storage capacity values of each water level layer by using the siltation storage capacity. The application completes siltation calculation and water level-storage capacity relationship correction by fully utilizing initial storage capacity data, daily running data and water level data of the reservoir, without complex external terrain measurement and data measurement at a specific time or scene, and the calculation process is simple. The calculated value is a proportional value based on an initial value, conforms to an existing economic dispatching system of a hydropower station, and has a low water level-storage capacity relationship correction cost.

Description

Technical Field

[0001] This invention belongs to the field of data processing technology, specifically relating to a method, apparatus, and equipment for correcting the relationship between water level and reservoir capacity based on reservoir siltation. Background Technology

[0002] With the operation of hydropower stations, in engineering practice, natural river cross-sections are affected by topography, geological conditions, flood changes, upstream sediment flushing, rainfall erosion, debris flows, landslides, and domestic or natural waste along the riverbanks, leading to reservoir siltation and a certain degree of reservoir capacity loss. This reduces the adjustable reservoir capacity, especially for small-capacity hydropower stations that do not have or do not fully have daily regulation capacity. The reservoir capacity loss caused by siltation directly changes the relationship between water level and reservoir capacity, thereby affecting emergency response time, water consumption rate, and economic benefits.

[0003] In existing technologies, the calculation of reservoir siltation is mainly carried out through methods such as ultrasonic measurement, underwater topographic measurement, and total station data acquisition. However, these methods are greatly affected by the precision of the instruments and equipment, the selection of the measurement section, the measurement method, the measurement location and density, etc., and the calculation process is complex. Summary of the Invention

[0004] The purpose of this invention is to provide a method, apparatus, and equipment for correcting the water level-reservoir capacity relationship based on reservoir siltation, in order to solve the technical problems in the prior art that are greatly affected by the precision of the instruments and equipment, the selection of the measurement section, the measurement method, the measurement location and density, and the complex calculation process.

[0005] To achieve the above objectives, the present invention adopts the following technical solution:

[0006] The first aspect provides a method for correcting the water level-reservoir capacity relationship based on reservoir siltation, including:

[0007] Obtain the initial reservoir capacity V corresponding to each water level layer in the reservoir area. 0i Where i = 1, 2, ..., n, and n represents the number of water level layers;

[0008] Acquire data on the operational status changes of each water level and flow rate in the reservoir area over a period of time. ij Where j = 1, 2, ..., m, m represents the number of flow levels;

[0009] Based on the operational status change data A ij Compared with the initial running status data A ij0 The proportional relationship between them is used to calculate the relative siltation parameter P of the reservoir area;

[0010] Based on the initial storage capacity value V 0iThe siltation capacity V is calculated from the relative siltation parameter P of the reservoir area, and the actual reservoir capacity V' of each water level layer is then calculated using the siltation capacity V. He Make corrections.

[0011] In one possible design, the operating state change data A ij At least include the water level change time A1 ij Total station active power output value A2 ij and the opening degree A3 of all station gates. ij .

[0012] In one possible design, based on the operational state change data A ij Compared with the initial running status data A ij0 The proportional relationship between them is used to calculate the relative siltation parameter P of the reservoir area, including:

[0013] Based on the data A of each operating state change ij With the corresponding initial running status data A ij0 The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor K. g The calculation formula is as follows:

[0014]

[0015] Where g = 1, 2, ..., q, q represents the number of types of running status change data, and Δt represents the time period for acquiring running status change data;

[0016] Calculate the sedimentation parameter factor K for each sedimentation parameter separately. g The corresponding parameter weight M g The calculation formula is as follows:

[0017]

[0018] in,

[0019] Based on each sedimentation parameter factor K g and the weight M of each parameter g The relative siltation parameter P of the reservoir area is calculated using the following formula:

[0020]

[0021] Where T represents the transpose of the matrix.

[0022] In one possible design, after calculating the relative sedimentation parameter P of the reservoir area, the method further includes:

[0023] Calculate the relative sedimentation parameter P for different reservoir areas during L time periods. ll = 1, 2, ..., L;

[0024] Based on the relative siltation parameter P of multiple reservoir areas l The average siltation parameter P in the reservoir area was calculated. aver The calculation formula is as follows:

[0025]

[0026] In one possible design, based on the initial storage capacity value V 0i The siltation capacity V is calculated from the relative siltation parameter P of the reservoir area, including:

[0027] Based on the initial storage capacity value V 0i Based on the relative siltation parameter P of the reservoir area, the reservoir capacity ΔV after siltation is calculated using the following formula:

[0028] ΔV=PV 0i (5)

[0029] Based on the reservoir capacity ΔV after siltation, the silted-up reservoir capacity V is calculated using the following formula:

[0030] V = (1-P)V 0i (6).

[0031] In one possible design, the siltation capacity V is used to calculate the actual reservoir capacity V' of each water level layer. He Make corrections, including:

[0032] Using the siltation capacity V, the actual reservoir capacity V' of each water level layer was calculated. He The correction is made, and the correction formula is as follows:

[0033]

[0034] Where e represents the e-th water level in layer E, and E represents the number of water level layers in the reservoir area, from the dead water level to the upper limit water level, based on a preset unit step size; H e This represents the calibration value of the water level at layer e. This represents the original reservoir capacity at the e-th water level.

[0035] In one possible design, the actual reservoir capacity V' of each water level layer is calculated using the siltation reservoir capacity V. He After modification, the method further includes:

[0036] For each water level layer, the actual reservoir capacity V' He Verification is required.

[0037] The second aspect provides a device for correcting the water level-reservoir capacity relationship based on reservoir siltation, including:

[0038] The initial reservoir capacity acquisition module is used to obtain the initial reservoir capacity value V corresponding to each water level layer in the reservoir area. 0i Where i = 1, 2, ..., n, and n represents the number of water level layers;

[0039] The data acquisition module is used to acquire data on the changes in the operating status of each water level and flow rate in the reservoir over a period of time. ij Where j = 1, 2, ..., m, m represents the number of flow levels;

[0040] The siltation parameter calculation module is used to calculate data A based on changes in operating status. ij Compared with the initial running status data A ij0 The proportional relationship between them is used to calculate the relative siltation parameter P of the reservoir area;

[0041] The storage capacity correction module is used to adjust the storage capacity based on the initial storage capacity value V. 0i The siltation capacity V is calculated from the relative siltation parameter P of the reservoir area, and the actual reservoir capacity V' of each water level layer is then calculated using the siltation capacity V. He Make corrections.

[0042] In one possible design, the operating state change data A ij At least include the water level change time A1 ij Total station active power output value A2 ij and the opening degree A3 of all station gates. ij .

[0043] In one possible design, based on the changing data A of the operating state... ij Compared with the initial running status data A ij0 When calculating the relative siltation parameter P of the reservoir area based on the proportional relationship between the two, the data acquisition module is specifically used for the following purposes:

[0044] Based on the data A of each operating state change ij With the corresponding initial running status data A ij0 The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor K. g The calculation formula is as follows:

[0045]

[0046] Where g = 1, 2, ..., q, q represents the number of types of running status change data, and Δt represents the time period for acquiring running status change data;

[0047] Calculate the sedimentation parameter factor K for each sedimentation parameter separately. g The corresponding parameter weight M g The calculation formula is as follows:

[0048]

[0049] in,

[0050] Based on each sedimentation parameter factor K g and the weight M of each parameter g The relative siltation parameter P of the reservoir area is calculated using the following formula:

[0051]

[0052] Where T represents the transpose of the matrix.

[0053] In one possible design, the sedimentation parameter calculation module is also used for:

[0054] Calculate the relative sedimentation parameter P for different reservoir areas during L time periods. l l = 1, 2, ..., L;

[0055] Based on the relative siltation parameter P of multiple reservoir areas l The average siltation parameter P in the reservoir area was calculated. aver The calculation formula is as follows:

[0056]

[0057] In one possible design, based on the initial storage capacity value V 0i When calculating the siltation capacity V based on the relative siltation parameter P of the reservoir area, the reservoir capacity correction module is specifically used for:

[0058] Based on the initial storage capacity value V 0i Based on the relative siltation parameter P of the reservoir area, the reservoir capacity ΔV after siltation is calculated using the following formula:

[0059] ΔV=PV 0i (5)

[0060] Based on the reservoir capacity ΔV after siltation, the silted-up reservoir capacity V is calculated using the following formula:

[0061] V = (1-P)V 0i (6).

[0062] In one possible design, the actual reservoir capacity V' of each water level layer is calculated using the siltation reservoir capacity V. He When making corrections, the warehouse capacity correction module is specifically used for:

[0063] Using the siltation capacity V, the actual reservoir capacity V' of each water level layer was calculated. He The correction is made, and the correction formula is as follows:

[0064]

[0065] Where e represents the e-th water level in layer E, and E represents the number of water level layers in the reservoir area, from the dead water level to the upper limit water level, based on a preset unit step size; H e This represents the calibration value of the water level at layer e. This represents the original reservoir capacity at the e-th water level.

[0066] In one possible design, the device further includes:

[0067] The verification module is used to verify the actual reservoir capacity V' of each water level. He Verification is required.

[0068] Thirdly, the present invention provides a computer device comprising a memory, a processor, and a transceiver connected in sequence and communication, wherein the memory is used to store a computer program, the transceiver is used to send and receive messages, and the processor is used to read the computer program and execute the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any possible design of the first aspect.

[0069] Fourthly, the present invention provides a computer-readable storage medium storing instructions that, when executed on a computer, perform the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any possible design of the first aspect.

[0070] Fifthly, the present invention provides a computer program product containing instructions that, when the instructions are executed on a computer, cause the computer to perform the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any possible design of the first aspect.

[0071] The advantages of this invention compared to the prior art are:

[0072] This invention can calculate reservoir siltation and correct the water level-reservoir ratio by fully utilizing initial reservoir capacity data, daily reservoir operation data, and water level data. In other words, it can quickly calculate the degree of siltation based on changes in daily operation data, thus eliminating the need for complex external topographic measurements or data measurements at specific times or in specific scenarios. The calculation process is very simple. At the same time, since the calculated value is a proportional value based on the initial value, it is fully in line with the existing economic dispatch system of hydropower stations, so there is no need to change the existing economic dispatch method, and the cost of correcting the water level-reservoir ratio is low. Attached Figure Description

[0073] Figure 1 A flowchart illustrating the method for correcting the water level-reservoir capacity relationship based on reservoir siltation, as provided in an embodiment of the present invention. Detailed Implementation

[0074] To more clearly illustrate the technical solutions in the embodiments of the present invention or the prior art, the present invention will be briefly introduced below in conjunction with the accompanying drawings and descriptions of the embodiments or the prior art. Obviously, the following description of the structure of the accompanying drawings is only some embodiments of the present invention. For those skilled in the art, other drawings can be obtained based on these drawings without creative effort. It should be noted that the description of these embodiments is for the purpose of helping to understand the present invention, but does not constitute a limitation of the present invention.

[0075] Example

[0076] To address the technical problems in existing technologies, which are significantly affected by factors such as the precision of instruments and equipment, the selection of measurement sections, measurement methods, measurement locations, and density, resulting in complex calculation processes, this application proposes a method for correcting the water level-storage capacity relationship based on reservoir siltation. This method can complete reservoir siltation calculation and water level-storage capacity relationship correction by fully utilizing initial reservoir capacity data, daily reservoir operation data, and water level data. In other words, it can quickly calculate the degree of reservoir siltation based on changes in daily operation data, thus eliminating the need for complex external topographic surveys or data measurements at specific times or in specific scenarios, making the calculation process very simple. Furthermore, since the calculated value is a proportional value based on the initial value, it fully conforms to the existing economic dispatch system of hydropower stations, thus requiring no changes to existing economic dispatch methods, and resulting in low cost for correcting the water level-storage capacity relationship.

[0077] The method for correcting the water level-storage capacity relationship based on reservoir siltation, as provided in the embodiments of this application, will be described in detail below.

[0078] It should be noted that, since small-capacity hydropower stations do not have, or do not fully have, daily regulation capacity, the method for correcting the water level-capacity relationship based on reservoir siltation provided in this application is mainly applied to the calculation of reservoir siltation and the correction of the water level-capacity relationship for small-capacity hydropower stations. Of course, it is understood that the method in this application can also be applied to the calculation of reservoir siltation and the correction of the water level-capacity relationship for hydropower stations with conventional or large reservoir capacities; this is not limited here. The method in this application can be executed on any terminal device with any operating system. Terminal devices include, but are not limited to, industrial computers, desktop computers, smartphones, and mobile iPads; this is not limited here. For ease of description, unless otherwise specified, this application uses an industrial computer as the execution subject in its embodiments. It is understood that the execution subject does not constitute a limitation on the embodiments of this application; other types of terminal devices may be used as the execution subject in some other embodiments.

[0079] like Figure 1The diagram shown is a flowchart of a method for correcting the water level-storage capacity relationship based on reservoir siltation, provided in an embodiment of this application. The method for correcting the water level-storage capacity relationship based on reservoir siltation includes, but is not limited to, steps S1 to S4:

[0080] Step S1. Obtain the initial reservoir capacity value V corresponding to each water level layer in the reservoir area. 0i Where i = 1, 2, ..., n, and n represents the number of water level layers;

[0081] It should be noted that the initial reservoir capacity of each water level layer in this application embodiment was calibrated during the construction of the hydropower station, that is, the initial reservoir capacity of each water level layer in this application embodiment is a preset fixed value.

[0082] Step S2. Obtain data on the operational status changes of each water level and flow rate in the reservoir over a period of time. ij Where j = 1, 2, ..., m, m represents the number of flow levels;

[0083] It should be noted that the aforementioned operating status change data A ij At least include the water level change time A1 ij Total station active power output value A2 ij and the opening degree A3 of all station gates. ij Here, water level change time refers to the time elapsed from one water level to another in the reservoir area; total station active power output refers to the total active power output of the entire station when the unit gate opening is 0; and total station gate opening refers to the total station gate opening when the water level remains at the current level, either when the unit is shut down for peak shaving or when the total station active power is 0. Of course, it is understood that the operational status change data in this embodiment is not limited to the examples above, and can also be other operational status change data, such as power generation, the number of times the total station gates actuate, etc., which are not specifically limited here.

[0084] For example: Water level change time A1 ij It could be: the time T during which the water level changes from 553m to 554m when the flow rate is 3000m³ / s. 33 Or, when the water level is 552m and the flow rate is 2000m³ / s, the time T for the water level to change to 554m. 22 The total station's active power output is A2. ij This can be calculated when the water level is 553m and the flow rate is 2000m³ / s, showing the total active power output P of the entire station. 32 Or, when the water level is 553m and the flow rate is 3000m³ / s, the total active power output P of the station is... 32 The opening degree of all gates at the station is A3. ij This can be defined as follows: when the water level is 553m, the flow rate is 4000m³ / s, the total active power output of the station is 0, and the water level is maintained at 553m, the opening degree G of all the station's gates is considered.34 .

[0085] Step S3. Based on the running status change data A ij Compared with the initial running status data A ij0 The proportional relationship between them is used to calculate the relative siltation parameter P of the reservoir area;

[0086] In one specific implementation of step S3, based on the operating state change data A ij Compared with the initial running status data A ij0 The proportional relationship between them is used to calculate the relative siltation parameter P of the reservoir area, including:

[0087] Step S31. Based on the change data A of each operating state ij With the corresponding initial running status data A ij0 The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor K. g The calculation formula is as follows:

[0088]

[0089] Where g = 1, 2, ..., q, q represents the number of types of running status change data, and Δt represents the time period for acquiring running status change data;

[0090] Specifically, when the operating status changes data A ij At least include the water level change time A1 ij Total station active power output value A2 ij and the opening degree A3 of all station gates. ij At that time, the data A based on each operating state change ij With the corresponding initial running status data A ij0 The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor K. g ,include:

[0091] According to the water level change time A1 ij Time A1 of initial water level change ij0 The proportional relationship between them is based on the total active power output of the station, A2. ij Compared with the initial total active power output value A2 ij0 The proportional relationship, and based on the opening degree A3 of the entire station gate. ij With the initial total gate opening A3 ij Based on the proportional relationship, the first sedimentation parameter factor K1, the second sedimentation parameter factor K2, and the third sedimentation parameter factor K3 are calculated respectively, using the following specific calculation formulas:

[0092]

[0093] If there is siltation in the reservoir area, then A1 ij The value will decrease, A2 ij The value will decrease, A3 ij The value will increase, i.e., A1 ij -A1 ij0 <0, A2 ij -A2 ij0 <0,A3 ij -A3 ij0 If the value is greater than 0, then the trend of siltation in the reservoir area Y = |T| can be determined based on the changes in water level over time, the total active power output of the station, and the absolute values ​​of the gate openings relative to their initial values. mn -T mn0 |+P mn -P mn0 |+G mn -G mn0 The larger the Y value, the more severe the reservoir siltation (the initial value of Y is 0). The siltation parameter factor values ​​K1, K2, and K3 can be determined by the magnitude of the change.

[0094] Step S32. Calculate the sedimentation parameter factor K for each sedimentation parameter. g The corresponding parameter weight M g The calculation formula is as follows:

[0095]

[0096] in,

[0097] Specifically, the parameter weight M1 corresponding to the first sedimentation parameter factor K1 for water level change time, the parameter weight M2 corresponding to the second sedimentation parameter factor K2 for the total station active power output, and the total station gate opening A3 are calculated respectively. ij The third sedimentation parameter M2 factor K3 corresponds to the parameter weight M3, and M1+M2+M3=1.

[0098] Step S33. Based on each sedimentation parameter factor K g and the weight M of each parameter g The relative siltation parameter P of the reservoir area is calculated using the following formula:

[0099]

[0100] Where T represents the transpose of the matrix.

[0101] In one specific implementation, after calculating the relative siltation parameter P of the reservoir area, the method further includes:

[0102] Calculate the relative sedimentation parameter P for different reservoir areas during L time periods. l l = 1, 2, ..., L;

[0103] Based on the relative siltation parameter P of multiple reservoir areas l The average siltation parameter P in the reservoir area was calculated. aver The calculation formula is as follows:

[0104]

[0105] Based on the above, this embodiment calculates the relative siltation parameter P of the reservoir area for different time periods within L time intervals. l The average siltation parameter P in the reservoir area was calculated. aver This yields the geometric mean of the relative siltation parameters in the reservoir area, which is then used for subsequent siltation capacity calculations, further improving the accuracy of the results. Preferably, in this embodiment, each time layer of L can be randomly selected to ensure the authenticity of the time layer acquisition.

[0106] Step S4. Based on the initial storage capacity value V 0i The siltation capacity V is calculated from the relative siltation parameter P of the reservoir area, and the actual reservoir capacity V' of each water level layer is then calculated using the siltation capacity V. He Make corrections.

[0107] In one specific implementation of step S4, based on the initial storage capacity value V 0i The siltation capacity V is calculated from the relative siltation parameter P of the reservoir area, including:

[0108] Step S41. Based on the initial storage capacity value V 0i Based on the relative siltation parameter P of the reservoir area, the reservoir capacity ΔV after siltation is calculated using the following formula:

[0109] ΔV=PV 0i (5)

[0110] It should be noted that the reservoir capacity value ΔV after siltation refers to the remaining reservoir capacity after the original capacity has been affected by siltation.

[0111] Step S42. Calculate the silted-up reservoir capacity V based on the reservoir capacity ΔV after siltation. The calculation formula is as follows:

[0112] V = (1-P)V 0i (6).

[0113] It should be noted that the siltation capacity V refers to the volume of the siltation itself.

[0114] In one specific implementation, the actual reservoir capacity V' of each water level layer is calculated using the siltation reservoir capacity V. He Make corrections, including:

[0115] Using the siltation capacity V, the actual reservoir capacity V' of each water level layer was calculated. He The correction is made, and the correction formula is as follows:

[0116]

[0117] Where e represents the e-th water level in layer E, and E represents the number of water level layers in the reservoir area, from the dead water level to the upper limit water level, based on a preset unit step size; H e This represents the calibration value of the water level at layer e. This represents the original reservoir capacity at the e-th water level.

[0118] In one specific implementation, the actual reservoir capacity V' of each water level layer is calculated using the siltation reservoir capacity V. He After modification, the method further includes:

[0119] For each water level layer, the actual reservoir capacity V' H e Verification is required.

[0120] Specifically, the actual reservoir capacity value of each water level layer was verified using the following methods:

[0121]

[0122] Preferably, the verification is performed e-2 times, and e≥3.

[0123] Based on the above-disclosed content, the embodiments of this application can complete the calculation of reservoir siltation and the correction of the water level-reservoir capacity relationship by making full use of initial reservoir capacity data, daily operation data of the reservoir area, and water level data. That is, the degree of siltation in the reservoir area can be quickly calculated based on changes in daily operation data, so there is no need to carry out complex external topographic surveys or data measurements at specific times or in specific scenarios. The calculation process is very simple. At the same time, since the calculated value is a proportional value based on the initial value, it is fully in line with the existing economic dispatch system of hydropower stations, so there is no need to change the existing economic dispatch method, and the cost of correcting the water level-reservoir capacity relationship is low.

[0124] The second aspect provides a device for correcting the water level-reservoir capacity relationship based on reservoir siltation, including:

[0125] The initial reservoir capacity acquisition module is used to obtain the initial reservoir capacity value V corresponding to each water level layer in the reservoir area. 0i Where i = 1, 2, ..., n, and n represents the number of water level layers;

[0126] The data acquisition module is used to acquire data on the changes in the operating status of each water level and flow rate in the reservoir over a period of time. ij Where j = 1, 2, ..., m, m represents the number of flow levels;

[0127] The siltation parameter calculation module is used to calculate data A based on changes in operating status. ij Compared with the initial running status data A ij0 The proportional relationship between them is used to calculate the relative siltation parameter P of the reservoir area;

[0128] The storage capacity correction module is used to adjust the storage capacity based on the initial storage capacity value V. 0i The siltation capacity V is calculated from the relative siltation parameter P of the reservoir area, and the actual reservoir capacity V' of each water level layer is then calculated using the siltation capacity V. He Make corrections.

[0129] In one possible design, the operating state change data A ij At least include the water level change time A1 ij Total station active power output value A2 ij and the opening degree A3 of all station gates. ij .

[0130] In one possible design, based on the changing data A of the operating state... ij Compared with the initial running status data A ij0 When calculating the relative siltation parameter P of the reservoir area based on the proportional relationship between the two, the data acquisition module is specifically used for the following purposes:

[0131] Based on the data A of each operating state change ij With the corresponding initial running status data A ij0 The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor K. g The calculation formula is as follows:

[0132]

[0133] Where g = 1, 2, ..., q, q represents the number of types of running status change data, and Δt represents the time period for acquiring running status change data;

[0134] Calculate the sedimentation parameter factor K for each sedimentation parameter separately. g The corresponding parameter weight M g The calculation formula is as follows:

[0135]

[0136] in,

[0137] Based on each sedimentation parameter factor K g and the weight M of each parameter g The relative siltation parameter P of the reservoir area is calculated using the following formula:

[0138]

[0139] Where T represents the transpose of the matrix.

[0140] In one possible design, the sedimentation parameter calculation module is also used for:

[0141] Calculate the relative sedimentation parameter P for different reservoir areas during L time periods. l l = 1, 2, ..., L;

[0142] Based on the relative siltation parameter P of multiple reservoir areas l The average siltation parameter P in the reservoir area was calculated. aver The calculation formula is as follows:

[0143]

[0144] In one possible design, based on the initial storage capacity value V 0i When calculating the siltation capacity V based on the relative siltation parameter P of the reservoir area, the reservoir capacity correction module is specifically used for:

[0145] Based on the initial storage capacity value V 0i Based on the relative siltation parameter P of the reservoir area, the reservoir capacity ΔV after siltation is calculated using the following formula:

[0146] ΔV=PV 0i (5)

[0147] Based on the reservoir capacity ΔV after siltation, the silted-up reservoir capacity V is calculated using the following formula:

[0148] V = (1-P)V 0i (6).

[0149] In one possible design, the actual reservoir capacity V' of each water level layer is calculated using the siltation reservoir capacity V. He When making corrections, the warehouse capacity correction module is specifically used for:

[0150] Using the siltation capacity V, the actual reservoir capacity V' of each water level layer was calculated. He The correction is made, and the correction formula is as follows:

[0151]

[0152] Where e represents the e-th water level in layer E, and E represents the number of water level layers in the reservoir area, from the dead water level to the upper limit water level, based on a preset unit step size; H e This represents the calibration value of the water level at layer e. This represents the original reservoir capacity at the e-th water level.

[0153] In one possible design, the device further includes:

[0154] The verification module is used to verify the actual reservoir capacity V' of each water level. H e Verification is required.

[0155] The working process, working details and technical effects of the aforementioned device provided in the second aspect of this embodiment can be found in the method described in the first aspect or any possible design of the first aspect, and will not be repeated here.

[0156] Thirdly, the present invention provides a computer device comprising a memory, a processor, and a transceiver connected in sequence and communication, wherein the memory is used to store a computer program, the transceiver is used to send and receive messages, and the processor is used to read the computer program and execute the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any possible design of the first aspect.

[0157] Specifically, the memory may include, but is not limited to, Random-Access Memory (RAM), Read-Only Memory (ROM), Flash Memory, First-In-First-Out (FIFO) Memory, and / or First-In-Last-Out (FILO) Memory, etc.; the processor may not be limited to the STM32F105 series microprocessor; the transceiver may be, but is not limited to, a WiFi (Wireless Fidelity) wireless transceiver, a Bluetooth wireless transceiver, a GPRS (General Packet Radio Service) wireless transceiver, and / or a ZigBee (a low-power LAN protocol based on the IEEE 802.15.4 standard) wireless transceiver, etc. Furthermore, the computer device may also include, but is not limited to, a power module, a display screen, and other necessary components.

[0158] The working process, working details and technical effects of the aforementioned computer device provided in the third aspect of this embodiment can be found in the method described in the first aspect or any possible design of the first aspect, and will not be repeated here.

[0159] Fourthly, the present invention provides a computer-readable storage medium storing instructions that, when executed on a computer, perform the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any possible design of the first aspect.

[0160] The computer-readable storage medium refers to a carrier for storing data, which may include, but is not limited to, floppy disks, optical disks, hard disks, flash memory, USB flash drives and / or memory sticks, etc. The computer may be a general-purpose computer, a special-purpose computer, a computer network, or other programmable devices.

[0161] The working process, working details and technical effects of the aforementioned computer-readable storage medium provided in the fourth aspect of this embodiment can be found in the method described in the first aspect or any possible design of the first aspect, and will not be repeated here.

[0162] Fifthly, the present invention provides a computer program product containing instructions that, when the instructions are executed on a computer, cause the computer to perform the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any possible design of the first aspect.

[0163] The working process, working details and technical effects of the aforementioned computer program product containing instructions provided in the fifth aspect of this embodiment can be found in the method described in the first aspect or any possible design of the first aspect, and will not be repeated here.

[0164] Finally, it should be noted that the above description is merely a preferred embodiment of the present invention and is not intended to limit the scope of protection of the present invention. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of the present invention should be included within the scope of protection of the present invention.

Claims

1. A method for correcting the water level-reservoir capacity relationship based on reservoir siltation, characterized in that, include: Obtain the initial reservoir capacity corresponding to each water level layer in the reservoir area. Where i = 1, 2, ..., n, and n represents the number of water level layers; Acquire operational status change data of each water level and flow rate in the reservoir area over a period of time. Where j = 1, 2, ..., m, m represents the number of flow levels; Based on operational status change data Compared with initial running status data The proportional relationship between them is used to calculate the relative siltation parameters of the reservoir area. ; Based on the initial storage capacity Relative siltation parameters of the reservoir area The siltation capacity was calculated. and utilize the accumulated reservoir capacity The actual reservoir capacity values ​​for each water level were corrected separately. The data on changes in operating status At least include the time of water level change Total station contribution value and the opening degree of all station gates ; Based on operational status change data Compared with initial running status data The proportional relationship between them is used to calculate the relative siltation parameters of the reservoir area. ,include: Based on the data of each operating state change With the corresponding initial running status data The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor. The calculation formula is as follows: ；（1） Where g = 1, 2, ..., q, q represents the number of data types related to changes in running state. Indicates the time period for acquiring data on changes in operating status; Calculate each sedimentation parameter factor separately Corresponding parameter weights The calculation formula is as follows: ；（2） in, ; Based on each sedimentation parameter factor and the weight of each parameter The relative siltation parameters of the reservoir area were calculated. The calculation formula is as follows: ；（3） in, Represents the transpose of a matrix; Based on the initial storage capacity Relative siltation parameters of the reservoir area The siltation capacity was calculated. ,include: Based on the initial storage capacity Relative siltation parameters of the reservoir area Calculate the reservoir capacity after siltation. The calculation formula is as follows: ；（5） Based on the reservoir capacity after siltation The siltation capacity was calculated. The calculation formula is as follows: （6）。 2. The method for correcting the water level-reservoir capacity relationship based on reservoir siltation according to claim 1, characterized in that, The relative siltation parameters of the reservoir area were calculated. Subsequently, the method further includes: Calculate the relative siltation parameters of different reservoir areas during L time periods. , ; Based on relative siltation parameters of multiple reservoir areas The average siltation parameters of the reservoir area were calculated. The calculation formula is as follows: （4）。 3. The method for correcting the water level-reservoir capacity relationship based on reservoir siltation according to claim 1, characterized in that, Utilizing silt storage capacity The actual reservoir capacity values ​​for each water level layer are corrected, including: Utilizing silt storage capacity The actual reservoir capacity values ​​at each water level were respectively... The correction is made, and the correction formula is as follows: ；（7） in, This represents the e-th water level in the E-level water level, where E represents the number of water level layers in the reservoir area, from the dead water level to the upper limit water level, based on a preset unit step size. This represents the calibration value of the water level at layer e. This represents the original reservoir capacity at the e-th water level.

4. The method for correcting the water level-reservoir capacity relationship based on reservoir siltation according to claim 1, characterized in that, Utilizing the accumulated reservoir capacity After correcting the actual reservoir capacity values ​​for each water level layer, the method further includes: The actual reservoir capacity value for each water level layer. Verification is required.

5. A device for correcting the water level-reservoir capacity relationship based on reservoir siltation, characterized in that, include: The initial reservoir capacity acquisition module is used to obtain the initial reservoir capacity value corresponding to each water level layer in the reservoir area. Where i = 1, 2, ..., n, and n represents the number of water level layers; The data acquisition module is used to acquire data on the changes in the operating status of each water level and flow rate in the reservoir over a period of time. Where j = 1, 2, ..., m, m represents the number of flow levels; The siltation parameter calculation module is used to calculate data based on changes in operating status. Compared with initial running status data The proportional relationship between them is used to calculate the relative siltation parameters of the reservoir area. ; The storage capacity correction module is used to adjust the storage capacity based on the initial storage capacity value. Relative siltation parameters of the reservoir area The siltation capacity was calculated. and utilize the accumulated reservoir capacity The actual reservoir capacity values ​​at each water level were respectively... Make corrections; The data on changes in operating status At least include the time of water level change Total station contribution value and the opening degree of all station gates ; Based on operational status change data Compared with initial running status data The proportional relationship between them is used to calculate the relative siltation parameters of the reservoir area. ,include: Based on the data of each operating state change With the corresponding initial running status data The proportional relationship between them is used to calculate the corresponding sedimentation parameter factor. The calculation formula is as follows: ；（1） Where g = 1, 2, ..., q, q represents the number of data types related to changes in running state. Indicates the time period for acquiring data on changes in operating status; Calculate each sedimentation parameter factor separately Corresponding parameter weights The calculation formula is as follows: ；（2） in, ; Based on each sedimentation parameter factor and the weight of each parameter The relative siltation parameters of the reservoir area were calculated. The calculation formula is as follows: ；（3） in, Represents the transpose of a matrix; Based on the initial storage capacity Relative siltation parameters of the reservoir area The siltation capacity was calculated. ,include: Based on the initial storage capacity Relative siltation parameters of the reservoir area Calculate the reservoir capacity after siltation. The calculation formula is as follows: ；（5） Based on the reservoir capacity after siltation The siltation capacity was calculated. The calculation formula is as follows: （6）。 6. A storage medium, characterized in that, The storage medium stores instructions that, when executed on a computer, perform the method for correcting the water level-reservoir capacity relationship based on reservoir siltation as described in any one of claims 1 to 4.

Images