Method for determining river channel longitudinal profile based on river channel cross-section topography, related equipment

Abstract

The present application discloses a method for determining a river channel longitudinal profile based on the river channel cross-section topography and related devices. The method includes: obtaining the main channel cross-section data of multiple cross-sections in the target river reach, where the main channel cross-section data includes the main channel cross-section area of each cross-section at the corresponding cross-section floodplain water level; for each cross-section, among multiple preset water level values greater than the thalweg elevation of the cross-section, determining a target water level value, where the river channel range interval where the thalweg is located at the target water level value of the cross-section is the target interval, and the river channel cross-section area within the target interval is greater than or equal to the product of the main channel cross-section area and a preset ratio value; taking the target water level value as the deep trough elevation of the corresponding cross-section; using the deep trough elevations of multiple cross-sections in the target river reach to determine the river channel longitudinal profile of the target river reach. The present application can reduce the sawteeth in the river channel longitudinal profile, so as to facilitate accurately identifying the longitudinal adjustment law of the river channel and reduce the level requirements for relevant technical personnel.

Description

Technical Field

[0001] This application relates to the fields of river evolution simulation, river regulation and protection, and specifically relates to a method for determining a river longitudinal profile based on river cross-section topography and related equipment. Background Art

[0002] The longitudinal profile of an alluvial river is an important tool for studying the morphology of river thalwegs and their dynamic evolution, and is an important research content in river regulation and protection work. As an important index for judging the stability of a river channel and its longitudinal evolution trend, the morphological characteristics of the longitudinal profile of an alluvial river directly affect the hydrodynamic characteristics of the water flow, can effectively display the process of riverbed erosion and deposition, and mastering its dynamic change law has important application value for revealing the evolution of river thalweg morphology, river (navigation) channel design and maintenance, flood prediction, etc.

[0003] In related technologies, a thalweg longitudinal profile can be drawn through cross-section drawing software to reflect the longitudinal changes of a river channel. However, due to the large variation in the elevation of thalweg points in the cross-section, the obtained thalweg longitudinal profile is mostly significantly serrated, which interferes with accurately identifying the longitudinal adjustment law of the river channel and requires a relatively high level of relevant technical personnel. Summary of the Invention

[0004] The embodiments of this application provide a method for determining a river longitudinal profile based on river cross-section topography and related equipment, which can reduce the serrations in the river longitudinal profile, facilitate accurately identifying the longitudinal adjustment law of the river channel, and reduce the level requirement for relevant technical personnel.

[0005] On the one hand, this application provides a method for determining a river longitudinal profile based on river cross-section topography, and the method includes:

[0006] Obtain the main channel cross-section data of multiple cross-sections in the target river reach, where the main channel cross-section data includes the main channel cross-section area of each cross-section at the corresponding cross-section floodplain water level;

[0007] For each cross-section, among multiple preset water level values greater than the elevation of the thalweg point of the cross-section, determine the target water level value. Wherein, the river channel range interval where the thalweg point of the cross-section is located at the target water level value is the target interval, and the cross-section area of the river channel within the target interval is greater than or equal to the product of the main channel cross-section area and the preset proportional value;

[0008] Take the target water level value as the thalweg elevation of the corresponding cross-section;

[0009] Use the thalweg elevations of multiple cross-sections in the target river reach to determine the river longitudinal profile of the target river reach.

[0010] In some embodiments, the obtaining the main channel cross-section data of multiple cross-sections in the target river reach includes:

[0011] For each cross-section, obtain the range interval of the main channel of the cross-section;

[0012] Determine the area of the region of the cross-section within the corresponding main channel range interval and at the corresponding cross-section flood-level, and use it as the main channel cross-sectional area of the cross-section at the corresponding cross-section flood-level.

[0013] In some embodiments, the obtaining the range interval of the main channel of the cross-section includes:

[0014] Obtain the observation data of the cross-section, where the observation data includes the starting distances and riverbed elevations of multiple observation points in the cross-section;

[0015] Based on the riverbed elevation, determine the thalweg points among the multiple observation points;

[0016] Use the low-water level of the cross-section to divide the multiple observation points sorted in sequence according to the starting distances, and obtain multiple first observation point intervals;

[0017] Take the first observation point interval where the thalweg point is located as the range interval of the dry-season river channel of the cross-section;

[0018] Based on the range interval of the dry-season river channel, determine the range interval of the main channel of the cross-section.

[0019] In some embodiments, the determining the range interval of the main channel of the cross-section based on the range interval of the dry-season river channel includes:

[0020] Among the starting distances of the multiple observation points in the range interval of the dry-season river channel, determine the maximum starting distance and the minimum starting distance;

[0021] Among the multiple observation points sorted in sequence according to the starting distances, detect whether there is a target observation point that satisfies the first preset condition or the second preset condition. Among them, the first preset condition includes that the target observation point does not exist in the range interval of the dry-season river channel, the target observation point is adjacent to the observation point where the maximum starting distance is located among the multiple observation points sorted in sequence according to the starting distances, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the maximum starting distance is located satisfy the first preset magnitude relationship. The second preset condition includes that the target observation point does not exist in the range interval of the dry-season river channel, the target observation point is adjacent to the observation point where the minimum starting distance is located among the multiple observation points sorted in sequence according to the starting distances, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the minimum starting distance is located satisfy the second preset magnitude relationship;

[0022] If there is the target observation point, add the target observation point to the range interval of the dry-season river channel to obtain an expanded range interval of the dry-season river channel;

[0023] Based on the expanded range interval of the dry river channel, determine the main channel range interval of the cross-section.

[0024] In some embodiments, the observed data of the cross-section includes the starting distances and riverbed elevations of multiple observation points in the cross-section. Among the multiple observation points, there is a thalweg point, which is determined based on the riverbed elevation. Before determining the target water level value among multiple preset water level values greater than the elevation of the thalweg point of the cross-section, it further includes:

[0025] For each of the preset water level values, use the preset water level value to divide the multiple observation points sorted in sequence according to the starting distance, and obtain multiple second observation point intervals;

[0026] Take the second observation point interval where the thalweg point is located as the river channel range interval where the thalweg point is located under the preset water level value of the cross-section.

[0027] In some embodiments, multiple preset water level values greater than the elevation of the thalweg point of the cross-section are determined through the following steps:

[0028] Determine the sum of the elevation of the thalweg point and the preset elevation to obtain the target elevation;

[0029] Increase the target elevation sequentially according to the preset step size to obtain multiple of the preset water level values.

[0030] In some embodiments, the method of determining the longitudinal profile of the river channel using the thalweg elevations of multiple cross-sections in the target river reach includes:

[0031] Perform moving average processing on the thalweg elevations of multiple cross-sections in the target river reach in sequence according to the cross-section distance from the dam to obtain the processed thalweg elevations of multiple cross-sections;

[0032] Connect the processed thalweg elevations of multiple cross-sections according to the cross-section distance from the dam to obtain the longitudinal profile of the river channel of the target river reach.

[0033] In some embodiments, the cross-section floodplain water level is determined through the following steps:

[0034] Obtain the measured historical floodplain water levels of multiple target hydrological stations in the target river reach;

[0035] Perform linear interpolation processing on the measured historical floodplain water levels of multiple target hydrological stations to obtain the cross-section floodplain water level of each cross-section in the target river reach.

[0036] On the other hand, an embodiment of the present application provides a device for determining the longitudinal profile of a river channel based on the river channel cross-section topography, including:

[0037] A first acquisition module, configured to acquire the main channel cross-section data of multiple cross-sections in a target river reach, where the main channel cross-section data includes the main channel cross-section area of each cross-section at the corresponding cross-section floodplain water level;

[0038] A first determination module, configured to determine, for each cross-section, a target water level value from multiple preset water level values greater than the thalweg point elevation of the cross-section, where the river channel range interval where the thalweg point of the cross-section is located at the target water level value is a target interval, and the river channel cross-section area within the target interval is greater than or equal to the product of the main channel cross-section area and a preset ratio value;

[0039] A second determination module, configured to use the target water level value as the deep channel elevation of the corresponding cross-section;

[0040] A third determination module, configured to determine the river channel longitudinal profile of the target river reach by using the deep channel elevations of multiple cross-sections in the target river reach.

[0041] On the other hand, the present application also provides a computer device, where the computer device includes:

[0042] One or more processors;

[0043] A memory; and

[0044] One or more applications, where 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 river channel longitudinal profile based on the river channel cross-section topography.

[0045] On the other hand, the present application also 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 river channel longitudinal profile based on the river channel cross-section topography.

[0046] 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 river channel longitudinal profile based on the river channel cross-section topography as described in any one of the above.

[0047] The method and related equipment for determining the longitudinal profile of a river channel based on the topographic profile of the river channel cross-section provided by the embodiments of the present application obtain the main channel cross-section data of multiple cross-sections in the target river reach. The main channel cross-section data includes the main channel cross-section area of each cross-section at the corresponding cross-section floodplain water level. For each cross-section, among multiple preset water level values greater than the elevation of the thalweg point of the cross-section, a target water level value is determined. The thalweg range interval where the thalweg point of the cross-section is located at the target water level value is the target interval, and the cross-section area of the river channel within the target interval is greater than or equal to the product of the main channel cross-section area and the preset ratio value. The target water level value is used as the deep channel elevation of the corresponding cross-section. The longitudinal profile of the target river reach is determined by using the deep channel elevations of multiple cross-sections in the target river reach. The embodiments of the present application determine the target water level value with the corresponding cross-section area of the river channel based on the product of the main channel cross-section area of the cross-section at the corresponding cross-section floodplain water level and the preset ratio value, and use it as the deep channel elevation of the corresponding cross-section, and then determine the longitudinal profile of the target river reach. Compared with the thalweg longitudinal profile drawn based on the elevation of the thalweg point of the cross-section, the longitudinal profile of the river channel has fewer zigzags, which is convenient for accurately identifying the longitudinal adjustment law of the river channel and reducing the level requirements for relevant technical personnel. Description of the Drawings

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

[0049] Figure 1 It is a schematic diagram of an embodiment of the method for determining the longitudinal profile of a river channel based on the topographic profile of the river channel cross-section provided by the embodiments of the present application;

[0050] Figure 2 It is a schematic diagram of the thalweg point elevation and deep channel elevation of the cross-sections at different cross-section distances from the dam in the target river reach provided by the embodiments of the present application;

[0051] Figure 3 It is a schematic diagram of the longitudinal profile of the target river reach at different measurement times provided by the embodiments of the present application;

[0052] Figure 4 It is a schematic diagram of the structure of an embodiment of the device for determining the longitudinal profile of a river channel based on the topographic profile of the river channel cross-section provided by the embodiments of the present application;

[0053] Figure 5 It is a schematic diagram of the terminal structure of an embodiment of the computer device provided by the embodiments of the present application.

[0054] Among them, Figure 2 、 3 are color pictures to facilitate distinguishing different objects in the same picture. Detailed implementation manners

[0055] The following will clearly and completely describe the technical solutions in the embodiments of the present application with reference to the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part rather than all of the embodiments of the present application. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative efforts fall within the protection scope of the present application.

[0056] In the description of the present application, the terms "first" and "second" are only used for descriptive purposes and cannot 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 described features. In the description of the present application, "a plurality of" means two or more unless otherwise specifically defined.

[0057] In the present application, the phrase "in some embodiments" is used to mean "serving as an example, illustration, or description". 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. In order for any person skilled in the art to implement and use the present application, the following description is provided. Details are set forth for the purpose of explanation in the following description. It should be understood that those of ordinary skill in the art can recognize that the present application can be implemented without 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.

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

[0059] The embodiments of the present application provide a method and related equipment for determining the longitudinal profile of a river based on the cross-sectional topography of the river. The aim is to determine the target water level value corresponding to the river channel cross-sectional area based on the product of the main channel cross-sectional area at the corresponding cross-sectional floodplain water level and a preset ratio value, and use it as the deep channel elevation of the corresponding cross-section. The connection of the deep channel elevations of each cross-section along the target river section is used as the longitudinal profile line of the river channel, and the automatic recognition of the deep channel elevation of the cross-section and the automatic drawing of the longitudinal profile line are realized through programs such as MATLAB, which can efficiently, quickly, and simply identify the temporal variation law of the longitudinal profile of the alluvial river channel. The following will be described in detail respectively.

[0060] In one embodiment, referring to Figure 1 , the method for determining the longitudinal profile of a river based on the cross-sectional topography of the river includes:

[0061] 101. Obtain the main channel cross-sectional data of multiple cross-sections in the target river section, where the main channel cross-sectional data includes the main channel cross-sectional area of each cross-section at the corresponding cross-sectional floodplain water level.

[0062] In the embodiments of the present application, the target river section is a specified river section, and the target river section includes multiple cross-sections of the river. In the main channel cross-sectional data, each cross-section corresponds to a main channel cross-sectional area at the corresponding cross-sectional floodplain water level. Among them, since each cross-section has a river channel main channel, and the main channel has a certain width, the corresponding main channel cross-sectional area can be determined by combining the corresponding cross-sectional floodplain water level of the cross-section.

[0063] In some embodiments of the present application, step 101 may include: for each cross-section, obtain the main channel range interval of the cross-section. Among them, the cross-section may include multiple observation points, and each observation point has a corresponding starting distance and riverbed elevation. The maximum value of the main channel range interval may be the maximum value of the starting distances of multiple observation points in the main channel, and the minimum value of the main channel range interval may be the minimum value of the starting distances of multiple observation points in the main channel. At this time, the main channel range interval is the starting distance range of the main channel of the cross-section; determine the area of the region of the cross-section within the corresponding main channel range interval and at the corresponding cross-sectional floodplain water level, and use it as the main channel cross-sectional area of the cross-section at the corresponding cross-sectional floodplain water level. Among them, the main channel cross-sectional area can be calculated by the trapezoidal area summation method of adjacent observation points.

[0064] In some embodiments of the present application, obtaining the main channel range interval of a cross-section may include: obtaining the observation data of the cross-section, where the observation data includes the starting distances and riverbed elevations of multiple observation points in the cross-section, and the starting distances and riverbed elevations 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 cross-section to divide the multiple observation points sorted in ascending order of starting distance, obtaining multiple first observation point intervals, for example, the multiple observation points sorted in ascending order of 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 deeper than the low water level of the cross-section, and each subsequence is a first observation point interval; taking the first observation point interval where the thalweg point is located as the low water channel range interval of the cross-section, that is, the starting distance range of the low water channel of the cross-section; determining the main channel range interval of the cross-section based on the low water channel range interval.

[0065] 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 the 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 rapid.

[0066] In some embodiments of the present application, determining the main channel range interval of the cross-section based on the low water channel range interval may include: determining the maximum starting distance and the minimum starting distance among the starting distances of multiple observation points in the low water channel range interval; detecting whether there are target observation points that meet the first preset condition or the second preset condition among the multiple observation points sorted in ascending order of starting distance; if there are target observation points, adding the target observation points to the low water channel range interval to obtain an expanded low water channel range interval; determining the main channel range interval of the cross-section based on the expanded low water channel range interval, for example, the expanded low water channel range interval can be directly used as the main channel range interval of the cross-section, so that the determined main channel range interval is more accurate.

[0067] Among them, based on the characteristic that the riverbed elevations on both banks increase in the area outside the dry river channel range within the main channel range interval, the first preset condition is set as follows: The first preset condition includes that the target observation point does not exist in the dry river channel range interval, the target observation point is adjacent to the observation point with the largest starting distance among multiple observation points sorted in sequence according to the starting distance, and a first preset size relationship is satisfied between the riverbed elevation of the target observation point and the riverbed elevation of the observation point with the largest starting distance. The first preset size relationship can be, for example, that the riverbed elevation of the target observation point is greater than the riverbed elevation of the observation point with the largest starting distance, that is, it conforms to the above characteristics (within the main channel range interval, the riverbed elevation on the right bank of the area outside the dry river channel range interval increases). Therefore, the target observation point can be added to the dry river channel range interval to obtain a more accurate main channel range interval.

[0068] Similarly, the second preset condition includes that the target observation point does not exist in the dry river channel range interval, the target observation point is adjacent to the observation point with the smallest starting distance among multiple observation points sorted in sequence according to the starting distance, and a second preset size relationship is satisfied between the riverbed elevation of the target observation point and the riverbed elevation of the observation point with the smallest starting distance. The second preset size relationship can be, for example, that the riverbed elevation of the target observation point is greater than the riverbed elevation of the observation point with the smallest starting distance, that is, it conforms to the above characteristics (within the main channel range interval, the riverbed elevation on the left bank of the area outside the dry river channel range interval increases). Therefore, the target observation point can be added to the dry river channel range interval to obtain a more accurate main channel range interval.

[0069] After adding the target observation point to the dry river channel range interval, the starting distances of multiple observation points in the expanded dry river channel range interval still satisfy the ascending order.

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

[0071] In some embodiments of the present application, an example is given to illustrate the refinement content of step 101. Specifically, step 101 may include:

[0072] Step 1, standardize the long-sequence observation data of the measured large cross-section to form a data set including the starting distance X, the riverbed elevation Y, and the mileage L from the dam.

[0073] 1.1, Sort out the measured large cross-section observation data of the target river section for a long time series. For the data of a certain measurement, if the number of cross-sections is set to mm and the number of cross-section observation points is set to dd, then the river cross-section topography X and Y are two-dimensional arrays:

[0074] X(i,j), Y(i,j), L(i)

[0075] In the formula, 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; L is the distance of the measured large cross-section from the dam site (i.e., the mileage from the dam); i is the measurement cross-section number; j is the maximum value of the number of observation points for each cross-section. Among them, i = 1, 2, …, mm; j = 1, 2, …, dd.

[0076] 1.2, Standardize the measured large cross-section topographic data. Store the data of different cross-sections in the same measurement in the same Excel file in the format of starting distance X and riverbed elevation Y. At the same time, store the mileage L from the dam of different cross-sections in the same measurement in different worksheets (sheets) of the above Excel file, and name the data of different measurements by the measurement time and store them in different Excel files.

[0077] Step 2, Determine the typical water level data set of the measured large cross-section.

[0078] 2.1, Collect and sort out the measured historical low water level and floodplain water level data of the target hydrological station in the target river section;

[0079] 2.2, Assume that the water surface slope between hydrological stations remains unchanged, and use linear interpolation to obtain the characteristic water levels of the measured large cross-section, forming a data set including cross-section name, mileage L from the dam, low water level Z_low and floodplain water level Z_bankfull.

[0080] Step 3, Identify the main channel range of the measured large cross-section.

[0081] 3.1, Use the find statement in software such as MATLAB to determine the lowest riverbed elevation and location of the measured large cross-section. The lowest riverbed elevation is the thalweg elevation, denoted as Z. The cross-section node where the thalweg is located is denoted as d_thaw, that is, Z = min(Y), d_thaw = find(Y = min(Y)).

[0082] 3.2, Extract the area where the cross-section riverbed elevation is less than the low water level Z_low of the cross-section. According to the continuity of the cross-section nodes, divide the area where the cross-section riverbed elevation is less than the low water level Z_low into different sub-intervals, and then determine the interval where the thalweg is located as the low water channel range.

[0083] 3.3. Based on the characteristic that the elevations on both the left and right banks increase in the area outside the dry river channel within the main channel range, using the while conditional statement in software such as MATLAB, in the order of increasing starting distance, taking the increase in the elevations of adjacent nodes on both the left and right sides of the dry river channel as the standard, expand the dry river channel range determined in step 3.2 until the elevations on both the left and right banks decrease, and determine the range between the left and right nodes at this time as the main channel range.

[0084] 3.4. Calculate the cross-sectional area A of the main channel at the floodplain water level by using the trapezoidal area summation method of adjacent observation points.

[0085] 102. For each cross-section, among multiple preset water level values greater than the elevation of the thalweg point of the cross-section, determine the target water level value. Among them, the river channel range interval where the thalweg point of the cross-section is located at the target water level value is the target interval, and the cross-sectional area of the river channel within the target interval is greater than or equal to the product of the cross-sectional area of the main channel and the preset ratio value.

[0086] 103. Take the target water level value as the thalweg elevation of the corresponding cross-section.

[0087] In the embodiments of the present application, the thalweg elevation refers to the riverbed elevation of the thalweg point. Multiple preset water level values greater than the thalweg elevation of the cross-section can be set based on actual needs. In order to make the determined river channel longitudinal profile more accurate, it is necessary to determine the target water level value that can reflect the actual characteristics of the river channel longitudinal profile among multiple preset water level values greater than the thalweg elevation of the cross-section. For example, the river channel range interval where the thalweg point of the cross-section is located at the target water level value can be set as the target interval. At this time, the cross-sectional area of the river channel within the target interval needs to be greater than or equal to the product of the cross-sectional area of the main channel and the preset ratio value, so that the target water level value can better represent the actual characteristics of the river channel longitudinal profile. The preset ratio value can be, for example, 10%. In addition, if there are multiple preset water level values that meet the above conditions (the cross-sectional area of the river channel within the target interval is greater than or equal to the product of the cross-sectional area of the main channel and the preset ratio value), then take the minimum value of the multiple preset water level values that meet the above conditions as the target water level value.

[0088] In some embodiments of the present application, before determining the target water level value among multiple preset water level values greater than the thalweg elevation of the cross-section, it may further include: for each preset water level value, using the preset water level value to divide multiple observation points sorted in sequence according to the starting distance to obtain multiple second observation point intervals. The division method of the multiple second observation point intervals is similar to that of the multiple first observation point intervals, which will not be elaborated here; regarding the second observation point interval where the thalweg is located as the river channel range interval where the thalweg of the cross-section is located under the preset water level value, the determination method of the river channel range interval where the thalweg is located is similar to that of the low-flow river channel range interval, which will not be elaborated here. And the determination method of the river channel cross-sectional area within the river channel range interval where the thalweg of the cross-section is located under the preset water level value is similar to that of the main channel cross-sectional area of the cross-section at the corresponding cross-section floodplain water level, which will not be elaborated here.

[0089] In some embodiments of the present application, multiple preset water level values greater than the thalweg elevation of the cross-section can be determined through the following steps: determining the sum of the thalweg elevation and the preset elevation to obtain the target elevation. The preset elevation can be, for example, 0.5 meters; increasing the target elevation sequentially according to the preset step size to obtain multiple preset water level values.

[0090] In a further embodiment, increasing the target elevation sequentially according to the preset step size to obtain multiple preset water level values may include: since the number of nodes in the cross-section within the preset elevation range above the thalweg elevation is too small to effectively smooth the river channel longitudinal profile, first quickly search above the target elevation at the first preset step size (for example, 0.1 meter) for a new target elevation where the number of consecutive nodes near the cross-section where the thalweg is located is greater than the preset number (for example, 6); then, increasing the new target elevation sequentially at the second preset step size (for example, 0.02 meter) to obtain multiple preset water level values. Among them, the second preset step size is smaller than the first preset step size to facilitate accurately finding the target water level value that meets the area condition among the multiple preset water level values.

[0091] In some embodiments of the present application, an example is given to illustrate the refined content of steps 102 and 103. Specifically, steps 102 and 103 may include:

[0092] Step 4, identifying the elevation of the deep trough of the measured large cross-section.

[0093] 4.1. Using the while conditional statement in software such as MATLAB, the trapezoidal area summation method for adjacent observation points is adopted to calculate the cross-sectional area ans_A of the section where the thalweg point is located under different preset water levels ans_Z. Among them, the initial water level of ans_Z is 0.5 m above the thalweg point elevation, that is, ans_Z = Z + 0.5 m, and the step size is 0.02 m. It continues until the cross-sectional area of the section where the thalweg is located under ans_Z is greater than or equal to 10% of the main channel cross-sectional area A under the floodplain water level, that is, ans_A ≥ 0.1 A, and the preset water level at this time is determined as the thalweg elevation of the measured cross-section.

[0094] 4.2. Using the for loop statement in MATLAB software, calculate the thalweg elevations of different cross-sections in the same measurement.

[0095] 104. Using the thalweg elevations of multiple cross-sections in the target river reach to determine the longitudinal profile of the river channel in the target river reach.

[0096] In the embodiments of the present application, taking the target water level value as the thalweg elevation of the corresponding cross-section means taking the target water level value as the depth characteristic value of the main channel of the corresponding cross-section. Using the thalweg elevations of multiple cross-sections in the target river reach, the longitudinal profile of the river channel in the target river reach can be drawn.

[0097] In some embodiments of the present application, step 104 may include: performing a moving average process on the thalweg elevations of multiple cross-sections in the target river reach in sequence according to the cross-section distance from the dam mileage, so as to obtain the processed thalweg elevations of multiple cross-sections, making the transition between the thalweg elevations of multiple cross-sections smoother. Among them, the cross-section distance from the dam mileage refers to the distance of the cross-section relative to the target dam site. It can be understood that the distances of different cross-sections in the target river reach relative to the target dam site are different, so different cross-sections can be distinguished by the cross-section distance from the dam mileage; connecting the processed thalweg elevations of multiple cross-sections according to the cross-section distance from the dam mileage to obtain the longitudinal profile of the river channel in the target river reach, making the sawteeth in the longitudinal profile of the river channel fewer. Based on the longitudinal profile of the river channel, the trough depth conditions and trough depth change conditions of different cross-sections in the target river reach can be known, so as to efficiently, quickly and simply judge the time-series change law of the longitudinal profile of the alluvial river channel.

[0098] In some embodiments of the present application, an example is given to illustrate the refined content of step 104. Specifically, step 104 may include:

[0099] Step 5. Draw the longitudinal profile of the river channel in the same measurement.

[0100] 5.1. Based on the calculation results in step 4.2, use the moving average method to smooth the thalweg elevations of adjacent cross-sections along the way.

[0101] 5.2. Use the plot function in software such as MATLAB to plot the relationship diagram between the measured elevation of the deep trough of the large cross-section and the mileage L from the dam for the same measurement, that is, the longitudinal profile of the river channel.

[0102] In some embodiments of the present application, after step 104, it may further include:

[0103] Step 6. Determine the variation law of the longitudinal profile of the river channel.

[0104] 6.1. Use the for loop statement in software such as MATLAB to calculate the longitudinal profiles of the river channel for different measurements;

[0105] 6.2. Use the plot function in software such as MATLAB to plot the longitudinal profiles of the river channel for different measurements on the same graph to obtain the variation law of the longitudinal profile of the river channel.

[0106] In the technical solution disclosed in this embodiment, based on the product of the main channel cross-sectional area at the corresponding cross-sectional flat flood level and the preset proportional value, the target water level value with the corresponding river channel cross-sectional area is determined and used as the elevation of the deep trough of the corresponding cross-section, and then the longitudinal profile of the target river section is determined. Compared with the longitudinal profile of the thalweg drawn based on the elevation of the thalweg point of the cross-section, there are fewer sawtooth shapes in the longitudinal profile of the river channel, which is convenient for accurately identifying the longitudinal adjustment law of the river channel and reducing the level requirements for relevant technical personnel.

[0107] Next, with reference to Figure 1 , an example description of the method for determining the longitudinal profile of the river channel based on the river channel cross-sectional topography is given. Specifically, identifying the variation law of the longitudinal profile of the target river section includes the following steps:

[0108] Step 1. Standardize the long-term observation data of the measured large cross-section to form a data set including the starting distance X, the riverbed elevation Y, and the mileage L from the dam.

[0109] Collect and organize the long-term historical observation data of the measured large cross-section of the target river section, store the data of different cross-sections for the same measurement in the same Excel file in the format of the starting distance X and the elevation Y, as shown in Table 1 below. At the same time, sort out the cross-sectional data of all measurements, store the mileage L from the dam of different cross-sections for the same measurement in different sheets of the above Excel file, and store the data of different measurements named by the measurement time in different Excel files.

[0110] Table 1 Standardized data set of different cross-sections for the same measurement

[0111]

[0112] Step 2. Determine the typical water level data set of the measured large cross-section.

[0113] Collect and organize the data of the low water level and bankfull water level of typical hydrological stations or typical cross-sections in the target river reach. Assuming that the water surface slope between hydrological stations or typical cross-sections remains unchanged, use linear interpolation to obtain the characteristic water levels of the measured large cross-section, and form a data set including the cross-section name, the mileage L from the dam, the low water level Z_low, and the bankfull water level Z_bankfull, as shown in Table 2 below.

[0114] Table 2 Data set of characteristic parameters of measured large cross-section

[0115]

[0116] Step 3: Identification of the main channel range of the measured large cross-section.

[0117] Taking cross-section 66 in 2001 as an example, the identification process of the main channel range of the cross-section data is illustrated. Use the find statement in MATLAB to determine the elevation of the lowest point of the cross-section (i.e., the thalweg elevation Z = 28.1m), and the cross-section node number d_thaw where the thalweg is located is 48. Use the find statement in MATLAB to extract the area where the riverbed elevation of the cross-section is less than the low water level Z_low = 36.84m of the cross-section. The continuous node range of this area is [13, 50], including the cross-section node 48 where the thalweg is located, and then determine this area as the low water channel range.

[0118] Then, based on the characteristics of the increasing elevation of the left and right banks of the area outside the low water channel within the main channel range, use the while conditional statement in MATLAB. In the order of increasing starting distance, with the increase of the elevation of adjacent nodes on the left and right as the standard, expand the low water channel range until the elevation of the left and right banks decreases, and determine the range between the left and right nodes at this time as the main channel range. The continuous node range of the main channel area of this cross-section is [3, 53].

[0119] Step 4: Identification of the deep channel elevation of the measured large cross-section.

[0120] Adopt the trapezoidal area summation method of adjacent observation points to calculate the main channel cross-sectional area A within the node range [3, 53] below the bankfull water level Z_bankfull = 42.82m.

[0121] Considering that the elevation of the thalweg point is Z = 28.1 m, and 0.5 m above it, the number of nodes in the cross-section is too small to effectively smooth the longitudinal profile of the river channel. Therefore, first, using the while conditional statement in MATLAB, with a step size of 0.1 m, search for the preset water level (i.e., ans_Z1 = 32.3 m) where the number of consecutive nodes near cross-section node 48 where the thalweg point is located, above 0.5 m of the elevation of the thalweg point (i.e., above the elevation of ans_Z = 28.6 m), is greater than the preset number (e.g., 6) to facilitate subsequent smoothing of the longitudinal profile of the river channel. Then, using the while conditional statement in MATLAB with a step size of 0.02 m, calculate the cross-sectional area ans_A of the interval where the thalweg point is located at different preset water levels by summing the trapezoidal areas of adjacent observation points until the cross-sectional area of the interval where the thalweg point is located at the preset water level is greater than or equal to 10% of the cross-sectional area A of the main channel at the floodplain water level (i.e., ans_A ≥ 0.1 A). And determine the water level at this time as the thalweg elevation of the measured large cross-section, that is, the thalweg elevation of cross-section 66 in 2001 is 32.50 m.

[0122] Then, using the for loop statement in MATLAB, calculate the thalweg elevation and thalweg elevation of different cross-sections in the same measurement. The calculation results of the thalweg elevation and thalweg elevation of multiple cross-sections in the target reach in 2003 are as Figure 2 shown.

[0123] Step five, draw the longitudinal profile of the river channel for the same measurement.

[0124] To better display the variation law of the longitudinal profile of the river channel, based on the calculation results of the thalweg elevation of different cross-sections in the same measurement obtained in step four, use the moving average method to smooth the thalweg elevation of adjacent cross-sections along the way. Then, use the plot function in MATLAB to draw the relationship diagram between the thalweg elevation of the cross-section and the mileage L from the dam in the target reach for the same measurement, that is, the longitudinal profile of the river channel.

[0125] Step six, determine the variation law of the longitudinal profile of the river channel.

[0126] Using the for loop statement in MATLAB, adopt the above steps to calculate the longitudinal profiles of the river channel for different measurements in the target reach. Then use the plot function in MATLAB to draw the longitudinal profiles of the river channel for different measurements in the target reach on the same graph. The longitudinal profiles of the river channel for different measurements (e.g., the two measurements in 2001 and 2021) in the target reach are for example Figure 3 shown.

[0127] It can be seen that based on the measured large cross-section data of a long time series, the present invention uses the preset water level corresponding to the preset proportion value (such as 10%) of the floodplain area within the main channel of the river as the elevation of the deep trough of the cross-section, and takes the connection line of the deep trough elevations of each cross-section along the way as the longitudinal profile line of the river channel, which can effectively reduce the sawtooth-like fluctuations of the longitudinal profile of the thalweg of the river channel, and realizes the automatic identification of the deep trough elevation of the cross-section and the automatic drawing of the longitudinal profile line through software such as MATLAB. This technology is based on the measured large cross-section data, with relatively small data acquisition difficulty, low cost and high accuracy, and can efficiently, quickly and simply identify the time-series change law of the longitudinal profile of the alluvial river channel, which is of great significance for predicting the characteristics of river channel evolution, formulating appropriate river (navigation) regulation plans and flood prediction, etc.

[0128] In order to better implement the method for determining the longitudinal profile of a river channel based on the river channel cross-section topography in the embodiments of the present application, on the basis of the method for determining the longitudinal profile of a river channel based on the river channel cross-section topography, the embodiments of the present application also provide a device for determining the longitudinal profile of a river channel based on the river channel cross-section topography, as Figure 4 shown, the device 400 for determining the longitudinal profile of a river channel based on the river channel cross-section topography includes:

[0129] A first acquisition module 401, configured to acquire the main channel cross-section data of multiple cross-sections in the target river reach, where the main channel cross-section data includes the main channel cross-section area of each cross-section at the floodplain water level corresponding to the cross-section;

[0130] A first determination module 402, configured to determine a target water level value for each cross-section from multiple preset water level values greater than the elevation of the thalweg point of the cross-section, where the river channel range interval where the thalweg point of the cross-section is located at the target water level value is the target interval, and the river channel cross-section area within the target interval is greater than or equal to the product of the main channel cross-section area and the preset proportion value;

[0131] A second determination module 403, configured to use the target water level value as the elevation of the deep trough of the corresponding cross-section;

[0132] A third determination module 404, configured to determine the longitudinal profile of the target river reach by using the deep trough elevations of multiple cross-sections in the target river reach.

[0133] The embodiments of the present application also provide a computer device, which integrates any device for determining the longitudinal profile of a river channel based on the river channel cross-section topography provided in the embodiments of the present application. As Figure 5 shown, it shows the structural schematic diagram of the computer device involved in the embodiments of the present application. Specifically:

[0134] This computer device may include a processor 501 with one or more processing cores, a memory 502 with one or more computer-readable storage media, a power supply 503, an input unit 504 and other components. Those skilled in the art can understand, Figure 5The computer device structure shown in the figure does not constitute a limitation on the computer device. It may include more or fewer components than shown, or combine certain components, or have different component arrangements. Among them:

[0135] The processor 501 is the control center of the computer device. It connects various parts of the entire computer device through various interfaces and circuits. 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-mentioned modem processor may not be integrated into the processor 501 either.

[0136] 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 mainly includes 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 the sound playback function, image playback function, etc.); the data storage area can store the data created according to the use of the computer device. In addition, the memory 502 can include high-speed random access memory, and can also include non-volatile memory, such as at least one magnetic disk storage device, flash memory device, or other volatile solid-state storage devices. Correspondingly, the memory 502 can also include a memory controller to provide the processor 501 with access to the memory 502.

[0137] The computer device also includes a power supply 503 that powers each component. Preferably, the power supply 503 can be logically connected to the processor 501 through a power management system, so as to implement functions such as management of charging, discharging, and power consumption management through the power management system. The power supply 503 can also 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.

[0138] The computer device may also 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.

[0139] Although not shown, the computer device may further include a display unit and the like, 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 achieve various functions as follows:

[0140] Obtain the main channel cross-section data of multiple cross-sections in the target river section. The main channel cross-section data includes the main channel cross-section area of each cross-section at the corresponding cross-section floodplain water level; for each cross-section, determine the target water level value among multiple preset water level values greater than the thalweg elevation of the cross-section. Wherein, the river channel range interval where the thalweg is located at the target water level value of the cross-section is the target interval, and the river channel cross-section area within the target interval is greater than or equal to the product of the main channel cross-section area and the preset ratio value; use the target water level value as the deep channel elevation of the corresponding cross-section; use the deep channel elevations of multiple cross-sections in the target river section to determine the river channel longitudinal profile of the target river section.

[0141] Those of ordinary skill in the art can understand that all or part of the steps in the above-mentioned various methods can be completed by instructions or by controlling relevant hardware through instructions. The instructions can be stored in a computer-readable storage medium and loaded and executed by a processor.

[0142] For this reason, 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 one of the river channel longitudinal profile determination methods based on the river channel cross-section topography provided by the embodiments of the present application. For example, when the computer program is loaded by a processor, it can execute the following steps:

[0143] Obtain the main channel cross-section data of multiple cross-sections in the target river section. The main channel cross-section data includes the main channel cross-section area of each cross-section at the corresponding cross-section floodplain water level; for each cross-section, determine the target water level value among multiple preset water level values greater than the thalweg elevation of the cross-section. Wherein, the river channel range interval where the thalweg is located at the target water level value of the cross-section is the target interval, and the river channel cross-section area within the target interval is greater than or equal to the product of the main channel cross-section area and the preset ratio value; use the target water level value as the deep channel elevation of the corresponding cross-section; use the deep channel elevations of multiple cross-sections in the target river section to determine the river channel longitudinal profile of the target river section.

[0144] 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 longitudinal profile of a river channel based on the cross-sectional topography of the river channel as described in any one of the above, for example:

[0145] Obtain the main channel cross-section data of multiple cross-sections in the target river reach. The main channel cross-section data includes the main channel cross-sectional area of each cross-section at the corresponding cross-section floodplain water level; for each cross-section, among multiple preset water level values greater than the thalweg elevation of the cross-section, determine the target water level value. The river channel range interval where the thalweg is located at the target water level value of the cross-section is the target interval, and the cross-sectional area of the river channel within the target interval is greater than or equal to the product of the main channel cross-sectional area and the preset ratio value; use the target water level value as the deep channel elevation of the corresponding cross-section; use the deep channel elevations of multiple cross-sections in the target river reach to determine the longitudinal profile of the target river reach.

[0146] In the above embodiments, the descriptions of the respective embodiments have their own emphases. 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.

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

[0148] The specific implementation of each of the above operations can refer to the foregoing embodiments, and details will not be repeated here.

[0149] The above has introduced in detail a method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel 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 descriptions of the above embodiments are 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 longitudinal profile of a river channel based on the cross-sectional topography of the river channel, characterized in that The method includes: Obtaining the main channel cross-section data of multiple cross-sections in the target river reach, where the main channel cross-section data includes the main channel cross-sectional area of each cross-section at the corresponding cross-section floodplain water level; For each cross-section, determining a target water level value among multiple preset water level values greater than the thalweg elevation of the cross-section. Wherein, the river channel range interval where the thalweg is located at the target water level value of the cross-section is the target interval, and the river channel cross-sectional area within the target interval is greater than or equal to the product of the main channel cross-sectional area and the preset ratio value; Taking the target water level value as the thalweg elevation of the corresponding cross-section; Using the thalweg elevations of multiple cross-sections in the target river reach to determine the longitudinal profile of the river channel of the target river reach.

2. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 1, wherein The obtaining of the main channel cross-section data of multiple cross-sections in the target river reach includes: For each cross-section, obtaining the main channel range interval of the cross-section; Determining the area of the region of the cross-section within the corresponding main channel range interval and at the corresponding cross-section floodplain water level, and taking it as the main channel cross-sectional area of the cross-section at the corresponding cross-section floodplain water level.

3. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 2, wherein The obtaining of the main channel range interval of the cross-section includes: Obtaining the observation data of the cross-section, where the observation data includes the starting distances and riverbed elevations of multiple observation points in the cross-section; Based on the riverbed elevation, determining the thalweg among the multiple observation points; Using the low water level of the cross-section to divide the multiple observation points sorted in sequence according to the starting distances, obtaining multiple first observation point intervals; Taking the first observation point interval where the thalweg is located as the low-flow river channel range interval of the cross-section; Based on the low-flow river channel range interval, determining the main channel range interval of the cross-section.

4. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 3, wherein The determining of the main channel range interval of the cross-section based on the low-flow river channel range interval includes: Determining the maximum starting distance and the minimum starting distance among the starting distances of multiple observation points in the low-flow river channel range interval; Detecting whether there is a target observation point that satisfies the first preset condition or the second preset condition among the multiple observation points sorted in sequence according to the starting distances. Wherein, the first preset condition includes that the target observation point does not exist in the low-flow river channel range interval, the target observation point is adjacent to the observation point where the maximum starting distance is located among the multiple observation points sorted in sequence according to the starting distances, and there is a first preset size relationship between the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the maximum starting distance is located. The second preset condition includes that the target observation point does not exist in the low-flow river channel range interval, the target observation point is adjacent to the observation point where the minimum starting distance is located among the multiple observation points sorted in sequence according to the starting distances, and there is a second preset size relationship between the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the minimum starting distance is located; If there is the target observation point, adding the target observation point to the low-flow river channel range interval to obtain an expanded low-flow river channel range interval; Based on the expanded low-flow river channel range interval, determining the main channel range interval of the cross-section.

5. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 1, characterized in that, The observation data of the cross-section includes the starting distances and riverbed elevations of multiple observation points in the cross-section. Among the multiple observation points, there is a thalweg point, which is determined based on the riverbed elevation. Before determining the target water level value among multiple preset water level values greater than the elevation of the thalweg point of the cross-section, it further includes: For each of the preset water level values, using the preset water level value, divide the multiple observation points sorted in sequence according to the starting distance to obtain multiple second observation point intervals; Take the second observation point interval where the thalweg point is located as the river channel range interval where the thalweg point of the cross-section is located under the preset water level value.

6. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 1 or 5, characterized in that The multiple preset water level values greater than the elevation of the thalweg point of the cross-section are determined through the following steps: Determine the sum of the elevation of the thalweg point and the preset elevation to obtain the target elevation; Increase the target elevation sequentially according to the preset step length to obtain the multiple preset water level values.

7. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 1, characterized in that The method of determining the longitudinal profile of the river channel using the thalweg elevations of multiple cross-sections in the target river reach includes: Perform moving average processing on the thalweg elevations of multiple cross-sections in the target river reach in sequence according to the cross-section distance from the dam to obtain the processed thalweg elevations of multiple cross-sections; Connect the processed thalweg elevations of multiple cross-sections according to the cross-section distance from the dam to obtain the longitudinal profile of the river channel of the target river reach.

8. The method for determining the longitudinal profile of a river channel based on the cross-sectional topography of the river channel according to claim 1, wherein The flat flood water level of the cross-section is determined through the following steps: Obtain the measured historical flat flood water levels of multiple target hydrological stations in the target river reach; Perform linear interpolation processing on the measured historical flat flood water levels of multiple target hydrological stations to obtain the flat flood water level of each cross-section in the target river reach.

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 longitudinal profile of a river channel 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 longitudinal profile of a river channel based on the river channel cross-section topography according to any one of claims 1 to 8.

Images