River section terrain-based near-shore riverbed erosion and deposition amount determination method and related equipment

Abstract

The invention discloses a method for determining the scouring and silting amount of a near-shore riverbed based on the cross-section terrain of a river channel and related equipment. The method comprises the following steps: acquiring low-water river channel ranges of a plurality of sections in a target river reach at a plurality of measurement times; determining a near-shore river channel area of the corresponding section based on the low-water river channel range; aiming at each measurement time of each section, determining the section area under the near-shore river channel area; and according to the sectional area, determining the near-shore riverbed erosion and deposition amount of the target river reach. According to the method, the erosion and deposition amount of the near-shore riverbed can be determined more efficiently and accurately.

Description

Technical Field

[0001] This application relates to the fields of fluvial geomorphology and hydraulic engineering, and specifically relates to a method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of a river channel and related equipment. Background Art

[0002] The study of the evolution law of the nearshore riverbed in alluvial rivers is an important research content in the field of river channels. The scouring and silting volume of the nearshore riverbed can intuitively reflect the scouring and sedimentation processes of the nearshore riverbed, and reveal the interaction mechanism among water flow, sediment, and the riverbed. Mastering its dynamic change law has important application value for predicting the stability of riverbanks, designing and maintaining river (navigation) channels, and the safety management of coastal infrastructure.

[0003] In related technologies, relevant technicians can determine the scouring and silting volume of the nearshore riverbed based on past experience. However, this method has high requirements for the level of relevant technicians, and there is a certain degree of subjectivity and randomness, making it difficult to meet the efficient and accurate engineering requirements. Summary of the Invention

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

[0005] On the one hand, this application provides a method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of a river channel. The method includes: Obtain the range of the dry-season river channel at multiple cross-sections in a target river reach at multiple measurement times; Based on the range of the dry-season river channel, determine the nearshore river channel area of the corresponding cross-section; For each measurement time of each cross-section, determine the cross-sectional area under the nearshore river channel area; Based on the cross-sectional area, determine the scouring and silting volume of the nearshore riverbed of the target river reach.

[0006] In some embodiments, the step of determining the scouring and silting volume of the nearshore riverbed of the target river reach based on the cross-sectional area includes: Compare the cross-sectional areas of the same cross-section at adjacent measurement times to obtain the change value of the cross-sectional area; Based on the change value of the cross-sectional area of adjacent cross-sections at the adjacent measurement times, determine the scouring and silting volume of the nearshore riverbed of the adjacent cross-sections at the adjacent measurement times; Based on the scouring and silting volume of the nearshore riverbed of adjacent cross-sections at the adjacent measurement times, determine the scouring and silting volume of the nearshore riverbed of the target river reach.

[0007] In some embodiments, the step of determining the scouring and silting volume of the nearshore riverbed of adjacent cross-sections at the adjacent measurement times based on the change value of the cross-sectional area of adjacent cross-sections at the adjacent measurement times includes: Obtain the mileage from the dam for multiple cross-sections; Compare the mileage from the dam of adjacent cross-sections to obtain the change value of the mileage from the dam; Determine the near-shore riverbed scouring and silting volume of adjacent cross-sections in the adjacent measurement times according to the change value of the cross-sectional area and the change value of the mileage from the dam of adjacent cross-sections in the adjacent measurement times.

[0008] In some embodiments, the determining the near-shore riverbed scouring and silting volume of the target river reach according to the near-shore riverbed scouring and silting volume of adjacent cross-sections in the adjacent measurement times includes: For all cross-sections in the target river reach, determine the sum of the near-shore riverbed scouring and silting volumes of adjacent cross-sections in the adjacent measurement times to obtain the total near-shore riverbed scouring and silting volume of the target river reach in the adjacent measurement times; Determine the near-shore riverbed scouring and silting volume of the target river reach according to the total near-shore riverbed scouring and silting volume.

[0009] In some embodiments, the determining the cross-sectional area under the near-shore river channel region for each measurement time of each cross-section includes: For each measurement time of each cross-section, determine multiple target observation points under the near-shore river channel region among multiple observation points of the corresponding cross-section at the corresponding measurement time; Determine the cross-sectional geometric shape under the near-shore river channel region according to the multiple target observation points and the low water level of the corresponding cross-section; Based on the cross-sectional geometric shape, determine the cross-sectional area under the near-shore river channel region.

[0010] In some embodiments, the determining the near-shore river channel region of the corresponding cross-section based on the low-flow river channel range includes: Based on the low-flow river channel range, determine the low-flow river channel width of the corresponding cross-section at the corresponding measurement time; Take 0.2 times of the low-flow river channel width as the critical value of the near-shore river channel region; Determine the near-shore river channel region of the corresponding cross-section according to the boundary value of the low-flow river channel range and the critical value of the near-shore river channel region.

[0011] In some embodiments, the determining the near-shore river channel region of the corresponding cross-section according to the boundary value of the low-flow river channel range and the critical value of the near-shore river channel region includes: Determine the near-shore river channel region of the corresponding cross-section at the corresponding measurement time according to the boundary value of the low-flow river channel range and the critical value of the near-shore river channel region; Take the union of the near-shore river channel regions of the corresponding cross-section at multiple measurement times to obtain the near-shore river channel region of the corresponding cross-section.

[0012] In some embodiments, the obtaining the low-flow river channel ranges of multiple cross-sections in the target river reach at multiple measurement times includes: Obtain the observation data of each cross-section at each measurement time. The observation data includes the starting distance and riverbed elevation of multiple observation points in the corresponding cross-section at the corresponding measurement time; Based on the riverbed elevation, determine the thalweg points among the multiple observation points; Use the low water level of the corresponding cross-section to divide the multiple observation points sorted in sequence according to the starting distance, and obtain multiple observation point intervals; According to the observation point interval where the thalweg point is located, determine the low water channel range of the corresponding cross-section at the corresponding measurement time.

[0013] On the other hand, an embodiment of the present application provides a device for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography, including: A first acquisition module, configured to acquire the low water channel ranges of multiple cross-sections in a target river reach at multiple measurement times; A first determination module, configured to determine the nearshore channel area of the corresponding cross-section based on the low water channel range; A second determination module, configured to determine the cross-sectional area under the nearshore channel area for each measurement time of each cross-section; A third determination module, configured to determine the erosion and deposition volume of the nearshore riverbed of the target river reach according to the cross-sectional area.

[0014] On the other hand, the present application further provides a computer device, which includes: One or more processors; A memory; and One or more applications, wherein the one or more applications are stored in the memory and are configured to be executed by the processor to implement the steps in any one of the methods for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography.

[0015] On the other hand, the present application further provides a computer-readable storage medium, on which a computer program is stored, and the computer program is loaded by a processor to execute the steps in any one of the methods for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography.

[0016] On the other hand, an embodiment of the present application provides a computer program product, including a computer program or instruction, and the computer program or instruction is executed by a processor to implement the method for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography as described in any one of the above.

[0017] The method and related equipment for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography provided by the embodiments of the present application obtain the low-flow river channel ranges of multiple cross-sections in a target river reach at multiple measurement times. Based on the low-flow river channel ranges, the nearshore river channel areas of the corresponding cross-sections are determined. For each measurement time of each cross-section, the cross-sectional area under the nearshore river channel area is determined, and the erosion and deposition volume of the nearshore riverbed of the target river reach is determined according to the cross-sectional area. According to the low-flow river channel ranges of multiple cross-sections in the target river reach, the embodiments of the present application determine the nearshore river channel areas of the corresponding cross-sections, then determine the cross-sectional areas under the nearshore river channel areas, and determine the erosion and deposition volume of the nearshore riverbed of the target river reach according to the cross-sectional areas, realizing the calculation of the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography. Compared with the related technicians who determine the erosion and deposition volume of the nearshore riverbed according to past experience, the determination of the erosion and deposition volume of the nearshore riverbed is more efficient and accurate. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the drawings required for the description of the embodiments will be briefly introduced below. Obviously, the drawings in the following description 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 efforts.

[0019] Figure 1 It is a schematic diagram of an embodiment of the method for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography provided in the embodiments of the present application; Figure 2 It is a schematic diagram of the nearshore river channel area of a typical cross-section in the target river reach provided in the embodiments of the present application; Figure 3 It is a schematic diagram of the time variation process of the erosion and deposition volume of the nearshore riverbed of the target river reach provided in the embodiments of the present application; Figure 4 It is a schematic diagram of the structure of an embodiment of the device for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography provided in the embodiments of the present application; Figure 5 It 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 OF THE EMBODIMENTS

[0020] The technical solutions in the embodiments of the present application will be clearly and completely described below with reference to the drawings in the embodiments of the present application. Obviously, the described embodiments are only some embodiments of the present application, rather than all embodiments. Based on the embodiments of the present application, all other embodiments obtained by those skilled in the art without creative efforts belong to the scope of protection of the present application.

[0021] In the description of the present application, the terms "first" and "second" are used only for descriptive purposes and should not be construed as indicating or implying relative importance or implicitly specifying the quantity of the indicated technical features. Thus, the features defined with "first" and "second" may explicitly or implicitly include one or more of the said features. In the description of the present application, "a plurality of" means two or more, unless otherwise specifically defined.

[0022] In the present application, the phrase "in some embodiments" is used to mean "serving as an example, illustration, or explanation". Any embodiment described as "in some embodiments" in the present application is not necessarily construed as being more preferred or having more advantages than other embodiments. The following description is provided to enable any person skilled in the art to implement and use the present application. In the following description, details are set forth for purposes of explanation. It should be understood that those 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 to be accorded the widest scope consistent with the principles and features disclosed in the present application.

[0023] It should be noted that since the system of the embodiments 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, which is actually time information. It can be understood that in subsequent embodiments, if dimensions, quantities, positions, etc. are mentioned, they are all corresponding data existences for the computer device to process, and specific details are not elaborated here.

[0024] The embodiments of the present application provide a method and related devices for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography, aiming to automatically identify the nearshore river channel area and statistically analyze the erosion and deposition volume of the riverbed in this area by standardizing the long-term measured large cross-section data and combining the low-water level data set based on the long-term measured large cross-section data. The embodiments of the present application can efficiently, quickly, and simply identify and analyze the variation law of the erosion and deposition of the nearshore riverbed of an alluvial river, providing a scientific basis for river channel regulation and bank collapse prediction, and the following will be described in detail respectively.

[0025] In one embodiment, referring to Figure 1 , the method for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography includes: 101. Obtain the low-water river channel ranges of multiple cross-sections in a target river reach at multiple measurement times.

[0026] In an embodiment of the present application, the target river reach is a specified river reach, and multiple cross-sections of the river channel are included in the target river reach. A measurement instance refers to different measurements, that is, multiple measurements will be performed on multiple cross-sections in the target river reach to determine the low-flow channel ranges of multiple cross-sections in the target river reach at multiple measurement instances. The low-flow channel range refers to the channel range below the low water level of the cross-section, and this channel range can be calculated using the corresponding data obtained from the measurement of the cross-section.

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

[0028] In some embodiments of the present application, the low water level is determined through the following steps: obtaining the measured historical low water levels of multiple target hydrological stations in the target river reach; assuming that the water surface slope between hydrological stations remains unchanged, performing linear interpolation processing on the measured historical low water levels of multiple target hydrological stations to obtain the low water level of each cross-section in the target river reach, so as to make the determination of the low water level more convenient and fast.

[0029] In some embodiments of the present application, an example is given to illustrate the refined content of step 101. Specifically, step 101 may include: Step 1, standardize the long-sequence observation data of the measured large cross-section to form a data set including the starting distance X and the riverbed elevation Y.

[0030] 1.1, sort out the long-time sequence measured large cross-section observation data of the target river reach. For a specific cross-section, assuming the number of cross-sections is mm, the number of measurement times is nn, and the number of cross-section observation points is dd, then the river channel cross-section topography X, Y is a three-dimensional array: X(i, j, k), Y(i, j, k) Wherein, X is the distance of the observation point from the starting point of the cross-section (i.e., the starting distance); Y is the measured riverbed elevation of the observation point; i is the measurement cross-section number; j is the number of observations; k is the maximum number of observation points for each cross-section. Where i = 1, 2, …, mm; j = 1, 2, …, nn; k = 1, 2, …, dd.

[0031] 1.2, Standardize the measured large cross-section topographic data. Store the data of different cross-sections in the same measuring run in the same Excel file in the format of starting distance X and riverbed elevation Y. The data of different cross-sections are named after the cross-section name and stored in different Excel files.

[0032] Step 2, Determine the measured low-water level dataset of the large cross-section.

[0033] 2.1, Collect and organize the measured historical low-water levels of the target hydrological station in the target river reach. 2.2, Assume that the water surface slope between hydrological stations remains unchanged, and use linear interpolation to obtain the low-water levels of the measured large cross-section, forming a dataset including cross-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 cross-section number, i = 1, 2, …, mm.

[0034] Step 3, Identify the dry river channel area in the target river reach.

[0035] 3.1, Use the min function in software such as MATLAB to determine the lowest elevation (thalweg elevation Z) and its position (thalweg node d_thaw) of the measured large cross-section, i.e., Z = min(Y); d_thaw = find(Y == min(Y)).

[0036] 3.2, Extract the area where the measured riverbed elevation of the observation points in the measured large cross-section is less than the low-water level Z_low. According to the continuity of the cross-section nodes, determine the range where the thalweg is located as the dry river channel range, and record its left and right nodes as ans_node_left and ans_node_right. Take the difference between the starting distances of the left and right nodes as the dry river channel width ans_W, i.e., ans_W = X(node_right) - X(node_left).

[0037] 102. Based on the dry river channel range, determine the near-shore river channel area of the corresponding cross-section.

[0038] In the embodiments of the present application, the near-shore river channel refers to the part of the river channel close to the shore. The near-shore river channel area can be determined based on the dry river channel range of the corresponding cross-section. Since there are near-shore river channels on both the left and right banks of the river, the near-shore river channel area may include at least one of the left-bank near-shore river channel area and the right-bank near-shore river channel area.

[0039] In some embodiments of the present application, step 102 may include: based on the range of the dry river channel, determining the dry river channel width of the corresponding cross-section at the corresponding measurement time. For example, the span of the starting distance range of the dry river channel may be used as the corresponding dry river channel width; determining the critical value of the near-shore river channel area smaller than the dry river channel width. The critical value of the near-shore river channel area refers to the size of the range of the near-shore river channel within the dry river channel. For example, the product of the average value of the dry river channel widths of all cross-sections in all measurement times in the target river section and a preset ratio may be used as the critical value of the near-shore river channel area, and the preset ratio may be 20% for example; according to the boundary value of the dry river channel range and the critical value of the near-shore river channel area, determining the near-shore river channel area of the corresponding cross-section. For example, the sum value or difference value between the boundary value of the dry river channel range and the critical value of the near-shore river channel area, and the boundary value of the dry river channel range may be used as the near-shore river channel area of the corresponding cross-section.

[0040] In some embodiments of the present application, according to the boundary value of the dry river channel range and the critical value of the near-shore river channel area, determining the near-shore river channel area of the corresponding cross-section may include: according to the boundary value of the dry river channel range and the critical value of the near-shore river channel area, determining the near-shore river channel area of the corresponding cross-section at the corresponding measurement time. For example, the sum value or difference value between the boundary value of the dry river channel range and the critical value of the near-shore river channel area, and the boundary value of the dry river channel range may be used as the near-shore river channel area of the corresponding cross-section at the corresponding measurement time; taking the union of the near-shore river channel areas of the corresponding cross-section at multiple measurement times to obtain the near-shore river channel area of the corresponding cross-section, ensuring that the corresponding cross-section has a unique near-shore river channel area, so as to facilitate subsequent calculation of the cross-sectional area and the scouring and silting amount of the near-shore riverbed of each measurement time of the cross-section based on the same near-shore river channel area.

[0041] In some embodiments of the present application, an example is given to illustrate the refinement content of step 102. Specifically, step 102 may include: Step 4, determining the critical value of the near-shore river channel area of the target river section.

[0042] 4.1, using the for loop statement in MATLAB to identify the dry river channel areas of different measurement times of the same measured large cross-section and calculate the dry river channel width W(j).

[0043] 4.2, using the for loop statement in MATLAB to identify the dry river channel areas of different measurement times of each measured large cross-section and calculate the dry river channel width W(i).

[0044] 4.3, taking the average value W_Typ of the dry river channel widths of different measurement times of each measured large cross-section in the target river section as the representative river width of the target river section, and determining 0.2 times of this value as the critical value of the near-shore river channel area of the target river section, that is, WCV = 0.20 W_Typ.

[0045] Step 5, identify the nearshore river channel area of the target river reach.

[0046] 5.1, According to the left and right nodes ans_node_left and ans_node_right of the low-flow river channel determined in 3.2, determine the nearshore river channel range. Among them, the left-bank nearshore river channel area is [X(ans_node_left), X(ans_node_left)+WCV], and the right-bank nearshore river channel area is [X(ans_node_right)-WCV, X(ans_node_right)].

[0047] 5.2, Use the for loop statement in MATLAB to determine the nearshore river channel areas of different measurement times of the same measured cross-section, and take their union as the nearshore river channel area of the target cross-section.

[0048] 5.3, Use the for loop statement in MATLAB to determine the nearshore river channel areas of different cross-sections.

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

[0050] In the embodiments of the present application, since the cross-section is a vertical plane in the river channel, therefore, for each measurement time of each cross-section, the cross-sectional area under the nearshore river channel area of the cross-section can be determined.

[0051] In some embodiments of the present application, since the nearshore river channel area is for the low water level of the cross-section, therefore, step 103 may include: for each measurement time of each cross-section, among the multiple observation points of the corresponding cross-section at the corresponding measurement time, determine multiple target observation points under the nearshore river channel area. Among them, 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 river channel area; according to the multiple target observation points and the low water level of the corresponding cross-section, determine the cross-sectional geometric shape under the nearshore river channel area. For example, the multiple target observation points and the water level line of the low water level of the corresponding cross-section can be used as the boundary values of the cross-sectional geometric shape under the nearshore river channel area, so as to determine the cross-sectional geometric shape under the nearshore river channel area; based on the cross-sectional geometric shape, determine the cross-sectional area under the nearshore river channel area. For example, the cross-sectional area A under the nearshore river channel area can be calculated by the trapezoidal area summation method.

[0052] 104. According to the cross-sectional area, determine the nearshore riverbed scouring and silting volume of the target river reach.

[0053] In the embodiments of the present application, the nearshore riverbed scouring and silting volume of the target river reach refers to the scouring and silting volume within the nearshore riverbed range of the target river reach.

[0054] In some embodiments of the present application, step 104 may include: comparing the cross-sectional areas of the same cross-section at adjacent measurement times to obtain a cross-sectional area change value, which may be, for example, the difference between the cross-sectional areas of the cross-section at adjacent measurement times; determining the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times according to the cross-sectional area change value of adjacent cross-sections at adjacent measurement times, wherein the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times may have a positive correlation with the cross-sectional area change value of adjacent cross-sections at adjacent measurement times; determining the nearshore riverbed scouring and silting volume of the target river reach according to the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times.

[0055] In some embodiments of the present application, determining the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times according to the cross-sectional area change value of adjacent cross-sections at adjacent measurement times may include: obtaining the mileage from the dam of multiple cross-sections; comparing the mileage from the dam of adjacent cross-sections to obtain a mileage change value, which may be, for example, the difference between the mileage from the dam of adjacent cross-sections, that is, the distance between adjacent cross-sections; determining the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times according to the cross-sectional area change value of adjacent cross-sections at adjacent measurement times and the mileage change value of adjacent cross-sections. For example, multiplying the cross-sectional area change value of adjacent cross-sections at adjacent measurement times by the mileage change value of adjacent cross-sections to obtain the change amount of scouring and silting volume during the adjacent measurement times of adjacent cross-sections, and using it as the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times.

[0056] In some embodiments of the present application, determining the nearshore riverbed scouring and silting volume of the target river reach according to the nearshore riverbed scouring and silting volume of adjacent cross-sections at adjacent measurement times may include: for all cross-sections in the target river reach, determining the sum of the nearshore riverbed scouring and silting volumes of adjacent cross-sections at adjacent measurement times to obtain the total nearshore riverbed scouring and silting volume of the target river reach at adjacent measurement times. For example, summing up the change amounts of scouring and silting volume of all cross-sections in the target river reach during adjacent measurement times to obtain the total nearshore riverbed scouring and silting volume of the target river reach at adjacent measurement times; determining the nearshore riverbed scouring and silting volume of the target river reach according to the total nearshore riverbed scouring and silting volume. The nearshore riverbed scouring and silting volume of the target river reach may include the total nearshore riverbed scouring and silting volume of the target river reach at each measurement time, so as to analyze the change law of the nearshore riverbed scouring and silting volume of the target river reach based on the total nearshore riverbed scouring and silting volume of the target river reach at different measurement times, and provide a scientific basis for river channel regulation and bank collapse prediction.

[0057] In some embodiments of the present application, an example is given to illustrate the refined content of step 104. Specifically, step 104 may include: Step 6, calculating the scouring and silting volume of the nearshore riverbed area.

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

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

[0060] 6.3. Calculate the erosion and deposition volume of the nearshore riverbed between adjacent measured large cross-sections using the trapezoidal method, i.e., V i,j =(ΔA i,j + ΔA i+1,j) × (L i+1 - L i ) / 2.

[0061] 6.4. Calculate the erosion and deposition volume of the nearshore riverbed for adjacent measurement times in the target river reach, i.e., V j = ∑V i,j 。

[0062] Step 7. Plot the change process of the erosion and deposition volume of the nearshore riverbed in the target river reach.

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

[0064] In the technical solution disclosed in this embodiment, according to the low-flow river channel range of multiple cross-sections in the target river reach, the nearshore river channel area of the corresponding cross-section is determined, and then the cross-sectional area under the nearshore river channel area is determined. The erosion and deposition volume of the nearshore riverbed in the target river reach is determined based on the cross-sectional area, realizing the calculation of the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography. Compared with the prior art where relevant technicians determine the erosion and deposition volume of the nearshore riverbed based on past experience, the determination of the erosion and deposition volume of the nearshore riverbed is more efficient and accurate.

[0065] Moreover, by standardizing the measured data, constructing the low-water level data set, identifying the area, and calculating the erosion and deposition volume, the automatic identification of the alluvial river channel area and the automatic plotting of the change of the erosion and deposition characteristic quantity are realized. It solves the problems of strong subjectivity in defining the nearshore riverbed range and low calculation efficiency in the prior art, and has the advantages of data standardization, accurate calculation, and dynamic analysis.

[0066] Next, with reference to Figure 1 , an example description of the method for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography will be given. Specifically, to identify the change law of the erosion and deposition volume of the nearshore riverbed, the following steps are included: Step 1. Standardize the long-term observation data of the measured large cross-section.

[0067] Collect and organize the historical observation data of the long-term measured cross-sections of the target river reach. Store the data of different measurement times at the same cross-section in the same Excel file in the format of starting distance X and elevation Y, as shown in Table 1 below. The data of different cross-sections are named after the cross-section name and stored in different Excel files.

[0068] Table 1

[0069] Step 2: Construct a dataset of the low water levels of the measured cross-sections.

[0070] Collect and organize the low water level data of typical hydrological stations or typical cross-sections in the target river reach. Assume that the water surface slope between hydrological stations or typical cross-sections remains unchanged, and use linear interpolation to obtain the low water levels of the measured cross-sections, forming a dataset containing cross-section name Csname, distance from the dam L, and low water level Z_low, as shown in Table 2 below.

[0071] Table 2

[0072] Step 3: Identify the low-flow river channel area in the target river reach.

[0073] Taking the measured cross-section data of Jing 175 cross-section in 2003 as an example, the process of identifying the low-flow river channel area is described. First, use the min function in MATLAB to determine the lowest elevation of the measured cross-section, that is, the thalweg elevation Z = 1.40, and use the find statement to determine the cross-section node number d_thaw where the thalweg is located = 11; then, use the find statement in MATLAB to extract the area where the cross-section riverbed elevation is less than the low water level Z_low = 22.09 m of the cross-section. According to the continuity of the cross-section nodes, the interval where the thalweg is located is determined as the low-flow river channel range, and its left and right nodes ans_node_left = 10, ans_node_right = 43, including the cross-section node 11 where the thalweg is located; finally, take the difference between the starting distances of the left and right nodes as the low-flow river channel width ans_W, that is, ans_W = X(43) - X(10) = 1024.

[0074] Step 4: Determine the critical values of the near-shore river channel area in the target river reach.

[0075] First, use the for loop statement in MATLAB to identify the low-flow channel area of different measurement times during 2001 - 2021 at Section Jing 175 and calculate the low-flow channel width W(1,j); then, use the for loop statement in MATLAB to identify the low-flow channel areas of different measurement times of each measured large cross-section between Section Jing 175 and Section Jing 179+1 and calculate the low-flow channel width W(i,j); finally, take the average value W_Typ of the low-flow channel widths of different measurement times of each measured large cross-section from Section Jing 175 to Section Jing 179+1 as the representative channel width of the target reach, W_Typ = mean(W(i,j)) = 864.55, and determine 20% of this value as the critical value of the near-shore channel area of the target reach, that is, WCV = 0.20 W_Typ = 173.

[0076] Step 5: Identify the near-shore channel area of the target reach.

[0077] First, based on the left and right nodes of the low-flow channel of the cross-section data of Section Jing 175 in 2003 determined in Step 3, ans_node_left = 10 and ans_node_right = 43, determine the near-shore channel range, where the left-bank near-shore 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, use the for loop statement in MATLAB to determine the near-shore channel areas of different measurement times at Section Jing 175 and take their union as the near-shore channel area of the target cross-section, where the range of the left-side near-shore channel area is [45.48, 253, 28]; finally, use the for loop statement in MATLAB to determine the near-shore channel areas of different cross-sections. The near-shore channel areas of typical cross-sections in the target reach are as Figure 2 shown in

[0078] Step 6: Calculate the sediment erosion and deposition volume of the near-shore riverbed.

[0079] First, use the trapezoidal area method to calculate the cross-sectional area A of the near-shore channel area below the low-water level. Use the for loop statement in MATLAB to calculate the cross-sectional areas A of the near-shore channel areas of the same measured large cross-section at different measurement times j and their change values △A j = A j - A j+1 ; then, use the for loop statement in MATLAB to calculate the cross-sectional areas A of the near-shore channel areas of different measurement times of each measured cross-section i,j and their change values △A i,j; Furthermore, the trapezoidal method is used to calculate the scouring and silting volume of the nearshore riverbed between adjacent measured large cross-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 volumes of the nearshore riverbed for adjacent measurement times in the target reach is used as the scouring and silting volume of the nearshore riverbed in the reach, that is, V j =∑V i,j .

[0080] Step 7, draw the variation process of the scouring and silting volume of the nearshore riverbed in the target reach.

[0081] Based on the calculation results in Step 6, use the plot function in MATLAB to draw the time variation process of the scouring and silting volume of the nearshore riverbed. The variation process of the scouring and silting volume of the left nearshore riverbed from Section Jing 175 to Section Jing 179+1 from 2003 to 2021 is as Figure 3 shown.

[0082] Based on the measured large cross-section data with a long time series, taking the left and right bank boundaries of the dry-season river channel where the thalweg point is located as the boundaries of the region, taking 20% of the average width of the dry-season river channel in the target reach as the regional critical value, taking the range of the regional critical value in the river channel area near the left and right bank boundaries in the dry-season river channel as the region, and realizing the automatic recognition of the region and the automatic drawing of the variation of the nearshore scouring and silting characteristic quantities through the MATLAB program. This technology is based on the measured large cross-section data, with less difficulty in obtaining data, lower cost, and high accuracy. It can efficiently, quickly, and accurately identify the temporal variation law of the scouring and silting volume in the nearshore riverbed area, which is of great significance for predicting the characteristics of river channel evolution and formulating appropriate river (navigation) regulation plans, etc.

[0083] To better implement the method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography in the embodiments of the present application, on the basis of the method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography, an apparatus for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography is further provided in the embodiments of the present application, as Figure 4 shown. The apparatus 400 for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography includes: The first acquisition module 401 is used to acquire the dry-season river channel ranges of multiple cross-sections in the target reach at multiple measurement times; The first determination module 402 is used to determine the nearshore river channel area of the corresponding cross-section based on the dry-season river channel range; The second determination module 403 is used to determine the cross-sectional area under the nearshore river channel area for each measurement time of each cross-section; The third determination module 404 is used to determine the scouring and silting volume of the nearshore riverbed in the target reach according to the cross-sectional area.

[0084] The embodiment of the present application further provides a computer device, which integrates any of the devices for determining the erosion and deposition volume of the nearshore riverbed based on the river channel cross-section topography provided by the embodiment of the present application. As Figure 5 shown, it shows a schematic structural diagram of the computer device involved in the embodiment of the present application. Specifically: The computer device may include components such as a processor 501 with one or more processing cores, a memory 502 with one or more computer-readable storage media, a power supply 503, and an input unit 504. Those skilled in the art can understand that Figure 5 the computer device structure shown in does not limit the computer device. It may include more or fewer components than shown, or combine certain components, or have different component arrangements. Among them: The processor 501 is the control center of the computer device. It connects various parts of the entire computer device through various interfaces and lines. By running or executing software programs and / or modules stored in the memory 502, and by calling the data stored in the memory 502, it executes various functions of the computer device and processes data, thereby monitoring the computer device as a whole. Optionally, the processor 501 may include one or more processing cores; preferably, the processor 501 may integrate an application processor and a modem processor. Among them, the application processor mainly processes the operating system, user interface, application programs, etc., and the modem processor mainly processes wireless communication. It can be understood that the above modem processor may not be integrated into the processor 501.

[0085] The memory 502 can be used to store software programs and modules. The processor 501 executes various functional applications and data processing by running the software programs and modules stored in the memory 502. The memory 502 may mainly include a program storage area and a data storage area. Among them, the program storage area can store the operating system, application programs required for at least one function (such as a sound playback function, an image playback function, etc.); the data storage area can store data created according to the use of the computer device. In addition, the memory 502 may include a high-speed random access memory, and may also include a non-volatile memory, such as at least one magnetic disk storage device, a flash memory device, or other volatile solid-state storage devices. Correspondingly, the memory 502 may also include a memory controller to provide the processor 501 with access to the memory 502.

[0086] The computer device further includes a power supply 503 for powering each component. Preferably, the power supply 503 can be logically connected to the processor 501 through a power management system, so as to manage functions such as charging, discharging, and power consumption management through the power management system. The power supply 503 may further include any components such as one or more DC or AC power supplies, a recharge system, a power failure detection circuit, a power converter or inverter, and a power status indicator.

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

[0088] Although not shown, the computer device may further include a display unit, etc., which will not be elaborated here. Specifically, in this embodiment, the processor 501 in the computer device will load the executable files corresponding to the processes of one or more application programs into the memory 502 according to the following instructions, and the processor 501 will run the application programs stored in the memory 502 to implement various functions as follows: Obtain the low-flow river channel ranges of multiple cross-sections in a target river reach at multiple measurement times; based on the low-flow river channel ranges, determine the near-shore river channel areas of the corresponding cross-sections; for each measurement time of each cross-section, determine the cross-sectional area under the near-shore river channel area; according to the cross-sectional area, determine the scouring and silting volume of the near-shore riverbed in the target river reach.

[0089] Those of ordinary skill in the art can understand that all or part of the steps in the various methods of the above embodiments can be completed through instructions, or through instructions to control related hardware. The instructions can be stored in a computer-readable storage medium and loaded and executed by a processor.

[0090] Therefore, an embodiment of the present application provides a computer-readable storage medium, which may include: read-only memory (ROM, Read Only Memory), random access memory (RAM, Random Access Memory), a magnetic disk or an optical disc, etc. A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps in any method for determining the scouring and silting volume of the near-shore riverbed based on the river channel cross-section topography provided by the embodiments of the present application. For example, the computer program loaded by the processor can execute the following steps: Obtain the low-flow river channel ranges of multiple cross-sections in a target river reach at multiple measurement times; based on the low-flow river channel ranges, determine the near-shore river channel areas of the corresponding cross-sections; for each measurement time of each cross-section, determine the cross-sectional area under the near-shore river channel area; according to the cross-sectional area, determine the scouring and silting volume of the near-shore riverbed in the target river reach.

[0091] Embodiments of the present application also provide a computer program product or a computer program. The computer program product or the computer program includes computer instructions, and the computer instructions are stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium, and the processor executes the computer instructions, so that the electronic device executes to implement the method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography as described in any one of the above, for example: Obtain the dry-season river channel ranges of multiple cross-sections in a target river reach; based on the dry-season river channel ranges, determine the nearshore river channel areas corresponding to the cross-sections; for each survey of each cross-section, determine the cross-sectional area under the nearshore river channel area; and determine the scouring and silting volume of the nearshore riverbed of the target river reach according to the cross-sectional area.

[0092] In the above embodiments, the descriptions of the respective embodiments have their own focuses. For parts not detailed in a certain embodiment, reference may be made to the detailed descriptions of other embodiments above, and details will not be repeated here.

[0093] In specific implementation, the above units or structures can be implemented as independent entities, or can be combined arbitrarily to be implemented as the same or several entities. For the specific implementation of the above units or structures, reference may be made to the method embodiments above, and details will not be repeated here.

[0094] For the specific implementation of each of the above operations, reference may be made to the foregoing embodiments, and details will not be repeated here.

[0095] The above has introduced in detail a method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography and related devices provided by the embodiments of the present application. Specific examples are used in this article to elaborate on the principle and implementation manner of the present application. The description of the above embodiments is only used to help understand the method and its core idea of the present application; at the same time, for those skilled in the art, according to the idea of the present application, there will be changes in the specific implementation manner and application scope. In summary, the content of this specification should not be construed as a limitation to the present application.

Claims

1. A method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of the river channel, characterized in that The method includes: Obtaining the low-flow river channel ranges of multiple cross-sections in a target river reach at multiple measurement times; Based on the low-flow river channel ranges, determining the near-shore river channel areas of the corresponding cross-sections; For each measurement time of each cross-section, determining the cross-sectional area under the near-shore river channel area; Based on the cross-sectional area, determining the scouring and silting volume of the near-shore riverbed of the target river reach.

2. The method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography as described in claim 1, wherein, The determining the scouring and silting volume of the near-shore riverbed of the target river reach based on the cross-sectional area includes: Comparing the cross-sectional areas of the same cross-section at adjacent measurement times to obtain a cross-sectional area change value; Based on the cross-sectional area change values of adjacent cross-sections at the adjacent measurement times, determining the scouring and silting volume of the near-shore riverbed of the adjacent cross-sections at the adjacent measurement times; Based on the scouring and silting volume of the near-shore riverbed of the adjacent cross-sections at the adjacent measurement times, determining the scouring and silting volume of the near-shore riverbed of the target river reach.

3. The method for determining the erosion and deposition volume of the nearshore riverbed based on the cross-sectional topography of the river channel according to claim 2, characterized in that, The determining the scouring and silting volume of the near-shore riverbed of the adjacent cross-sections at the adjacent measurement times based on the cross-sectional area change values of the adjacent cross-sections at the adjacent measurement times includes: Obtaining the distances from the dams of multiple cross-sections; Comparing the distances from the dams of adjacent cross-sections to obtain a distance-from-dam change value; Based on the cross-sectional area change values of the adjacent cross-sections at the adjacent measurement times and the distance-from-dam change values of the adjacent cross-sections, determining the scouring and silting volume of the near-shore riverbed of the adjacent cross-sections at the adjacent measurement times.

4. The method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of the river channel according to claim 2, wherein, The determining the scouring and silting volume of the near-shore riverbed of the target river reach based on the scouring and silting volume of the near-shore riverbed of the adjacent cross-sections at the adjacent measurement times includes: For all cross-sections in the target river reach, determining the sum of the scouring and silting volume of the near-shore riverbed of the adjacent cross-sections at the adjacent measurement times to obtain the total scouring and silting volume of the near-shore riverbed of the target river reach at the adjacent measurement times; Based on the total scouring and silting volume of the near-shore riverbed, determining the scouring and silting volume of the near-shore riverbed of the target river reach.

5. The method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of the river channel according to claim 1, characterized in that, The determining the cross-sectional area under the near-shore river channel area for each measurement time of each cross-section includes: For each measurement time of each cross-section, among the multiple observation points of the corresponding cross-section at the corresponding measurement time, determining multiple target observation points under the near-shore river channel area; Based on the multiple target observation points and the low-water level of the corresponding cross-section, determining the cross-sectional geometric shape under the near-shore river channel area; Based on the cross-sectional geometric shape, determining the cross-sectional area under the near-shore river channel area.

6. The method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of the river channel according to claim 1, wherein The determining the near-shore river channel area of the corresponding cross-section based on the low-flow river channel range includes: Based on the low-flow river channel range, determining the low-flow river channel width of the corresponding cross-section at the corresponding measurement time; Taking 0.2 times of the low-flow river channel width as the critical value of the near-shore river channel area; Based on the boundary values of the low-flow river channel range and the critical value of the near-shore river channel area, determining the near-shore river channel area of the corresponding cross-section.

7. The method for determining the scouring and silting volume of the nearshore riverbed based on the cross-sectional topography of the river channel according to claim 6, wherein The determining the near-shore river channel area of the corresponding cross-section based on the boundary values of the low-flow river channel range and the critical value of the near-shore river channel area includes: Based on the boundary values of the low-flow river channel range and the critical value of the near-shore river channel area, determining the near-shore river channel area of the corresponding cross-section at the corresponding measurement time; Taking the union of the near-shore river channel areas of the corresponding cross-section at multiple measurement times to obtain the near-shore river channel area of the corresponding cross-section.

8. The method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography as claimed in claim 1, wherein The obtaining the low-flow river channel ranges of multiple cross-sections in a target river reach at multiple measurement times includes: Obtain the observation data of each cross-section at each measurement time, where the observation data includes the starting distance and riverbed elevation of multiple observation points in the corresponding cross-section at the corresponding measurement time; Based on the riverbed elevation, determine the thalweg points among the multiple observation points; Use the low water level of the corresponding cross-section to divide the multiple observation points sorted in sequence according to the starting distance, and obtain multiple observation point intervals; According to the observation point interval where the thalweg point is located, determine the low water channel range of the corresponding cross-section at the corresponding measurement time.

9. A computer device, characterized in that, The computer device includes: One or more processors; A memory; And one or more application programs, where the one or more application programs are stored in the memory and are configured to be executed by the processor to implement the steps in the method for determining the scouring and silting volume of the nearshore riverbed based on the river channel cross-section topography according to any one of claims 1 to 8.

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

Images