A dynamic adjustment-based iterative calculation method for water balance in storage belts

Abstract

The present invention provides an iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment, which relates to the technical field of water resource management of regulating storage belts / reservoirs and ponds. The method comprises: constructing a basic model framework for calculating water balance of a regulating storage belt, setting input time series data and parameters, and then performing a first-level sequential dynamic simulation of water volume of the regulating storage belt, and processing the obtained time series data through a low-pass filtering process; then performing a second-level sequential dynamic simulation of water volume of the regulating storage belt, calculating the relative difference value of the total inflow of the regulating storage belt before and after the adjustment to perform a convergence check, performing iterative calculation until the relative difference value is less than a set threshold or reaches a maximum number of iterations, and outputting the iterative calculation result of water balance of the regulating storage belt. The present invention introduces inflow uncertainty adjustment and additional outflow adjustment to express the uncertainty caused by manual control of gates and incomplete monitoring data, significantly improving the effectiveness of the calculation results, and can provide reliable technical support for the scientific management of water resources in the regulating storage belt.

Description

Technical Field

[0001] The invention relates to the technical field of water resource management of a regulating storage belt / reservoir, and in particular to an iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment. Background Art

[0002] The storage belt is a comprehensive utilization project that integrates sewage interception, storage, purification, reuse and ecological landscape. It intercepts water resources such as initial rainwater, farmland tailwater, and leaked sewage for recycling and reduces pollution load emissions. The storage belt can not only effectively control the entry of pollutants into downstream water bodies, but also improve water quality through purification and treatment, and ultimately achieve sustainable use of water resources. In order to scientifically manage the water resources in the storage belt, it is necessary to accurately calculate the water volume in the storage belt.

[0003] In actual operation, there are multiple challenges in accurately calculating the amount of water in the storage belt. First, the shape of the storage belt is irregular, and there are many inlets and outlets. When monitoring the amount of water in and out, it is impossible to cover all the inlets and outlets in real time and comprehensively. Therefore, it is difficult to fully calculate the water balance in the storage belt based on the data obtained by monitoring. Secondly, in general, the main inlets and outlets of the storage belt are controlled by gates, and the degree of opening and closing of the gates is mostly adjusted by on-site operators based on personal experience. Therefore, the inflow ratio of upstream water into the storage belt has a certain degree of randomness. In the traditional method of calculating the amount of water in the storage belt, it is usually simplified to simplify the inflow and outflow of the storage belt, such as setting only one inlet and one outlet, and using a static calculation of the inflow ratio of the storage belt, that is, to calculate the change of water in the storage belt in the form of a fixed inflow ratio. Since this calculation method is inconsistent with the actual operation of the storage belt, it will lead to a huge deviation between the calculation result and the actual water volume of the storage belt, which will affect the management and scheduling of the storage belt and bring hidden dangers and risks to the management and protection of water volume and quality of downstream water bodies. Summary of the invention

[0004] Purpose of the invention: The purpose of the present invention is to provide an iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment, which can reflect the inflow and outflow of multiple channels of the regulating storage belt and dynamically adjust the inflow and outflow ratio of the regulating storage belt.

[0005] Technical solution: A method for iterative calculation of water balance in a storage belt based on dynamic adjustment, comprising the following steps:

[0006] S1. Construct the basic model framework for water balance calculation of the storage belt: collect the time series data of the target storage belt, including the time series data of potential inflow, known outflow, measured water storage, measured rainfall, and measured evaporation in the same simulation period; set the model parameters, including overflow threshold, dead storage capacity, inflow uncertainty adjustment factor, additional outflow ratio tolerance adjustment factor, maximum number of iterations, adjustment time window, and simulation period;

[0007] S2, using the basic model framework for water balance calculation of the regulating storage belt constructed in step S1 to perform a first-level sequential dynamic simulation of the regulating storage belt water volume, read the time series data within the simulation period, update the simulated water storage volume based on the principle of conservation of matter, and record the first-level sequential dynamic simulation results of the regulating storage belt water volume, including the time series of the calculated values ​​of the simulated water storage volume at the end of each day within the simulation period, the time series of the calculated values ​​of the overflow volume each day within the simulation period, and the time series of the calculated values ​​of the additional inflow volume each day within the simulation period;

[0008] S3, based on the first-level sequential dynamic simulation result of the water volume of the regulating and storage belt obtained in step S2, a low-pass filtering algorithm is used to adjust the inflow time series data and the additional outflow time series data of the regulating and storage belt;

[0009] S4, based on the inflow time series data and the additional outflow time series data obtained in step S3, perform a secondary sequential dynamic simulation of the storage belt, calculate the simulated water storage capacity based on the principle of conservation of matter, and record the secondary sequential dynamic simulation results of the water volume in the storage belt, including the time series of the calculated values ​​of the simulated water storage capacity at the end of each day in the simulation period, the time series of the calculated values ​​of the overflow volume each day in the simulation period, and the time series of the calculated values ​​of the additional outflow volume each day in the simulation period;

[0010] S5. Calculate the relative difference of the total inflow of the storage zone before and after the adjustment;

[0011] S6. If the relative difference value is less than the set threshold, stop the iteration and output the calculation results of the water balance of the regulating storage belt, including the time series data of the water storage at the end of each day, the adjusted inflow, the additional outflow and the overflow; if the relative difference value is greater than or equal to the set threshold, repeat steps S2 to S5 to continue the next round of iteration until the relative difference value is less than the set threshold or the maximum number of iterations is reached, and output the dynamically adjusted iterative calculation results of the water balance of the regulating storage belt, including the time series data of the water storage at the end of each day, the adjusted inflow, the additional outflow and the overflow.

[0012] Specifically, step S2 includes: performing water balance calculation according to the simulation time step of day scale:

[0013] Update the simulated water storage volume, and read the measured water storage volume on the first day of the simulation period from the measured water storage volume time series data as the initial value of the water storage volume in the regulation storage zone , based on the principle of conservation of matter, update the simulated water storage:

[0014] ,

[0015] Where: is the water storage capacity of the storage zone at the beginning of day t, is the simulated water storage capacity of the storage belt at the end of day t, is the potential inflow of the storage zone on day t, is the known outflow of the storage belt on day t, is the rainfall in the storage zone on day t, is the evaporation of the storage zone on day t;

[0016] If the above calculation If it is less than the dead storage capacity, the simulated water storage at the end of day t is updated to:

[0017] ,

[0018] If the simulated water storage at the end of day t is Greater than the measured water storage on day t+1 , then calculate and record the overflow on day t :

[0019] ,

[0020] If the simulated water storage at the end of day t is Less than the measured water storage volume on day t+1 , then calculate and record the additional inflow on day t :

[0021] ,

[0022] After completing the above adjustment calculations, update the simulated water storage at the end of day t :

[0023] ,

[0024] If overflow or additional inflow does not exist, the corresponding overflow or additional inflow is set to 0; the simulated water storage at the end of day t is The initial value of water storage calculated as the next time step , repeat the above steps to obtain the time series data of simulated water storage at the end of each day, daily overflow, and daily additional inflow during the simulation period.

[0025] Specifically, step S3 includes:

[0026] S31, based on the adjustment time window, dividing the simulation period into a plurality of non-overlapping time windows;

[0027] S32, based on the first-level sequential dynamic simulation results of the storage belt water volume obtained in step S2, calculate the cumulative inflow, cumulative outflow, cumulative overflow and cumulative additional inflow in each time window;

[0028] S33, calculating the adjustment coefficient of each time window, and using the adjustment coefficient to calculate and obtain the adjusted inflow time series data and the additional outflow time series data.

[0029] Specifically, in step S32, the accumulated inflow includes the potential inflow and the sum of rainfall in all time steps in the time window to be calculated; the accumulated outflow includes the known outflow and the sum of evaporation in all time steps in the time window to be calculated.

[0030] Specifically, in step S33, the adjustment coefficient includes an inflow adjustment coefficient and an additional outflow adjustment coefficient, and the calculation formula is:

[0031] ,

[0032] ,

[0033] ,

[0034] Where: is the flow adjustment requirement of the tcth time window, is the cumulative inflow in the tcth time window, is the cumulative overflow in the tcth time window, is the accumulated additional inflow in the tcth time window, is the inflow adjustment coefficient of the tcth time window, is the additional outflow adjustment coefficient of the tcth time window.

[0035] Preferably, step S33 further includes:

[0036] Constrain the inflow adjustment coefficient and additional outflow adjustment coefficient of the tcth time window:

[0037] ,

[0038] Where: is the inflow adjustment coefficient of the tcth time window, is the additional outflow adjustment coefficient of the tcth time window, X is the inflow uncertainty adjustment factor, , T is the additional outflow ratio tolerance adjustment factor.

[0039] Specifically, step S33 includes:

[0040] like ,and , then calculate the inflow adjustment coefficient of the tcth time window and correct the inflow time series in this time window:

[0041] ,

[0042] ,

[0043] Where: is the adjusted rainfall in the storage zone in the tcth time window, is the adjusted potential inflow of the storage zone in the tcth time window;

[0044] like and , then calculate the additional outflow adjustment coefficient of the tcth time window and increase the additional outflow time series in this time window:

[0045] ,

[0046] Where: is the additional outflow of the storage belt in the tcth time window;

[0047] like , no adjustment is made.

[0048] Specifically, step S4 includes: performing water balance calculation of the storage belt according to the simulation time step of day scale:

[0049] Update the simulated water storage volume, and read the measured water storage volume on the first day of the simulation period from the measured water storage volume time series data as the initial value of the water storage volume in the regulation storage zone , based on the principle of conservation of matter, update the simulated water storage:

[0050] ,

[0051] Where: is the water storage capacity of the storage zone at the beginning of day t, is the simulated water storage capacity of the storage belt at the end of day t, is the adjusted potential inflow of the storage zone on day t, is the known outflow of the storage belt on day t, is the adjusted rainfall in the storage zone on day t, is the evaporation of the storage zone on day t;

[0052] If the above calculation If it is less than the dead storage capacity, the simulated water storage at the end of day t is updated to:

[0053] ,

[0054] If the simulated water storage at the end of day t is Greater than the measured water storage on day t+1 , the flow rate is calculated and can be divided into the following three cases:

[0055] If the measured water storage volume on day t+1 is Greater than or equal to the overflow threshold :

[0056] ,

[0057] If the measured water storage volume on day t+1 is Less than overflow threshold , and is less than the end of the tth day :

[0058] ,

[0059] ,

[0060] limit Greater than or equal to 0;

[0061] If the simulated water storage at the end of day t is Less than or equal to the overflow threshold :

[0062] ,

[0063] ,

[0064] limit Greater than or equal to 0;

[0065] Record the overflow volume on day t calculated in this round , and record the updated additional outflow on day t :

[0066] ,

[0067] Update the simulated water storage at the end of day t again :

[0068] ,

[0069] If the overflow or additional outflow on the tth day does not exist, the corresponding overflow or additional outflow is set to 0; the simulated water storage at the end of the tth day is set to The initial value of water storage calculated as the next time step , repeat the above steps to obtain the time series data of water storage, adjusted inflow, additional outflow and overflow at the end of each day during the simulation period.

[0070] Specifically, in step S5, the relative difference value of the total inflow of the regulating storage zone before and after the adjustment is calculated as follows:

[0071] ,

[0072] Where: C is the relative difference value, is the total inflow of the storage zone before adjustment, is the adjusted total inflow into the storage zone.

[0073] Specifically, in step S6, repeating steps S2 to S5 to continue the next round of iteration includes: updating the adjustment time window, adding 1 day to the original adjustment time window, using the time series data updated in step S5 to execute step S2 of the next round of iteration, using the updated adjustment time window to execute step S3, and continuing to execute steps S4 and S5 to complete the next round of iteration.

[0074] Beneficial effects: Compared with the prior art, the significant effects of the present invention are: the present invention improves the traditional method of calculating the water balance of the regulating and storage belt by statistically calculating the information of the regulating and storage belt, introduces the uncertainty adjustment of the inflow and the additional outflow adjustment to express the uncertainty caused by the manual operation of the gate and the incomplete monitoring network, significantly improves the effectiveness of the calculation results, can reflect the multi-channel outflow of the regulating and storage belt as a whole, and can realize dynamic adjustment calculation at the same time, and adopts the iterative algorithm of low-pass filtering to improve the reliability of the calculation results of the water volume of the regulating and storage belt. The calculation method provided by the present invention can provide reliable data support for the scientific management of the water volume of the regulating and storage belt. BRIEF DESCRIPTION OF THE DRAWINGS

[0075] Figure 1 It is a flow chart of the method of the present invention.

[0076] Figure 2 It is a schematic diagram of the overflow of the regulating storage belt in Example 1 of the present invention.

[0077] Figure 3 Schematic diagram of potential inflow of the storage belt in Example 1 of the present invention.

[0078] Figure 4 It is a schematic diagram of the outflow flow of the storage device according to the first embodiment of the present invention.

[0079] Figure 5It is a schematic diagram of rainfall in the storage zone according to Example 1 of the present invention.

[0080] Figure 6 It is a schematic diagram comparing the simulated water storage capacity and the actually measured water storage capacity of the storage belt in Example 1 of the present invention. DETAILED DESCRIPTION

[0081] The present invention is further described below in conjunction with the accompanying drawings.

[0082] Example 1

[0083] See also Figure 1 As shown, this embodiment provides an iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment, comprising the following steps:

[0084] S1. Construct the basic model framework for calculating water balance in the storage belt. First, determine the input time series data, including:

[0085] Potential inflow: Potential inflow time series data, containing two columns: timestamp (day-scale interval) and inflow (m³ / s).

[0086] Known outflow: Time series data of known outflow, containing two columns: timestamp (day-scale interval) and outflow (m³ / s).

[0087] Measured water storage: time series data of measured water storage in the storage zone, including two columns: timestamp (day scale interval) and storage capacity (m³).

[0088] Rainfall: Measured rainfall time series, containing two columns: timestamp (hourly interval) and rainfall (mm / h).

[0089] Evaporation: Measured evaporation time series, including two columns: timestamp (hourly interval) and evaporation (mm / h).

[0090] Collect the above time series data in the target storage zone and store it in CSV format. After determining the input time series data, determine the basic parameters of the model, including:

[0091] Overflow threshold: The storage volume threshold that triggers overflow, set according to the actual design capacity of the storage tank.

[0092] Dead storage capacity: the minimum storage capacity of the regulating reservoir, which is set according to the actual design capacity of the regulating reservoir.

[0093] Inflow uncertainty adjustment factor: In view of the uncertainty of inflow into the storage zone, the maximum increase ratio of additional inflow relative to potential inflow is set.

[0094] Additional outflow ratio tolerance adjustment factor: In view of the uncertainty of storage belt outflow, the maximum ratio of additional outflow to inflow is set.

[0095] Maximum number of iterations: The maximum number of iterations for the model to simulate.

[0096] Adjust the time window: the initial value of the time window parameter of the low-pass filter algorithm in the model.

[0097] Simulation period: the set simulation start time and simulation end time.

[0098] After completing the above model settings, align the potential inflow, known outflow, measured water storage, rainfall and evaporation data to ensure that the above data cover the same time range and that the same time range includes the simulation period.

[0099] To facilitate calculation, the time series data is converted to a time step:

[0100] The units of the potential inflow and known outflow time series are converted from m³ / s to m³ / d using the following formula:

[0101]

[0102] Where: is the converted volume of potential inflow or known outflow on day t, in m³ / d, is the flow rate before conversion on day t, in m³ / s, and 86400 is the total number of seconds in a day.

[0103] The units of rainfall and evaporation time series are converted from mm / h to m³ / d, and the calculation formula is:

[0104]

[0105] Where: is the volume of rainfall or evaporation on day t after conversion, in m³ / d. is the surface area of ​​the storage zone or reservoir in the study area, in m 2 , It is the rainfall or evaporation at the i-th hour on the t-th day before conversion, in mm / h, 24 is the total number of hours in a day, and k is the unit conversion coefficient.

[0106] S2. Use the basic model framework for calculating the water balance in the regulating storage belt constructed in step S1 to perform a first-level sequential dynamic simulation of the water volume in the regulating storage belt.

[0107] S21. Calculate water balance using the day scale as the simulation time step:

[0108] Update the simulated water storage volume, and read the measured water storage volume on the first day of the simulation period from the measured water storage volume time series data as the initial value of the water storage volume in the regulation storage zone , based on the principle of conservation of matter, update the simulated water storage:

[0109] ,

[0110] Where: is the water storage capacity of the storage zone at the beginning of day t, is the simulated water storage capacity of the storage belt at the end of day t, is the potential inflow of the storage zone on day t, is the known outflow of the storage belt on day t, is the rainfall in the storage zone on day t, is the evaporation of the storage zone on day t;

[0111] If the above calculation If it is less than the dead storage capacity, the simulated water storage at the end of day t is updated to:

[0112] ,

[0113] If the simulated water storage at the end of day t is Greater than the measured water storage on day t+1 , then calculate and record the overflow on day t :

[0114] ,

[0115] If the simulated water storage at the end of day t is Less than the measured water storage volume on day t+1 , then calculate and record the additional inflow on day t :

[0116] ,

[0117] After completing the above adjustment calculations, update the simulated water storage at the end of day t :

[0118] ,

[0119] If overflow or additional inflow does not exist, the corresponding overflow or additional inflow is set to 0; the simulated water storage at the end of day t is The initial value of water storage calculated as the next time step .

[0120] S22, repeat the above steps and record the first-level sequential dynamic simulation results of the water volume in the storage belt, including:

[0121] Time series of calculated values ​​of simulated water storage at the end of each day during the simulation period, time series of calculated values ​​of overflow volume per day during the simulation period, and time series of calculated values ​​of additional inflow volume per day during the simulation period.

[0122] S3, adjust the inflow and outflow of the storage and regulation based on the low-pass filtering algorithm:

[0123] S31. Divide the simulation period into a plurality of non-overlapping time windows according to the adjustment time window.

[0124] S32, based on the first-level sequential dynamic simulation results of the water volume in the storage belt obtained in step S2, calculate the cumulative volume in each time window:

[0125] Cumulative inflow: The potential inflow and rainfall for all time steps within this time window.

[0126] Cumulative outflow: The known outflow and evaporation at all time steps within this time window.

[0127] Cumulative overflow volume: the total overflow volume within the time window.

[0128] Cumulative additional inflow: the sum of additional inflows in this time window.

[0129] S33, determine the adjustment coefficient of each time window and the adjusted inflow and outflow:

[0130] First, calculate the traffic adjustment requirements for each time window :

[0131] ,

[0132] Where: is the flow adjustment requirement of the tcth time window, is the cumulative inflow in the tcth time window, is the cumulative overflow in the tcth time window, is the accumulated additional inflow in the tcth time window.

[0133] like ,and , then calculate the inflow adjustment coefficient and correct the inflow time series within the time window:

[0134]

[0135] Where: is the inflow adjustment coefficient of the tcth time window.

[0136] Inflow adjustment factor To constrain:

[0137]

[0138] Where: is the inflow uncertainty adjustment factor, .

[0139] Using the Inflow Adjustment Factor Correct the inflow time series:

[0140] ,

[0141] ,

[0142] Where: is the adjusted rainfall in the storage zone in the tcth time window, is the adjusted potential inflow of the storage zone in the tcth time window;

[0143] like and , then calculate the additional outflow adjustment coefficient and increase the additional outflow time series within the window:

[0144]

[0145] Where: is the additional outflow adjustment coefficient of the tcth time window.

[0146] Since the additional outflow and inflow are highly correlated, in this step, the additional outflow time series is constructed based on the inflow time series instead of modifying the known outflow time series, because the known outflow time series usually records the process controlled by humans, while the additional outflow time series records the natural process that is not controlled by humans.

[0147] Adjustment factor for additional outflow To constrain:

[0148]

[0149] Where: It is the tolerance adjustment factor of the additional outflow ratio.

[0150] Using the additional outflow adjustment factor Generate additional outflow time series:

[0151] ,

[0152] Where: is the additional outflow of the storage zone in the tcth time window.

[0153] like , no adjustment is made.

[0154] S34. Record the adjusted inflow and outflow time series, including:

[0155] The time series of adjusted rainfall during the simulation period, the time series of adjusted potential inflow during the simulation period, and the time series of additional outflow during the simulation period.

[0156] S4. Use the inflow and outflow time series adjusted in step S3 to perform secondary sequential dynamic simulation of the storage zone.

[0157] S41. Carry out water balance calculations in sequence according to the simulation time step of day scale:

[0158] Similar to step S2, the water storage is first updated, and the measured water storage on the first day of the simulation period is read from the measured water storage time series data as the initial value of the water storage in the regulation and storage belt. , based on the principle of conservation of matter, update the simulated water storage:

[0159] ,

[0160] Where: is the water storage capacity of the storage zone at the beginning of day t, is the simulated water storage capacity of the storage belt at the end of day t, is the adjusted potential inflow of the storage zone on day t, is the known outflow of the storage belt on day t, is the adjusted rainfall in the storage zone on day t, is the evaporation of the storage zone on day t;

[0161] If the above calculation If it is less than the dead storage capacity, the simulated water storage at the end of day t is updated to:

[0162] ,

[0163] If the simulated water storage at the end of day t is Greater than the measured water storage on day t+1 , the flow rate is calculated and can be divided into the following three cases:

[0164] If the measured water storage volume on day t+1 is Greater than or equal to the overflow threshold :

[0165] ,

[0166] If the measured water storage volume on day t+1 is Less than overflow threshold , and is less than the end of the tth day :

[0167] ,

[0168] ,

[0169] limit Greater than or equal to 0;

[0170] If the simulated water storage at the end of day t is Less than or equal to the overflow threshold :

[0171] ,

[0172] ,

[0173] limit Greater than or equal to 0.

[0174] Record the overflow volume on day t calculated in this round , and record the updated additional outflow on day t :

[0175] ,

[0176] Update the simulated water storage at the end of day t again :

[0177] ,

[0178] If the overflow or additional outflow on the tth day does not exist, the corresponding overflow or additional outflow is set to 0; the simulated water storage at the end of the tth day is set to The initial value of water storage calculated as the next time step .

[0179] S42, repeat the above steps, and record the secondary sequential dynamic simulation results of the water volume in the storage belt, including:

[0180] Time series of calculated values ​​of simulated water storage at the end of each day during the simulation period, time series of calculated values ​​of overflow volume per day during the simulation period, and time series of calculated values ​​of additional outflow volume per day during the simulation period.

[0181] S5. Check the convergence of the above calculations and calculate the relative difference in total inflow before and after the adjustment:

[0182]

[0183] Where: C is the relative difference value, is the total inflow of the storage zone before adjustment, is the adjusted total inflow into the storage zone.

[0184] S6. If the relative difference value , the iteration stops and the time series data of the water storage volume, adjusted inflow, additional outflow and overflow at the end of each day in the storage zone are returned. , it has not converged, and the next round of iteration is continued, and steps S2 to S5 are repeated, specifically:

[0185] Update and adjust the time window setting to add 1 day to the original time window setting; use the time series data updated in step S34 to perform the calculation in step S2; use the updated time window setting to perform the calculation in step S3, and continue to execute steps S4 and S5, repeat the iterations until the results converge or the set maximum number of iterations is reached, and obtain the dynamically adjusted iterative calculation results of the water balance in the regulating storage belt, including the time series data of the water storage volume at the end of each day in the regulating storage belt, the adjusted inflow, the additional outflow and the overflow.

[0186] In this embodiment, the data of a certain storage belt in 2023 is substituted into the above calculation method for statistics and calculation. For key data, please refer to Figures 2 to 5 As shown, Figure 2 This is a schematic diagram of the overflow of the storage belt, and the overflow is concentrated in the period from May to December; Figure 3 This is a schematic diagram of the potential inflow of the storage belt, and the potential inflow is also concentrated in the period from May to December; Figure 4 This is a schematic diagram of the outflow of the storage belt, including outflow data of 4 different pumping stations, evaporation data and additional outflow data; Figure 5 This is a schematic diagram of rainfall in the storage zone.

[0187] After calculation and statistics, the potential inflow of the regulating belt in 2023 is 8.8 million cubic meters, the total water pumping volume of the pumping station (known outflow) is 2.44 million cubic meters, the evaporation is 100,000 cubic meters, the rainfall is 120,000 cubic meters, the overflow threshold of the regulating belt is 180,000 cubic meters, and the dead storage capacity of the regulating belt is 40,000 cubic meters. The above data and the time series data of the water storage volume of the regulating belt are used as input for iterative calculation. After iteration, the water balance data of the regulating belt is obtained. In 2023, the actual inflow of the regulating belt is 2.62 million cubic meters, the additional outflow is 120,000 cubic meters, and the overflow is 60,000 cubic meters.

[0188] Please refer to Figure 6 As shown in the figure, it is a comparison chart between the simulated water storage capacity of the regulating storage belt and the measured water storage capacity of the regulating storage belt. Figure 6 It can be seen intuitively that in the statistical period from May to December, the simulated water storage in the storage belt is highly consistent with the measured water storage in the storage belt, which confirms the effectiveness of this method.

Claims

1. A method for iterative calculation of water balance in a storage belt based on dynamic adjustment, characterized in that: The following steps are involved: S1. Construct the basic model framework for water balance calculation of the storage belt: collect the time series data of the target storage belt, including the time series data of potential inflow, known outflow, measured water storage, measured rainfall, and measured evaporation in the same simulation period; set the model parameters, including overflow threshold, dead storage capacity, inflow uncertainty adjustment factor, additional outflow ratio tolerance adjustment factor, maximum number of iterations, adjustment time window, and simulation period; S2, using the basic model framework for water balance calculation of the regulating storage belt constructed in step S1 to perform a first-level sequential dynamic simulation of the regulating storage belt water volume, read the time series data within the simulation period, update the simulated water storage volume based on the principle of conservation of matter, and record the first-level sequential dynamic simulation results of the regulating storage belt water volume, including the time series of the calculated values ​​of the simulated water storage volume at the end of each day within the simulation period, the time series of the calculated values ​​of the overflow volume each day within the simulation period, and the time series of the calculated values ​​of the additional inflow volume each day within the simulation period; S3, based on the first-level sequential dynamic simulation result of the water volume of the regulating and storage belt obtained in step S2, a low-pass filtering algorithm is used to adjust the inflow time series data and the additional outflow time series data of the regulating and storage belt; S4, based on the inflow time series data and the additional outflow time series data obtained in step S3, perform a secondary sequential dynamic simulation of the storage belt, update the simulated water storage capacity based on the principle of conservation of matter, and record the secondary sequential dynamic simulation results of the water volume in the storage belt, including the time series of the calculated values ​​of the simulated water storage capacity at the end of each day in the simulation period, the time series of the calculated values ​​of the overflow volume each day in the simulation period, and the time series of the calculated values ​​of the additional outflow volume each day in the simulation period; S5, calculate the relative difference of the total inflow of the storage belt before and after the adjustment; S6. If the relative difference value is less than the set threshold, stop the iteration and output the calculation results of the water balance of the regulating storage belt, including the time series data of the water storage at the end of each day, the adjusted inflow, the additional outflow and the overflow; if the relative difference value is greater than or equal to the set threshold, repeat steps S2 to S5 to continue the next round of iteration until the relative difference value is less than the set threshold or the maximum number of iterations is reached, and output the iterative calculation results of the water balance of the dynamically adjusted regulating storage belt, including the time series data of the water storage at the end of each day, the adjusted inflow, the additional outflow and the overflow.

2. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 1 is characterized in that: The step S2 comprises: performing water balance calculation according to the simulation time step on a daily scale: Update the simulated water storage volume, and read the measured water storage volume on the first day of the simulation period from the measured water storage volume time series data as the initial value of the water storage volume in the regulation storage zone , based on the principle of conservation of matter, update the simulated water storage: ， Where: is the water storage capacity of the storage zone at the beginning of day t, is the simulated water storage capacity of the storage belt at the end of day t, is the potential inflow of the storage zone on day t, is the known outflow of the storage belt on day t, is the rainfall in the storage zone on day t, is the evaporation of the storage zone on day t; If the above calculation If it is less than the dead storage capacity, the simulated water storage at the end of day t is updated to: ， If the simulated water storage at the end of day t is Greater than the measured water storage on day t+1 , then calculate and record the overflow on day t : ， If the simulated water storage at the end of day t is Less than the measured water storage volume on day t+1 , then calculate and record the additional inflow on day t : ， After completing the above calculations, update the simulated water storage at the end of day t : ， If overflow or additional inflow does not exist, the corresponding overflow or additional inflow is set to 0; the simulated water storage at the end of day t is The initial value of water storage calculated as the next time step , repeat the above steps to obtain the time series data of simulated water storage at the end of each day, daily overflow, and daily additional inflow during the simulation period.

3. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 1 is characterized in that: The step S3 comprises: S31, based on the adjustment time window, dividing the simulation period into a plurality of non-overlapping time windows; S32, based on the first-level sequential dynamic simulation results of the storage belt water volume obtained in step S2, calculate the cumulative inflow, cumulative outflow, cumulative overflow and cumulative additional inflow in each time window; S33, calculating the adjustment coefficient of each time window, and using the adjustment coefficient to calculate and obtain the adjusted inflow time series data and the additional outflow time series data.

4. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 3 is characterized in that: In step S32, the accumulated inflow includes the potential inflow and the sum of rainfall in all time steps in the time window to be calculated; the accumulated outflow includes the known outflow and the sum of evaporation in all time steps in the time window to be calculated.

5. The iterative calculation method for water balance of regulating storage belt based on dynamic adjustment according to claim 3 is characterized by: In step S33, the adjustment coefficient includes an inflow adjustment coefficient and an additional outflow adjustment coefficient, and the calculation formula is: ， ， ， Where: is the flow adjustment requirement of the tcth time window, is the cumulative inflow in the tcth time window, is the cumulative overflow in the tcth time window, is the accumulated additional inflow in the tcth time window, is the inflow adjustment coefficient of the tcth time window, is the additional outflow adjustment coefficient of the tcth time window.

6. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 5 is characterized in that: The step S33 further comprises: Constrain the inflow adjustment coefficient and additional outflow adjustment coefficient of the tcth time window: ， Where: is the inflow adjustment coefficient of the tcth time window, is the additional outflow adjustment coefficient of the tcth time window, X is the inflow uncertainty adjustment factor, , T is the additional outflow ratio tolerance adjustment factor.

7. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 6 is characterized by: The step S33 further comprises: like ,and , then calculate the inflow adjustment coefficient of the tcth time window and correct the inflow time series in this time window: ， ， Where: is the adjusted rainfall in the storage zone in the tcth time window, is the adjusted potential inflow of the storage zone in the tcth time window; like and , then calculate the additional outflow adjustment coefficient of the tcth time window and increase the additional outflow time series in this time window: ， Where: is the additional outflow of the storage belt in the tcth time window; like , no adjustment is made.

8. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 1 is characterized in that: The step S4 comprises: performing water balance calculation of the storage belt according to the simulation time step of day scale: Update the simulated water storage volume, and read the measured water storage volume on the first day of the simulation period from the measured water storage volume time series data as the initial value of the water storage volume in the regulation storage zone , based on the principle of conservation of matter, update the simulated water storage: ， Where: is the water storage capacity of the storage zone at the beginning of day t, is the simulated water storage capacity of the storage belt at the end of day t, is the adjusted potential inflow of the storage zone on day t, is the known outflow of the storage belt on day t, is the adjusted rainfall in the storage zone on day t, is the evaporation of the storage zone on day t; If the above calculation If it is less than the dead storage capacity, the simulated water storage at the end of day t is updated to: ， If the simulated water storage at the end of day t is Greater than the measured water storage on day t+1 , the flow rate is calculated and can be divided into the following three cases: If the measured water storage volume on day t+1 is Greater than or equal to the overflow threshold : ， If the measured water storage volume on day t+1 is Less than overflow threshold , and is less than the end of the tth day : ， ， limit Greater than or equal to 0; If the simulated water storage at the end of day t is Less than or equal to the overflow threshold : ， ， limit Greater than or equal to 0; Record the overflow volume on day t calculated in this round , and record the updated additional outflow on day t : ， Update the simulated water storage at the end of day t again : ， If the overflow or additional outflow on the tth day does not exist, the corresponding overflow or additional outflow is set to 0; the simulated water storage at the end of the tth day is set to The initial value of water storage calculated as the next time step , repeat the above steps to obtain the time series data of water storage, adjusted inflow, additional outflow and overflow at the end of each day during the simulation period.

9. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 1 is characterized in that: In step S5, the relative difference value of the total inflow of the regulating and storage belt before and after the adjustment is calculated as follows: ， Where: C is the relative difference value, is the total inflow of the storage zone before adjustment, is the adjusted total inflow into the storage zone.

10. The iterative calculation method for water balance of a regulating storage belt based on dynamic adjustment according to claim 1 is characterized in that: In step S6, repeating steps S2 to S5 to continue the next round of iteration includes: updating the time window, adding 1 day to the original time window, using the time series data updated in step S5 to execute step S2 of the next round of iteration, using the updated adjusted time window to execute step S3, and continuing to execute steps S4 and S5 to complete the next round of iteration.

Images