River type reservoir main stream reservoir inflow correction method, device, equipment and medium

Abstract

The invention relates to the technical field of reservoir dispatching, and discloses a riverway type reservoir main stream reservoir inflow correction method, device, equipment and medium, and the riverway type reservoir main stream reservoir inflow correction method comprises the following steps: constructing a one-dimensional hydrodynamic model based on actually measured hydrological data and topographic data of a riverway type reservoir; the initial roughness value is adjusted, the calculation error of the one-dimensional hydrodynamic model is reduced, and a first roughness value is obtained; the calculation time period is prolonged, and based on the first roughness value, a first calculation water level is obtained through calculation by means of a one-dimensional hydrodynamic model; the first roughness value is adjusted again, the calculation error of the one-dimensional hydrodynamic model is reduced, and a second roughness value is obtained; based on the second roughness value, a second calculated water level is obtained through calculation by means of the one-dimensional hydrodynamic model; determining a difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level and the actually measured water level; and correcting the main stream storage flow in the difference time period. The calculation precision of the main stream storage flow can be improved.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of reservoir operation, in particular to a river type reservoir main stream inflow correction method, device, equipment and medium. BACKGROUND

[0002] The inflow data of a reservoir is basic data for reservoir operation and management, hydrological forecasting and reservoir scheduling.

[0003] At present, the most commonly used inflow calculation method is the water balance method, which inversely calculates the inflow according to the reservoir water level and outflow. However, the current method does not consider the influence of the inflow propagation time. When the time interval is small, or the inflow process changes greatly, such as heavy floods or heavy rains, the current method will have a time and space mismatch in the correction results.

[0004] For a river type reservoir, it has the dual characteristics of river type reservoir and lake type reservoir. The relationship between reservoir water level and reservoir storage cannot be simply represented by the water level-storage curve. The traditional water balance method cannot be applied to river type reservoirs, and cannot provide an accurate inflow calculation method. SUMMARY

[0005] Therefore, the present application provides a river type reservoir main stream inflow correction method, device, equipment and medium to improve the calculation accuracy of the main stream inflow and make full use of the measured data to meet the business needs of fine scientific scheduling.

[0006] In a first aspect, the present application provides a river type reservoir main stream inflow correction method, which comprises: based on the measured hydrological data and topographic data of the river type reservoir, a one-dimensional hydrodynamic model of the river type reservoir is constructed, the measured hydrological data includes measured flow data and measured water level data, the measured flow data includes main stream inflow, branch stream inflow and along-path monitoring flow, and the measured water level data includes dam front water level and along-path monitoring water level; adjusting the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model to obtain the first roughness value of each river section; extending the calculation time period, and based on the first roughness value, the first calculation water level is calculated by using the one-dimensional hydrodynamic model; adjusting the first roughness value of each river section again to reduce the calculation error of the one-dimensional hydrodynamic model to obtain the second roughness value of each river section; based on the second roughness value, the second calculation water level is calculated by using the one-dimensional hydrodynamic model; the difference time period is determined in combination with the first roughness value, the second roughness value, the first calculation water level, the second calculation water level and the measured water level; and the main stream inflow in the difference time period is corrected.

[0007] In this implementation, in order to solve the problem that it is difficult to accurately obtain the main inflow reservoir flow under large flow fluctuations, a one-dimensional hydrodynamic model is constructed. The deviation of the main inflow reservoir flow is analyzed by adjusting the roughness and measurement time, and the main inflow reservoir flow is corrected based on the roughness. This fills the technical deficiency of judging the main inflow reservoir flow by the water level-flow relationship curve, solves the problem of inaccurate judgment of the main inflow reservoir flow when the water level in front of the dam or the main inflow reservoir flow changes significantly, and provides important technical support for the scientific scheduling and safe operation of river-type reservoirs.

[0008] In an optional embodiment, the initial roughness value of each river section is adjusted to reduce the calculation error of the one-dimensional hydrodynamic model, and a first roughness value of each river section is obtained; the calculation time period is extended, and based on the first roughness value, a first calculated water level is calculated using the one-dimensional hydrodynamic model, including: for a first preset time period, the initial roughness value of each river section is adjusted until the error between the calculated water level calculated by the one-dimensional hydrodynamic model of the channel-type reservoir and the first measured water level is less than an error threshold, and the first roughness value of each river section is obtained, and the first measured water level is the measured water level of the first preset time period; for a second preset time period, based on the first roughness value of each river section, a first calculated water level is calculated using the one-dimensional hydrodynamic model of the channel-type reservoir; the second preset time period includes the first preset time period and its preset extended time period.

[0009] In an optional embodiment, the first roughness value of each river section is adjusted again to reduce the calculation error of the one-dimensional hydrodynamic model, so as to obtain a second roughness value of each river section; based on the second roughness value, a second calculated water level is calculated using the one-dimensional hydrodynamic model, including: adjusting the first roughness value of each river section for a second preset time period until the error between the calculated water level calculated by the one-dimensional hydrodynamic model of the river-type reservoir and the second measured water level is less than an error threshold, so as to obtain the second roughness value of each river section; the second measured water level is the measured water level in the second preset time period; based on the second roughness value of each river section, a second calculated water level is calculated using the one-dimensional hydrodynamic model.

[0010] In this implementation, the first calculated water level is determined by extending the measurement time on the basis of determining the first roughness value in the first preset time period; and the second calculated water level is determined on the basis of determining the second roughness value in the second preset time period. The measurement time is fully considered, and the influence of the measurement time on the calculation of the roughness value and the influence of the change in the roughness value on the water level calculation can be analyzed, which facilitates the subsequent analysis of abnormal dry inflow reservoir flow.

[0011] In an optional embodiment, adjusting the initial roughness value of each river section includes: adjusting the initial roughness value of each river section section by section according to a first calibration order, the first calibration order being the order from downstream to upstream; adjusting the first roughness value of each river section includes: adjusting the first roughness value of each river section section by section according to a second calibration order, the second calibration order being the order from upstream to downstream.

[0012] In this implementation, different calibration orders are set for the two roughness value calibrations, which can improve the accuracy of roughness value determination.

[0013] In an optional embodiment, the initial roughness value of each river section is adjusted to reduce the calculation error of the one-dimensional hydrodynamic model, and before obtaining the first roughness value of each river section, the method includes: obtaining the along-line observation station information of the river-type reservoir; dividing the river channel of the river-type reservoir into multiple river sections based on the along-line observation station information; and assigning fixed values ​​to the initial roughness values ​​of the multiple river sections.

[0014] In this implementation method, the river section division is based on the observation stations along the river and is divided according to the hydraulic conditions of the river section, which can improve the rationality of the river section division and provide a basis for the subsequent roughness value calibration.

[0015] In an optional embodiment, the difference time period is determined in combination with the first roughness value, the second roughness value, the first calculated water level, the second calculated water level and the second measured water level, including: determining the difference river section based on the difference between the first roughness value and the second roughness value; for the difference river section, determining the difference time period within a second preset time period based on the difference between the first calculated water level, the second calculated water level and the second measured water level.

[0016] In this implementation, the impact of measurement time on roughness value calculation is analyzed to determine abnormal river sections, the impact of roughness value changes on water level calculation is analyzed, and abnormal time periods are determined to facilitate subsequent accurate correction of main inflow reservoir flows.

[0017] In an optional embodiment, the main inflow reservoir flow rate within the difference time period is corrected, including: calculating the roughness change rate from the first roughness value to the second roughness value; correcting the main inflow reservoir flow rate based on the flow correction coefficient, and the sum of the flow correction coefficient and the roughness change rate is one.

[0018] In this implementation, when the roughness changes, the flow rate will change accordingly if the water level does not change. Based on the change in roughness, the change in the flow rate into the main inflow reservoir is determined and corrected, and the correction calculation efficiency is high.

[0019] In the second aspect, the present invention provides a device for correcting the flow of the main inflow into the reservoir of a river-type reservoir, and the device for correcting the flow of the main inflow into the reservoir of a river-type reservoir comprises: a construction module for constructing a one-dimensional hydrodynamic model of the river-type reservoir based on the measured hydrological data and terrain data of the river-type reservoir, the measured hydrological data including measured flow data and measured water level data, the measured flow data including the main inflow into the reservoir flow, the tributary inflow into the reservoir flow and the flow monitored along the way, and the measured water level data including the water level in front of the dam and the water level monitored along the way; a first calculation module for adjusting the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model, and obtain to the first roughness value of each river section; extend the calculation time period, and obtain the first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model; the second calculation module is used to adjust the first roughness value of each river section again, reduce the calculation error of the one-dimensional hydrodynamic model, and obtain the second roughness value of each river section; based on the second roughness value, obtain the second calculated water level using the one-dimensional hydrodynamic model; the determination module is used to determine the difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level and the measured water level; the correction module is used to correct the main inflow reservoir flow within the difference time period.

[0020] In a third aspect, the present invention provides a computer device comprising: a memory and a processor, the memory and the processor being communicatively connected to each other, the memory storing computer instructions, and the processor executing the method for correcting the flow rate of a main inflow into a river-type reservoir according to the first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0021] In a fourth aspect, the present invention provides a computer-readable storage medium having computer instructions stored thereon, the computer instructions being used to enable a computer to execute the method for correcting the flow rate of a main inflow into a river-type reservoir according to the first aspect or any corresponding embodiment thereof.

[0022] In a fifth aspect, the present invention provides a computer program product comprising computer instructions for causing a computer to execute the method for correcting the flow rate of a main inflow into a river-type reservoir according to the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

[0023] In order to more clearly illustrate the specific embodiments of the present invention or the technical solutions in the prior art, the following briefly introduces the drawings required for use in the specific embodiments or the description of the prior art. Obviously, the drawings described below are some embodiments of the present invention. For ordinary technicians in this field, other drawings can be obtained based on these drawings without paying any creative work.

[0024] Figure 1is a flow chart of a river type reservoir main stream inflow correction method according to an embodiment of the present application;

[0025] Figure 2 is a flow chart of another river type reservoir main stream inflow correction method according to an embodiment of the present application;

[0026] Figure 3 is a water level change schematic diagram according to an embodiment of the present application;

[0027] Figure 4 is a structural block diagram of a river type reservoir main stream inflow correction device according to an embodiment of the present application;

[0028] Figure 5 is a hardware structure schematic diagram of a computer device according to an embodiment of the present application. DETAILED DESCRIPTION

[0029] In order to make the objects, technical solutions and advantages of embodiments of the present application clearer, the technical solutions in the embodiments of the present application will be described clearly and completely below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are some but not all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative work fall within the protection scope of the present application.

[0030] At present, the error in the existing inflow calculation method mainly comes from the magnification effect of reservoir water level record error. When the reservoir capacity is large and the backstepping period is short, even a small reservoir water level record error can also cause a large inflow calculation error, so as to cause the inflow process to appear jagged fluctuations, and even negative values and other unreasonable phenomena. In order to overcome the jagged fluctuation phenomenon of the inflow, the traditional method is to select a longer backstepping period to reduce the magnification effect of the reservoir water level record error, but the period is too long to homogenize the flood process, which causes the homogenization and flat head phenomenon of the flood peak, and it is difficult to reflect the flood peak information accurately. In addition, according to the relationship among the storage error, the inflow allowable error and the optimized backstepping period, a backstepping period dynamic adjustment method is constructed based on the flow test error specification and the reservoir water level record error, so as to automatically calculate the backstepping period and the calculated inflow at the corresponding forecast time which meets the real-time error control requirements.

[0031] As for river-type reservoirs, since their length is often longer, their inflow is generally composed of main inflow flow, tributary inflow flow, and uncontrolled interval flow. For the main inflow flow, the water level is generally obtained by the hydrological station through a water gauge, and then it is calculated with the help of the water level-flow relationship. When the water level in front of the dam or the inflow flow changes relatively slowly, the main inflow flow is generally relatively reliable; when the water level in front of the dam or the inflow flow changes greatly, there will be a certain error in the main inflow flow. However, for reservoir scheduling business, the more the water level in front of the dam or the inflow flow changes greatly, the more accurate main inflow flow data support is needed. Therefore, this application proposes a method for correcting the main inflow flow of river-type reservoirs to correct the main inflow flow under large fluctuations, improve the flow accuracy, and meet the business needs of refined and scientific scheduling.

[0032] According to an embodiment of the present invention, an embodiment of a method for correcting the flow of the main inflow into a river-type reservoir is provided. It should be noted that the steps shown in the flowchart of the accompanying drawings can be executed in a computer system such as a set of computer-executable instructions, and although a logical order is shown in the flowchart, in some cases, the steps shown or described can be executed in an order different from that shown here.

[0033] In this embodiment, a method for correcting the flow rate of the main inflow of a river-type reservoir is provided. Figure 1 This is a flow chart of a method for correcting the flow rate of a main inflow of a river-type reservoir according to an embodiment of the present invention. It should be noted that if there are substantially the same results, this embodiment does not use Figure 1 The process sequence shown is limited. Figure 1 As shown, the process includes the following steps:

[0034] Step S101 : constructing a one-dimensional hydrodynamic model of the river-type reservoir based on the measured hydrological data and topographic data of the river-type reservoir.

[0035] For a river-type reservoir, in order to ensure the accuracy of the constructed one-dimensional hydrodynamic model, corresponding reliable data are collected, including measured hydrological data and topographic data.

[0036] Specifically, the measured hydrological data include measured flow data and measured water level data. The measured flow data are the actually measured flow of the main stream and the tributary into the reservoir, and the flow monitored along the way by the hydrological stations along the way. The measured water level data include the actually measured water level in front of the dam, and the water level monitored along the way by the hydrological stations or water level stations along the way.

[0037] It can be understood that the main inflow reservoir flow and tributary inflow reservoir flow in this application are flow data that can be directly measured, not inflow flow data. In addition to the main inflow reservoir flow and tributary inflow reservoir flow, the inflow flow data also includes uncontrolled interval flow data, and the uncontrolled interval flow data cannot be directly measured. Therefore, in order to ensure the accuracy of model construction, this application does not consider uncontrolled interval flow data.

[0038] In other implementations, the measured hydrological data also includes rainfall data monitored by meteorological stations, etc.

[0039] Specifically, the topographic data includes the topographic cross-section information of the found river-type reservoir. It can be understood that the topographic data does not include underwater topographic information.

[0040] In other implementations, the terrain data also includes information such as the dam site location, inlet and outlet locations, etc.

[0041] Furthermore, a one-dimensional hydrodynamic model of the river-type reservoir is constructed by combining the above-mentioned measured hydrological data and topographic data.

[0042] Specifically, the one-dimensional hydrodynamic model is based on the Saint-Venant equations, which are the fundamental equations that describe the motion of water in rivers, channels, or pipes. The Saint-Venant equations consist of the continuity equation and the momentum equation.

[0043] Step S102: Adjust the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model and obtain the first roughness value of each river section; extend the calculation time period, and calculate the first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model.

[0044] Among them, roughness is a dimensionless number that measures the effect of the river boundary on the flow resistance of water. It can reflect the friction and resistance of water when flowing in the river. The greater the roughness, the stronger the flow capacity of the river. The roughness of a river-type reservoir is not static, but will change with the changes in the flow rate of the main inflow into the reservoir and the water level in front of the dam. Its change pattern is very complex. If the roughness value is inaccurate, the water levels of each section calculated by the one-dimensional hydrodynamic model cannot be consistent with the measured along-the-line monitoring water level. At this time, the main inflow into the reservoir calculated by the model cannot reflect the actual situation. Therefore, this application reduces the calculation error of the one-dimensional hydrodynamic model by finding a suitable roughness of the river-type reservoir.

[0045] Specifically, the river channel is divided into multiple sections according to the differences in river channels of river-type reservoirs.

[0046] In one implementation, a river is divided into multiple sections based on changes in hydraulic characteristics and river topography. Specifically, when the river topography changes, such as narrowing, widening, deepening, or shallowing, the river is divided into two sections based on the changed location.

[0047] In another implementation, the river is divided into multiple sections based on the distribution of observation stations along the river-type reservoir. Specifically, the river is divided into two sections based on the locations of the observation stations along the river.

[0048] An initial roughness value is set for each river section. Based on the initial roughness value, the water level along the reservoir in the time period to be corrected is calculated using a one-dimensional hydrodynamic model, and the calculated water level along the reservoir is compared with the measured water level in the time period to be corrected. When the initial roughness value is inaccurate, the water level comparison error is large. By continuously adjusting the initial roughness value, the water level comparison error is reduced, so that the calculated water level result is basically consistent with the measured value, and the corresponding initial roughness value is used as the first roughness value.

[0049] Furthermore, the time period to be corrected is extended, and based on the first roughness value, the water level along the reservoir in the extended time period is calculated using a one-dimensional hydrodynamic model as the first calculated water level.

[0050] Step S103 , adjusting the first roughness value of each river section again to reduce the calculation error of the one-dimensional hydrodynamic model, and obtaining the second roughness value of each river section; and calculating the second calculated water level using the one-dimensional hydrodynamic model based on the second roughness value.

[0051] Since the first roughness value is calibrated during the time period to be corrected, the first calculated water level is compared with the measured water level during the extended time period. When the first roughness value is inaccurate, the roughness value needs to be adjusted again for the extended time period.

[0052] Specifically, based on the first roughness value, a one-dimensional hydrodynamic model is used to calculate the water level along the reservoir during the extended time period, and the calculated water level along the reservoir is compared with the measured water level during the extended time period. By continuously adjusting the initial roughness value, the water level comparison error is reduced so that the calculated water level result is basically consistent with the measured value, and the corresponding first roughness value is used as the second roughness value.

[0053] Furthermore, based on the second roughness value, a one-dimensional hydrodynamic model is used to calculate the water level along the reservoir during the extended time period as the second calculated water level.

[0054] Step S104 : determining a difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level, and the measured water level.

[0055] During the extended time period, find the time period when the difference between the first roughness value and the second roughness value is large, and the difference between the first calculated water level and / or the second calculated water level and the measured water level is large, and determine this time period as the difference time period, that is, the main inflow reservoir flow measured during the difference time period has an error and needs to be corrected.

[0056] Step S105, correcting the main inflow reservoir flow rate within the difference time period.

[0057] The method for correcting the main inflow flow of a channel-type reservoir provided in this embodiment addresses the problem that it is difficult to accurately obtain the main inflow flow under conditions of large flow fluctuations. A one-dimensional hydrodynamic model is constructed, and the deviation of the main inflow flow is analyzed by adjusting the roughness and measurement time. The main inflow flow is corrected based on the roughness. This fills the gap in the existing technology for determining the main inflow flow through the water level-flow relationship curve, solves the problem of inaccurate determination of the main inflow flow when the water level in front of the dam or the main inflow flow changes significantly, and provides important technical support for the scientific scheduling and safe operation of channel-type reservoirs.

[0058] In this embodiment, another method for correcting the flow rate of the main inflow of a river-type reservoir is provided. Figure 2 is a flow chart of another method for correcting the flow rate of a main inflow of a river-type reservoir according to an embodiment of the present invention. It should be noted that if there are substantially the same results, this embodiment does not use Figure 2 The process sequence shown is limited. Figure 2 As shown, the process includes the following steps:

[0059] Step S201 : constructing a one-dimensional hydrodynamic model of the river-type reservoir based on the measured hydrological data and topographic data of the river-type reservoir.

[0060] For details, please see Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0061] Step S202: Adjust the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model and obtain the first roughness value of each river section; extend the calculation time period, and calculate the first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model.

[0062] Generally, observation stations along a river-type reservoir are set up based on the river topography. In this implementation, the river is divided into multiple sections based on the distribution of observation stations along the river-type reservoir.

[0063] Specifically, the information of the observation stations along the river-type reservoir is obtained, including the location information; based on the information of the observation stations along the river, the river channel of the river-type reservoir is divided into multiple river sections from upstream to downstream with the observation stations along the river as the division nodes, and the initial roughness values ​​of the multiple river sections are assigned fixed values.

[0064] Taking the observation stations along the river as the basis for river section division, and dividing according to the river hydraulic conditions, the rationality of river section division can be improved, and a basis for subsequent roughness value calibration is provided.

[0065] Exemplarily, refer to Table 1, which is a specific river channel and its roughness statistics table. As shown in Table 1, there are 14 water level observation stations along the river from the tail of the reservoir to the front of the dam (from upstream to downstream), and the river roughness can be divided into 14 sections A-N. The range of each river section is shown in Table 1, and the initial roughness value of each river section is set to 20.

[0066] Table 1: River channel and its roughness statistics table

[0067]

[0068]

[0069] Further, the above step S202 includes:

[0070] Step S2021, for the first preset time period, adjusting the initial roughness value of each river section until the calculation water level calculated by the one-dimensional hydrodynamic model of the river-type reservoir and the first measured water level have an error less than an error threshold, and obtaining the first roughness value of each river section.

[0071] The first measured water level is a measured water level in the first preset time period.

[0072] In an implementation manner, the initial roughness value of each river section is adjusted in sections according to a first calibration order, and the first calibration order is an order from downstream to upstream.

[0073] Specifically, for the first preset time period to be corrected, based on the initial roughness value of each river section, the one-dimensional hydrodynamic model is used to calculate the calculation water level along the river of the reservoir, and the calculation water level along the river of the reservoir is compared with the first measured water level. When the initial roughness value is inaccurate, the water level comparison error is larger, the initial roughness value corresponding to the most downstream river section is adjusted first, until the calculation water level along the river of the reservoir corresponding to the river section and the first measured water level have an error less than an error threshold, and then the initial roughness value of the next river section is adjusted. According to the order from downstream to upstream, the initial roughness value of each river section is adjusted in sections until the calculation water level calculated by the one-dimensional hydrodynamic model of the river-type reservoir and the first measured water level have an error less than an error threshold, and the adjusted initial roughness value of each river section is taken as the first roughness value.

[0074] Exemplarily, the first preset time period of the dry flow inflow to be corrected is from September 14 to September 30, and the roughness values of the segments are sequentially adjusted from the dam to the upstream, so that the water level change process calculated by the one-dimensional water dynamic model is basically consistent with the measured value, and the first roughness value is denoted as K1. For details, refer to Table 1.

[0075] In step S2022, the first calculated water level is calculated by using the one-dimensional water dynamic model of the river type reservoir based on the first roughness value of each river segment for the second preset time period.

[0076] The second preset time period includes the first preset time period and a preset extended time period of the first preset time period.

[0077] In an implementation manner, the second preset time period is a time period twice as long as the first preset time period.

[0078] Exemplarily, the second preset time period of the dry flow inflow to be corrected is from September 14 to October 30, and the first calculated water level is calculated by using the one-dimensional water dynamic model based on the first roughness value K1, and is denoted as L1.

[0079] In step S203, the first roughness value of each river segment is adjusted again to reduce the calculation error of the one-dimensional water dynamic model, and the second roughness value of each river segment is obtained, and the second calculated water level is calculated by using the one-dimensional water dynamic model based on the second roughness value.

[0080] Specifically, the above step S203 includes:

[0081] In step S2031, the first roughness value of each river segment is adjusted for the second preset time period until the calculation water level calculated by the one-dimensional water dynamic model of the river type reservoir is less than the error threshold value from the second measured water level, and the second roughness value of each river segment is obtained.

[0082] The second measured water level is the measured water level of the second preset time period.

[0083] In an implementation manner, the first roughness value of each river segment is adjusted in a segment-by-segment manner according to a second rating sequence, and the second rating sequence is a sequence from the upstream to the downstream.

[0084] Specifically, for the second preset time period, the first roughness value of each river section is used to calculate the reservoir water level along the way by using a one-dimensional hydrodynamic model, and the reservoir water level along the way is compared with the second measured water level. When the first roughness value is inaccurate, the water level comparison error is larger. The first roughness value of the uppermost river section is adjusted first until the error between the reservoir water level along the way of the river section and the second measured water level is less than the error threshold. Then the first roughness value of the next river section is adjusted. The initial roughness value of each river section is adjusted in sequence from upstream to downstream until the error between the calculated water level calculated by the one-dimensional hydrodynamic model of the river type reservoir and the second measured water level is less than the error threshold. The adjusted first roughness value of each river section is taken as the second roughness value.

[0085] For example, for the second preset time period of the main stream inflow into the reservoir, i.e. from September 14 to October 30, the roughness value of each section is adjusted in sequence from the reservoir tail to the downstream, so that the one-dimensional hydrodynamic model can basically match the measured value in the process of calculating the reservoir water level along the way. The second roughness value is denoted as K2. For details, please refer to Table 1.

[0086] In step S2032, the second calculated water level is calculated by using a one-dimensional hydrodynamic model based on the second roughness value of each river section.

[0087] For example, the second preset time period of the main stream inflow into the reservoir is from September 14 to October 30. The second calculated water level calculated by using a one-dimensional hydrodynamic model based on the second roughness value K2 is denoted as L2.

[0088] In this implementation, the first calculated water level is determined by extending the measurement time based on the first roughness value determined in the first preset time period. The second calculated water level is determined based on the second roughness value determined in the second preset time period. The influence of the measurement time on the calculation of the roughness value and the influence of the change of the roughness value on the calculation of the water level are fully considered, which facilitates subsequent analysis of abnormal main stream inflow into the reservoir.

[0089] In step S204, the difference time period is determined in combination with the first roughness value, the second roughness value, the first calculated water level, the second calculated water level and the measured water level.

[0090] Specifically, the above step S204 includes:

[0091] In step S2041, the difference river section is determined based on the difference between the first roughness value and the second roughness value.

[0092] The first roughness value and the second roughness value of each river section are compared. When the difference between the first roughness value and the second roughness value is large, it indicates that there is an error in the water level and flow monitoring result corresponding to the river section.

[0093] Specifically, the roughness difference values ​​of multiple river sections are calculated, and the river section with the largest difference value is selected as the difference river section.

[0094] For example, as shown in Table 1, the first roughness values ​​K1 and the second roughness values ​​K2 of 14 river sections are compared, and the river section B with the largest roughness difference is selected as the difference river section.

[0095] Step S2042: for the differential river section, based on the difference among the first calculated water level, the second calculated water level and the second measured water level, a differential time period is determined within a second preset time period.

[0096] For different river sections, compare the differences between the corresponding calculated water levels and the measured water levels.

[0097] Specifically, the difference between the first calculated water level and the second measured water level is compared, and the difference between the second calculated water level and the second measured water level is compared, and the time period in which the difference between the first calculated water level and the second measured water level is large and / or the difference between the second calculated water level and the second measured water level is large is taken as the difference time period.

[0098] For example, see Figure 3 , Figure 3 2 is a schematic diagram of water level changes according to an embodiment of the present invention. Figure 3 The middle shows the water level changes of the first calculated water level L1, the second calculated water level L2 and the measured water level in the river section B during the second preset time period from September 14 to October 30. Figure 3 As shown, during the revised period from September 14th to September 30th, the first calculated water level L1 agrees well with the measured water level, while the second calculated water level L2 deviates from the measured water level. During the extended period from October 1st to October 30th, both the first and second calculated water levels L1 and L2 agree well with the measured water levels. Therefore, it can be determined that the main inflow flow from September 14th to September 30th deviates.

[0099] In this implementation, the impact of measurement time on roughness value calculation is analyzed to determine abnormal river sections, the impact of roughness value changes on water level calculation is analyzed, and abnormal time periods are determined to facilitate subsequent accurate correction of main inflow reservoir flows.

[0100] Step S205, correcting the main inflow reservoir flow rate within the difference time period.

[0101] The river flow is related to the roughness value. The calculation formula for river flow is:

[0102]

[0103] Where Q is the flow rate at a certain section in the river channel, m3 / s; K is the roughness value of the river channel, A is the water flow area at a certain section in the river channel, m2; R is the hydraulic radius, m; i is the slope of the river channel bottom.

[0104] According to the river flow calculation formula, when the main inflow flow and water level are the same, and the terrain does not change, the corresponding roughness remains unchanged. When the roughness changes, the flow will change accordingly if the water level does not change. Therefore, this application corrects the main inflow flow based on the change in roughness of the river section where the difference occurs.

[0105] Specifically, the roughness change rate α from the first roughness value to the second roughness value is calculated as:

[0106]

[0107] Among them, when α is a positive value, it means that the flow rate of the main inflow into the reservoir is larger than the actual value; when α is a negative value, it means that the flow rate of the main inflow into the reservoir is smaller than the actual value.

[0108] Combined with the roughness change rate α, the flow correction coefficient is calculated as 1-α. Based on the flow correction coefficient, the dry inflow flow is corrected. The corrected dry inflow flow Q fix for:

[0109] Q fix =Q(1-α).

[0110] For example, referring to Table 1, the first roughness value of the river section B is 22.22, the second roughness value is 23.53, and the roughness change rate is:

[0111]

[0112] That is, the main inflow flow rate is 5.56% less than the actual value, so the flow correction coefficient is 1.0556. From September 14 to September 30, the corrected main inflow flow rate Qfix is:

[0113] Q fix =1.0556Q.

[0114] The method for correcting the main inflow flow of a river-type reservoir provided in this embodiment addresses the existing difficulty in accurately determining the main inflow flow under conditions of large flow fluctuations. By constructing a one-dimensional hydrodynamic model and utilizing the model's roughness variation to determine the deviation and correction results of the main inflow flow, the method addresses the existing technical deficiencies in determining the main inflow flow using water level-flow curves. It solves the problem of inaccurate determination of the main inflow flow under conditions of significant fluctuations in the dam's water level or main inflow flow. The entire correction process is computationally efficient, free of human intervention, and possesses distinct scientific merit, providing important technical support for the scientific scheduling and safe operation of river-type reservoirs.

[0115] In this embodiment, a device for correcting the flow rate of a main inflow of a river-type reservoir is also provided. The device is used to implement the above-mentioned embodiments and preferred embodiments, and the details that have been described will not be repeated. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the devices described in the following embodiments are preferably implemented in software, implementation using hardware, or a combination of software and hardware, is also possible and contemplated.

[0116] This embodiment provides a device for correcting the flow rate of a main inflow into a river-type reservoir. Figure 4 is a structural block diagram of a device for correcting flow rate of a main inflow into a river-type reservoir according to an embodiment of the present invention. Figure 4 As shown in the figure, the flow correction device for the main inflow of the river type reservoir includes:

[0117] Construction module 401 is used to construct a one-dimensional hydrodynamic model of a river-type reservoir based on the measured hydrological data and topographic data of the river-type reservoir. The measured hydrological data includes measured flow data and measured water level data. The measured flow data includes the main inflow reservoir flow, tributary inflow reservoir flow and along-line monitoring flow. The measured water level data includes the water level in front of the dam and the water level monitored along the dam.

[0118] The first calculation module 402 is used to adjust the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model and obtain the first roughness value of each river section; extend the calculation time period, and calculate the first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model.

[0119] The second calculation module 403 is used to adjust the first roughness value of each river section again to reduce the calculation error of the one-dimensional hydrodynamic model and obtain the second roughness value of each river section; based on the second roughness value, the one-dimensional hydrodynamic model is used to calculate the second calculated water level.

[0120] The determination module 404 is configured to determine a difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level, and the measured water level.

[0121] The correction module 405 is used to correct the main inflow reservoir flow within the difference time period.

[0122] In some optional implementations, the first calculation module 402 includes:

[0123] The first calculation unit is used to adjust the initial roughness value of each river section for a first preset time period until the error between the calculated water level calculated by the one-dimensional hydrodynamic model of the river-type reservoir and the first measured water level is less than an error threshold, thereby obtaining the first roughness value of each river section, where the first measured water level is the measured water level for the first preset time period.

[0124] The second calculation unit is used to calculate the first calculated water level based on the first roughness value of each river section using a one-dimensional hydrodynamic model of the river-type reservoir for a second preset time period; the second preset time period includes the first preset time period and its preset extended time period.

[0125] In some optional implementations, the second calculation module 403 includes:

[0126] The third calculation unit is used to adjust the first roughness value of each river section for a second preset time period until the error between the calculated water level calculated by the one-dimensional hydrodynamic model of the river-type reservoir and the second measured water level is less than an error threshold, thereby obtaining the second roughness value of each river section; the second measured water level is the measured water level for the second preset time period.

[0127] The fourth calculation unit is configured to calculate a second calculated water level based on the second roughness value of each river section using a one-dimensional hydrodynamic model.

[0128] In some optional implementations, the first computing unit includes:

[0129] The first calibration subunit is used to adjust the initial roughness value of each river section section by section according to a first calibration sequence, and the first calibration sequence is the sequence from downstream to upstream.

[0130] In some optional implementations, the third computing unit includes:

[0131] The second calibration subunit is used to adjust the first roughness value of each river section section by section according to the second calibration sequence, and the second calibration sequence is from upstream to downstream.

[0132] In some optional implementations, the first calculation module 402 includes:

[0133] The acquisition unit is used to obtain the information of observation stations along the river type reservoir.

[0134] The division unit is used to divide the river channel of the river-type reservoir into multiple river sections based on the information of the observation stations along the river.

[0135] The assignment unit is used to assign fixed values ​​to the initial roughness values ​​of multiple river sections.

[0136] In some optional implementations, the determining module 404 includes:

[0137] The first determining unit is configured to determine a difference river section based on a difference between the first roughness value and the second roughness value.

[0138] The second determining unit is configured to determine, for the different river section, a difference time period within a second preset time period based on the difference between the first calculated water level, the second calculated water level and the second measured water level.

[0139] In some optional implementations, the correction module 405 includes:

[0140] The fifth calculation unit is configured to calculate a roughness change rate when the first roughness value changes to a second roughness value.

[0141] The correction unit is used to correct the flow of the main inflow reservoir based on the flow correction coefficient, and the sum of the flow correction coefficient and the roughness change rate is one.

[0142] The further functional description of each of the above modules and units is the same as that of the above corresponding embodiments and will not be repeated here.

[0143] The flow correction device for the main inflow of a river-type reservoir in this embodiment is presented in the form of a functional unit, where the unit refers to an ASIC (Application Specific Integrated Circuit) circuit, a processor and memory that executes one or more software or fixed programs, and / or other devices that can provide the above functions.

[0144] The embodiment of the present invention also provides a computer device having the above Figure 4 The shown figure shows the flow correction device for the main inflow of a river-type reservoir.

[0145] See also Figure 5 , Figure 5 is a structural diagram of a computer device provided by an optional embodiment of the present invention, such as Figure 5 As shown, the computer device includes: one or more processors 10, memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components utilize different buses to communicate with each other and can be installed on a common mainboard or installed in other ways as needed. The processor can process the instructions executed in the computer device, including instructions stored in the memory or on the memory to display the graphical information of the GUI on an external input / output device (such as, a display device coupled to the interface). In some optional embodiments, if necessary, multiple processors and / or multiple buses can be used together with multiple memories and multiple memories. Equally, multiple computer devices can be connected, and each device provides part of the necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 5 A processor 10 is taken as an example.

[0146] The processor 10 may be a central processing unit, a network processor, or a combination thereof. The processor 10 may further include a hardware chip. The hardware chip may be an application-specific integrated circuit, a programmable logic device, or a combination thereof. The programmable logic device may be a complex programmable logic device, a field programmable gate array, a general purpose array logic, or any combination thereof.

[0147] The memory 20 stores instructions that can be executed by at least one processor 10, so as to enable at least one processor 10 to execute the method shown in the above embodiment.

[0148] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system and application programs required for at least one function; the data storage area may store data created based on the use of the computer device, etc. In addition, the memory 20 may include a high-speed random access memory, and may also include a non-transient memory, such as at least one disk storage device, a flash memory device, or other non-transient solid-state storage device. In some optional embodiments, the memory 20 may optionally include a memory remotely located relative to the processor 10, and these remote memories may be connected to the computer device via a network. Examples of the above-mentioned network include, but are not limited to, the Internet, an intranet, a local area network, a mobile communication network, and combinations thereof.

[0149] The memory 20 may include a volatile memory, such as a random access memory; the memory may also include a non-volatile memory, such as a flash memory, a hard disk or a solid-state drive; the memory 20 may also include a combination of the above types of memory.

[0150] The computer device further includes an input device 30 and an output device 40. The processor 10, the memory 20, the input device 30 and the output device 40 may be connected via a bus or other means. Figure 5 The bus connection is taken as an example.

[0151] The input device 30 can receive input digital or character information and generate key signal input related to user settings and function control of the computer device, such as a touch screen, a keypad, a mouse, a trackpad, a touch pad, an indicator stick, one or more mouse buttons, a trackball, a joystick, etc. The output device 40 can include a display device, an auxiliary lighting device (e.g., an LED), and a tactile feedback device (e.g., a vibration motor). The above-mentioned display device includes but is not limited to a liquid crystal display, a light emitting diode, a display, and a plasma display. In some optional embodiments, the display device can be a touch screen.

[0152] The embodiments of the present application further provide a computer readable storage medium, and the method according to the embodiments of the present application can be implemented in hardware, firmware, or recorded in a storage medium, or stored in a remote storage medium or a non-transitory machine readable storage medium and downloaded to a local storage medium through network, so that the method described herein can be processed by such software on a storage medium using a general purpose computer, a special purpose processor, or programmable or special hardware. The storage medium can be a magnetic disk, an optical disk, a read-only memory, a random access memory, a flash memory, a hard disk, or a solid state disk, etc. Further, the storage medium can also include a combination of the above-mentioned memories. It can be understood that the computer, the processor, the microprocessor controller, or the programmable hardware includes a storage component that can store or receive software or computer code, when the software or computer code is accessed and executed by the computer, the processor, or the hardware, the method shown in the above embodiments is implemented.

[0153] Part of the present application can be applied as a computer program product, for example, computer program instructions, when executed by a computer, the operation of the computer can invoke or provide the method and / or technical solutions according to the present application. Those skilled in the art should understand that the form of computer program instructions in computer readable medium includes but is not limited to source file, executable file, installation package file, etc. Correspondingly, the way of computer program instructions executed by computer includes but is not limited to: the computer directly executes the instructions, or the computer compiles the instructions and then executes the corresponding compiled program, or the computer reads and executes the instructions, or the computer reads and installs the instructions and then executes the corresponding installed program. Here, the computer readable medium can be any available computer readable storage medium or communication medium accessible to the computer.

[0154] Although the embodiments of the present application are described in conjunction with the accompanying drawings, various modifications and changes can be made by those skilled in the art without departing from the spirit and scope of the present application, and such modifications and changes fall within the scope defined by the appended claims.

Claims

1. A method for correcting the flow rate of a main inflow of a river-type reservoir, characterized in that: The method comprises: Based on the measured hydrological data and topographic data of the river-type reservoir, a one-dimensional hydrodynamic model of the river-type reservoir is constructed, wherein the measured hydrological data includes measured flow data and measured water level data, wherein the measured flow data includes the main inflow reservoir flow, the tributary inflow reservoir flow, and the flow monitored along the channel, and the measured water level data includes the water level in front of the dam and the water level monitored along the channel; Adjusting the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model and obtain a first roughness value of each river section; extending the calculation time period, and calculating a first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model; Adjusting the first roughness value of each river section again to reduce the calculation error of the one-dimensional hydrodynamic model, thereby obtaining a second roughness value of each river section; and calculating a second calculated water level using the one-dimensional hydrodynamic model based on the second roughness value. determining a difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level, and the measured water level; The main inflow reservoir flow rate within the difference time period is corrected.

2. The method for correcting the flow rate of the main inflow of a river-type reservoir according to claim 1 is characterized in that: The initial roughness value of each river section is adjusted to reduce the calculation error of the one-dimensional hydrodynamic model, thereby obtaining a first roughness value of each river section; Extending the calculation time period, and calculating a first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model, includes: For a first preset time period, adjusting the initial roughness value of each river section until the error between a calculated water level calculated by the one-dimensional hydrodynamic model of the channel-type reservoir and a first measured water level is less than an error threshold, thereby obtaining a first roughness value for each river section, where the first measured water level is the measured water level during the first preset time period; For a second preset time period, based on the first roughness value of each river section, a first calculated water level is calculated using the one-dimensional hydrodynamic model of the river-type reservoir; the second preset time period includes the first preset time period and its preset extended time period.

3. The method for correcting the flow rate of the main inflow of a river-type reservoir according to claim 2 is characterized in that: said adjusting the first roughness value of each river section again to reduce the calculation error of the one-dimensional hydrodynamic model, and obtaining a second roughness value of each river section; The method of calculating a second calculated water level based on the second roughness value using the one-dimensional hydrodynamic model includes: For the second preset time period, adjusting the first roughness value of each river section until the error between the calculated water level calculated by the one-dimensional hydrodynamic model of the channel-type reservoir and the second measured water level is less than an error threshold, thereby obtaining a second roughness value for each river section; the second measured water level is the measured water level during the second preset time period; Based on the second roughness value of each river section, a second calculated water level is calculated using the one-dimensional hydrodynamic model.

4. The method for correcting the flow rate of the main inflow of a river-type reservoir according to any one of claims 1 to 3, characterized in that: The adjustment of the initial roughness value of each river section includes: Adjusting the initial roughness value of each river section section by section according to a first calibration order, wherein the first calibration order is from downstream to upstream; The adjusting the first roughness value of each river section includes: The first roughness value of each river section is adjusted section by section according to a second calibration order, where the second calibration order is from upstream to downstream.

5. The method for correcting the flow rate of the main inflow of a river-type reservoir according to claim 1 is characterized in that: The adjusting of the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model and obtain the first roughness value of each river section includes: Obtaining information of observation stations along the river-type reservoir; Dividing the river channel of the river-type reservoir into multiple river sections based on the information of the observation stations along the river; The initial roughness values ​​of the plurality of river sections are assigned fixed values.

6. The method for correcting the flow rate of the main inflow of a river-type reservoir according to claim 1, characterized in that: The determining of the difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level, and the second measured water level includes: determining a difference river section based on a difference between the first roughness value and the second roughness value; For the differential river section, a differential time period is determined within the second preset time period based on the difference between the first calculated water level, the second calculated water level and the second measured water level.

7. The method for correcting the flow rate of the main inflow of a river-type reservoir according to claim 1, characterized in that: The correction of the main inflow to the reservoir flow rate within the difference time period includes: Calculating a roughness change rate from the first roughness value to the second roughness value; The flow rate of the main inflow reservoir is corrected based on a flow correction coefficient, and the sum of the flow correction coefficient and the roughness change rate is one.

8. A flow correction device for the main inflow of a river-type reservoir, characterized in that: The device comprises: a construction module for constructing a one-dimensional hydrodynamic model of a river-type reservoir based on measured hydrological data and topographic data of the river-type reservoir, wherein the measured hydrological data includes measured flow data and measured water level data, wherein the measured flow data includes main inflow reservoir flow, tributary inflow reservoir flow, and flow monitored along the reservoir, and the measured water level data includes water level in front of the dam and water level monitored along the reservoir; a first calculation module, configured to adjust the initial roughness value of each river section to reduce the calculation error of the one-dimensional hydrodynamic model and obtain a first roughness value of each river section; extend the calculation time period, and calculate a first calculated water level based on the first roughness value using the one-dimensional hydrodynamic model; a second calculation module, configured to adjust the first roughness value of each river section again to reduce the calculation error of the one-dimensional hydrodynamic model, thereby obtaining a second roughness value of each river section; and calculate a second calculated water level based on the second roughness value using the one-dimensional hydrodynamic model; a determination module, configured to determine a difference time period by combining the first roughness value, the second roughness value, the first calculated water level, the second calculated water level, and the measured water level; A correction module is used to correct the main inflow reservoir flow within the difference time period.

9. A computer device, characterized in that: include: A memory and a processor, wherein the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the method for correcting the flow rate of the main inflow of a river-type reservoir according to any one of claims 1 to 7 by executing the computer instructions.

10. A computer-readable storage medium, characterized in that The computer-readable storage medium stores computer instructions, and the computer instructions are used to enable a computer to execute the method for correcting the flow rate of the main inflow of a river-type reservoir according to any one of claims 1 to 7.

Images