A method, device, computer equipment and storage medium for inbound flow reverse estimation

Abstract

The present invention relates to the technical field of hydrological simulation, and discloses a method, device, computer equipment and storage medium for inverse calculation of inflow flow. The method comprises: obtaining water level data in front of the dam, a reservoir capacity relationship curve and outflow flow data; performing low-pass filtering on the water level data in front of the dam to generate water level data after low-pass filtering; and inversely calculating the inflow flow according to the water level data, the reservoir capacity relationship curve and the outflow flow data after low-pass filtering. The present invention realizes noise reduction of the water level data in front of the dam by performing low-pass filtering on the water level data in front of the dam, thereby reducing the influence of the water level fluctuation in front of the dam caused by external interference and the high-frequency fluctuation of the water level on the accuracy of the calculation of the inflow flow. The smoothing processing of the water level data in front of the dam is realized, and the accuracy of inverse calculation of the inflow flow according to the water level data is improved. The problem of low calculation accuracy when using the smoothing method to process the inflow flow in the related art is solved.

Description

Technical Field

[0001] The present invention relates to the technical field of hydrological simulation, and in particular to a method, device, computer equipment and storage medium for inverse estimation of inflow. Background Art

[0002] Most of the inflow of reservoirs is obtained by reverse calculation from the water level in front of the dam and the outflow using the principle of water balance, but the calculated flow process generally has obvious sawtooth fluctuations, and even negative values ​​appear. At present, in order to solve this problem, the calculation results of the inflow are usually smoothed to eliminate abnormal sawtooth in the data. Commonly used correction methods include sliding average method and five-point cubic method. The calculation of sliding average is based on the inflow data of the past period of time, but these data cannot accurately reflect the current inflow situation. Therefore, the smoothed flow data will lag, resulting in deviations in the peak time. And the degree of smoothing of the sliding average is affected by the window size. If the window is too large, the flow process will be over-flattened, causing the peak flow of the flood process to deviate. If the window is too small, the smoothing effect is poor. The five-point cubic method is to smooth the data of five points by fitting a cubic polynomial, but there may be distortion after smoothing the steeper inflow curve. Therefore, there is a problem in the prior art that the calculation accuracy is not high when using the smoothing method to process the inflow. Summary of the invention

[0003] In view of this, the present invention provides an inflow flow reverse estimation method, device, computer equipment and storage medium to solve the problem of low calculation accuracy when using conventional smoothing methods to process inflow flow.

[0004] In a first aspect, the present invention provides a method for reverse estimating inflow flow, comprising: obtaining water level data in front of the dam, a reservoir capacity relationship curve and outflow flow data; performing low-pass filtering on the water level data in front of the dam to generate water level data after low-pass filtering; and reverse estimating the inflow flow based on the water level data after low-pass filtering, the reservoir capacity relationship curve and the outflow flow data.

[0005] In the embodiment of the present invention, by low-pass filtering the water level data in front of the dam, the noise of the water level data in front of the dam is reduced, thereby reducing the influence of the water level fluctuation in front of the dam caused by external interference and the high-frequency fluctuation of the water level on the accuracy of the calculation of the inflow flow. The accuracy of the inflow flow inversely inferred from the water level data, the reservoir capacity relationship curve and the outflow flow data is improved. The problem of low calculation accuracy when using the smoothing method to process the inflow flow in the related art is solved.

[0006] In an optional embodiment, the water level data in front of the dam is low-pass filtered to generate water level data after low-pass filtering, including: using a low-pass filter to low-pass filter the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to generate water level data after low-pass filtering.

[0007] In the embodiment of the present invention, the water level data in front of the dam is denoised by low-pass filtering, thereby reducing the impact of water level fluctuations in front of the dam caused by external interferences such as wind direction, ships, unit power generation, and gate opening and closing.

[0008] In an optional embodiment, a low-pass filter is used to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to generate water level data after low-pass filtering, including: using a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to obtain water level data after a primary filtering; performing a secondary filtering on the water level data after the primary filtering to generate water level data after low-pass filtering.

[0009] In the embodiment of the present invention, the water level data after the primary filtering is subjected to secondary filtering, which solves the problem of phase delay in the water level data after the primary filtering and further improves the accuracy of the inflow flow calculation.

[0010] In an optional embodiment, the water level data after the first filtering is filtered for the second time to generate the water level data after the low-pass filtering, including: flipping the water level data after the first filtering to obtain the flipped water level data; using a low-pass filter to generate the water level data after the second filtering according to the preset sampling frequency, the preset cutoff frequency and the flipped water level data; flipping the water level data after the second filtering to obtain the water level data after the low-pass filtering.

[0011] In an embodiment of the present invention, by flipping the water level data after the first filtering, performing a second filtering on the flipped water level data, and flipping the second filtered water level data again, the purpose of offsetting the phase offset of the water level data is achieved, and the phase lag problem existing in the traditional method is solved.

[0012] In an optional embodiment, after inferring the inflow flow based on the low-pass filtered water level data, the reservoir capacity relationship curve and the outflow flow data, the method further includes: correcting the inflow flow based on a preset minimum flow and a preset lag time step.

[0013] In the embodiment of the present invention, by correcting the inflow flow, the influence of the start and stop of the unit on the water level observation and the inflow flow estimation is further solved, and the accuracy of the inflow flow calculation is improved.

[0014] In an optional implementation, the outflow flow data includes: downstream flow and other water intake.

[0015] In the embodiment of the present invention, the purpose of improving the accuracy of inbound flow reverse estimation is achieved by improving the comprehensiveness of outbound flow data.

[0016] In an optional embodiment, the inflow flow is reversed based on the water level data after low-pass filtering, the reservoir capacity relationship curve and the outflow flow data, including: using the water balance method to calculate the inflow flow based on the sum of the outflow flow, other water withdrawals and the reservoir tolerance, and the reservoir tolerance is determined by the reservoir capacity relationship curve and the water level data after low-pass filtering.

[0017] In the embodiment of the present invention, the water level data after low-pass filtering is used to calculate the inflow flow, thereby achieving the effect of eliminating abnormal jagged edges in the data. Under the premise of ensuring water balance, the accuracy of the inflow flow reverse calculation is improved to the maximum extent.

[0018] In an optional embodiment, the water level data after low-pass filtering is data of multiple hours, the reservoir tolerance is the difference in reservoir capacity corresponding to the water level data after low-pass filtering of any two adjacent hours, and the reservoir capacity corresponding to the water level data after low-pass filtering is determined by the water level data after low-pass filtering and the water level and reservoir capacity curve.

[0019] In the embodiment of the present invention, the water level data after low-pass filtering uses multiple hours of data, and multiple hours of data of the inflow flow can be obtained, which provides a calculation basis for the correction of the inflow flow.

[0020] In an optional embodiment, the incoming flow is corrected according to a preset minimum flow and a preset lag time step, including: obtaining the incoming flow for multiple hours; taking the maximum value of the average of the incoming flow of the current hour and the incoming flow of the previous hour and the preset minimum flow value as the corrected incoming flow for the current hour; taking the difference between the incoming flow of the current hour and the corrected incoming flow of the current hour as the corrected flow difference for the current hour; and dividing the corrected flow difference into the corrected inflows of the preset lag time step after the current hour.

[0021] In the embodiment of the present invention, low-pass filtering has eliminated most of the abnormal values, but there is still a small amount of flow that is affected by the start-up and shutdown of the unit on the water level and cannot be completely solved by low-pass filtering. By correcting the incoming flow, the abnormal flow caused by the start-up and shutdown of the unit is further solved, thereby improving the accuracy of the incoming flow calculation.

[0022] In an optional implementation, the corrected flow difference is divided into preset lag time steps of corrected inflow flows after the current hour, including: dividing the corrected flow difference into preset lag time steps of corrected inflow flows after the current hour in a linear proportion.

[0023] In an embodiment of the present invention, the flow difference caused by the flow correction is linearly allocated to the preset lag time step after the current hour to correct the incoming flow, thereby achieving the effect of avoiding abnormal negative values ​​of the incoming flow and further improving the accuracy of the inward flow reverse calculation.

[0024] In an optional implementation, the preset sampling frequency is 24 Hz and the preset cutoff frequency is 3 Hz.

[0025] In the embodiment of the present invention, preferred values ​​of the preset sampling frequency and the preset cutoff frequency are provided so as to more accurately calculate the inflow flow.

[0026] In an optional implementation, the preset lag time step is 3.

[0027] In the embodiment of the present invention, an optimal value of the preset lag time step is given to more accurately calculate the inflow flow.

[0028] In the second aspect, the present invention provides an inflow flow reverse estimation device, including: an acquisition module, used to obtain the water level data in front of the dam, the reservoir capacity relationship curve and the outflow flow data; a low-pass filtering module, used to perform low-pass filtering on the water level data in front of the dam to generate the water level data after low-pass filtering; an inflow flow reverse estimation module, used to reverse the inflow flow based on the water level data after low-pass filtering, the reservoir capacity relationship curve and the outflow flow data.

[0029] In an optional embodiment, the low-pass filtering module includes: a data generating unit, which is used to use a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to generate water level data after low-pass filtering.

[0030] In an optional embodiment, the data generation unit includes: a primary filtering subunit, which uses a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to obtain water level data after the primary filtering; a secondary filtering subunit, which performs secondary filtering on the water level data after the primary filtering to generate water level data after low-pass filtering.

[0031] In an optional embodiment, the secondary filtering sub-unit includes: a first flipping sub-module, used to flip the water level data after the primary filtering to obtain the flipped water level data; a data generation sub-module, used to generate the water level data after the secondary filtering according to the preset sampling frequency, the preset cutoff frequency and the flipped water level data by using a low-pass filter; a second flipping sub-module, used to flip the water level data after the secondary filtering to obtain the water level data after the low-pass filtering.

[0032] In an optional embodiment, the device further includes: a correction module, configured to correct the incoming flow rate according to a preset minimum flow rate and a preset lag time step.

[0033] In an optional implementation, the outflow flow data includes: downstream flow and other water intake.

[0034] In an optional embodiment, the inflow flow reverse calculation module includes: an inflow flow calculation unit, which is used to calculate the inflow flow according to the sum of the downstream flow, other water withdrawals and reservoir tolerance using a water balance method, and the reservoir tolerance is determined by a reservoir capacity relationship curve and water level data after low-pass filtering.

[0035] In an optional embodiment, the water level data after low-pass filtering is data of multiple hours, the reservoir tolerance is the difference in reservoir capacity corresponding to the water level data after low-pass filtering of any two adjacent hours, and the reservoir capacity corresponding to the water level data after low-pass filtering is determined by the water level data after low-pass filtering and the water level and reservoir capacity curve.

[0036] In an optional embodiment, the correction module includes: an inflow flow acquisition unit, used to obtain the inflow flow of multiple hours; a selection unit, used to take the maximum value of the average value of the inflow flow of the current hour and the inflow flow of the previous hour and the preset minimum flow value as the corrected inflow flow of the current hour; a corrected flow difference calculation unit, used to take the difference between the inflow flow of the current hour and the corrected inflow flow of the current hour as the corrected flow difference of the current hour; and a division unit, used to divide the corrected flow difference into preset lag time steps after the current hour.

[0037] In an optional implementation, the division unit includes: a linear division subunit, which is used to divide the corrected flow difference into preset lag time steps of corrected inflow flows after the current hour in a linear proportion.

[0038] In an optional implementation, the preset sampling frequency is 24 Hz and the preset cutoff frequency is 3 Hz.

[0039] In an optional implementation, the preset lag time step is 3.

[0040] In a third aspect, the present invention provides a computer device, comprising: a memory and a processor, the memory and the processor are communicatively connected to each other, the memory stores computer instructions, and the processor executes the inflow flow reverse estimation method of the above-mentioned first aspect or any corresponding embodiment thereof by executing the computer instructions.

[0041] 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 inflow flow reverse estimation method of the first aspect or any corresponding embodiment thereof. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0043] Figure 1 is a schematic diagram of a flow chart of an inbound flow reverse estimation method according to an embodiment of the present invention;

[0044] Figure 2 is a flow chart of another method for inverse estimation of inflow flow according to an embodiment of the present invention;

[0045] Figure 3 is a comparison diagram of the original water level and the secondary filtered water level data according to an embodiment of the present invention;

[0046] Figure 4 is a comparison diagram of the original inflow flow and the inflow flow after secondary filtering of the water level according to an embodiment of the present invention;

[0047] Figure 5 is a comparison diagram of the inflow flow before and after correction according to an embodiment of the present invention;

[0048] Figure 6 is a structural block diagram of an inflow flow reverse estimation device according to an embodiment of the present invention;

[0049] Figure 7 It is a schematic diagram of the hardware structure of a computer device according to an embodiment of the present invention. DETAILED DESCRIPTION

[0050] In order to make the purpose, technical solution and advantages of the embodiments of the present invention clearer, the technical solution in the embodiments of the present invention will be clearly and completely described below in conjunction with the drawings in the embodiments of the present invention. Obviously, the described embodiments are part of the embodiments of the present invention, not all of the embodiments. Based on the embodiments of the present invention, all other embodiments obtained by those skilled in the art without creative work are within the scope of protection of the present invention.

[0051] It should be noted that, in the description of the present invention, the terms "first", "second", etc. are used to distinguish similar objects, and are not necessarily used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged where appropriate, so that the embodiments of the present invention described herein can be implemented in an order other than those illustrated or described herein. In addition, the terms "including" and "having" and any of their variations are intended to cover non-exclusive inclusions, for example, a process, method, system, product or device that includes a series of steps or units is not necessarily limited to those steps or units clearly listed, but may include other steps or units that are not clearly listed or inherent to these processes, methods, products or devices. The terms "installed", "connected", and "connected" should be understood in a broad sense, for example, it can be a fixed connection, a detachable connection, or an integral connection; it can be a mechanical connection or an electrical connection; it can be directly connected, it can be indirectly connected through an intermediate medium, it can also be the internal communication of two elements, it can be a wireless connection, or it can be a wired connection. For those of ordinary skill in the art, the specific meanings of the above terms in the present invention can be understood in specific circumstances.

[0052] Most of the inflow of reservoirs is obtained by reverse calculation from the water level in front of the dam and the outflow using the principle of water balance. However, the calculated flow process generally has obvious sawtooth fluctuations, and even negative values ​​appear. This has a serious impact on a series of work such as hydrological analysis, flood control and power generation scheduling research. At present, this problem is usually dealt with by smoothing the calculation results of the inflow flow to eliminate abnormal sawtooth in the data and smooth the flow process. Common correction methods include sliding average method and five-point cubic method. Sliding average is a commonly used time series smoothing method. Its basic principle is to take the average value of the data within a period of time as the value at the current moment, then slide the time window backward by one unit, continue to calculate the average value, and so on. Different sliding average methods can use different weights, such as simple sliding average, exponential sliding average, weighted sliding average, etc. The difference between these methods lies in the different weights used when calculating the average value. The five-point cubic method is also a data smoothing algorithm. It uses a window containing five data points for moving smoothing. The basic principle is to use a cubic polynomial function to fit the data points in the window to obtain a smoothed data point.

[0053] These two methods are widely used because of their simplicity, but the inflow flow obtained by these two methods is not very accurate. The calculation of the sliding average is based on the inflow flow data of the past period of time, and these data cannot accurately reflect the current inflow flow situation. Therefore, the smoothed flow data will lag, causing the peak time to deviate. And the degree of smoothing of the sliding average is affected by the window size. If the selected window is too large, the flow process will be over-flattened, causing the peak flow of the flood process to deviate. If the window is too small, the smoothing effect will be poor. The five-point cubic smoothing method smoothes the data of five points through cubic polynomial fitting, so there may be distortion after smoothing for very steep inflow flow curves. In addition, since the two methods smooth the inflow flow process, the water balance of the corrected value is not considered during the processing process, which will lead to large errors in the calculated water volume of the corrected value and the original value, resulting in water imbalance.

[0054] According to an embodiment of the present invention, a method for inbound flow inference 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.

[0055] In this embodiment, a method for inferring inbound traffic is provided, which can be used in the above-mentioned mobile terminal, such as a central processing unit or a server, Figure 1 is a flow chart of the method for inverse estimation of inbound flow according to an embodiment of the present invention. Figure 1 As shown, the process includes the following steps:

[0056] Step S101, obtaining the water level data in front of the dam, the reservoir capacity relationship curve and the outflow flow data.

[0057] Optionally, the water level in front of the dam is the water level upstream of the dam, and the water level in front of the dam can be water level data at multiple times, such as water level data for multiple days, multiple hours, etc. The reservoir capacity relationship curve can be expressed as: V = f (y), where y is the water level data in front of the dam (water level data), V is the water level reservoir capacity (reservoir capacity), and f is the function of V with respect to y. The outflow flow data is the flow data flowing out of the reservoir.

[0058] Step S102, low-pass filtering is performed on the water level data in front of the dam to generate water level data after low-pass filtering.

[0059] Optionally, the water level data in front of the dam obtained in step S101 is low-pass filtered, that is, low-frequency signals are allowed to pass normally, and high-frequency signals exceeding a set critical value are blocked and weakened, and finally the water level data after low-pass filtering is obtained.

[0060] Step S103, reversely inferring the inflow flow based on the water level data after low-pass filtering, the reservoir capacity relationship curve and the outflow flow data.

[0061] Optionally, the low-pass filtered water level data obtained in step S102 is used to calculate (reverse-calculate) the inflow flow rate based on the reservoir capacity relationship curve and outflow flow data obtained in step S101 using the water balance principle.

[0062] In the embodiment of the present invention, by low-pass filtering the water level data in front of the dam, the noise of the water level data in front of the dam is reduced, thereby reducing the influence of the water level fluctuation in front of the dam caused by external interference and the high-frequency fluctuation of the water level on the accuracy of the calculation of the inflow flow. The accuracy of the inflow flow inversely inferred from the water level data, the reservoir capacity relationship curve and the outflow flow data is improved. The problem of low calculation accuracy when using the smoothing method to process the inflow flow in the related art is solved.

[0063] In this embodiment, a method for inferring inbound traffic is provided, which can be used in the above-mentioned mobile terminal, such as a central processing unit or a server, Figure 2 FIG. 1 is a flow chart of another method for inferring inflow flow according to an embodiment of the present invention. Figure 2 As shown, the process includes the following steps:

[0064] Step S201, obtain the water level data before the dam, the reservoir capacity relationship curve and the outflow flow data. Figure 1 Step S101 of the illustrated embodiment will not be described in detail here.

[0065] Step S202: low-pass filter the water level data in front of the dam according to the preset sampling frequency and the preset cut-off frequency to generate water level data after low-pass filtering. s ,f c ,k) is the transfer function of the low-pass filter, where f s is the sampling frequency, f c is the cut-off frequency, and k is the order of the low-pass filter. It should be noted that the higher the order of the low-pass filter, the better the filtering effect and the higher the cost, because the higher the order, the more complex the structure of the low-pass filter. Therefore, the order k and f of the low-pass filter are s 、f c As a preferred implementation, in the embodiment of the present invention, the order of the preset low-pass filter is 3, the preset sampling frequency is 24 Hz, and the preset cutoff frequency is 3 Hz, that is, k=3, f s =24hz, f c =3hz. Specifically, the above step S202 includes:

[0066] Step S2021: low-pass filter the water level data in front of the dam according to the preset sampling frequency and the preset cutoff frequency to obtain the water level data after the primary filtering. Optionally, the water level data in front of the dam is filtered once to eliminate the high-frequency oscillation effect of external interference such as wind and waves on the water level. The water level data after the primary filtering The following formula can be used for calculation:

[0067]

[0068] Where, Y is the water level data in front of the dam, H(f s ,f c ,k) is the transfer function of the k-order low-pass filter, which is used to describe the attenuation degree of the low-pass filter for signals of different frequencies. When the frequency of the water level data signal in front of the dam is greater than the cut-off frequency f c When , the signal passes through the filter with a greater attenuation and the amplitude is reduced.

[0069] Step S2022, performing secondary filtering on the water level data after the primary filtering to generate water level data after low-pass filtering. Optionally, the water level data after the primary filtering is filtered again using a low-pass filter to obtain water level data after low-pass filtering. In an embodiment of the present invention, performing secondary filtering on the water level data after the primary filtering solves the problem of phase delay in the water level data after the primary filtering, and further improves the accuracy of the calculation of the inflow flow.

[0070] In some optional implementations, the above step S2022 includes:

[0071] Step a1, flip the filtered water level data to obtain flipped water level data. Specifically, the filtered water level data obtained in step S2021 is flipped Flip it to get the water level data after flipping

[0072] Step a2, using a low-pass filter to generate secondary filtered water level data according to a preset sampling frequency, a preset cutoff frequency and the flipped water level data. Specifically, the preset sampling frequency and the preset cutoff frequency determined in step S2021 are used to filter the flipped water level data. Perform secondary filtering to offset the phase offset and obtain the water level data after secondary filtering. The calculation formula is as follows:

[0073]

[0074] Step a3, flip the water level data after secondary filtering to obtain the water level data after low-pass filtering. Specifically, the water level data after secondary filtering obtained in step a2 is Flip and restore to obtain the low-pass filtered water level data (zero phase delay dam front water level low-pass filtered data) Y lp .

[0075] In the embodiment of the present invention, the purpose of offsetting the phase offset of the water level data is achieved by flipping the water level data after the primary filtering, performing secondary filtering on the flipped water level data, and flipping the water level data after the secondary filtering again, thereby solving the phase lag problem existing in the traditional method. Original water level data (water level data before the dam) Y and secondary filtered water level data (water level data after low-pass filtering) Y lp For example Figure 3 As shown in the figure, it can be seen that the original water level data Y has more sawtooth, and the secondary filtered water level data Y lp There is no obvious sawtooth, thus achieving the effect of improving the accuracy of reverse estimation of inflow flow.

[0076] Step S203, inversely infer the inflow flow based on the low-pass filtered water level data, the reservoir capacity relationship curve and the outflow flow data. Optionally, the outflow flow data includes: outflow flow and other water withdrawals, and other water withdrawals may include evaporation, emanation, infiltration and irrigation water withdrawal.

[0077] In some optional embodiments, step S203 includes: using a water balance method to calculate the inflow according to the sum of the outflow, other water withdrawals and the reservoir tolerance, and the reservoir tolerance is determined by the reservoir capacity relationship curve and the water level data after low-pass filtering. In some optional embodiments, the water level data after low-pass filtering is data for multiple hours, and the reservoir tolerance is the difference in the reservoir capacity corresponding to the low-pass filtered water level data of any two adjacent hours, and the reservoir capacity corresponding to the low-pass filtered water level data is determined by the low-pass filtered water level data and the water level and reservoir capacity curve. It can be understood that when the water level data Y before the dam: {y1, y2...y n} is hourly series data (subscript n represents the nth hour), the water level data after one filtering Water level data after flipping Water level data after secondary filtering And the water level data Y after low-pass filtering lp : The data are for multiple hours.

[0078] Using the water balance method (the difference between the water input and the water output in the system is the water storage capacity in the system), the inflow flow QI is calculated according to the following formula:

[0079] QI=QO+QW+ΔV(Y lp )

[0080] Where, QI: {qi1, qi2...qin} is the inflow flow, QO: {qo1, qo2...qo n} is the discharge flow, QW: {qw1, qw2...qw n} is other water consumption, ΔV(Y lp ) is the water level data Y after filtering lp The reservoir tolerance is the water level data after low-pass filtering in any two adjacent hours. Substitute into the storage capacity curve V=f(y) and obtain the difference in storage capacity.

[0081] Step S204, correct the inflow flow according to the preset minimum flow and the preset lag time step. Optionally, after step S203, the effect of reducing sawtooth and eliminating most abnormal values ​​has been achieved. Figure 4 As shown in the figure, it can be seen that the calculated inflow QI (inflow after secondary filtering of water level) after two low-pass filtering has significantly less sawtooth and abnormal values ​​(flow less than 0) than the original inflow Q. However, there is still a small amount of flow that cannot be completely resolved by low-pass filtering due to the influence of the start and stop of the unit on the water level, and there are abnormal negative values. Step S204 is based on the preset minimum flow q min The preset lag time step m is used to correct the inflow flow QI, and the corrected inflow flow QIC is obtained. min It can be set according to the actual operation of the reservoir. In the embodiment of the present invention, the minimum flow rate q is preset min =3m 3 / s, the preset lag time step m can be selected according to the calculation requirements. As a preferred implementation, the preset lag time step m=3.

[0082] In the embodiment of the present invention, by correcting the inflow flow, the influence of the start and stop of the unit on the water level observation and the inflow flow calculation is further solved, and the accuracy of the inflow flow calculation is improved. Figure 5 As shown, the secondary corrected inflow QIC (corrected inflow) solves the anomaly of flow less than 0 compared with the inflow QI (inflow before correction) after secondary filtering of the water level.

[0083] In some optional implementations, the above step S204 includes:

[0084] Step b1, obtain the inflow flow for multiple hours. Obtain the inflow flow QI for multiple hours calculated in step S203: {qi1, qi2...qi n}.

[0085] Step b2: The maximum value among the average value of the current hour's inflow flow, the previous hour's inflow flow, and the preset minimum flow value is used as the corrected inflow flow for the current hour. t Less than the preset minimum flow rate q min When the inflow flow for that hour is t The corrected inflow flow rate for any hour is t Calculated according to the following formula:

[0086] qq t =max{(qi t +qi t-1 ) / 2,q min}

[0087] In the formula, qi t is the inbound flow of the current hour, qi t-1 is the inflow flow in the previous hour, q min To preset the minimum flow rate, t It is the corrected inbound flow rate of the current hour, that is, the corrected inbound flow rate.

[0088] Step b3, taking the difference between the current hour's inflow and the corrected inflow of the current hour as the corrected flow difference of the current hour. Optionally, the flow correction will result in a flow difference, that is, the corrected flow difference is the difference between the corrected inflow of the current hour and the inflow of the current hour (qi t -qic t ).

[0089] Step b4, divide the corrected flow difference into the preset lag time step corrected inflow flow after the current hour. Optionally, divide the corrected flow difference obtained in step b3 into the preset lag time step corrected inflow flow after the current hour, thereby reducing the impact of flow correction and further improving the accuracy of inflow flow calculation.

[0090] In some optional implementations, the above step b4 includes: dividing the corrected flow difference into the preset lag time step corrected inflow flow after the current hour in linear proportion. Specifically, the corrected flow difference can be divided using the following formula:

[0091]

[0092] In the formula, qic t+i is the corrected inflow flow in any hour after the current hour, qi t+i is the inflow flow in any hour after the current hour, qi t is the inflow flow of the current hour, qic tis the corrected inflow flow of the current hour, (qi t -qic t ) is the corrected flow difference of the current hour, and m is the preset lag time step. Thus, the corrected inflow flow QIC can be obtained: {qic1, qic2...qic n}.

[0093] The inflow reverse estimation method provided in this embodiment performs low-pass filtering and noise reduction on the water level data in front of the dam to reduce the influence of water level fluctuations in front of the dam caused by external interference such as wind direction, ships, unit power generation, gate opening and closing, and solves the influence of high-frequency fluctuations in the water level that has the greatest impact on the inflow from the source. The sawtooth of the inflow process calculated by the water level data after low-pass filtering is significantly reduced, and the inflow process is restored to the maximum extent under the premise of ensuring water balance. At the same time, by correcting the inflow, the influence of the opening and closing of the unit on water level observation and inflow calculation is further solved. Compared with the existing methods, this method overcomes the drawbacks of the traditional method of calculating the inflow, such as hysteresis and water imbalance. It achieves the avoidance of phase lag caused by traditional methods, and can significantly improve the sawtooth and negative value problems of the inflow, and the method is simpler and more general.

[0094] In this embodiment, a device for inflow flow reverse prediction is also provided, which is used to implement the above-mentioned embodiments and preferred implementation modes, and will not be repeated hereafter. As used below, the term "module" can be a combination of software and / or hardware that implements a predetermined function. Although the device described in the following embodiments is preferably implemented in software, the implementation of hardware, or a combination of software and hardware, is also possible and conceivable.

[0095] This embodiment provides an inflow flow reverse estimation device, such as Figure 6 As shown, it includes: an acquisition module 601, which is used to obtain the water level data in front of the dam, the reservoir capacity relationship curve and the outflow flow data; a low-pass filtering module 602, which is used to perform low-pass filtering on the water level data in front of the dam to generate the water level data after low-pass filtering; an inflow flow reverse estimation module 603, which is used to reverse the inflow flow according to the water level data after low-pass filtering, the reservoir capacity relationship curve and the outflow flow data.

[0096] Through the above module, the water level data in front of the dam is low-pass filtered to achieve noise reduction of the water level data in front of the dam, thereby reducing the influence of the water level fluctuation in front of the dam caused by external interference and the high-frequency fluctuation of the water level on the accuracy of the calculation of the inflow flow. The purpose of smoothing the water level data in front of the dam is achieved, and the accuracy of the inflow flow inversely deduced based on the water level data, the reservoir capacity relationship curve and the outflow flow data is improved. The problem of low calculation accuracy when using the smoothing method to process the inflow flow in the related technology is solved.

[0097] In an optional embodiment, the low-pass filtering module includes: a data generating unit, which is used to use a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to generate water level data after low-pass filtering.

[0098] In an optional embodiment, the data generation unit includes: a primary filtering subunit, which uses a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to obtain water level data after the primary filtering; a secondary filtering subunit, which performs secondary filtering on the water level data after the primary filtering to generate water level data after low-pass filtering.

[0099] In an optional embodiment, the secondary filtering sub-unit includes: a first flipping sub-module, used to flip the water level data after the primary filtering to obtain the flipped water level data; a data generation sub-module, used to generate the water level data after the secondary filtering according to the preset sampling frequency, the preset cutoff frequency and the flipped water level data by using a low-pass filter; a second flipping sub-module, used to flip the water level data after the secondary filtering to obtain the water level data after the low-pass filtering.

[0100] In an optional embodiment, the device further includes: a correction module, configured to correct the incoming flow rate according to a preset minimum flow rate and a preset lag time step.

[0101] In an optional implementation, the outflow flow data includes: downstream flow and other water intake.

[0102] In an optional embodiment, the inflow flow reverse calculation module includes: an inflow flow calculation unit, which is used to calculate the inflow flow according to the sum of the downstream flow, other water withdrawals and reservoir tolerance using a water balance method, and the reservoir tolerance is determined by a reservoir capacity relationship curve and water level data after low-pass filtering.

[0103] In an optional embodiment, the water level data after low-pass filtering is data of multiple hours, the reservoir tolerance is the difference in reservoir capacity corresponding to the water level data after low-pass filtering of any two adjacent hours, and the reservoir capacity corresponding to the water level data after low-pass filtering is determined by the water level data after low-pass filtering and the water level and reservoir capacity curve.

[0104] In an optional embodiment, the correction module includes: an inflow flow acquisition unit, used to obtain the inflow flow of multiple hours; a selection unit, used to take the maximum value of the average value of the inflow flow of the current hour and the inflow flow of the previous hour and the preset minimum flow value as the corrected inflow flow of the current hour; a corrected flow difference calculation unit, used to take the difference between the inflow flow of the current hour and the corrected inflow flow of the current hour as the corrected flow difference of the current hour; and a division unit, used to divide the corrected flow difference into preset lag time steps after the current hour.

[0105] In an optional implementation, the division unit includes: a linear division subunit, which is used to divide the corrected flow difference into preset lag time steps of corrected inflow flows after the current hour in a linear proportion.

[0106] In an optional implementation, the preset sampling frequency is 24 Hz and the preset cutoff frequency is 3 Hz.

[0107] In an optional implementation, the preset lag time step is 3.

[0108] 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.

[0109] The inflow flow reverse estimation device 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.

[0110] The embodiment of the present invention also provides a computer device having the above Figure 6 The inflow flow reverser is shown.

[0111] See also Figure 7 , Figure 7 is a schematic diagram of the structure of a computer device provided by an optional embodiment of the present invention, such as Figure 7As shown, the computer device includes: one or more processors 10, a memory 20, and interfaces for connecting various components, including high-speed interfaces and low-speed interfaces. Various components are connected to each other using different buses for communication, 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 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. Similarly, multiple computer devices can be connected, and each device provides some necessary operations (for example, as a server array, a group of blade servers, or a multi-processor system). Figure 7 A processor 10 is taken as an example.

[0112] 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 a dedicated 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.

[0113] The memory 20 stores instructions executable by at least one processor 10, so that at least one processor 10 executes the method shown in the above embodiment.

[0114] The memory 20 may include a program storage area and a data storage area, wherein the program storage area may store an operating system, an application required for at least one function; the data storage area may store data created according to 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 arranged 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.

[0115] 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.

[0116] The computer device also 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 7 The example of connecting through bus is taken in the following.

[0117] The input device 30 can receive input digital or character information, and generate key signal input related to the user settings and function control of the computer device, such as a touch screen, a keypad, a mouse, a track pad, a touch pad, an indicator bar, one or more mouse buttons, a trackball, a joystick, etc. The output device 40 may include a display device, an auxiliary lighting device (e.g., an LED) and a tactile feedback device (e.g., a vibration motor), etc. 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.

[0118] The embodiment of the present invention also provides a computer-readable storage medium. The method according to the embodiment of the present invention can be implemented in hardware, firmware, or can be implemented as a computer code that can be recorded in a storage medium, or can be implemented as a computer code that is originally stored in a remote storage medium or a non-temporary machine-readable storage medium and will be stored in a local storage medium through a network download, so that the method described herein can be stored in such software processing on a storage medium using a general-purpose computer, a dedicated processor, or programmable or dedicated hardware. Among them, the storage medium can be a magnetic disk, an optical disk, a read-only storage memory, a random access memory, a flash memory, a hard disk or a solid-state hard disk, etc.; further, the storage medium can also include a combination of the above types of memories. It can be understood that a computer, a processor, a microprocessor controller, or 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 a computer, a processor, or hardware, the method shown in the above embodiment is implemented.

[0119] Although the embodiments of the present invention have been described in conjunction with the accompanying drawings, those skilled in the art may make various modifications and variations without departing from the spirit and scope of the present invention, and such modifications and variations are all within the scope defined by the appended claims.

Claims

1. A method for inverse estimation of inbound flow, characterized in that: The method comprises: Obtain dam front water level data, reservoir capacity relationship curve and outflow flow data; Performing low-pass filtering on the water level data in front of the dam to generate low-pass filtered water level data, thereby eliminating the high-frequency oscillation effect of wind and waves on the water level; Reversely estimate the inflow flow rate based on the water level data after low-pass filtering, the reservoir capacity relationship curve and the outflow flow data; The inflow flow is corrected according to the preset minimum flow and the preset lag time step, which solves the problem that the start and stop of the unit affect the water level observation and the inflow flow calculation; The step of performing low-pass filtering on the water level data in front of the dam to generate water level data after low-pass filtering comprises: Using a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to obtain the water level data after primary filtering; The water level data after the primary filtering is subjected to secondary filtering to generate water level data after low-pass filtering.

2. The method for inward flow estimation according to claim 1, characterized in that: The performing secondary filtering on the water level data after the primary filtering to generate water level data after low-pass filtering comprises: Flipping the filtered water level data to obtain flipped water level data; Using a low-pass filter to generate secondary filtered water level data according to a preset sampling frequency, a preset cutoff frequency and the flipped water level data; The water level data after the secondary filtering is flipped to obtain the water level data after the low-pass filtering.

3. The method for inward flow estimation according to claim 1, characterized in that: The outflow flow data includes: downstream flow and other water intake.

4. The method for inward flow estimation according to claim 3 is characterized in that: The reverse calculation of the inflow flow according to the low-pass filtered water level data, the reservoir capacity relationship curve and the outflow flow data includes: The water balance method is used to calculate the inflow flow based on the sum of the outflow flow, other water withdrawals and the reservoir tolerance, and the reservoir tolerance is determined by the reservoir capacity relationship curve and the water level data after low-pass filtering.

5. The method for inward flow estimation according to claim 4, characterized in that: The water level data after low-pass filtering is data of multiple hours, and the reservoir tolerance is the difference in reservoir capacity corresponding to the water level data after low-pass filtering of any two adjacent hours. The reservoir capacity corresponding to the water level data after low-pass filtering is determined by the water level data after low-pass filtering and the water level and reservoir capacity curve.

6. The method for inward flow estimation according to claim 5, characterized in that: The method of correcting the inflow flow according to the preset minimum flow and the preset lag time step includes: Get multiple hours of inbound traffic; The maximum value among the average value of the inflow flow rate of the current hour, the inflow flow rate of the previous hour and the preset minimum flow value is taken as the corrected inflow flow rate of the current hour; The difference between the inflow flow rate of the current hour and the corrected inflow flow rate of the current hour is used as the corrected flow difference of the current hour; The corrected flow difference is divided into preset lag time steps after the current hour to correct the incoming flow.

7. The method for inward flow estimation according to claim 6, characterized in that: The step of dividing the corrected flow difference into a preset lag time step number of corrected inflow flows after the current hour includes: The corrected flow difference is divided into the preset lag time step corrected inflow flows after the current hour in linear proportion.

8. The method for inflow flow inversion according to claim 1, characterized in that: The preset sampling frequency is 24 Hz, and the preset cutoff frequency is 3 Hz.

9. The method for inward flow estimation according to claim 1, characterized in that: The preset lag time step is 3.

10. A device for inflow flow reverse estimation, characterized in that: The device comprises: The acquisition module is used to obtain the water level data in front of the dam, the reservoir capacity relationship curve and the outflow flow data; A low-pass filtering module is used to perform low-pass filtering on the water level data in front of the dam to generate water level data after low-pass filtering, thereby eliminating the high-frequency oscillation effect of wind and waves on the water level; An inflow reverse estimation module, used for reversely estimating the inflow according to the water level data after the low-pass filtering, the reservoir capacity relationship curve and the outflow data; The inflow flow is corrected according to the preset minimum flow and the preset lag time step, which solves the problem that the start and stop of the unit affect the water level observation and the inflow flow calculation; The step of performing low-pass filtering on the water level data in front of the dam to generate water level data after low-pass filtering comprises: Using a low-pass filter to perform low-pass filtering on the water level data in front of the dam according to a preset sampling frequency and a preset cutoff frequency to obtain the water level data after primary filtering; The water level data after the primary filtering is subjected to secondary filtering to generate water level data after low-pass filtering.

11. 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 inflow flow reverse estimation method according to any one of claims 1 to 9 by executing the computer instructions.

12. 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 inflow flow reverse estimation method according to any one of claims 1 to 9.

Images