Method and related equipment for determining nearshore riverbed erosion and deposition based on river cross-sectional topography

Abstract

This application discloses a method and related equipment for determining the amount of erosion and deposition in the nearshore riverbed based on river cross-sectional topography. The method includes obtaining the low-flow channel ranges of multiple sections within a target river section during multiple measurements; determining the nearshore channel area of ​​the corresponding sections based on the low-flow channel ranges; determining the cross-sectional area underlying the nearshore channel area for each section at each measurement; and determining the amount of erosion and deposition in the nearshore riverbed of the target river section based on the cross-sectional area. This application can more efficiently and accurately determine the amount of erosion and deposition in the nearshore riverbed.

Description

Technical Field

[0001] The present application relates to the fields of river geomorphology and water conservancy engineering, and specifically to a method and related equipment for determining the amount of scouring and silting of a nearshore riverbed based on the cross-sectional topography of a river channel. Background Art

[0002] Studying the evolution of the nearshore riverbed of alluvial rivers is a crucial area of ​​research within the river channel field. The amount of erosion and sedimentation in the nearshore riverbed provides a direct reflection of the erosion and deposition processes, revealing the interaction between water flow, sediment, and the riverbed. Understanding its dynamic evolution has important applications in bank stability prediction, river and waterway engineering design and maintenance, and the safety management of coastal infrastructure.

[0003] In related technologies, technicians can determine the amount of erosion and siltation in the nearshore riverbed based on past experience. However, this method requires a high level of technical expertise and is subject to certain subjectivity and randomness, making it difficult to meet the requirements of efficient and accurate engineering. Summary of the Invention

[0004] The embodiments of the present application provide a method and related equipment for determining the amount of scouring and silting of the nearshore riverbed based on the cross-sectional topography of the river channel, aiming to determine the amount of scouring and silting of the nearshore riverbed more efficiently and accurately.

[0005] In one aspect, the present application provides a method for determining the amount of scouring and silting of a nearshore riverbed based on the topography of a river cross section, the method comprising:

[0006] Obtain the low-water channel range of multiple sections in the target river section at multiple measurement times;

[0007] Based on the dry river channel range, determine the nearshore river channel area of ​​the corresponding section;

[0008] For each measurement of each cross section, determine the cross-sectional area beneath the nearshore channel region;

[0009] The amount of scouring and silting of the near-shore riverbed of the target river section is determined based on the cross-sectional area.

[0010] In some embodiments, determining the amount of erosion and deposition of the nearshore riverbed of the target river section based on the cross-sectional area includes:

[0011] Comparing the cross-sectional areas of the same cross-section at adjacent measurements to obtain a cross-sectional area change value;

[0012] Determine the amount of scouring and sedimentation of the nearshore riverbed of the adjacent section at the adjacent measurement times according to the change in cross-sectional area of ​​the adjacent section at the adjacent measurement times;

[0013] The amount of scouring and silting of the near-shore riverbed of the target river section is determined based on the amount of scouring and silting of the near-shore riverbed of adjacent sections at the adjacent measurements.

[0014] In some embodiments, determining the amount of nearshore riverbed erosion and deposition of the adjacent section at the adjacent measurement times based on the change in cross-sectional area of ​​the adjacent section at the adjacent measurement times includes:

[0015] Obtain the distance from the dam to multiple sections;

[0016] Comparing the distances from the dam of adjacent sections to obtain a change in distance from the dam;

[0017] The amount of scouring and silting of the nearshore riverbed of the adjacent section at the adjacent measurement times is determined based on the change in the cross-sectional area of ​​the adjacent section at the adjacent measurement times and the change in the mileage of the adjacent section from the dam.

[0018] In some embodiments, determining the amount of erosion and deposition of the nearshore riverbed of the target river section based on the amount of erosion and deposition of the nearshore riverbed of adjacent sections at adjacent measurements includes:

[0019] For all sections in the target river section, the sum of the erosion and deposition of the near-shore riverbed of adjacent sections during the adjacent measurements is determined to obtain the total erosion and deposition of the near-shore riverbed of the target river section during the adjacent measurements;

[0020] The amount of scouring and silting of the nearshore riverbed of the target river section is determined based on the total scouring and silting amount of the nearshore riverbed.

[0021] In some embodiments, determining the cross-sectional area under the nearshore channel region for each measurement of each cross-section includes:

[0022] For each measurement of each section, determining a plurality of target observation points under the nearshore channel area from a plurality of observation points of the corresponding section at the corresponding measurement;

[0023] Determining the cross-sectional geometry below the nearshore channel region based on the plurality of target observation points and the low water levels of the corresponding cross-sections;

[0024] Based on the cross-sectional geometry, the cross-sectional area beneath the nearshore channel region is determined.

[0025] In some embodiments, determining the nearshore channel area of ​​the corresponding section based on the low-flow channel range includes:

[0026] Based on the dry river channel range, determining the dry river channel width of the corresponding section at the corresponding measurement time;

[0027] 0.2 times the width of the dry river channel is used as the critical value of the nearshore river channel area;

[0028] The nearshore river channel area of ​​the corresponding section is determined according to the boundary value of the low-water river channel range and the critical value of the nearshore river channel area.

[0029] In some embodiments, determining the nearshore channel area of ​​the corresponding section based on the boundary value of the low-flow channel range and the critical value of the nearshore channel area includes:

[0030] Determine the nearshore channel area of ​​the corresponding section at the corresponding measurement time according to the boundary value of the dry river channel range and the critical value of the nearshore channel area;

[0031] The nearshore channel area of ​​the corresponding section at multiple measurements is taken as the union to obtain the nearshore channel area of ​​the corresponding section.

[0032] In some embodiments, obtaining the low-flow river channel ranges of multiple sections in the target river section at multiple measurements includes:

[0033] Obtaining observation data of each section at each measurement, wherein the observation data includes starting point distances and riverbed elevations of a plurality of observation points of the corresponding section at the corresponding measurement;

[0034] Determining a deep point among the plurality of observation points based on the riverbed elevation;

[0035] Using the low water level of the corresponding section, the plurality of observation points sequentially sorted according to the starting point distances are divided to obtain a plurality of observation point intervals;

[0036] According to the observation point interval where the deep point is located, the dry river channel range of the corresponding section in the corresponding measurement time is determined.

[0037] On the other hand, an embodiment of the present application provides a device for determining the amount of scouring and silting of a nearshore riverbed based on a river cross-sectional topography, comprising:

[0038] The first acquisition module is used to obtain the low-flow river channel range of multiple sections in the target river section at multiple measurement times;

[0039] A first determining module is used to determine the nearshore river channel area of ​​the corresponding section based on the dry river channel range;

[0040] A second determination module is configured to determine the cross-sectional area under the nearshore channel region for each measurement of each cross section;

[0041] The third determination module is used to determine the amount of erosion and deposition of the near-shore riverbed of the target river section based on the cross-sectional area.

[0042] On the other hand, the present application further provides a computer device, comprising:

[0043] one or more processors;

[0044] Memory; and

[0045] One or more applications, wherein the one or more applications are stored in the memory and configured to be executed by the processor to implement the steps of any one of the methods for determining the amount of scouring and deposition of a nearshore riverbed based on river cross-section topography.

[0046] On the other hand, the present application also provides a computer-readable storage medium having a computer program stored thereon, which is loaded by a processor to execute the steps of any of the methods for determining the amount of scouring and deposition of the nearshore riverbed based on the river cross-sectional topography.

[0047] On the other hand, an embodiment of the present application provides a computer program product, including a computer program or instructions, which are executed by a processor to implement the method for determining the amount of nearshore riverbed erosion and deposition based on river cross-sectional topography as described in any of the above items.

[0048] The embodiment of the present application provides a method and related equipment for determining the amount of scouring and silting of the nearshore riverbed based on the cross-sectional topography of the river channel, obtains the low-flow river channel range of multiple sections in the target river section at multiple measurements, determines the nearshore river channel area of ​​the corresponding section based on the low-flow river channel range, determines the cross-sectional area under the nearshore river channel area for each measurement of each section, and determines the amount of scouring and silting of the nearshore riverbed of the target river section based on the cross-sectional area. The embodiment of the present application determines the nearshore river channel area of ​​the corresponding section based on the low-flow river channel range of multiple sections in the target river section at multiple measurements, then determines the cross-sectional area under the nearshore river channel area, and determines the amount of scouring and silting of the nearshore riverbed of the target river section based on the cross-sectional area, thereby realizing the calculation of the amount of scouring and silting of the nearshore riverbed based on the cross-sectional topography of the river channel. Compared with the determination of the amount of scouring and silting of the nearshore riverbed by relevant technical personnel based on previous experience, the determination of the amount of scouring and silting of the nearshore riverbed is more efficient and accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0049] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following briefly introduces the drawings required for use in the description of the embodiments. Obviously, the drawings described below are only some embodiments of the present application. For those skilled in the art, other drawings can be obtained based on these drawings without creative work.

[0050] Figure 1 This is a schematic diagram of an embodiment of a method for determining the amount of scouring and deposition of a near-shore riverbed based on river cross-sectional topography provided in an embodiment of the present application;

[0051] Figure 2 Schematic diagram of a nearshore channel area of ​​a typical section in a target river section provided in an embodiment of the present application;

[0052] Figure 3 1 is a schematic diagram of the time variation process of the scouring and silting amount of the near-shore riverbed of the target river section provided in the embodiment of the present application;

[0053] Figure 4 This is a schematic structural diagram of an embodiment of a device for determining the amount of scouring and sedimentation of a near-shore riverbed based on river cross-sectional topography provided in an embodiment of the present application;

[0054] Figure 5 This is a schematic diagram of the terminal structure of an embodiment of the computer device provided in the embodiments of the present application. DETAILED DESCRIPTION

[0055] The following will be combined with the drawings in the embodiments of this application to clearly and completely describe the technical solutions in the embodiments of this application. Obviously, the embodiments described are only part of the embodiments of this application, not all of the embodiments. Based on the embodiments in this application, all other embodiments obtained by those skilled in the art without making any creative efforts shall fall within the scope of protection of this application.

[0056] In the description of this application, the terms "first" and "second" are used for descriptive purposes only and should not be understood to indicate or imply relative importance or implicitly specify the number of the technical features indicated. Therefore, features specified as "first" or "second" may explicitly or implicitly include one or more of the described features. In the description of this application, "plurality" means two or more, unless otherwise specifically specified.

[0057] In this application, the phrase "in some embodiments" is used to mean "serving as an example, illustration or illustration". Any embodiment described in this application as "in some embodiments" is not necessarily to be construed as being preferred or advantageous over other embodiments. The following description is given to enable any person skilled in the art to implement and use the present application. In the following description, details are listed for the purpose of explanation. It should be understood that a person of ordinary skill in the art can recognize that the present application can be implemented without using these specific details. In other instances, well-known structures and processes are not elaborated in detail to avoid obscuring the description of the present application with unnecessary details. Therefore, the present application is not intended to be limited to the embodiments shown, but is consistent with the widest scope consistent with the principles and features disclosed in this application.

[0058] It should be noted that since the system of the embodiment of the present application is executed in a computer device, the processing objects of each computer device exist in the form of data or information. For example, time is actually time information. It can be understood that if size, quantity, position, etc. are mentioned in subsequent embodiments, the corresponding data exists for the computer device to process. The details will not be repeated here.

[0059] The present application provides a method and related equipment for determining nearshore riverbed erosion and deposition based on river channel cross-sectional topography. These methods aim to automatically identify nearshore channel areas and calculate and analyze the riverbed erosion and deposition within these areas based on a long series of measured large-section data, through standardized processing of this long series of measured large-section data and combined with a low-water level dataset. The present application provides efficient, rapid, and simple identification and analysis of the changing patterns of nearshore riverbed erosion and deposition in alluvial river channels, providing a scientific basis for river channel management and bank collapse prediction. These methods are described in detail below.

[0060] In one embodiment, referring to Figure 1 The methods for determining the amount of erosion and deposition of the nearshore riverbed based on the river cross-sectional topography include:

[0061] 101. Obtain the low-flow river channel range of multiple sections in the target river section at multiple measurement times.

[0062] In the embodiments of the present application, the target river section is a designated river section, which includes multiple sections of the river channel. Measurements refer to different measurements, i.e., multiple measurements are performed on multiple sections of the target river section to determine the low-flow channel range of the multiple sections of the target river section during the multiple measurements. The low-flow channel range refers to the channel range of the section below the low-flow level, which can be calculated using the corresponding data obtained from the cross-section measurements.

[0063] In some embodiments of the present application, step 101 may include: obtaining observation data of each section in each measurement, the observation data including the starting point distance and riverbed elevation of multiple observation points of the corresponding section in the corresponding measurement, the starting point distance and riverbed elevation can be obtained based on the actual measurement of the observation points; based on the riverbed elevation, determining the deepest point among the multiple observation points, for example, the observation point with the deepest riverbed elevation among the multiple observation points can be used as the deepest point; using the low water level of the corresponding section, dividing the multiple observation points sorted in order according to the starting point distance, and obtaining For example, multiple observation point intervals can be used as an observation point sequence, where multiple observation points are sorted in ascending order according to the distance from the starting point. The riverbed elevations of the multiple observation points in the observation point sequence are then compared with the low water level of the section in turn, thereby obtaining multiple subsequences. The riverbed elevation of each observation point in each subsequence is lower than the low water level of the section, and each subsequence is an observation point interval. Based on the observation point interval where the deep-sea point is located, the low-water river channel range of the corresponding section in the corresponding measurement (that is, the starting point distance range of the section in the low-water river channel in the corresponding measurement) is determined.

[0064] In some embodiments of the present application, the low water level is determined by the following steps: obtaining the measured historical low water levels of multiple target hydrological stations in the target river section; assuming that the water surface gradient between the hydrological stations is unchanged, linear interpolation processing is performed on the measured historical low water levels of multiple target hydrological stations to obtain the low water level of each section in the target river section, so as to make the determination of the low water level more convenient and quick.

[0065] In some embodiments of the present application, the detailed content of step 101 is illustrated. Specifically, step 101 may include:

[0066] Step 1: Standardize the measured long-sequence observation data of large sections to form a data set including starting point distance X and riverbed elevation Y.

[0067] 1.1, sort out the long-term series of large-section observation data of the target river section. For a specific section, set the section number as mm, the number of observations as nn, and the number of cross-section observation points as dd, then the river section topography X and Y are three-dimensional arrays:

[0068] X(i, j, k) , Y(i, j, k)

[0069] Where X is the distance from the observation point to the starting point of the section (i.e., the starting distance); Y is the measured riverbed elevation at the observation point; i is the measurement section number; j is the number of observations; and k is the maximum number of observation points for each section. Here, i = 1, 2, …, mm; j = 1, 2, …, nn; and k = 1, 2, …, dd.

[0070] 1.2. Standardization of measured large-section terrain data. The data of different sections of the same measurement are stored in the same Excel file in the format of starting point distance X and riverbed elevation Y. The data of different sections are named with the section name and stored in different Excel files.

[0071] Step 2: Determine the measured large-section low water level dataset.

[0072] 2.1. Collect and organize the measured historical low water levels of the target hydrological stations in the target river sections;

[0073] 2.2, assuming that the water surface gradient between hydrological stations remains unchanged, linear interpolation is used to obtain the low water level of the measured large section, forming a data set containing the section name Csname, mileage from the dam L, and low water level Z_low, i.e., Csname(i), L(i), Z_low(i), where i is the measurement section number, i = 1, 2, …, mm.

[0074] Step 3: Identify the low-flow channel area of ​​the target river section.

[0075] 3.1. Use the min function in MATLAB or other software to determine the lowest point elevation (thaw point elevation Z) and its position (thaw point node d_thaw) of the measured large section, that is, Z = min(Y); d_thaw = find(Y==min(Y)).

[0076] 3.2. Extract the area where the measured riverbed elevation of the large-section observation point is less than the low-water level Z_low. According to the continuity of the section nodes, the interval where the deep point is located is determined as the low-water river channel range. Its left and right nodes are recorded as ans_node_left and ans_node_right. The difference between the starting points of the left and right nodes is taken as the low-water river channel width ans_W, that is, ans_W=X(node_right)-X(node_left).

[0077] 102. Based on the dry river channel range, determine the nearshore river channel area of ​​the corresponding section.

[0078] In the embodiments of the present application, the nearshore channel refers to the portion of the channel close to the bank in the river. The nearshore channel area can be determined based on the low-water channel range of the corresponding section. Since the river channel has nearshore channels on both the left and right banks, the nearshore channel area can include at least one of the left bank nearshore channel area and the right bank nearshore channel area.

[0079] In some embodiments of the present application, step 102 may include: determining the dry river channel width of the corresponding section at the corresponding measurement time based on the dry river channel range, for example, the span from the starting point of the dry river channel to the range can be used as the corresponding dry river channel width; determining a nearshore river channel area critical value that is smaller than the dry river channel width, the nearshore river channel area critical value refers to the range size of the nearshore river channel within the dry river channel, for example, the product of the average value of the dry river channel width of all sections in the target river section at all measurements and a preset ratio can be used as the nearshore river channel area critical value, the preset ratio can be, for example, 20%; determining the nearshore river channel area of ​​the corresponding section according to the boundary value of the dry river channel range and the nearshore river channel area critical value, for example, the sum or difference between the boundary value of the dry river channel range and the nearshore river channel area critical value, and the boundary value of the dry river channel range can be used as the nearshore river channel area of ​​the corresponding section.

[0080] In some embodiments of the present application, determining the nearshore channel area of ​​the corresponding section based on the boundary value of the dry river channel range and the critical value of the nearshore channel area can include: determining the nearshore channel area of ​​the corresponding section in the corresponding measurement time based on the boundary value of the dry river channel range and the critical value of the nearshore channel area, for example, the sum or difference between the boundary value of the dry river channel range and the critical value of the nearshore channel area, and the boundary value of the dry river channel range can be used as the nearshore channel area of ​​the corresponding section in the corresponding measurement time; taking the union of the nearshore channel areas of the corresponding section in multiple measurements to obtain the nearshore channel area of ​​the corresponding section, ensuring that the corresponding section has a unique nearshore channel area, so as to facilitate the subsequent calculation of the cross-sectional area and the nearshore riverbed scouring and deposition amount of the section under the nearshore channel area for each measurement time based on the same nearshore channel area.

[0081] In some embodiments of the present application, the detailed content of step 102 is illustrated. Specifically, step 102 may include:

[0082] Step 4: Determine the critical value of the nearshore channel area of ​​the target river section.

[0083] 4.1. Using the for loop statement in MATLAB, identify the low-flow river channel area at different measurement times of the same large-scale cross-section and calculate the low-flow river channel width W(j).

[0084] 4.2. Using the for loop statement in MATLAB, identify the low-flow river channel area at different measurement times of each measured large section and calculate the low-flow river channel width W(i).

[0085] 4.3. The average value of the low-water channel width of each large section of the target river section at different times is taken as the representative river width of the target river section, and 0.2 times of this value is determined as the critical value of the nearshore channel area of ​​the target river section, that is, WCV=0.20 W_Typ.

[0086] Step 5: Identify the nearshore channel area of ​​the target river section.

[0087] 5.1, according to the left and right nodes ans_node_left and ans_node_right of the dry river channel determined in 3.2, determine the nearshore channel range, where the left bank nearshore channel area is [X(ans_node_left), X(ans_node_left)+WCV], and the right bank nearshore channel area is [X(ans_node_right)-WCV, X(ans_node_right)].

[0088] 5.2. Use the for loop statement in MATLAB to determine the nearshore channel area of ​​the same measured large section at different times, and use their union as the nearshore channel area of ​​the target section.

[0089] 5.3, use the for loop statement in MATLAB to determine the nearshore channel area at different sections.

[0090] 103. For each measurement of each section, determine the cross-sectional area under the nearshore channel area.

[0091] In the embodiment of the present application, since the cross section is a vertical plane in the river channel, the cross-sectional area of ​​the cross section under the nearshore channel region can be determined for each measurement of each cross section.

[0092] In some embodiments of the present application, since the nearshore channel area refers to the low water level of the section, step 103 may include: for each measurement of each section, determining multiple target observation points under the nearshore channel area among multiple observation points of the corresponding section in the corresponding measurement, wherein each observation point has a corresponding starting distance and riverbed elevation, and the starting distance corresponding to the target observation point is within the starting distance range of the nearshore channel area; determining the cross-sectional geometry under the nearshore channel area based on the multiple target observation points and the low water level of the corresponding section, for example, the water level lines of the multiple target observation points and the low water level of the corresponding section can be used as the boundary values ​​of the cross-sectional geometry under the nearshore channel area, thereby determining the cross-sectional geometry under the nearshore channel area; determining the cross-sectional area under the nearshore channel area based on the cross-sectional geometry, for example, the cross-sectional area A under the nearshore channel area can be calculated by the trapezoidal area summation method.

[0093] 104. Determine the amount of scouring and silting of the nearshore riverbed of the target river section based on the cross-sectional area.

[0094] In the embodiment of the present application, the amount of erosion and deposition in the near-shore riverbed of the target river section refers to the amount of erosion and deposition within the near-shore riverbed of the target river section.

[0095] In some embodiments of the present application, step 104 may include: comparing the cross-sectional areas of the same section in adjacent measurements to obtain a cross-sectional area change value, where the cross-sectional area change value may be, for example, the difference in the cross-sectional areas of the section in adjacent measurements; determining the amount of scouring and silting of the nearshore riverbed of the adjacent section in the adjacent measurements based on the cross-sectional area change values ​​of the adjacent sections in the adjacent measurements, wherein the amount of scouring and silting of the nearshore riverbed of the adjacent section in the adjacent measurements may be positively correlated with the cross-sectional area change values ​​of the adjacent section in the adjacent measurements; determining the amount of scouring and silting of the nearshore riverbed of the target river section based on the amount of scouring and silting of the nearshore riverbed of the adjacent section in the adjacent measurements.

[0096] In some embodiments of the present application, determining the amount of scouring and silting of the nearshore riverbed of adjacent sections in adjacent measurements based on the change values ​​of the cross-sectional areas of adjacent sections in adjacent measurements can include: obtaining the mileage from the dam of multiple sections; comparing the mileage from the dam of adjacent sections to obtain the change value of the mileage from the dam, which can be, for example, the difference between the mileages from the dam of adjacent sections, that is, the distance between adjacent sections; determining the amount of scouring and silting of the nearshore riverbed of adjacent sections in adjacent measurements based on the change values ​​of the cross-sectional areas of adjacent sections in adjacent measurements and the change values ​​of the mileage from the dam of adjacent sections, for example, the change values ​​of the cross-sectional areas of adjacent sections in adjacent measurements can be multiplied by the change values ​​of the mileage from the dam of adjacent sections to obtain the change amount of scouring and silting of the adjacent sections during the adjacent measurements, and use it as the amount of scouring and silting of the nearshore riverbed of the adjacent sections in the adjacent measurements.

[0097] In some embodiments of the present application, determining the nearshore riverbed erosion and deposition volume of a target river section based on the nearshore riverbed erosion and deposition volume of adjacent sections in adjacent measurements may include: determining the sum of the nearshore riverbed erosion and deposition volume of adjacent sections in adjacent measurements for all sections in the target river section to obtain the total nearshore riverbed erosion and deposition volume of the target river section in adjacent measurements. For example, the changes in erosion and deposition volume of all sections in the target river section during adjacent measurements may be summed to obtain the total nearshore riverbed erosion and deposition volume of the target river section in adjacent measurements; and determining the nearshore riverbed erosion and deposition volume of the target river section based on the total nearshore riverbed erosion and deposition volume. The nearshore riverbed erosion and deposition volume of the target river section may include the total nearshore riverbed erosion and deposition volume of the target river section in each measurement, so that based on the total nearshore riverbed erosion and deposition volume of the target river section in each different measurement, the changing pattern of the nearshore riverbed erosion and deposition volume of the target river section can be analyzed, thereby providing a scientific basis for river channel management and bank collapse prediction.

[0098] In some embodiments of the present application, the detailed content of step 104 is illustrated. Specifically, step 104 may include:

[0099] Step 6: Calculate the amount of scouring and silting in the nearshore riverbed area.

[0100] 6.1, using the for loop statement in MATLAB, calculate the cross-sectional area A of the nearshore channel at different times of measurement of the same large cross-section j and its change value △A j =A j -A j+1 .

[0101] 6.2, using the for loop statement in MATLAB, calculate the cross-sectional area A of the nearshore channel at different measurement times for each measured section i,j and its change value △A i,j .

[0102] 6.3, the trapezoidal method is used to calculate the amount of erosion and deposition of the nearshore riverbed between adjacent measured large sections, that is, V i,j =(△A i,j +△A i+1,j) ×(L i+1 -L i ) / 2.

[0103] 6.4, calculate the amount of erosion and sedimentation of the near-shore riverbed at adjacent measurements in the target river section, i.e. V j =∑V i,j .

[0104] Step 7: Draw the change process of riverbed erosion and sedimentation near the bank of the target river section.

[0105] 7.1 Based on the calculation results in step 6.4, use the plot function in MATLAB to plot the time variation process of the nearshore riverbed erosion and deposition volume and analyze the variation pattern of the nearshore riverbed erosion and deposition volume.

[0106] In the technical solution disclosed in this embodiment, the nearshore channel area of ​​the corresponding section is determined based on the low-water river channel range of multiple sections in the target river section at multiple measurements, and then the cross-sectional area under the nearshore channel area is determined. The nearshore riverbed scouring and silting amount of the target river section is determined based on the cross-sectional area, thereby realizing the calculation of the nearshore riverbed scouring and silting amount based on the river section topography. Compared with the determination of the nearshore riverbed scouring and silting amount based on previous experience by relevant technical personnel, the determination of the nearshore riverbed scouring and silting amount is more efficient and accurate.

[0107] Furthermore, by standardizing measured data, constructing a low-water level dataset, identifying regions, and calculating erosion and deposition, the system achieves automated identification of alluvial river channel areas and automated mapping of changes in erosion and deposition characteristics. This addresses the existing issues of subjective nearshore riverbed definition and low computational efficiency, offering advantages such as data standardization, precise calculations, and dynamic analysis.

[0108] Below, refer to Figure 1 , an example is given to illustrate the method for determining the amount of erosion and deposition of the near-shore riverbed based on the river cross-sectional topography. Specifically, the identification of the variation pattern of the amount of erosion and deposition of the near-shore riverbed includes the following steps:

[0109] Step 1: Standardization of long-sequence observation data of large cross-sections.

[0110] Collect and organize historical observation data from long sequences of large sections of the target river reach. Data from different measurements of the same section are stored in the same Excel file in the format of starting point distance X and elevation Y, as shown in Table 1. Data from different sections are named after the section and stored in different Excel files.

[0111] Table 1

[0112]

[0113] Step 2: Construct a dataset of measured low water levels in large sections.

[0114] The low water level data of typical hydrological stations or typical sections of the target river section are collected and collated. Assuming that the water surface gradient between hydrological stations or typical sections remains unchanged, linear interpolation is used to obtain the low water level of the measured large section, forming a data set containing the section name Csname, mileage from the dam L, and low water level Z_low, as shown in Table 2 below.

[0115] Table 2

[0116]

[0117] Step three: Identify the dry river channel area of ​​the target river section.

[0118] Taking the measured large-section data of Jing 175 section in 2003 as an example, the process of identifying the low-flow channel area is explained. First, the min function in MATLAB is used to determine the lowest point elevation of the measured large section, namely the deep-water point elevation Z=1.40, and the find statement is used to determine the section node number d_thaw=11 where the deep-water point is located. Then, the find statement in MATLAB is used to extract the area where the riverbed elevation of the section is less than the low-water level Z_low=22.09m of the section. According to the continuity of the section nodes, the interval where the deep-water point is located is determined as the low-flow channel range, and its left and right nodes ans_node_left=10 and ans_node_right=43 include the section node 11 where the deep-water point is located. Finally, the difference between the starting point distances of the left and right nodes is used as the dry-flow channel width ans_W, that is, ans_W= X(43)-X(10) =1024.

[0119] Step 4: Determine the critical value of the nearshore channel area of ​​the target river section.

[0120] First, the for loop statement in MATLAB was used to identify the low-flow river channel area of ​​Jing 175 section at different measurement times from 2001 to 2021 and calculate the low-flow river channel width W(1,j); then, the for loop statement in MATLAB was used to identify the low-flow river channel area of ​​each measured large section between Jing 175 section and Jing 179+1 section at different measurement times and calculate the low-flow river channel width W(i,j); finally, the average value W_Typ of the low-flow river channel width of each measured large section between Jing 175 section and Jing 179+1 section at different measurement times was used as the representative river width of the target river section, W_Typ=mean(W(i,j))=864.55, and 20% of this value was determined as the critical value of the nearshore river channel area of ​​the target river section, that is, WCV=0.20 W_Typ=173.

[0121] Step 5: Identify the nearshore channel area of ​​the target river section.

[0122] First, according to the low-water channel data of section Jing 175 in 2003, the left and right nodes ans_node_left=10 and ans_node_right=43 are determined in step 3 to determine the nearshore channel range. The left bank nearshore channel area is [X(ans_node_left)=72,X(ans_node_left)+WCV=245], and the right bank channel area is [X(ans_node_right)-WCV=923,X(ans_node_right)=1096]. Then, the nearshore channel area of ​​section Jing 175 at different times is determined by the for loop statement in MATLAB, and the union of them is used as the nearshore channel area of ​​the target section. The left bank nearshore channel area ranges from [45.48,253,28]. Finally, the nearshore channel area of ​​different sections is determined by the for loop statement in MATLAB. The nearshore channel area of ​​a typical section in the target river section is as follows: Figure 2 As shown in .

[0123] Step six: Calculate the amount of erosion and sedimentation in the nearshore riverbed.

[0124] First, the trapezoidal area method is used to calculate the cross-sectional area A of the nearshore channel area below the low water level. The for loop statement in MATLAB is used to calculate the cross-sectional area A of the nearshore channel area at different measurement times of the same measured large section. j and its change value △A j =A j -A j+1 Then, use the for loop statement in MATLAB to calculate the cross-sectional area A of the nearshore channel area at different times of each measured section. i,j and its change value △A i,j ; Then, the trapezoidal method is used to calculate the amount of erosion and deposition of the near-shore riverbed between adjacent measured large sections, that is, V i,j =(△A i,j +△A i+1,j )×(L i+1 -L i ) / 2; Finally, the sum of the scouring and silting volume of the near-shore riverbed of the target river section in adjacent measurements is taken as the scouring and silting volume of the near-shore riverbed of the river section, that is, V j =∑V i,j .

[0125] Step seven: Draw the changes in the amount of erosion and sedimentation in the riverbed near the bank of the target river section.

[0126] Based on the calculation results in step 6, use the plot function in MATLAB to plot the time variation process of the nearshore riverbed erosion and deposition. The variation process of the nearshore riverbed erosion and deposition on the left side of the Jing 175 section to the Jing 179+1 section from 2003 to 2021 is shown as follows: Figure 3 shown.

[0127] The present invention is based on long-term series of measured large-section data. The boundaries of the left and right banks of the dry river channel where the fen are located are used as the boundaries of the region, 20% of the average width of the dry river channel in the target river section is used as the regional critical value, and the range of the critical values ​​of the river channel area near the boundaries of the left and right banks in the dry river channel is used as the region. The MATLAB program is used to realize the automatic identification of the region and the automatic mapping of the changes in the nearshore scouring and silting characteristic quantities. This technology is based on measured large-section data, has low difficulty and low cost in data acquisition, high accuracy, and can efficiently, quickly and accurately identify the temporal changes in the scouring and silting amount in the nearshore riverbed area. It is of great significance for predicting the evolution characteristics of the river channel and formulating appropriate river (navigation) regulation plans.

[0128] In order to better implement the method for determining the amount of scouring and silting of the near-shore riverbed based on the river cross-sectional topography in the embodiment of the present application, on the basis of the method for determining the amount of scouring and silting of the near-shore riverbed based on the river cross-sectional topography, the embodiment of the present application also provides a device for determining the amount of scouring and silting of the near-shore riverbed based on the river cross-sectional topography, such as Figure 4 As shown, the device 400 for determining the amount of erosion and deposition of the near-shore riverbed based on the river cross-sectional topography includes:

[0129] The first acquisition module 401 is used to obtain the low-flow river channel range of multiple sections in the target river section at multiple measurement times;

[0130] The first determining module 402 is used to determine the nearshore channel area of ​​the corresponding section based on the low-water channel range;

[0131] The second determination module 403 is used to determine the cross-sectional area under the nearshore channel region for each measurement of each cross section;

[0132] The third determining module 404 is used to determine the amount of erosion and deposition of the near-shore riverbed of the target river section based on the cross-sectional area.

[0133] The present application also provides a computer device that integrates any of the apparatuses for determining the amount of erosion and sedimentation of a near-shore riverbed based on the river cross-sectional topography provided in the present application. Figure 5 , which shows a schematic diagram of the structure of the computer device involved in the embodiment of the present application, specifically:

[0134] The computer device may include one or more processing core processors 501, one or more computer readable storage media memories 502, a power supply 503, an input unit 504 and other components. Those skilled in the art will understand that Figure 5 The computer device structure shown in the figure is not intended to limit the computer device, and may include more or fewer components than shown in the figure, or combine certain components, or arrange components differently.

[0135] Processor 501 is the control center of the computer device. It connects the various components of the entire computer device using various interfaces and lines. By running or executing software programs and / or modules stored in memory 502 and accessing data stored in memory 502, it performs various functions of the computer device and processes data, thereby providing overall monitoring of the computer device. Optionally, processor 501 may include one or more processing cores; preferably, processor 501 may integrate an application processor and a modem processor, wherein the application processor primarily processes the operating system, user interface, and application programs, while the modem processor primarily handles wireless communications. It is understood that the modem processor may not be integrated into processor 501.

[0136] Memory 502 can be used to store software programs and modules. Processor 501 executes various functional applications and data processing by running the software programs and modules stored in memory 502. Memory 502 may primarily include a program storage area and a data storage area. The program storage area may store an operating system and at least one application required for a function (such as sound playback or image playback); the data storage area may store data generated based on the use of the computer device. Furthermore, memory 502 may include high-speed random access memory (RAM) and non-volatile memory, such as at least one disk storage device, flash memory device, or other volatile solid-state storage device. Accordingly, memory 502 may also include a memory controller to provide processor 501 with access to memory 502.

[0137] The computer device also includes a power supply 503 for supplying power to various components. Preferably, the power supply 503 can be logically connected to the processor 501 via a power management system, thereby enabling the power management system to manage charging, discharging, and power consumption. The power supply 503 can also include one or more DC or AC power supplies, a recharging system, a power failure detection circuit, a power converter or inverter, a power status indicator, and other arbitrary components.

[0138] The computer device may further include an input unit 504 , which may be configured to receive input digital or character information and generate keyboard, mouse, joystick, optical or trackball signal input related to user settings and function control.

[0139] Although not shown, the computer device may further include a display unit, etc., which will not be described in detail here. Specifically, in this embodiment, the processor 501 in the computer device will load the executable files corresponding to one or more application processes into the memory 502 according to the following instructions, and the processor 501 will run the application stored in the memory 502 to implement various functions as follows:

[0140] Obtain the low-flow channel range of multiple sections in the target river section at multiple measurements; determine the nearshore channel area of ​​the corresponding section based on the low-flow channel range; determine the cross-sectional area under the nearshore channel area for each measurement of each section; and determine the nearshore riverbed erosion and deposition volume of the target river section based on the cross-sectional area.

[0141] Those skilled in the art will appreciate that all or part of the steps in the various methods of the above embodiments may be accomplished by instructions, or by controlling related hardware through instructions. The instructions may be stored in a computer-readable storage medium and loaded and executed by a processor.

[0142] To this end, embodiments of the present application provide a computer-readable storage medium, which may include a read-only memory (ROM), random access memory (RAM), a disk, or an optical disk. A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps of any of the methods for determining nearshore riverbed erosion and deposition based on river channel cross-sectional topography provided in embodiments of the present application. For example, the computer program, when loaded by the processor, may execute the following steps:

[0143] Obtain the low-flow channel range of multiple sections in the target river section at multiple measurements; determine the nearshore channel area of ​​the corresponding section based on the low-flow channel range; determine the cross-sectional area under the nearshore channel area for each measurement of each section; and determine the nearshore riverbed erosion and deposition volume of the target river section based on the cross-sectional area.

[0144] The embodiments of the present application further provide a computer program product or computer program, which includes computer instructions stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium and executes the computer instructions, so that the electronic device executes the method for determining the amount of erosion and deposition of a nearshore riverbed based on a river channel cross-sectional topography as described above, for example:

[0145] Obtain the low-flow channel range of multiple sections in the target river section at multiple measurements; determine the nearshore channel area of ​​the corresponding section based on the low-flow channel range; determine the cross-sectional area under the nearshore channel area for each measurement of each section; and determine the nearshore riverbed erosion and deposition volume of the target river section based on the cross-sectional area.

[0146] In the above embodiments, the description of each embodiment has its own focus. For parts that are not described in detail in a certain embodiment, please refer to the detailed description of other embodiments above and will not be repeated here.

[0147] In specific implementation, the above units or structures can be implemented as independent entities, or can be arbitrarily combined to implement as the same or several entities. The specific implementation of the above units or structures can refer to the previous method embodiments and will not be repeated here.

[0148] The specific implementation of the above operations can be found in the previous embodiments and will not be repeated here.

[0149] The above is a detailed introduction to a method for determining the amount of nearshore riverbed scouring and deposition based on the river cross-sectional topography and related equipment provided in an embodiment of the present application. Specific examples are used herein to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea; at the same time, for technical personnel in this field, based on the ideas of the present application, there will be changes in the specific implementation methods and application scope. In summary, the content of this specification should not be understood as a limitation on the present application.

Claims

1. A method for determining the amount of scouring and silting of a nearshore riverbed based on river cross-sectional topography, characterized in that: The method comprises: Obtain the low-water channel range of multiple sections in the target river section at multiple measurement times; Based on the dry river channel range, determine the nearshore river channel area of ​​the corresponding section; For each measurement of each cross section, determine the cross-sectional area beneath the nearshore channel region; Determining the amount of scouring and sedimentation of the nearshore riverbed of the target river section based on the cross-sectional area; Determining the nearshore channel area of ​​the corresponding section based on the dry channel range includes: determining the dry channel width of the corresponding section at the corresponding measurement time based on the dry channel range; using 0.2 times the dry channel width as a nearshore channel area critical value; and determining the nearshore channel area of ​​the corresponding section based on the boundary value of the dry channel range and the nearshore channel area critical value; The method of determining the nearshore channel area of ​​the corresponding section according to the boundary value of the low-water channel range and the critical value of the nearshore channel area includes: determining the nearshore channel area of ​​the corresponding section in the corresponding measurement times according to the boundary value of the low-water channel range and the critical value of the nearshore channel area; and taking the union of the nearshore channel areas of the corresponding section in multiple measurements to obtain the nearshore channel area of ​​the corresponding section.

2. The method for determining the amount of scouring and sedimentation of the near-shore riverbed based on the river cross-sectional topography according to claim 1, characterized in that: Determining the amount of scouring and silting of the nearshore riverbed of the target river section based on the cross-sectional area includes: Comparing the cross-sectional areas of the same cross-section at adjacent measurements to obtain a cross-sectional area change value; Determine the amount of scouring and sedimentation of the nearshore riverbed of the adjacent section at the adjacent measurement times according to the change in cross-sectional area of ​​the adjacent section at the adjacent measurement times; The amount of scouring and silting of the near-shore riverbed of the target river section is determined based on the amount of scouring and silting of the near-shore riverbed of adjacent sections at the adjacent measurements.

3. The method for determining the amount of scouring and sedimentation of the near-shore riverbed based on the river cross-sectional topography according to claim 2, characterized in that: Determining the amount of scouring and silting of the nearshore riverbed of the adjacent section at the adjacent measurement times based on the change in cross-sectional area of ​​the adjacent section at the adjacent measurement times includes: Obtain the distance from the dam to multiple sections; Comparing the distances from the dam of adjacent sections to obtain a change in distance from the dam; The amount of scouring and silting of the nearshore riverbed of the adjacent section at the adjacent measurement times is determined based on the change in the cross-sectional area of ​​the adjacent section at the adjacent measurement times and the change in the mileage of the adjacent section from the dam.

4. The method for determining the amount of scouring and sedimentation of the near-shore riverbed based on the river cross-sectional topography according to claim 2, characterized in that: Determining the amount of scouring and silting of the near-shore riverbed of the target river section according to the amount of scouring and silting of the near-shore riverbed of the adjacent sections at the adjacent measurements includes: For all sections in the target river section, the sum of the erosion and deposition of the near-shore riverbed of adjacent sections during the adjacent measurements is determined to obtain the total erosion and deposition of the near-shore riverbed of the target river section during the adjacent measurements; The amount of scouring and silting of the nearshore riverbed of the target river section is determined based on the total scouring and silting amount of the nearshore riverbed.

5. The method for determining the amount of scouring and sedimentation of the near-shore riverbed based on the river cross-sectional topography according to claim 1, characterized in that: Determining the cross-sectional area below the nearshore channel region for each measurement of each cross-section includes: For each measurement of each section, determining a plurality of target observation points under the nearshore channel area from a plurality of observation points of the corresponding section at the corresponding measurement; Determining the cross-sectional geometry below the nearshore channel region based on the plurality of target observation points and the low water levels of the corresponding cross-sections; Based on the cross-sectional geometry, the cross-sectional area beneath the nearshore channel region is determined.

6. The method for determining the amount of scouring and sedimentation of the near-shore riverbed based on the river cross-sectional topography according to claim 1, characterized in that: The method of obtaining the low-flow river channel ranges of multiple sections in the target river section at multiple measurement times includes: Obtaining observation data of each section at each measurement, wherein the observation data includes starting point distances and riverbed elevations of a plurality of observation points of the corresponding section at the corresponding measurement; Determining a deep point among the plurality of observation points based on the riverbed elevation; Using the low water level of the corresponding section, the plurality of observation points sequentially sorted according to the starting point distances are divided to obtain a plurality of observation point intervals; According to the observation point interval where the deep point is located, the dry river channel range of the corresponding section in the corresponding measurement time is determined.

7. A computer device, characterized in that: The computer device comprises: one or more processors; Memory; and one or more applications, wherein the one or more applications are stored in the memory and configured to be executed by the processor to implement the steps of the method for determining the amount of nearshore riverbed erosion and deposition based on river channel cross-sectional topography as described in any one of claims 1 to 6.

8. A computer-readable storage medium, characterized in that A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps in the method for determining the amount of scouring and deposition of a nearshore riverbed based on the river cross-sectional topography according to any one of claims 1 to 6.

Images