A rapid calculation method for sediment content in river sections

Abstract

The application discloses a method for quickly measuring the sediment content of a river section, and belongs to the technical field of water conservancy engineering. The method comprises: obtaining a historically measured water level Z1, flow rate Q and sediment content S data series of a target river section; determining different preset water levels Z2 and their associated water depths H and water-passing area A data according to the topographic data of the target river section to generate a Z2~A database and an H=F(Z) relationship; determining the section average flow velocity U data respectively associated with each Z1 data according to the Z1 data series, the Q data series and the Z2~A database to generate a Z1~U database; calibrating the undetermined coefficients in the to-be-solved relationship between U and Z through the Z1~U database to obtain the U=G(Z) relationship; calibrating the undetermined coefficients in the to-be-solved relationship between S and Z through the historically measured Z1 data series and the S data series to obtain the S=W(Z) relationship; substituting the measured water level Z0 data collected at any time into the S=W(Z) relationship to obtain the sediment content S0 data associated with the Z0 data.

Description

Technical Field

[0001] The present invention belongs to the technical field of water conservancy engineering, and specifically relates to a method for quickly calculating the sediment content of a river section. Background Art

[0002] Sediment content, that is, the volume of sediment carried by a unit of water, is related to the evolution and development of river channels and the design, construction and operation management of water-related projects. It is a key focus of hydrological monitoring in the field of water conservancy engineering technology. The sediment content of the river channel has a certain spatial variability along the depth and width of the river, so it is necessary to determine a sediment content value that can represent the overall level of the section, that is, the cross-section sediment content.

[0003] At present, the measurement methods of cross-section sediment content mainly include: direct measurement method, collecting water samples at different locations of the cross-section and measuring the sediment content, and taking the average value of the relevant results as the cross-section sediment content. Its limitations are reflected in the long cycle, high cost, and the measurement accuracy is affected by personnel operation and equipment; it is difficult to implement under extreme weather conditions and there are safety hazards. The analytical model method combines the sediment content measurement data at the designated location on site with the theoretical model to calculate the cross-section sediment content, which also faces the problems of difficulty in measurement and safety under extreme weather. The statistical fitting method establishes a relationship model between the cross-section sediment content and hydrological factors (such as flow) through historical data fitting. Its limitations are reflected in the over-reliance on historical data, and the applicability and accuracy of the model are limited by the richness of historical data.

[0004] Therefore, the current method for measuring sediment content in river sections has certain shortcomings in terms of measurement speed or accuracy. Summary of the invention

[0005] The purpose of the embodiments of the present application is to provide a method for quickly calculating the sediment content of a river section, which can solve the problem that the current method for calculating the sediment content of a river section is insufficient in terms of calculation speed or calculation accuracy.

[0006] In a first aspect, an embodiment of the present application provides a method for quickly calculating the sediment content of a river section, the method comprising:

[0007] Obtain the topographic data of the target river section and the data sequence of the historically measured water level Z1, cross-section flow Q and cross-section sediment content S, and set several preset water level Z2 data;

[0008] According to the topographic data of the target river section, the water depth H data and the water area A data associated with different preset water level Z2 data are determined to generate a Z2~A database and a H=F(Z) relationship;

[0009] According to the Z1 data sequence, the Q data sequence and the Z2~A database, the cross-sectional average flow velocity U data respectively associated with each Z1 data is determined to generate the Z1~U database;

[0010] According to the river dynamics analytical formula U=C(H) and H=F(Z), the undetermined relationship between U and Z is established, and the undetermined coefficients in the undetermined relationship between U and Z are calibrated through the Z1~U database to obtain the relationship U=G(Z);

[0011] According to the analytical formula of water-sediment dynamics S=D(U,H) and the relationship between H=F(Z) and U=G(Z), the relationship between S and Z to be analyzed is established, and the coefficients to be determined in the relationship between S and Z to be analyzed are calibrated through the historical measured Z1 data series and S data series to obtain the relationship S=W(Z);

[0012] For the measured water level Z0 data collected at any time, substitute it into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data.

[0013] In a second aspect, an embodiment of the present application provides a device for quickly measuring the sediment content of a river section, the device comprising:

[0014] The acquisition module is used to obtain the topographic data of the target river section and the data sequence of the historically measured water level Z1, cross-section flow Q and cross-section sediment content S, and set a number of preset water level Z2 data;

[0015] The first determination module is used to determine the water depth H data and the water area A data associated with different preset water level Z2 data according to the topographic data of the target river section, so as to generate a Z2~A database and a H=F(Z) relationship;

[0016] The second determination module is used to determine the cross-sectional average flow velocity U data respectively associated with each Z1 data according to the Z1 data sequence, the Q data sequence and the Z2~A database to generate the Z1~U database;

[0017] The first analytical module is used to establish the relationship between U and Z according to the river dynamics analytical formula U=C(H) and H=F(Z), and calibrate the undetermined coefficients in the relationship between U and Z through the Z1~U database to obtain the relationship U=G(Z);

[0018] The second analytical module is used to establish the undetermined relationship between S and Z according to the analytical formula of water-sediment dynamics S=D(U,H) and the relationship between H=F(Z) and U=G(Z), and calibrate the undetermined coefficients in the undetermined relationship between S and Z through the historical measured Z1 data sequence and S data sequence to obtain the relationship S=W(Z);

[0019] The third determination module is used to substitute the measured water level Z0 data collected at any time into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data.

[0020] In a third aspect, an embodiment of the present application provides an electronic device, which includes a processor and a memory, wherein the memory stores programs or instructions that can be run on the processor, and when the program or instructions are executed by the processor, the steps of the method described in the first aspect are implemented.

[0021] In a fourth aspect, an embodiment of the present application provides a readable storage medium, on which a program or instruction is stored, and when the program or instruction is executed by a processor, the steps of the method described in the first aspect are implemented.

[0022] In a fifth aspect, an embodiment of the present application provides a chip, comprising a processor and a communication interface, wherein the communication interface is coupled to the processor, and the processor is used to run a program or instruction to implement the method described in the first aspect.

[0023] In a sixth aspect, an embodiment of the present application provides a computer program product, which is stored in a storage medium and is executed by at least one processor to implement the method described in the first aspect.

[0024] In an embodiment of the present application, by acquiring the topographic data of the target river section and the data sequences of the historically measured water level Z1, cross-sectional flow Q and cross-sectional sediment content S, and setting a number of preset water level Z2 data; according to the topographic data of the target river section, the water depth H data and the water flow area A data associated with different preset water level Z2 data are determined to generate a Z2~A database and the relationship H=F(Z); according to the Z1 data sequence, the Q data sequence and the Z2~A database, the cross-sectional average flow velocity U data respectively associated with each Z1 data is determined to generate a Z1~U database; the topographic data, flow data and water level information are effectively integrated to improve data utilization. According to the river dynamics analytical formula U=C(H) and H=F(Z), the relationship between U and Z is established, and the coefficients to be determined in the relationship between U and Z are calibrated through the Z1~U database to obtain the relationship between U=G(Z); according to the water-sediment dynamics analytical formula S=D(U,H) and the relationship between H=F(Z) and U=G(Z), the relationship between S and Z is established, and the coefficients to be determined in the relationship between S and Z are calibrated through the historical measured Z1 data sequence and S data sequence to obtain the relationship between S=W(Z); using historical measured data to analyze and calibrate the formula based on river dynamics and water-sediment dynamics can improve the accuracy of prediction and ensure the effectiveness of the analyzed relationship in actual scenarios. For the measured water level Z0 data collected at any time, substitute it into the S=W(Z) relationship to obtain the sediment content S0 data associated with the Z0 data. The real-time updated measured water level Z0 data can quickly reflect the changes in the river channel sediment content S0 data, which is convenient for reducing the time and cost of measurement and responding to emergencies. The embodiments of the present application, by gradually establishing the relationship between the water level and related hydrological characteristics, achieve efficient use of data and intelligent calculation, reduce dependence on a large amount of on-site monitoring and experimental testing work, and can improve the calculation speed and accuracy of the sediment content of the river section. BRIEF DESCRIPTION OF THE DRAWINGS

[0025] Figure 1 It is a flow chart of a method for quickly calculating the sediment content of a river section provided in an embodiment of the present application;

[0026] Figure 2 is a schematic diagram of a cross section of a target river provided in an embodiment of the present application;

[0027] Figure 3 It is a schematic diagram of the Z2~A database provided in the embodiment of the present application;

[0028] Figure 4 It is a schematic diagram of the Z2~H database and the H= F(Z) relationship provided in the embodiment of the present application;

[0029] Figure 5It is a schematic diagram of the Z1~U database and the U=G(Z) relationship provided in the embodiment of the present application;

[0030] Figure 6 is a comparative schematic diagram of the sand content measurement results provided in the embodiments of the present application;

[0031] Figure 7 It is a structural diagram of a device for quickly measuring the sediment content of a river section provided in an embodiment of the present application;

[0032] Figure 8 It is a schematic diagram of the hardware structure of the electronic device according to an embodiment of the present application. DETAILED DESCRIPTION

[0033] The following will be combined with the drawings of the embodiments of the present application to clearly describe the technical solutions of the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all of the embodiments. All other embodiments obtained by ordinary technicians in this field based on the embodiments of the present application belong to the scope of protection of this application.

[0034] The terms "first", "second", etc. in the specification and claims of this application are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the data used in this way can be interchangeable under appropriate circumstances, so that the embodiments of the present application can be implemented in an order other than those illustrated or described here, and the objects distinguished by "first", "second", etc. are generally of one type, and the number of objects is not limited. For example, the first object can be one or more. In addition, "and / or" in the specification and claims represents at least one of the connected objects, and the character " / " generally indicates that the objects associated with each other are in an "or" relationship.

[0035] In response to the problems arising from the related technologies, the embodiments of the present application provide a method for quickly calculating the sediment content of a river section, which can solve the problem of insufficient speed or accuracy in calculating the sediment content of a river section in the related technologies.

[0036] The following is a detailed description of the method for quickly calculating the sediment content in a river section provided by the embodiment of the present application through specific embodiments and their application scenarios in conjunction with the accompanying drawings.

[0037] Figure 1 A flow chart of a method for rapidly calculating the sediment content in a river section provided in an embodiment of the present application.

[0038] like Figure 1 As shown, the method for quickly calculating the sediment content of a river section may include steps 110 to 160, and the method is applied to a device for quickly calculating the sediment content of a river section, as shown below:

[0039] Step 110, obtaining the topographic data of the target river section and the data sequence of the historically measured water level Z1, the cross-sectional flow Q and the cross-sectional sediment content S, and setting a number of preset water level Z2 data;

[0040] Step 120, determining the water depth H data and the water area A data associated with different preset water level Z2 data according to the topographic data of the target river section, so as to generate a Z2~A database and a relational expression H=F(Z);

[0041] Step 130, determining the cross-sectional average flow velocity U data respectively associated with each Z1 data according to the Z1 data sequence, the Q data sequence and the Z2~A database, so as to generate the Z1~U database;

[0042] Step 140, according to the river dynamics analytical formula U=C(H) and H=F(Z), establish the to-be-solved relation between U and Z, and calibrate the to-be-determined coefficients in the to-be-solved relation between U and Z through the Z1~U database to obtain the relation U=G(Z);

[0043] Step 150, according to the analytical formula of water-sediment dynamics S=D(U,H) and the relationship formulas H=F(Z) and U=G(Z), the relationship formula to be analyzed between S and Z is established, and the undetermined coefficients in the relationship formula to be analyzed between S and Z are calibrated through the historical measured Z1 data sequence and S data sequence to obtain the relationship formula S=W(Z);

[0044] Step 160, for the measured water level Z0 data collected at any time, substitute it into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data.

[0045] In an embodiment of the present application, by acquiring the topographic data of the target river section and the data sequences of the historically measured water level Z1, cross-sectional flow Q and cross-sectional sediment content S, and setting a number of preset water level Z2 data; according to the topographic data of the target river section, the water depth H data and the water flow area A data associated with different preset water level Z2 data are determined to generate a Z2~A database and the relationship H=F(Z); according to the Z1 data sequence, the Q data sequence and the Z2~A database, the cross-sectional average flow velocity U data respectively associated with each Z1 data is determined to generate a Z1~U database; the topographic data, flow data and water level information are effectively integrated to improve data utilization. According to the river dynamics analytical formula U=C(H) and H=F(Z), the relationship between U and Z is established, and the coefficients to be determined in the relationship between U and Z are calibrated through the Z1~U database to obtain the relationship between U=G(Z); according to the water-sediment dynamics analytical formula S=D(U,H) and the relationship between H=F(Z) and U=G(Z), the relationship between S and Z is established, and the coefficients to be determined in the relationship between S and Z are calibrated through the historical measured Z1 data sequence and S data sequence to obtain the relationship between S=W(Z); using historical measured data to analyze and calibrate the formula based on river dynamics and water-sediment dynamics can improve the accuracy of prediction and ensure the effectiveness of the analyzed relationship in actual scenarios. For the measured water level Z0 data collected at any time, substitute it into the S=W(Z) relationship to obtain the sediment content S0 data associated with the Z0 data. The real-time updated measured water level Z0 data can quickly reflect the changes in the river channel sediment content S0 data, which is convenient for reducing the time and cost of measurement and responding to emergencies. The embodiments of the present application, by gradually establishing the relationship between the water level and related hydrological characteristics, achieve efficient use of data and intelligent calculation, reduce dependence on a large amount of on-site monitoring and experimental testing work, and can improve the calculation speed and accuracy of the sediment content of the river section.

[0046] In a possible embodiment, step 120 may specifically include the following steps:

[0047] S1. Obtain the terrain data set corresponding to the target river section from the terrain data of the target river section (x i , Δ i ), x i is the lateral distance from the river bank, Δ i The horizontal distance from the river bank is x i The target river section has no severe scouring or siltation during the preset period of time.

[0048] For any Z2 data in the preset water level, execute S2-S4:

[0049] S2. According to the terrain data set, the water-passing section of the target river section is divided into several sub-sections, and the areas of each sub-section are accumulated to obtain the water-passing area A data of the entire water-passing section;

[0050] S3, calculating the wetted circumference data of each sub-section, and accumulating the wetted circumference data of each sub-section to obtain the wetted circumference χ data of the entire section;

[0051] S4, determining the ratio of the A data to the χ data as the water depth H data;

[0052] S5, repeat the above steps S2 to S4 to obtain the water depth H data and the water area A data corresponding to different preset water levels Z2, and form a Z2~A database;

[0053] S6. Perform linear fitting on the Z2 data and the H data to establish the relationship H = F(Z).

[0054] Involving terrain datasets (x i , Δ i ):

[0055] x i : Refers to the horizontal distance from the river bank edge, indicating the distance from the current position to the river bank (usually a positive value).

[0056] Δ i : corresponds to x i The riverbed elevation represents the vertical elevation value of the point.

[0057] Water-passing section: The cross-section of water flowing through the river channel, including the enclosed space formed above the riverbed and below the water surface.

[0058] Sub-section: Divide the water-passing section into several small parts to facilitate accurate calculation of the overall area of ​​the water-passing section and the corresponding wetted perimeter.

[0059] Wetted perimeter χ: the length of the river channel in contact with water, involving the effective length of the riverbed and the bank.

[0060] Specifically, S2-S4 are described below respectively:

[0061] According to the terrain dataset (x i ,Δ i ), divide the water-passing section into multiple sub-sections.

[0062] Exemplarily, the water-passing section is divided into five sub-sections, and the lateral distance and elevation of each sub-section are: (x1, Δ1), (x2, Δ2), (x3, Δ3), (x4, Δ4), and (x5, Δ5).

[0063] Calculate the area of ​​each sub-section and add these areas together to obtain the water flow area A data of the entire water flow section.

[0064] Calculate the wetted perimeter data of each sub-section in turn, accumulate all the wetted perimeter data, and obtain the wetted perimeter χ data of the entire section.

[0065] Exemplarily, the wet perimeter of the sub-section is equal to the length of the hypotenuse of a right triangle, and the lengths of the two right-angled sides of the right triangle are respectively the length along the section and the elevation difference.

[0066] The obtained water area A data is calculated by ratio with the wetted perimeter χ data to obtain the water depth H data:

[0067] The following uses a natural river as an example to illustrate S1-S6:

[0068] S1. Select a wide and shallow section of a stable river section and obtain a dataset of the section topography (x i , Δ i ), that is, the lateral distance from the river bank is x i The riverbed elevation Δ i , and draw Figure 2 Cross-sectional topographic map;

[0069] S2. According to the terrain data set, the water-passing section is divided into several sub-sections. By accumulating the areas of each sub-section, the water-passing area A of the entire section is obtained, as shown in formula (1):

[0070]

[0071] Among them, Z is the water level, N is the total number of sub-sections;

[0072] S3. Calculate the wetted perimeter data of each sub-section and obtain the wetted perimeter χ data of the entire section by accumulation, as shown in formula (2):

[0073]

[0074] S4, the ratio of A data to χ data, is the water depth H data;

[0075] S5. Repeat the above steps S2 to S4 to obtain the water depth H data and water area A data corresponding to different preset water level Z2 data, thereby forming a Z2~A database, such as Figure 3 As shown in the figure, and the Z2~H database, such as Figure 4 As shown;

[0076] S6. Based on the Π theorem, linear fitting is performed on Z2 and H data ( Figure 4 ), thus establishing the relationship H = F(Z), that is, H = 0.8304Z-156.88, and the corresponding correlation coefficient R 2 = 0.9997.

[0077] In this way, the water depth and water flow area of ​​the river at different water levels can be accurately assessed, providing an important basis for river protection and water conservancy project operation and management.

[0078] In a possible embodiment, step 130 may specifically include the following steps:

[0079] S7, determining the A data respectively associated with the Z1 data according to the Z1 data sequence and the Z2~A database;

[0080] S8. Determine U data associated with Z1 data according to A data associated with Z1 data and cross-sectional flow Q data corresponding to different Z1 data measured in history, wherein U=Q / A;

[0081] S9. According to the U data associated with different Z1 data, a Z1~U database is obtained.

[0082] The following are the instructions:

[0083] S7. Using the established Z2~A database, search for the water-passing area A data corresponding to the historical water level Z1 data. Specifically, the A data can be obtained by interpolation or direct search method.

[0084] S8. Calculate the average flow velocity U data based on the water flow area A data associated with the Z1 data and the historical measured flow rate Q data, where the Q data is the historical measured cross-sectional flow rate data corresponding to the Z1 data, and the A data is the water flow area data obtained in S7.

[0085] S9, establish a Z1~U database, record different historical water level Z1 data and corresponding calculated flow rate U data into a database, and form a Z1~U database. Figure 5 shown.

[0086] Therefore, by obtaining the water flow area A data associated with the historical water level Z1 data, the accurate reflection of the actual hydrological conditions is ensured. The actual flow Q data is used for calculation, so that the evaluation of the flow velocity U data is more in line with the actual situation and provides more reliable data support. The construction of the Z1~U database realizes the effective connection between water level, water flow area and flow velocity, which is convenient for subsequent data analysis and prediction.

[0087] In a possible embodiment, step 140 may specifically include the following steps:

[0088] S10. According to the river dynamics analytical formula U=C(H), determine the average flow velocity calculation formula:

[0089] U=J 1 / 2 / n · H 2 / 3

[0090] in, J For slope descent, J is a constant; n is the Manning coefficient, n is a constant;

[0091] S11. Combine the average velocity calculation formula with the relationship H = F(Z) to obtain the relationship between U and Z to be analyzed. J 1 / 2 / n is the coefficient to be determined;

[0092] S12, calibrate the undetermined coefficients in the undetermined relation between U and Z according to the Z1~U database, and determine J 1 / 2 / n to obtain the relationship U=G(Z).

[0093] The following is an explanation based on the formula:

[0094]

[0095] S10. According to river dynamics theory, the calculation formula of U is shown in formula (3), where J is the slope drop, which is approximately a constant, i.e. the bottom slope of the river channel; c is the roughness, c=H 1 / 6 / n; n is the Manning coefficient, which is approximately a constant. Thus, we can get formula (4);

[0096] S11. Combining the above formula (4) with the relationship H=0.8304Z-156.88 in S6, we can get formula (5); where: J 1 / 2 / n is the coefficient to be determined;

[0097] S12, use the Z1~U database described in S9 to calibrate, and obtain J 1 / 2 / n=0.4847(Correlation coefficient R 2 =0.9901), so U=0.4847·(0.8304Z-156.88) 2 / 3 ,like Figure 5 shown.

[0098] In a possible embodiment, step 150 may specifically include the following steps:

[0099] S13. According to the analytical formula of water-sand dynamics S=D(U,H), determine the first calculation formula of S:

[0100]

[0101] Among them, k1, k2, m are constants; S * is the force of water flow entraining sand; g is the acceleration of gravity, g=9.8 m / s 2 ;ω is the sediment settling velocity, ω is a constant;

[0102] S14. Define K = k1·k2 / (g·ω) m , K is a constant, and the second calculation formula of S is obtained according to the first calculation formula of S:

[0103] ;

[0104] S15. Combine the second calculation formula of S with the relationship formula H = F(Z) and the relationship formula U = G(Z) to obtain the relationship formula to be resolved between S and Z, where K and m are the coefficients to be determined;

[0105] S16. Calibrate the unknown coefficients in the unsolved relationship between S and Z through the historical measured Z1 data sequence and S data sequence, determine the K value and m value, and obtain the relationship S=W(Z).

[0106] The following is an explanation based on the formula:

[0107]

[0108] S13. According to the theory of water-sediment dynamics, the calculation formula of the cross-section sediment content S is shown in formula (6), where k1, k2, and m are constants; S * is the force of water flow entraining sand; g is the acceleration of gravity, g=9.8 m / s 2 ; ω is the sediment settling velocity. For stable river sections, ω can be approximated as a constant.

[0109] S14. Define K = k1·k2 / (g·ω) m , then K is a constant, so the calculation formula of S is as follows:

[0110] As shown;

[0111] S15. Combine formula (7) with H = 0.8304Z - 156.88, U = 0.4847·( 0.8304Z - 156.88) 2 / 3 By combining them, we can get the following relationship between S and Z, as shown in formula (8).

[0112] Among them, K and m are unknown coefficients;

[0113] S16. Obtain the historically measured Z1 data series and S data series, calibrate the undetermined coefficients in the unsolved relationship between S and Z, and obtain K=0.0059, m=9.88 (correlation coefficient R 2=0.9768), then we have the relationship S=W(Z):

[0114] S=2.8×10 -12 × (0.8304Z-156.88) 9.88 .

[0115] The effect of using the S=W(Z) relationship is as follows: Figure 6 shown.

[0116] In a possible embodiment, step 160 may specifically include the following steps:

[0117] S17, using water level detection equipment to collect measured water level Z0 data of the river section at any time;

[0118] S18. Substitute the Z0 data into the relationship S=W(Z) to calculate the corresponding sand content S0 data at any time.

[0119] For example, using a water gauge or other equipment, the cross-sectional water level Z0 = 203.602m is measured at a certain moment;

[0120] Substitute Z0=203.602m into S=2.8×10 -12 × (0.8304Z - 156.88) 9.88 The relationship is used to calculate and obtain the sand content S0 at the corresponding time = 0.1518 kg / m 3 .

[0121] In a possible embodiment, it further includes:

[0122] When the Z0 data is measured, the cross-sectional flow velocity U0 data is calculated by the relationship U=G(Z);

[0123] Through the Z2~A database, determine the water area A0 data associated with the Z0 data;

[0124] According to the U0 data and A0 data, the cross-sectional flow Q0 data is determined, where Q0 = U0·A0.

[0125] At a certain moment, the water level Z0 data is monitored. Use the relationship U = G(Z) to substitute the Z0 data to obtain the flow velocity U0 data. Through the Z2~A database, find the water flow area A0 data associated with the Z0 data.

[0126] Therefore, by measuring the water level Z0 data in real time, the current flow velocity and flow can be quickly evaluated, providing instant data support for river hydrological monitoring.

[0127] In a possible embodiment, it further includes:

[0128] When the cross-sectional flow Q0 data is measured, the Z0 data associated with the Q0 data is determined according to the Z2~A database and the U=G(Z) relationship;

[0129] According to the Z0 data and the relationship S= W(Z), the cross-section sediment content S0 data associated with the Q0 data is determined.

[0130] Use the Z2~A database and the velocity relationship U = G(Z) to find the Z0 data associated with the Q0 data, and use the established sediment content relationship S = W(Z) to substitute the Z0 data to determine the S0 data.

[0131] Therefore, by accurately calculating the flow and sediment content related to the water level Z0 data, it can help managers make scientific decisions, such as the formulation of relevant strategies such as flood prevention and water conservancy project operation.

[0132] In a possible embodiment, for the measured water level Z0 acquired at any time, after substituting into the relationship S=W(Z) to obtain S0 associated with Z0, the method further includes:

[0133] According to the S0 data, the cross-section sediment transport rate T0 data is determined, T0= Q0·S0.

[0134] The cross-section sediment transport rate is the amount of sediment transported through the river section per unit time (such as per second).

[0135] In an embodiment of the present application, by acquiring the topographic data of the target river section and the data sequences of the historically measured water level Z1, cross-sectional flow Q and cross-sectional sediment content S, and setting a number of preset water level Z2 data; according to the topographic data of the target river section, the water depth H data and the water flow area A data associated with different preset water level Z2 data are determined to generate a Z2~A database and the relationship H=F(Z); according to the Z1 data sequence, the Q data sequence and the Z2~A database, the cross-sectional average flow velocity U data respectively associated with each Z1 data is determined to generate a Z1~U database; the topographic data, flow data and water level information are effectively integrated to improve data utilization. According to the river dynamics analytical formula U=C(H) and H=F(Z), the relationship between U and Z is established, and the coefficients to be determined in the relationship between U and Z are calibrated through the Z1~U database to obtain the relationship between U=G(Z); according to the water-sediment dynamics analytical formula S=D(U,H) and the relationship between H=F(Z) and U=G(Z), the relationship between S and Z is established, and the coefficients to be determined in the relationship between S and Z are calibrated through the historical measured Z1 data sequence and S data sequence to obtain the relationship between S=W(Z); using historical measured data to analyze and calibrate the formula based on river dynamics and water-sediment dynamics can improve the accuracy of prediction and ensure the effectiveness of the analyzed relationship in actual scenarios. For the measured water level Z0 data collected at any time, substitute it into the S=W(Z) relationship to obtain the sediment content S0 data associated with the Z0 data. The real-time updated measured water level Z0 data can quickly reflect the changes in the river channel sediment content S0 data, which is convenient for reducing the time and cost of measurement and responding to emergencies. The embodiments of the present application, by gradually establishing the relationship between the water level and related hydrological characteristics, achieve efficient use of data and intelligent calculation, reduce dependence on a large amount of on-site monitoring and experimental testing work, and can improve the calculation speed and accuracy of the sediment content of the river section.

[0136] The method for rapidly calculating the sediment content of a river section provided in the embodiment of the present application can be executed by a rapid calculating device for calculating the sediment content of a river section. In the embodiment of the present application, the rapid calculating device for calculating the sediment content of a river section is used as an example to illustrate the method for rapidly calculating the sediment content of a river section provided in the embodiment of the present application.

[0137] Figure 7 700 is a block diagram of a device for quickly calculating the sediment content of a river section provided in an embodiment of the present application. The device 700 includes:

[0138] The acquisition module 710 is used to acquire the topographic data of the target river section and the data sequence of the historically measured water level Z1, the cross-section flow Q and the cross-section sediment content S, and set a number of preset water level Z2 data;

[0139] The first determination module 720 is used to determine the water depth H data and the water area A data associated with different preset water level Z2 data according to the topographic data of the target river section, so as to generate a Z2~A database and a H=F(Z) relationship;

[0140] The second determination module 730 is used to determine the cross-sectional average flow velocity U data respectively associated with each Z1 data according to the Z1 data sequence, the Q data sequence and the Z2~A database to generate the Z1~U database;

[0141] The first analytical module 740 is used to establish the relationship between U and Z to be analyzed according to the river dynamics analytical formula U=C(H) and H=F(Z), and calibrate the undetermined coefficients in the relationship between U and Z to be analyzed through the Z1~U database to obtain the relationship between U=G(Z);

[0142] The second analytical module 750 is used to establish the to-be-analyzed relation between S and Z according to the analytical formula of water-sediment dynamics S=D(U,H) and the relational expressions H=F(Z) and U=G(Z), and to calibrate the to-be-determined coefficients in the to-be-analyzed relation between S and Z through the historical measured Z1 data sequence and S data sequence, so as to obtain the relational expression S=W(Z);

[0143] The third determination module 760 is used to substitute the measured water level Z0 data collected at any time into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data.

[0144] In a possible embodiment, the first determining module 720 is specifically configured to execute:

[0145] S1. Obtain the terrain data set corresponding to the target river section from the terrain data of the target river section (x i , Δ i ), x i is the lateral distance from the river bank, Δ i The horizontal distance from the river bank is x i The target stable river section is the river section that has not experienced severe scouring or siltation within the preset period;

[0146] For any item Z2 of the preset water level, the first determination module 720 is specifically configured to execute:

[0147] S2. According to the terrain data set, the water-passing section of the target river section is divided into several sub-sections, and the areas of each sub-section are accumulated to obtain the water-passing area A data of the entire water-passing section;

[0148] S3, calculating the wetted circumference data of each sub-section, and accumulating the wetted circumference data of each sub-section to obtain the wetted circumference χ data of the entire section;

[0149] S4, determining the ratio of the A data to the χ data as the water depth H data;

[0150] S5, repeat the above steps S2 to S4 to obtain the water depth H data and the water area A data corresponding to different preset water levels Z2, and form a Z2~A database;

[0151] S6. Perform linear fitting on the Z2 data and the H data to establish the relationship H = F(Z).

[0152] In a possible embodiment, the second determining module 730 is specifically configured to execute:

[0153] S7, determining the A data respectively associated with the Z1 data according to the Z1 data sequence and the Z2~A database;

[0154] S8. Determine U data associated with Z1 data according to A data associated with Z1 data and cross-sectional flow Q data corresponding to different Z1 data measured in history, wherein U=Q / A;

[0155] S9. According to the U data associated with different Z1 data, a Z1~U database is obtained.

[0156] In a possible embodiment, the first parsing module 740 is specifically configured to execute:

[0157] S10. According to the river dynamics analytical formula U=C(H), determine the average flow velocity calculation formula:

[0158] U=J 1 / 2 / n · H 2 / 3

[0159] in, J For slope descent, J is a constant; n is the Manning coefficient, n is a constant;

[0160] S11. Combine the average velocity calculation formula with the relationship H = F(Z) to obtain the relationship between U and Z to be analyzed. J 1 / 2 / n is the coefficient to be determined;

[0161] S12, calibrate the undetermined coefficients in the undetermined relation between U and Z according to the Z1~U database, and determine J 1 / 2 / n to obtain the relationship U=G(Z).

[0162] In a possible embodiment, the second parsing module 750 is specifically configured to execute:

[0163] S13. According to the analytical formula of water-sand dynamics S=D(U,H), determine the first calculation formula of S:

[0164]

[0165] Among them, k1, k2, m are constants; S * is the force of water flow entraining sand; g is the acceleration of gravity, g=9.8 m / s 2 ;ω is the sediment settling velocity, ω is a constant;

[0166] S14. Define K = k1·k2 / (g·ω) m , K is a constant, and the second calculation formula of S is obtained according to the first calculation formula of S:

[0167] ;

[0168] S15. Combine the second calculation formula of S with the relationship formula H = F(Z) and the relationship formula U = G(Z) to obtain the relationship formula to be resolved between S and Z, where K and m are the coefficients to be determined;

[0169] S16. Calibrate the unknown coefficients in the unsolved relationship between S and Z through the historical measured Z1 data sequence and S data sequence, determine the K value and m value, and obtain the relationship S=W(Z).

[0170] In a possible embodiment, the third determining module 760 is specifically configured to:

[0171] S17, using water level detection equipment to collect measured water level Z0 data of the river section at any time;

[0172] S18. Substitute the Z0 data into the relationship S=W(Z) to calculate the corresponding sand content S0 data at any time.

[0173] In a possible embodiment, the apparatus 700 may further include a fourth determining module, specifically configured to:

[0174] When the Z0 data is measured, the cross-sectional flow velocity U0 data is calculated by the relationship U=G(Z);

[0175] Through the Z2~A database, determine the water area A0 data associated with the Z0 data;

[0176] According to the U0 data and A0 data, the cross-sectional flow Q0 data is determined, where Q0 = U0·A0.

[0177] In a possible embodiment, the apparatus 700 may further include a fifth determining module, specifically configured to:

[0178] When the cross-sectional flow Q0 data is measured, the Z0 data associated with the Q0 data is determined according to the Z2~A database and the U=G(Z) relationship;

[0179] According to the Z0 data and the relationship S= W(Z), the cross-section sediment content S0 data associated with the Q0 data is determined.

[0180] In a possible embodiment, the apparatus 700 may further include a sixth determining module, specifically configured to:

[0181] According to the S0 data, the cross-section sediment transport rate T0 data is determined, T0= Q0·S0.

[0182] In an embodiment of the present application, by acquiring the topographic data of the target river section and the data sequences of the historically measured water level Z1, cross-sectional flow Q and cross-sectional sediment content S, and setting a number of preset water level Z2 data; according to the topographic data of the target river section, the water depth H data and the water flow area A data associated with different preset water level Z2 data are determined to generate a Z2~A database and the relationship H=F(Z); according to the Z1 data sequence, the Q data sequence and the Z2~A database, the cross-sectional average flow velocity U data respectively associated with each Z1 data is determined to generate a Z1~U database; the topographic data, flow data and water level information are effectively integrated to improve data utilization. According to the river dynamics analytical formula U=C(H) and H=F(Z), the relationship between U and Z is established, and the coefficients to be determined in the relationship between U and Z are calibrated through the Z1~U database to obtain the relationship between U=G(Z); according to the water-sediment dynamics analytical formula S=D(U,H) and the relationship between H=F(Z) and U=G(Z), the relationship between S and Z is established, and the coefficients to be determined in the relationship between S and Z are calibrated through the historical measured Z1 data sequence and S data sequence to obtain the relationship between S=W(Z); using historical measured data to analyze and calibrate the formula based on river dynamics and water-sediment dynamics can improve the accuracy of prediction and ensure the effectiveness of the analyzed relationship in actual scenarios. For the measured water level Z0 data collected at any time, substitute it into the S=W(Z) relationship to obtain the sediment content S0 data associated with the Z0 data. The real-time updated measured water level Z0 data can quickly reflect the changes in the river channel sediment content S0 data, which is convenient for reducing the time and cost of measurement and responding to emergencies. The embodiments of the present application, by gradually establishing the relationship between the water level and related hydrological characteristics, achieve efficient use of data and intelligent calculation, reduce dependence on a large amount of on-site monitoring and experimental testing work, and can improve the calculation speed and accuracy of the sediment content of the river section.

[0183] The device provided in the embodiment of the present application can implement each process implemented in the above method embodiment. To avoid repetition, it will not be described again here.

[0184] Optionally, Figure 8 A schematic diagram of the hardware structure of an electronic device provided in an embodiment of the present application is shown.

[0185] The electronic device may include a processor 401 and a memory 402 storing computer program instructions.

[0186] Specifically, the processor 401 may include a central processing unit (CPU), or an application specific integrated circuit (ASIC), or may be configured to implement one or more integrated circuits of the embodiments of the present application.

[0187] The memory 402 may include a large capacity memory for data or instructions. By way of example and not limitation, the memory 402 may include a hard disk drive (HDD), a floppy disk drive, a flash memory, an optical disk, a magneto-optical disk, a tape, or a universal serial bus (USB) drive or a combination of two or more of these. Where appropriate, the memory 402 may include a removable or non-removable (or fixed) medium. Where appropriate, the memory 402 may be inside or outside the integrated gateway disaster recovery device. In a specific embodiment, the memory 402 is a non-volatile solid-state memory. In a specific embodiment, the memory 402 includes a read-only memory (ROM). Where appropriate, the ROM may be a mask-programmed ROM, a programmable ROM (PROM), an erasable PROM (EPROM), an electrically erasable PROM (EEPROM), an electrically rewritable ROM (EAROM) or a flash memory or a combination of two or more of these.

[0188] The processor 401 implements any one of the methods in the embodiments shown in the figure by reading and executing the computer program instructions stored in the memory 402.

[0189] In one example, the electronic device may further include a communication interface 404 and a bus 410. Figure 8 As shown, the processor 401, the memory 402, and the communication interface 404 are connected via a bus 410 and communicate with each other.

[0190] The communication interface 404 is mainly used to implement communication between various modules, devices, units and / or equipment in the embodiments of the present application.

[0191] Bus 410 includes hardware, software or both, and couples the components of the electronic device to each other. For example, but not limitation, the bus may include an accelerated graphics port (AGP) or other graphics bus, an enhanced industrial standard architecture (EISA) bus, a front-side bus (FSB), a hypertransport (HT) interconnect, an industrial standard architecture (ISA) bus, an infinite bandwidth interconnect, a low pin count (LPC) bus, a memory bus, a microchannel architecture (MCA) bus, a peripheral component interconnect (PCI) bus, a PCI-Express (PCI-X) bus, a serial advanced technology attachment (SATA) bus, a video electronics standard association local (VLB) bus or other suitable bus or a combination of two or more of these. Where appropriate, bus 410 may include one or more buses. Although the present application embodiment describes and shows a specific bus, the present application considers any suitable bus or interconnect.

[0192] The electronic device can execute the method in the embodiment of the present application, thereby realizing the combination Figure 1 Describe the method.

[0193] In addition, in combination with the method in the above embodiment, the embodiment of the present application can provide a computer-readable storage medium to implement. The computer-readable storage medium stores computer program instructions; when the computer program instructions are executed by the processor, the Figure 1 method.

[0194] It should be clear that the present application is not limited to the specific configuration and processing described above and shown in the figures. For the sake of simplicity, a detailed description of the known method is omitted here. In the above embodiments, several specific steps are described and shown as examples. However, the method process of the present application is not limited to the specific steps described and shown, and those skilled in the art can make various changes, modifications and additions, or change the order between the steps after understanding the spirit of the present application.

[0195] The functional blocks shown in the structural block diagram described above can be implemented as hardware, software, firmware or a combination thereof. When implemented in hardware, it can be, for example, an electronic circuit, an application-specific integrated circuit (ASIC), appropriate firmware, a plug-in, a function card, etc. When implemented in software, the elements of the present application are programs or code segments used to perform the required tasks. The program or code segment can be stored in a machine-readable medium, or transmitted on a transmission medium or a communication link by a data signal carried in a carrier. "Machine-readable medium" may include any medium capable of storing or transmitting information. Examples of machine-readable media include electronic circuits, semiconductor memory devices, ROMs, flash memories, erasable ROMs (EROMs), floppy disks, CD-ROMs, optical disks, hard disks, optical fiber media, radio frequency (RF) links, etc. The code segment can be downloaded via a computer network such as the Internet, an intranet, etc.

[0196] It should also be noted that the exemplary embodiments mentioned in this application describe some methods or systems based on a series of steps or devices. However, this application is not limited to the order of the above steps, that is, the steps can be performed in the order mentioned in the embodiment, or in a different order from the embodiment, or several steps can be performed simultaneously.

[0197] The above is only a specific implementation of the present application. Those skilled in the art can clearly understand that for the convenience and simplicity of description, the specific working processes of the systems, modules and units described above can refer to the corresponding processes in the aforementioned method embodiments, and will not be repeated here. It should be understood that the protection scope of the present application is not limited to this. Any technician familiar with the technical field can easily think of various equivalent modifications or replacements within the technical scope disclosed in this application, and these modifications or replacements should be included in the protection scope of this application.

Claims

1. A rapid calculation method for the sediment content of a river section, characterized in that: The method comprises: Obtain the topographic data of the target river section and the data sequence of the historically measured water level Z1, cross-section flow Q and cross-section sediment content S, and set several preset water level Z2 data; According to the topographic data of the target river section, the water depth H data and the water area A data associated with different preset water level Z2 data are determined to generate a Z2~A database and a H=F(Z) relationship; According to the Z1 data sequence, the Q data sequence and the Z2~A database, the cross-sectional average flow velocity U data respectively associated with each Z1 data is determined to generate the Z1~U database; According to the river dynamics analytical formula U=C(H) and H=F(Z), the undetermined relationship between U and Z is established, and the undetermined coefficients in the undetermined relationship between U and Z are calibrated through the Z1~U database to obtain the relationship U=G(Z); According to the analytical formula of water-sediment dynamics S=D(U,H) and the relationship between H=F(Z) and U=G(Z), the relationship between S and Z to be analyzed is established, and the coefficients to be determined in the relationship between S and Z to be analyzed are calibrated through the historical measured Z1 data series and S data series to obtain the relationship S=W(Z); For the measured water level Z0 data collected at any time, substitute it into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data; The method is to establish the undetermined relationship between S and Z according to the analytical formula of water-sediment dynamics S=D(U,H), H=F(Z) and U=G(Z), and calibrate the undetermined coefficients in the undetermined relationship between S and Z through the historical measured Z1 data sequence and S data sequence to obtain the relationship S=W(Z), including: S13. According to the analytical formula of water-sand dynamics S=D(U,H), determine the first calculation formula of S: Among them, k1, k2, m are constants; S * is the force of water flow entraining sand; g is the acceleration of gravity, g=9.8 m / s 2 ;ω is the sediment settling velocity, ω is a constant; S14. Define K = k1·k2 / (g·ω) m , K is a constant, and the second calculation formula of S is obtained according to the first calculation formula of S: ； S15. Combine the second calculation formula of S with the relationship formula H = F(Z) and the relationship formula U = G(Z) to obtain the relationship formula to be resolved between S and Z, where K and m are the coefficients to be determined; S16. Calibrate the unknown coefficients in the unsolved relationship between S and Z through the historical measured Z1 data sequence and S data sequence, determine the K value and m value, and obtain the relationship S=W(Z).

2. The method according to claim 1, characterized in that The method of determining the water depth H data and the water area A data associated with different preset water level Z2 data according to the topographic data of the target river section to generate a Z2~A database and a H=F(Z) relationship includes: S1. Obtain the terrain data set corresponding to the target river section from the terrain data of the target river section (x i , Δ i ), x i is the lateral distance from the river bank, Δ i The horizontal distance from the river bank is x i The target river section has no severe scouring or siltation during the preset period of time. For any Z2 data in the preset water level, execute S2-S4: S2. According to the terrain data set, the water-passing section of the target river section is divided into several sub-sections, and the areas of each sub-section are accumulated to obtain the water-passing area A data of the entire water-passing section; S3, calculating the wetted circumference data of each sub-section, and accumulating the wetted circumference data of each sub-section to obtain the wetted circumference χ data of the entire section; S4, determining the ratio of the A data to the χ data as the water depth H data; S5, repeat the above steps S2 to S4 to obtain the water depth H data and the water area A data corresponding to different preset water levels Z2, and form a Z2~A database; S6. Perform linear fitting on the Z2 data and the H data to establish the relationship H = F(Z).

3. The method according to claim 1, characterized in that Determining the cross-sectional average flow velocity U data respectively associated with each Z1 data according to the Z1 data sequence, the Q data sequence and the Z2~A database to generate the Z1~U database includes: S7, determining the A data respectively associated with the Z1 data according to the Z1 data sequence and the Z2~A database; S8. Determine U data associated with Z1 data according to A data associated with Z1 data and cross-sectional flow Q data corresponding to different Z1 data measured in history, wherein U=Q / A; S9. According to the U data associated with different Z1 data, a Z1~U database is obtained.

4. The method according to claim 1, characterized in that The method of establishing the relationship between U and Z based on the river dynamics analytical formula U=C(H) and H=F(Z) and calibrating the coefficients in the relationship between U and Z through the Z1~U database to obtain the relationship between U=G(Z) includes: S10. According to the river dynamics analytical formula U=C(H), determine the average flow velocity calculation formula: U=J 1 / 2 / n · H 2 / 3 in, J For slope descent, J is a constant; n is the Manning coefficient, n is a constant; S11. Combine the average velocity calculation formula with the relationship H = F(Z) to obtain the relationship between U and Z to be analyzed. J 1 / 2 / n is the coefficient to be determined; S12, calibrate the undetermined coefficients in the undetermined relation between U and Z according to the Z1~U database, and determine J 1 / 2 / n value to obtain the relationship U=G(Z).

5. The method according to claim 1, characterized in that Substituting the measured water level Z0 data collected at any time into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data includes: S17, using water level detection equipment to collect measured water level Z0 data of the river section at any time; S18. Substitute the Z0 data into the relationship S=W(Z) to calculate the corresponding sand content S0 data at any time.

6. The method according to claim 4, characterized in that The method further comprises: When the Z0 data is measured, the cross-sectional flow velocity U0 data is calculated by the relationship U=G(Z); Through the Z2~A database, determine the water area A0 data associated with the Z0 data; According to the U0 data and A0 data, the cross-sectional flow Q0 data is determined, where Q0 = U0·A0.

7. The method according to claim 1, characterized in that The method further comprises: When the cross-sectional flow Q0 data is measured, the Z0 data associated with the Q0 data is determined according to the Z2~A database and the U=G(Z) relationship; According to the Z0 data and the relationship S= W(Z), the cross-section sediment content S0 data associated with the Q0 data is determined.

8. The method according to claim 1, characterized in that After substituting the measured water level Z0 acquired at any time into the relationship S=W(Z) to obtain S0 associated with Z0, the method further includes: According to the S0 data, the cross-section sediment transport rate T0 data is determined, T0= Q0·S0.

9. A rapid measurement device for the sediment content of a river section, characterized in that: The device comprises: The acquisition module is used to obtain the topographic data of the target river section and the data sequence of the historically measured water level Z1, cross-section flow Q and cross-section sediment content S, and set a number of preset water level Z2 data; The first determination module is used to determine the water depth H data and the water area A data associated with different preset water level Z2 data according to the topographic data of the target river section, so as to generate a Z2~A database and a H=F(Z) relationship; The second determination module is used to determine the cross-sectional average flow velocity U data respectively associated with each Z1 data according to the Z1 data sequence, the Q data sequence and the Z2~A database to generate the Z1~U database; The first analytical module is used to establish the relationship between U and Z according to the river dynamics analytical formula U=C(H) and H=F(Z), and calibrate the undetermined coefficients in the relationship between U and Z through the Z1~U database to obtain the relationship U=G(Z); The second analytical module is used to establish the undetermined relationship between S and Z according to the analytical formula of water-sediment dynamics S=D(U,H) and the relationship between H=F(Z) and U=G(Z), and calibrate the undetermined coefficients in the undetermined relationship between S and Z through the historical measured Z1 data sequence and S data sequence to obtain the relationship S=W(Z); The third determination module is used to substitute the measured water level Z0 data collected at any time into the relationship S=W(Z) to obtain the sediment content S0 data associated with the Z0 data; The second parsing module is specifically used to execute: S13. According to the analytical formula of water-sand dynamics S=D(U,H), determine the first calculation formula of S: Among them, k1, k2, m are constants; S * is the force of water flow entraining sand; g is the acceleration of gravity, g=9.8 m / s 2 ;ω is the sediment settling velocity, ω is a constant; S14. Define K = k1·k2 / (g·ω) m , K is a constant, and the second calculation formula of S is obtained according to the first calculation formula of S: ； S15. Combine the second calculation formula of S with the relationship formula H = F(Z) and the relationship formula U = G(Z) to obtain the relationship formula to be resolved between S and Z, where K and m are the coefficients to be determined; S16. Calibrate the unknown coefficients in the unsolved relationship between S and Z through the historical measured Z1 data sequence and S data sequence, determine the K value and m value, and obtain the relationship S=W(Z).

10. The device according to claim 9, characterized in that The first determining module is specifically used to execute: S1. Obtain the terrain data set corresponding to the target river section from the terrain data of the target river section (x i , Δ i ), x i is the lateral distance from the river bank, Δ i The horizontal distance from the river bank is x i The target stable river section is the river section that has not experienced severe scouring or siltation within the preset period; For any item Z2 of the preset water level, the first determination module is specifically configured to execute: S2. According to the terrain data set, the water-passing section of the target river section is divided into several sub-sections, and the areas of each sub-section are accumulated to obtain the water-passing area A data of the entire water-passing section; S3, calculating the wetted circumference data of each sub-section, and accumulating the wetted circumference data of each sub-section to obtain the wetted circumference χ data of the entire section; S4, determining the ratio of the A data to the χ data as the water depth H data; S5, repeat the above steps S2 to S4 to obtain the water depth H data and the water area A data corresponding to different preset water levels Z2, and form a Z2~A database; S6. Perform linear fitting on the Z2 data and the H data to establish the relationship H = F(Z).

11. The device according to claim 9, characterized in that The second determination module is specifically used to execute: S7, determining the A data respectively associated with the Z1 data according to the Z1 data sequence and the Z2~A database; S8. Determine U data associated with Z1 data according to A data associated with Z1 data and cross-sectional flow Q data corresponding to different Z1 data measured in history, wherein U=Q / A; S9. According to the U data associated with different Z1 data, a Z1~U database is obtained.

12. The device according to claim 9, characterized in that The first parsing module is specifically used to execute: S10. According to the river dynamics analytical formula U=C(H), determine the average flow velocity calculation formula: U=J 1 / 2 / n · H 2 / 3 in, J For slope descent, J is a constant; n is the Manning coefficient, n is a constant; S11. Combine the average velocity calculation formula with the relationship H = F(Z) to obtain the relationship between U and Z to be analyzed. J 1 / 2 / n is the coefficient to be determined; S12, calibrate the undetermined coefficients in the undetermined relation between U and Z according to the Z1~U database, and determine J 1 / 2 / n to obtain the relationship U=G(Z).

13. The device according to claim 9, characterized in that The third determination module is specifically used to: S17, using water level detection equipment to collect measured water level Z0 data of the river section at any time; S18. Substitute the Z0 data into the relationship S=W(Z) to calculate the corresponding sand content S0 data at any time.

Images