River channel evolution scale prediction method and related equipment based on cross-sectional topography

Abstract

The present application discloses a method and related equipment for predicting the scale of river channel evolution based on cross-sectional topography. The method comprises: obtaining observation data of a target river section at multiple measurement times, as well as the operating years of the target river section dam at the corresponding measurement times; determining the flat-shoal river channel morphological parameters of each section in the target river section at the corresponding measurement times based on the observation data; determining a first calculation formula for the flat-shoal river channel morphological parameters, the first calculation formula including a balance value of the flat-shoal river channel morphological parameters at the river section scale and an attenuation coefficient of the flat-shoal river channel morphological parameters at the river section scale; determining the flat-shoal river channel morphological adjustment mode and time scale of the target river section based on the flat-shoal river channel morphological parameters at the river section scale, the balance value of the flat-shoal river channel morphological parameters at the river section scale, and the attenuation coefficient of the flat-shoal river channel morphological parameters at the river section scale, and obtaining a prediction result for the scale of the flat-shoal river channel morphological evolution. The present application can reduce the cost and difficulty of predicting the river channel morphological evolution law of the river channel downstream of a reservoir.

Description

Technical Field

[0001] The present application relates to the fields of river channel evolution simulation, river channel management and protection, and specifically to a river channel evolution scale prediction method based on cross-sectional topography and related equipment. Background Art

[0002] After the construction and operation of dams on alluvial rivers, the sediment retention function of reservoirs radically alters the flow and sediment regime downstream. Sediment concentrations in the outflowing water are significantly reduced, resulting in insufficient sediment replenishment, leading to a prolonged period of non-equilibrium scour and adjustment in the downstream channel. However, intense scour in the downstream channel can damage embankment structures, threatening flood control safety and ecological health. Therefore, predicting the channel morphology of downstream channels is an urgent issue.

[0003] In related technologies, channel morphology prediction for rivers downstream of reservoirs can be achieved through methods such as physical model testing. However, physical model testing is limited by site size, is costly, and difficult to sustain. Therefore, the current cost and difficulty of predicting channel morphology for rivers downstream of reservoirs are high. Summary of the Invention

[0004] The embodiments of the present application provide a river channel evolution scale prediction method and related equipment based on cross-sectional topography, aiming to reduce the prediction cost and difficulty of the river channel morphology evolution law of the river channel downstream of the reservoir.

[0005] In one aspect, the present application provides a method for predicting river channel evolution scale based on cross-sectional topography, the method comprising:

[0006] Obtaining observation data of a target river section at multiple measurement times and the operating years of the dam in the target river section at the corresponding measurement times, wherein the observation data includes starting point distances and riverbed elevations of multiple observation points at multiple sections in the target river section at the corresponding measurement times;

[0007] Determine the flat river channel morphological parameters of each section in the target river section at corresponding measurement times based on the observation data;

[0008] Determining a first calculation formula for flat-shoal river channel morphological parameters, wherein the first calculation formula includes a reach-scale flat-shoal river channel morphological parameter equilibrium value and a reach-scale flat-shoal river channel morphological parameter attenuation coefficient, both of which are calibrated based on the flat-shoal river channel morphological parameters of all sections in the target river reach at all measurements and the operating life of the dam;

[0009] Based on the flat river channel morphological parameters of the target river section at the river section scale, the balance value of the flat river channel morphological parameters at the river section scale, and the attenuation coefficient of the flat river channel morphological parameters at the river section scale, the flat river channel morphological adjustment mode and time scale of the target river section are determined, and the flat river channel morphological evolution scale prediction result is obtained.

[0010] In some embodiments, the river section-scale flat-shoal river channel morphological parameters include the river section-scale flat-shoal cross-sectional area, flat-shoal water depth, and flat-shoal river width, and the river section-scale flat-shoal river channel morphological parameter balance values ​​include the flat-shoal cross-sectional area balance value, the flat-shoal water depth balance value, and the flat-shoal river width balance value;

[0011] The flat river channel morphology adjustment mode of the target river section is determined by the following steps:

[0012] Determining the relative adjustment intensity of the flat-shoal river channel morphological parameters at the river reach scale based on the equilibrium value of the flat-shoal river channel morphological parameters at the river reach scale and the initial value of the flat-shoal river channel morphological parameters at the river reach scale, wherein the relative adjustment intensity includes the relative adjustment intensity of the flat-shoal cross-sectional area, the relative adjustment intensity of the flat-shoal water depth, and the relative adjustment intensity of the flat-shoal river width;

[0013] The flat beach river channel morphology adjustment mode is determined based on the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width.

[0014] In some embodiments, determining the flat beach river channel morphology adjustment mode based on the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width includes:

[0015] Determine the radar chart between the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width;

[0016] Coordinates are converted to the vertices in the radar map to respectively determine the relative adjustment intensity coordinate point of the flat beach water depth, the relative adjustment intensity coordinate point of the flat beach river width, and the vector between the center point of the radar map and the coordinate origin;

[0017] Based on the relative adjustment intensity coordinate points of the flat beach water depth, the relative adjustment intensity coordinate points of the flat beach river width, and the vector between the center point of the radar chart and the coordinate origin, the flat beach river channel morphology adjustment mode is determined, and the flat beach river channel morphology adjustment mode includes at least one of the downcutting type adjustment, the widening type adjustment, and the comprehensive type adjustment.

[0018] In some embodiments, determining the flat beach river channel morphology adjustment mode based on the relative adjustment intensity coordinate point of the flat beach water depth, the relative adjustment intensity coordinate point of the flat beach river width, and the vector between the center point of the radar chart and the coordinate origin includes:

[0019] Determine the first cross product value between the vector of the relative adjustment intensity coordinate point of the flat water depth and the coordinate origin, and the vector of the center point of the radar chart and the coordinate origin;

[0020] Determine the second cross product value between the vector of the relative adjustment intensity coordinate point of Pingtan River width and the coordinate origin, and the vector of the center point of the radar chart and the coordinate origin;

[0021] Determine a first ratio between the modulus of the vector of the relative adjustment intensity coordinate point of the flat beach water depth and the coordinate origin, and the modulus of the vector of the relative adjustment intensity coordinate point of the flat beach river width and the coordinate origin;

[0022] The flat river channel morphology adjustment mode is determined based on the positive and negative values ​​of the first cross product value and the second cross product value and the size of the first ratio.

[0023] In some embodiments, the time scale for adjusting the flat river channel morphology of the target river section is determined by the following steps:

[0024] For the target river section, determining a quantitative relationship between an attenuation coefficient of a river section-scale flat river channel morphological parameter and an adjustment equilibrium time of the river section-scale flat river channel morphological parameter;

[0025] The attenuation coefficient of the flat beach river channel morphological parameters at the river section scale is input into the quantitative relationship to determine the adjustment equilibrium time of the flat beach river channel morphological parameters at the river section scale. The flat beach river channel morphological adjustment time scale includes the adjustment equilibrium time of the flat beach river channel morphological parameters at the river section scale. The adjustment equilibrium time of the flat beach river channel morphological parameters at the river section scale includes the adjustment equilibrium time of the flat beach cross-sectional area, flat beach water depth and flat beach river width at the river section scale.

[0026] In some embodiments, determining, for the target river section, the quantitative relationship between the attenuation coefficient of the river section-scale flat river channel morphological parameter and the adjustment equilibrium time of the river section-scale flat river channel morphological parameter includes:

[0027] Determining, for each section in the target river reach, a second calculation formula for the flat-shoal river channel morphological parameters, wherein the second calculation formula includes a cross-sectional-scale flat-shoal river channel morphological parameter equilibrium value and a cross-sectional-scale flat-shoal river channel morphological parameter attenuation coefficient, both of which are calibrated based on the cross-sectional-scale flat-shoal river channel morphological parameters of the corresponding section in the target river reach and the operating life of the dam;

[0028] For each section in the target river section, determining the time at which the adjustment amount of the flat river channel morphological parameter of the corresponding section in the target river section at the cross-sectional scale reaches the target adjustment amount, and using this time as the initial proposed adjustment equilibrium time, wherein the target adjustment amount is determined based on the difference between the equilibrium value of the flat river channel morphological parameter of the corresponding section in the target river section at the cross-sectional scale and the initial value of the flat river channel morphological parameter of the corresponding section in the target river section at the cross-sectional scale;

[0029] The quantitative relationship is determined based on the attenuation coefficient of the flat river channel morphological parameters at the cross-sectional scale of all sections in the target river section and the initial adjustment equilibrium time.

[0030] In some embodiments, determining the quantitative relationship based on the attenuation coefficient of the flat river channel morphological parameters at the cross-sectional scale of all sections in the target river section and the initially proposed adjustment equilibrium time includes:

[0031] Generate a scatter plot of the relationship between the attenuation coefficient of the flat river channel morphological parameter and the initial adjustment equilibrium time at the cross-sectional scale for all sections in the target river section;

[0032] Determining upper and lower envelopes in the scatter plot;

[0033] The data of the upper and lower envelopes are averaged and then curve fitting is performed to obtain the quantitative relationship.

[0034] In some embodiments, the first calculation formula is determined by the following formula:

[0035]

[0036] in, is the value of the flat river channel morphological parameter at (x, t), x is the distance from the dam to the corresponding section in the target river section, t is the operating life of the dam, and Y e is the corresponding equilibrium value of the flat river channel morphological parameters, Φ(x, t) is the flat river channel morphological boundary function, Φ e is the boundary equilibrium value function of the flat river channel morphology, e is a natural constant, α is the attenuation coefficient of the corresponding flat river channel morphology parameter, v is the spatial propagation rate of the change of the corresponding flat river channel morphology parameter, and t0 is the initial value of the operating life of the dam.

[0037] On the other hand, an embodiment of the present application provides a device for predicting river channel evolution scale based on cross-sectional topography, comprising:

[0038] A first acquisition module is configured to acquire observation data of a target river section at multiple measurement times and the operating years of the dam in the target river section at the corresponding measurement times, wherein the observation data includes starting point distances and riverbed elevations of multiple observation points at multiple sections in the target river section at the corresponding measurement times;

[0039] A first determination module is used to determine the flat river channel morphological parameters of each section in the target river section at the corresponding measurement times based on the observation data;

[0040] a second determination module, configured to determine a first calculation formula for flat river channel morphological parameters, wherein the first calculation formula includes a balance value of the flat river channel morphological parameters at the river section scale and an attenuation coefficient of the flat river channel morphological parameters at the river section scale, both of which are calibrated based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating life of the dam;

[0041] The third determination module is used to determine the flat river channel morphology adjustment mode and time scale of the target river section based on the flat river channel morphology parameters of the target river section at the river section scale, the balance value of the flat river channel morphology parameters at the river section scale, and the attenuation coefficient of the flat river channel morphology parameters at the river section scale, so as to obtain the flat river channel morphology evolution scale prediction result.

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

[0043] one or more processors;

[0044] Memory; and

[0045] One or more applications, wherein the one or more applications are stored in the memory and configured to be executed by the processor to implement the steps in any one of the methods for predicting river channel evolution scale based on cross-sectional topography.

[0046] On the other hand, the present application also provides a computer-readable storage medium having a computer program stored thereon, which is loaded by a processor to execute the steps in any one of the methods for predicting river channel evolution scale based on cross-sectional topography.

[0047] On the other hand, an embodiment of the present application provides a computer program product, including a computer program or instructions, which are executed by a processor to implement the river channel evolution scale prediction method based on cross-sectional topography as described in any of the above items.

[0048] The embodiment of the present application provides a river channel evolution scale prediction method and related equipment based on cross-sectional topography, which obtains observation data of a target river section at multiple measurements and the operating years of the target river section dam at corresponding measurements; determines the flat river channel morphological parameters of each section in the target river section at corresponding measurements based on the observation data; determines a first calculation formula for the flat river channel morphological parameters, wherein the first calculation formula includes a balance value of the flat river channel morphological parameters at the river section scale and an attenuation coefficient of the flat river channel morphological parameters at the river section scale, both of which are calibrated based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating years of the dam; determines the flat river channel morphological adjustment mode and time scale of the target river section based on the flat river channel morphological parameters at the river section scale, and obtains a flat river channel morphological evolution scale prediction result. The embodiment of the present application can accurately predict the evolution law of the river channel morphology of the river channel downstream of the reservoir through the observation data of the target river section at multiple measurements and the years of operation of the dam. Compared with predicting the evolution law of the river channel morphology of the river channel downstream of the reservoir based on methods such as physical model experiments, the prediction cost and difficulty are lower. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0050] Figure 1 Schematic diagram of an embodiment of a method for predicting river channel evolution scale based on cross-sectional topography provided in an embodiment of the present application;

[0051] Figure 2 Schematic diagram of the evolution of the morphological parameters of each flat river channel at the river section a at multiple measurement times provided in the embodiment of the present application;

[0052] Figure 3 Schematic diagram of the evolution of the morphological parameters of each flat river channel at the river section scale of multiple measurements in the river section b provided in the embodiment of the present application;

[0053] Figure 4 : is a radar chart of the relative adjustment intensity of each flat river channel morphological parameter of the river section a in the embodiment of the present application (unit: %);

[0054] Figure 5 : is a radar chart of the relative adjustment intensity of each flat river channel morphological parameter of the river section b in the embodiment of the present application (unit: %);

[0055] Figure 61 is a schematic structural diagram of an embodiment of a device for predicting river channel evolution scale based on cross-sectional topography provided in an embodiment of the present application;

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

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

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

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

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

[0061] The embodiment of the present application provides a river channel evolution scale prediction method and related equipment based on cross-sectional topography. Based on a long sequence of measured large-section data, the method realizes the automatic calculation of the morphological adjustment mode and time scale of the flat-flat river channel of non-equilibrium alluvial rivers through the standardized processing of the measured large-section observation data, the construction of the target river section flat-flat river channel morphological parameter data set, and the calculation of the flat-flat river channel morphological evolution characteristic values. The results of the embodiment of the present application can be widely used for the rapid and accurate identification of river channel evolution trends in multiple fields such as riverbed evolution prediction, river channel management and protection. The following are detailed descriptions.

[0062] In one embodiment, referring to Figure 1 , the river channel evolution scale prediction methods based on cross-section topography include:

[0063] 101. Obtain observation data of the target river section at multiple survey times, as well as the operating years of the dam in the target river section at the corresponding survey times. The observation data include the starting distance and riverbed elevation of multiple observation points at multiple sections in the target river section at the corresponding survey times.

[0064] In the embodiments of the present application, the target river section is usually located downstream of the dam, that is, under the influence of the dam's reservoir sediment retention operation, the water and sediment situation entering the downstream target river section will be changed. Therefore, it is necessary to predict the evolution law of the flat river channel morphology of the target river section.

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

[0066] Step 1: Standardization of long-sequence large-section observation data.

[0067] 1.1. Sort out the long-sequence large-section observation data of the target river section. For the data of a certain measurement, set the number of sections to mm, the number of cross-section observation points to nn, and the number of measurements to yy, then the terrain data can be expressed as:

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

[0069] Where X is the distance from the observation point to the starting point of the section, also known as the starting distance; Y is the measured riverbed elevation at the observation point; T is the dam's operating age corresponding to each measurement; i is the measurement section number; j is the maximum number of observation points in each section; and k is the number of measurements. Here, i = 1, 2, …, mm; j = 1, 2, …, nn; and k = 1, 2, …, yy.

[0070] 1.2. Standardization of long-sequence field observation data. Topographic data of different field observations of large sections in the same survey are stored in the same Excel file in the format of starting point distance X and elevation Y. Data of different surveys are stored in different Excel files named by survey sequence number. The section observation time T is stored in an Excel file named by time.

[0071] 102. Based on the observation data, determine the flat river channel morphological parameters of each section in the target river section at the corresponding measurement times.

[0072] In the embodiments of the present application, the flat-beach channel morphological parameters of each section in the target river section refer to parameters used to characterize the flat-beach channel morphology of the corresponding section in the target river section. The flat-beach channel morphological parameters Y of each section in the target river section may, for example, include the flat-beach cross-sectional area A, flat-beach water depth H, and flat-beach river width B of the corresponding section in the target river section. Each section in the target river section corresponds to a set of flat-beach channel morphological parameters at each measurement.

[0073] In some embodiments of the present application, determining the flat-shoal channel morphological parameters of each section in the target river section at a corresponding measurement based on observational data may include: obtaining the flat-shoal water level at each section in the target river section at the corresponding measurement; and determining the flat-shoal channel morphological parameters of each section in the target river section at the corresponding measurement based on the flat-shoal water level and the observational data. Taking the flat-shoal river width of a section as an example, the geometric shape of the section can be determined for each section using the starting distance and riverbed elevation of multiple observation points in the section during the measurement, and the width of the geometric shape of the section at the flat-shoal water level is used as the flat-shoal river width. Taking the flat-shoal cross-sectional area of ​​a section as an example, the area of ​​the portion of the geometric shape of the section below the flat-shoal water level can be used as the flat-shoal cross-sectional area (this flat-shoal cross-sectional area can be calculated using the starting distance and riverbed elevation of multiple observation points in the section during the measurement using the trapezoidal area summation method of adjacent observation points). Taking the flat water depth of a section as an example, the ratio of the flat cross-sectional area of ​​the section to the flat river width of the section can be used as the flat water depth of the section.

[0074] In some embodiments of the present application, after step 102, it may also include: for each measurement, according to the flat river channel morphological parameters of multiple sections in the target river section in the corresponding measurement (for example, it may include the flat river section area, flat river depth and flat river width of each section in the target river section in the measurement), a spatial weighted average algorithm is used to calculate to obtain the flat river channel morphological parameters of the target river section at the river section scale at the corresponding measurement. For example, the change in mileage from the dam of adjacent sections in the target river section (that is, the distance between adjacent sections) can be used as a weight to perform weighted averaging on the flat river channel morphological parameters of multiple sections at the measurement to obtain the flat river channel morphological parameters of the target river section at the river section scale at the measurement.

[0075] The reach-scale flat-shoal channel morphological parameters of the target river section refer to parameters used to characterize the reach-scale flat-shoal channel morphology of the target river section. For example, the reach-scale flat-shoal channel morphological parameters Y of the target river section may include the reach-scale flat-shoal cross-sectional area A, flat-shoal water depth H, and flat-shoal river width B. Each measurement of the target river section corresponds to a set of reach-scale flat-shoal channel morphological parameters.

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

[0077] Step 2: Construct a dataset of flat river channel morphological parameters for all sections in the target river section.

[0078] 2.1. Identify the flat water level Z and flat river width B at each measured large section of the target river section at each measurement time. Use the for loop statement in MATLAB and the trapezoidal area summation method of adjacent observation points to calculate the main channel cross-sectional area under the flat water level, that is, the flat cross-sectional area A, and calculate the flat water depth H = A / B.

[0079] After step 102, the following steps may also be included:

[0080] 2.2. For the flat river channel morphological parameters of multiple sections in the target river section at the corresponding measurement times, the spatial weighted average method is used to calculate the flat river channel morphological parameter data at the river section scale at the corresponding measurement times.

[0081] 2.3. The flat cross-sectional area A, flat water depth H, and flat river width B of the target river section at the first measurement of the river section are recorded as the initial values ​​A0, H0, and B0 of the flat river channel morphological parameters at the river section scale.

[0082] 2.4. Construct a river channel morphological parameter dataset at the river section scale, including the flat cross-sectional area A, flat water depth H, flat river section B, and dam operation years T.

[0083] 103. Determine the first calculation formula for the flat river channel morphological parameters, wherein the first calculation formula includes the equilibrium value of the flat river channel morphological parameters at the river section scale and the attenuation coefficient of the flat river channel morphological parameters at the river section scale, both of which are calibrated based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating life of the dam.

[0084] In the embodiments of the present application, the first calculation formula for the flat-flat channel morphological parameters of a target river section characterizes the correlation between the flat-flat channel morphological parameters of the target river section and the operating age of the dam, that is, the relationship between the flat-flat channel morphological parameters of the target river section and the operating age of the dam. Because the flat-flat channel morphological parameters of the target river section tend to stabilize and reach an equilibrium state after a certain period of change, the first calculation formula includes a section-scale equilibrium value of the flat-flat channel morphological parameters and a section-scale attenuation coefficient of the flat-flat channel morphological parameters. The section-scale equilibrium value of the flat-flat channel morphological parameters refers to the flat-flat channel morphological parameters of the target river section when the flat-flat channel reaches equilibrium. The section-scale attenuation coefficient of the flat-flat channel morphological parameters characterizes the attenuation of the flat-flat channel morphological parameters of the target river section. Both the section-scale equilibrium value of the flat-flat channel morphological parameters and the section-scale attenuation coefficient of the flat-flat channel morphological parameters can be calibrated based on the flat-flat channel morphological parameters of all sections in the target river section at all measurements and the operating age of the dam.

[0085] In some embodiments of the present application, the first calculation formula for the flat river channel morphological parameters of the target river section is determined, for example, by the following formula:

[0086]

[0087] in, is the value of the flat river channel morphological parameters at (x, t), x is the distance from the dam to the corresponding section in the target river section, t is the operating life of the dam, and Y e is the corresponding equilibrium value of the flat river channel morphological parameters, Φ(x, t) is the flat river channel morphological boundary function, Φ e is the equilibrium value of the flat river channel morphological boundary function, e is a natural constant, α is the corresponding flat river channel morphological parameter attenuation coefficient, v is the spatial propagation rate of the change of the corresponding flat river channel morphological parameter, and t0 is the initial value of the dam's operating life. e The calibration can be performed based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating years of the dam at multiple measurements. The calibration method can be, for example, the least squares method.

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

[0089] Step 3: Calculate the scale characteristic value of the flat river channel morphological evolution of the target river section.

[0090] 3.1, construct a calculation method for the evolution scale of the horizontal, vertical and cross-sectional dimensions of a flat river channel.

[0091] The control equation for scale prediction of flat river channel morphological evolution is as follows:

[0092]

[0093] Where, is the value of the flat river channel morphological parameters (including the flat cross-sectional area A, flat water depth H, and flat river width B of the corresponding section in the target river section) at (x, t); x is the distance from the corresponding section in the target river section to the dam, and t is the operating life of the dam; v is the spatial propagation rate of the change of the flat river channel morphological parameters; α is the corresponding flat river channel morphological parameter attenuation coefficient; Y e is the corresponding flat river channel morphological parameter balance value (including the flat cross-sectional area balance value A of the target river section at the river section scale) e , Flat water depth balance value H e , Flat river width balance value B e ).

[0094] The above control equation is non-dimensionalized and recorded as:

[0095] , , ,

[0096] In the formula, the variable T e 、L e 、 are the equilibrium time, equilibrium distance, and characteristic propagation speed of the flat river channel morphology evolution; τ and ξ are dimensionless parameters.

[0097] Furthermore, we can get:

[0098]

[0099] The solution is:

[0100]

[0101] The first calculation formula for predicting the scale of flat river channel morphological evolution by variable reduction is as follows:

[0102]

[0103] Where Φ(x, t) is the flat river channel morphological boundary function. Considering x as a constant, this formula can predict the evolution process of the flat river channel morphological parameters Y in various dimensions (including the flat cross-sectional area A, flat water depth H, and flat river width B of the target river section at the river section scale), denoted by A e 、H e 、B eare the balanced values ​​of the flat cross-sectional area, flat water depth and flat river width of the target river section at the river section scale; A0, H0 and B0 are the initial values ​​of the flat cross-sectional area, flat water depth and flat river width of the target river section at the river section scale; α1, α2 and α3 are the attenuation coefficients of the flat cross-sectional area, flat water depth and flat river width of the target river section at the river section scale. It can be seen that the above formula takes into account the evolution of the flat river channel morphological parameters of the target river section in the time dimension and the spatial dimension, thereby making the prediction results more accurate. 3.2 Using the for loop statement in MATLAB, the least squares method is used to determine the balanced value A of the flat cross-sectional area, flat water depth and flat river width of the target river section at the river section scale. e 、H e 、B e and the values ​​of attenuation coefficients α1, α2, and α3 (wherein, for the attenuation coefficient α, the attenuation coefficients of the corresponding flat river channel morphological parameters at the cross-sectional scale can also be determined in turn, see the following content for details).

[0104] 104. Based on the flat river channel morphological parameters at the river section scale, the equilibrium value of the flat river channel morphological parameters at the river section scale, and the attenuation coefficient of the flat river channel morphological parameters at the river section scale, the flat river channel morphological adjustment mode and time scale of the target river section are determined to obtain the flat river channel morphological evolution scale prediction results.

[0105] In an embodiment of the present application, the flat-shoal river channel morphology evolution scale prediction results include a flat-shoal river channel morphology adjustment mode for the target river section and a flat-shoal river channel morphology adjustment time scale for the target river section. The flat-shoal river channel morphology adjustment mode may, for example, include at least one of a downcutting adjustment, a widening adjustment, and a combined adjustment. The flat-shoal river channel morphology adjustment time scale may, for example, include an adjustment equilibrium time for flat-shoal river channel morphology parameters.

[0106] In some embodiments of the present application, the flat river channel morphology adjustment mode of the target river section can be determined by the following steps: determining the relative adjustment intensity of the flat river channel morphology parameters at the river section scale based on the balance value of the flat river channel morphology parameters at the river section scale and the initial value of the flat river channel morphology parameters at the river section scale. For example, the difference between the balance value of the flat river channel morphology parameters at the river section scale and the initial value of the flat river channel morphology parameters at the river section scale can be determined, and then the ratio between the difference and the initial value of the flat river channel morphology parameters at the river section scale is determined and used as the relative adjustment intensity of the flat river channel morphology parameters; based on the relative adjustment intensity, the flat river channel morphology adjustment mode of the target river section is determined.

[0107] In some embodiments of the present application, the relative adjustment intensity includes the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width as an example for description. Based on the relative adjustment intensity, determining the flat beach river channel morphology adjustment mode of the target river section may include: determining the flat beach river channel morphology adjustment mode based on the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width. It can be seen that the embodiments of the present application respectively characterize the evolution of the flat beach river channel of the target river section in the cross-sectional dimension, longitudinal direction, and transverse direction by the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width, thereby more accurately predicting the flat beach river channel morphology adjustment mode of the target river section.

[0108] In some embodiments of the present application, determining the flat beach channel morphology adjustment mode based on the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width may include: determining a radar chart (e.g., Figure 4 or Figure 5 The vertices in the radar chart are coordinate-converted to determine the relative adjustment intensity coordinate points for the flat-shoal water depth, the relative adjustment intensity coordinate points for the flat-shoal river width, and the vector between the center point of the radar chart and the coordinate origin. The flat-shoal river channel morphological adjustment pattern is determined based on the relative adjustment intensity coordinate points for the flat-shoal water depth, the relative adjustment intensity coordinate points for the flat-shoal river width, and the vector between the center point of the radar chart and the coordinate origin. For example, the flat-shoal river channel morphological adjustment pattern can be determined based on the ratio between the relative adjustment intensity of the flat-shoal water depth and the relative adjustment intensity of the flat-shoal river width. The relative adjustment intensity of the flat-shoal cross-sectional area can be displayed in the radar chart to assist in determining the flat-shoal river channel morphological adjustment pattern, more intuitively helping to determine whether the change in the flat-shoal cross-sectional area is primarily contributed by the flat-shoal river width or the flat-shoal water depth, and thus determining the flat-shoal river channel morphological adjustment pattern.

[0109] In some embodiments of the present application, determining the flat beach river channel morphology adjustment mode based on the relative adjustment intensity coordinate point of the flat beach water depth, the relative adjustment intensity coordinate point of the flat beach river width, and the vector between the center point of the radar chart and the coordinate origin can include: determining a first cross product value between the relative adjustment intensity coordinate point of the flat beach water depth and the vector of the coordinate origin, and the vector between the center point of the radar chart and the coordinate origin; determining a second cross product value between the relative adjustment intensity coordinate point of the flat beach river width and the vector of the coordinate origin, and the vector between the center point of the radar chart and the coordinate origin; determining a first ratio between the modulus of the relative adjustment intensity coordinate point of the flat beach water depth and the vector of the coordinate origin, and the modulus of the relative adjustment intensity coordinate point of the flat beach river width and the vector of the coordinate origin; determining the flat beach river channel morphology adjustment mode based on the positive and negative signs of the first cross product value and the second cross product value, and the size of the first ratio.

[0110] In some embodiments of the present application, the time scale for adjusting the flat-shoal river channel morphology of a target river section can be determined by the following steps: determining, for the target river section, a quantitative relationship between a river section-scale flat-shoal river channel morphology parameter attenuation coefficient and a river section-scale adjustment equilibrium time for the flat-shoal river channel morphology parameter; and inputting the river section-scale flat-shoal river channel morphology parameter attenuation coefficient into the quantitative relationship to determine the river section-scale adjustment equilibrium time for the flat-shoal river channel morphology parameter. The flat-shoal river channel morphology adjustment time scale includes the river section-scale adjustment equilibrium time for the flat-shoal river channel morphology parameter, and the river section-scale adjustment equilibrium time for the flat-shoal river channel morphology parameter includes the river section-scale adjustment equilibrium time for the flat-shoal cross-sectional area, flat-shoal water depth, and flat-shoal river width.

[0111] In some embodiments of the present application, for a target river section, determining the quantitative relationship between the attenuation coefficient of the flat river channel morphological parameters at the river section scale and the adjustment equilibrium time of the flat river channel morphological parameters at the river section scale may include: determining a second calculation formula for the flat river channel morphological parameters for each section in the target river section, wherein the second calculation formula includes the equilibrium value of the flat river channel morphological parameters at the section scale and the attenuation coefficient of the flat river channel morphological parameters at the section scale, both of which are calibrated based on the flat river channel morphological parameters at the section scale of the corresponding section in the target river section and the operating years of the dam. The second calculation formula may refer to the first calculation formula and will not be elaborated here; for each section in the target river section, determining the time corresponding to the target adjustment amount of the flat river channel morphological parameters at the section scale of the corresponding section in the target river section, and using this as the initial proposed adjustment equilibrium time, The target adjustment amount is determined based on the difference between the equilibrium value of the flat-flat channel morphological parameters at the cross-sectional scale of the corresponding section in the target river reach and the initial value of the flat-flat channel morphological parameters at the cross-sectional scale of the corresponding section in the target river reach. For example, the proposed initial adjustment equilibrium time can be the time required for the adjustment amount of the flat-flat channel morphological parameters at the cross-sectional scale of the corresponding section in the target river reach to reach 90% of the difference between the equilibrium value of the flat-flat channel morphological parameters at the cross-sectional scale of the corresponding section in the target river reach and the initial value of the flat-flat channel morphological parameters at the cross-sectional scale of the corresponding section in the target river reach. A quantitative relationship is determined based on the attenuation coefficient of the flat-flat channel morphological parameters at the cross-sectional scale of all sections in the target river reach and the proposed initial adjustment equilibrium time. For example, the above quantitative relationship can be obtained by curve fitting based on the attenuation coefficient of the flat-flat channel morphological parameters at the cross-sectional scale of all sections in the target river reach and the proposed initial adjustment equilibrium time. The flat-flat channel morphological parameters at the cross-sectional scale of a single section in the target river reach refer to the flat-flat channel morphological parameters of that section in the target river reach at all measurements.

[0112] In some embodiments of the present application, determining a quantitative relationship based on the attenuation coefficient of the flat river channel morphological parameters of all sections in the target river section at the cross-sectional scale and the initial planned adjustment balance time can include: generating a scatter plot of the attenuation coefficient of the flat river channel morphological parameters of all sections in the target river section at the cross-sectional scale and the initial planned adjustment balance time; determining the upper and lower envelope lines in the scatter plot; taking the average of the data of the upper and lower envelope lines, and then performing curve fitting to obtain a more accurate quantitative relationship mentioned above.

[0113] In some embodiments of the present application, the attenuation coefficient of the flat beach river channel morphological parameters may, for example, include the attenuation coefficient of the flat beach cross-sectional area, the attenuation coefficient of the flat beach water depth, and the attenuation coefficient of the flat beach river width. Accordingly, the adjustment balance time of the flat beach river channel morphological parameters may include the adjustment balance time of the flat beach cross-sectional area, the adjustment balance time of the flat beach water depth, and the adjustment balance time of the flat beach river width.

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

[0115] Step 4: Analyze the morphological adjustment pattern of the flat river channel in the target river section.

[0116] 4.1, define “relative adjustment intensity D s " is the ratio of the difference between the equilibrium value and the initial value of the flat river channel morphological parameters to the initial value, which is used to characterize the relative magnitude of the adjustment amplitude in different dimensions. The for loop statement in MATLAB is used to calculate the relative adjustment intensity A of the flat cross-sectional area A, flat water depth H, and flat river width B of the target river section at the river section scale. s 、H s 、B s .

[0117] 4.2. Use the plot function in MATLAB to plot the relative adjustment intensity A of the flat cross-sectional area, flat water depth, and flat river width of the target river section at the river section scale. s 、H s 、B s Radar chart.

[0118] 4.3. Coordinate the vertices of the radar chart and calculate the coordinates of its center point (denoted as point D) using the MATLAB program. Construct the relative adjustment intensity coordinate points of the flat beach water depth, the relative adjustment intensity coordinate points of the flat beach river width, and the vector expressions of the center point of the radar chart and the coordinate origin (denoted as point O), respectively, and denote them as the relative adjustment intensity vectors of the flat beach water depth. , Pingtan River Width Relative Adjustment Intensity Vector and the center vector , calculate the relative adjustment intensity vector of the flat river width , Flat water depth relative adjustment intensity vector and the center vector The cross product value of the width of the flat river and the relative adjustment intensity vector modulus , Flat water depth relative adjustment intensity vector model The value of .

[0119] 4.4, the flat river channel morphology adjustment mode of the target river section is determined based on the positive and negative cross product values ​​of the flat river depth relative adjustment intensity vector and the central vector, the positive and negative cross product values ​​of the flat river width relative adjustment intensity vector and the central vector, and the relative magnitude relationship between the flat river depth relative adjustment intensity vector modulus and the flat river width relative adjustment intensity vector modulus. For example, if and The cross product of is positive, and , it is determined to be a down-cut adjustment; if and The cross product of is positive and , it is determined to be a widening type adjustment; if it cannot be determined to be a down-cut type adjustment or a widening type adjustment, it is determined to be a comprehensive type adjustment.

[0120] Step 5: Predict the time scale of the flat river channel morphology adjustment in the target river section.

[0121] 5.1. Using the second calculation formula for flat river channel morphological parameters, the equilibrium value and attenuation coefficient α' of the flat river channel morphological parameters in each dimension of each observation section of the target river section are calculated. Then, the time corresponding to the adjustment amount of the flat river channel morphological parameters in each dimension of each observation section of the target river section reaches 90% of the difference between its equilibrium value and the initial value is determined, which is defined as the initial adjustment equilibrium time T'.

[0122] 5.2. Use the plot function of MATLAB to draw the scatter plot of α'~T' for all observation sections in the target river section. After removing the abnormal points, determine the upper and lower envelopes of the corresponding scatter plot.

[0123] 5.3, the upper and lower envelope grid data are exported and averaged, and then the power curve is used to fit the target river section to obtain the quantitative relationship between α and T at the river section scale (i.e. T = aα b Taking into account the overall evolution characteristics of the target river section, the statistical laws of the measured sequence and the ease of operation, a can be approximately taken as 2.5 and b can be approximately taken as -1, that is, T=2.5α -1 , then the adjustment balance time of the flat cross-sectional area A of the target river section at the river section scale is T1=2.5 / α1, the adjustment balance time of the flat water depth H at the river section scale is T2=2.5 / α2, and the adjustment balance time of the flat river width B at the river section scale is T3=2.5 / α3.

[0124] 5.4. Use the plot function in MATLAB to plot the comparative relationship between the flat cross-sectional area A, flat water depth H, flat river width B, and equilibrium time T1, T2, and T3 of the target river section at the river section scale, and analyze the time scale characteristics of the cross-sectional morphological evolution of the target river section.

[0125] In the technical solution disclosed in this embodiment, by obtaining the observation data of the target river section at multiple measurements, the flat river channel morphological parameters of each section in the target river section at the corresponding measurements are determined based on the observation data, and then based on the flat river channel morphological parameters and the operating years of the dam, the river section scale flat river channel morphological parameter balance value and the river section scale flat river channel morphological parameter attenuation coefficient in the first calculation formula of the flat river channel morphological parameters of the target river section are determined. According to the river section scale flat river channel morphological parameters, the river section scale flat river channel morphological parameter balance value and the river section scale flat river channel morphological parameter attenuation coefficient, the flat river channel morphological adjustment mode and time scale of the target river section are determined, and the flat river channel morphological evolution scale prediction result is obtained. The embodiment of the present application can accurately predict the river channel morphological evolution law of the river channel downstream of the reservoir by using the observation data of the target river section at multiple measurements and the operating years of the dam. Compared with predicting the river channel morphological evolution law of the river channel downstream of the reservoir based on methods such as physical model experiments, the prediction cost and difficulty are lower.

[0126] Below, refer to Figure 1 This paper illustrates a method for predicting river channel evolution scale based on cross-sectional topography. Specifically, two typical river sections in the middle reaches of the Yangtze River (referred to as Section A and Section B) were selected as target river sections. The morphological adjustment patterns of the flat river channel in these target river sections were analyzed. The following steps were included:

[0127] Step 1: Standardize the long-sequence large-section observation data to form a data set including starting point distance X, riverbed elevation Y and observation time T.

[0128] Organize the long sequence of large-section observation data for the target river reach. Store the topographic data for different large-sections within the same survey in the same Excel file using the format of starting point distance (X) and elevation (Y), as shown in Table 1 below. Also, store the large-section topographic data for different surveys in separate Excel files named after the survey sequence, and store the observation time (T) in an Excel file named after the time.

[0129] Table 1 Standardized datasets of different measured large-section topography for a typical river section in the same survey

[0130]

[0131] Step 2: Construct a dataset of morphological parameters of the flat river channel in the target river section.

[0132] The measured flat water level Z and flat river width B of the target river section were identified. Using a for loop in MATLAB, the trapezoidal area summation method for adjacent observation points was used to sequentially calculate the main channel cross-sectional area at the flat water level Z, namely the flat cross-sectional area A. The flat water depth H was also calculated as A / B. Using spatial weighted averaging, the reach-scale flat channel morphological data were calculated (the flat cross-sectional area, flat water depth, and flat river width at the first measurement were taken as A0, H0, and B0, respectively). This generated a flat channel morphological parameter dataset consisting of the flat cross-sectional area A, flat water depth H, flat river section B, and dam operating life T, as shown in Table 2 below.

[0133] Table 2 Dataset of morphological parameters of flat river channel in typical river section a

[0134]

[0135] Step three: calculate the scale characteristic values ​​of the flat river channel morphology evolution in the target river section.

[0136] Taking two typical river sections (Section A and Section B) as examples, the calculation process of the scale characteristic value of the flat river channel morphological evolution is explained. Using the for statement in MATLAB, the least squares method is used to determine the flat cross-sectional area A, flat water depth H, and flat river width B of section a at the river section scale. e 、H e 、B e The values ​​of the attenuation coefficients α1, α2, and α3 of the time adjustment process of the flat cross-sectional area A, flat water depth H, and flat river width B at the river section scale are 0.06, 0.05, and 0.40, respectively. The equilibrium value A of the flat cross-sectional area A, flat water depth H, and flat river width B at the river section scale of river section b is e 、H e 、B e The values ​​of the attenuation coefficients α1, α2, and α3 for the time adjustment process of the flat cross-sectional area, flat water depth, and flat river width at the river section scale are 0.25, 0.20, and 0.60, respectively. Figure 2 and Figure 3 shown.

[0137] Step 4: Analyze the morphological adjustment pattern of the flat river channel in the target river section.

[0138] Calculate the relative adjustment intensity A of the flat cross-sectional area A, flat water depth H, and flat river width B of typical river section a and typical river section b at the river section scale s 、H s 、B sTake the values ​​and use the plot function in MATLAB to draw the relative adjustment intensity A of the flat cross-sectional area A, flat water depth H, and flat river width B of the typical river section a and the typical river section b at the river section scale. s 、H s 、B s The radar chart is as follows Figure 4 and Figure 5 shown.

[0139] The morphological adjustment pattern of the flat-shoal channel in the target river section was determined based on the positive and negative cross products of the relative adjustment intensity vectors of the flat-shoal depth and the central vector, the positive and negative cross products of the relative adjustment intensity vectors of the flat-shoal width and the central vector, and the relative magnitude relationships between the relative adjustment intensity vectors of the flat-shoal depth and the relative adjustment intensity vectors of the flat-shoal width. The results showed that the adjustment patterns of the flat-shoal channel in typical river sections a and b were, respectively, downcutting and comprehensive.

[0140] Step 5: Predict the time scale for adjusting the flat river channel morphology of the target river section.

[0141] In order to better show the temporal variation of the flat river channel morphology adjustment of typical river sections a and b, based on the calculation results of the attenuation coefficients α1, α2, and α3 of the flat cross-sectional area, flat water depth, and flat river width obtained in step 3, the quantitative relationship α~T of the target river section (i.e., T=aα b ), the adjustment balance time of the flat cross-sectional area A of the target river section at the river section scale is T1=2.5 / α1, the adjustment balance time of the flat water depth H is T2=2.5 / α2, and the adjustment balance time of the flat river width B is T3=2.5 / α3. The adjustment balance time of the flat cross-sectional area A, flat water depth H, and flat river width B of the target river section at the river section scale are calculated. The values ​​of T1, T2, and T3 of the typical river section a are 42, 50, and 6 respectively, and the values ​​of T1, T2, and T3 of the typical river section b are 10, 13, and 4 respectively.

[0142] By comparing the relationship between the adjustment equilibrium times T1, T2, and T3 of the flat cross-sectional area A, flat water depth H, and flat river width B of typical river sections a and b, the time scale characteristics of the cross-sectional morphology evolution of the target river section can be clarified.

[0143] It can be seen that the river channel evolution scale prediction method based on cross-sectional topography provided by the embodiment of the present application is based on long-sequence measured large-section data, and the data acquisition is relatively easy, low-cost, and highly accurate. Through the standardized processing of measured large-section observation data, the construction of a data set of flat-flat river channel morphological parameters of the target river section, and the calculation of the characteristic values ​​of the flat-flat river channel morphological evolution, the automatic calculation of the flat-flat river channel morphological adjustment mode and time scale of the non-equilibrium alluvial river can be realized, which has the advantages of accurate calculation, simple method, and fast operation. The embodiment of the present application is of great significance for quickly predicting the evolution trend of river channels, formulating corresponding management and protection plans, etc.

[0144] In order to better implement the river channel evolution scale prediction method based on cross-sectional topography in the embodiment of the present application, on the basis of the river channel evolution scale prediction method based on cross-sectional topography, the embodiment of the present application also provides a river channel evolution scale prediction device based on cross-sectional topography, such as Figure 6 As shown, the river channel evolution scale prediction device 600 based on cross-sectional topography includes:

[0145] The first acquisition module 601 is used to obtain observation data of a target river section at multiple measurement times and the operating years of the dam in the target river section at the corresponding measurement times, wherein the observation data includes the starting point distance and riverbed elevation of multiple observation points at multiple sections in the target river section at the corresponding measurement times;

[0146] The first determination module 602 is used to determine the flat river channel morphological parameters of each section in the target river section at the corresponding measurement time based on the observation data;

[0147] A second determination module 603 is configured to determine a first calculation formula for flat-shoal river channel morphological parameters, wherein the first calculation formula includes a reach-scale flat-shoal river channel morphological parameter equilibrium value and a reach-scale flat-shoal river channel morphological parameter attenuation coefficient, both of which are calibrated based on the flat-shoal river channel morphological parameters of all sections in the target river reach at all measurements and the operating life of the dam;

[0148] The third determination module 604 is used to determine the flat river channel morphology adjustment mode and time scale of the target river section based on the flat river channel morphology parameters of the target river section at the river section scale, the balance value of the flat river channel morphology parameters at the river section scale, and the attenuation coefficient of the flat river channel morphology parameters at the river section scale, so as to obtain the flat river channel morphology evolution scale prediction result.

[0149] The embodiment of the present application also provides a computer device that integrates any of the river channel evolution scale prediction devices based on cross-sectional topography provided in the embodiment of the present application. Figure 7 , which shows a schematic diagram of the structure of the computer device involved in the embodiment of the present application, specifically:

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

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

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

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

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

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

[0156] Obtain observation data of the target river section at multiple measurements and the operating years of the dam in the target river section at the corresponding measurements, the observation data including the starting point distances and riverbed elevations of multiple observation points in multiple sections in the target river section at the corresponding measurements; determine the flat river channel morphological parameters of each section in the target river section at the corresponding measurements based on the observation data; determine a first calculation formula for the flat river channel morphological parameters, wherein the first calculation formula includes a balance value of the flat river channel morphological parameters at the river section scale and an attenuation coefficient of the flat river channel morphological parameters at the river section scale, both of which are calibrated based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating years of the dam; determine the flat river channel morphological adjustment mode and time scale of the target river section based on the flat river channel morphological parameters at the river section scale, the balance value of the flat river channel morphological parameters at the river section scale and the attenuation coefficient of the flat river channel morphological parameters at the river section scale, and obtain a scale prediction result of the flat river channel morphological evolution.

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

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

[0159] Obtain observation data of the target river section at multiple measurements and the operating years of the dam in the target river section at the corresponding measurements, the observation data including the starting point distances and riverbed elevations of multiple observation points in multiple sections in the target river section at the corresponding measurements; determine the flat river channel morphological parameters of each section in the target river section at the corresponding measurements based on the observation data; determine a first calculation formula for the flat river channel morphological parameters, wherein the first calculation formula includes a balance value of the flat river channel morphological parameters at the river section scale and an attenuation coefficient of the flat river channel morphological parameters at the river section scale, both of which are calibrated based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating years of the dam; determine the flat river channel morphological adjustment mode and time scale of the target river section based on the flat river channel morphological parameters at the river section scale, the balance value of the flat river channel morphological parameters at the river section scale and the attenuation coefficient of the flat river channel morphological parameters at the river section scale, and obtain a scale prediction result of the flat river channel morphological evolution.

[0160] The embodiments of the present application further provide a computer program product or computer program, which includes computer instructions stored in a computer-readable storage medium. A processor of an electronic device reads the computer instructions from the computer-readable storage medium and executes the computer instructions, causing the electronic device to implement any of the above-described methods for predicting river channel evolution scale based on cross-sectional topography, for example:

[0161] Obtain observation data of the target river section at multiple measurements and the operating years of the dam in the target river section at the corresponding measurements, the observation data including the starting point distances and riverbed elevations of multiple observation points in multiple sections in the target river section at the corresponding measurements; determine the flat river channel morphological parameters of each section in the target river section at the corresponding measurements based on the observation data; determine a first calculation formula for the flat river channel morphological parameters, wherein the first calculation formula includes a balance value of the flat river channel morphological parameters at the river section scale and an attenuation coefficient of the flat river channel morphological parameters at the river section scale, both of which are calibrated based on the flat river channel morphological parameters of all sections in the target river section at all measurements and the operating years of the dam; determine the flat river channel morphological adjustment mode and time scale of the target river section based on the flat river channel morphological parameters at the river section scale, the balance value of the flat river channel morphological parameters at the river section scale and the attenuation coefficient of the flat river channel morphological parameters at the river section scale, and obtain a scale prediction result of the flat river channel morphological evolution.

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

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

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

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

Claims

1. A river channel evolution scale prediction method based on cross-sectional topography, characterized in that: The method comprises: Obtaining observation data of a target river section at multiple measurement times and the operating years of the dam in the target river section at the corresponding measurement times, wherein the observation data includes starting point distances and riverbed elevations of multiple observation points at multiple sections in the target river section at the corresponding measurement times; Determine the flat river channel morphological parameters of each section in the target river section at corresponding measurement times based on the observation data; Determining a first calculation formula for flat-shoal river channel morphological parameters, wherein the first calculation formula includes a reach-scale flat-shoal river channel morphological parameter equilibrium value and a reach-scale flat-shoal river channel morphological parameter attenuation coefficient, both of which are calibrated based on the flat-shoal river channel morphological parameters of all sections in the target river reach at all measurements and the operating life of the dam; Based on the flat-shoal river channel morphological parameters of the target river section at the river section scale, the equilibrium value of the flat-shoal river channel morphological parameters at the river section scale, and the attenuation coefficient of the flat-shoal river channel morphological parameters at the river section scale, the flat-shoal river channel morphological adjustment mode and time scale of the target river section are determined to obtain a flat-shoal river channel morphological evolution scale prediction result; The time scale for adjusting the flat river channel morphology of the target river section is determined by the following steps: For the target river section, determining a quantitative relationship between an attenuation coefficient of a river section-scale flat river channel morphological parameter and an adjustment equilibrium time of the river section-scale flat river channel morphological parameter; The attenuation coefficient of the flat beach river channel morphological parameters at the river section scale is input into the quantitative relationship to determine the adjustment equilibrium time of the flat beach river channel morphological parameters at the river section scale. The flat beach river channel morphological adjustment time scale includes the adjustment equilibrium time of the flat beach river channel morphological parameters at the river section scale. The adjustment equilibrium time of the flat beach river channel morphological parameters at the river section scale includes the adjustment equilibrium time of the flat beach cross-sectional area, flat beach water depth and flat beach river width at the river section scale.

2. The method for predicting river channel evolution scale based on cross-sectional topography according to claim 1, characterized in that: The reach-scale flat river channel morphological parameters include the reach-scale flat cross-sectional area, flat water depth, and flat river width; the reach-scale flat river channel morphological parameter balance values ​​include the flat cross-sectional area balance value, flat water depth balance value, and flat river width balance value; The flat river channel morphology adjustment mode of the target river section is determined by the following steps: Determining the relative adjustment intensity of the flat-shoal river channel morphological parameters at the river reach scale based on the equilibrium value of the flat-shoal river channel morphological parameters at the river reach scale and the initial value of the flat-shoal river channel morphological parameters at the river reach scale, wherein the relative adjustment intensity includes the relative adjustment intensity of the flat-shoal cross-sectional area, the relative adjustment intensity of the flat-shoal water depth, and the relative adjustment intensity of the flat-shoal river width; The flat beach river channel morphology adjustment mode is determined based on the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width.

3. The method for predicting river channel evolution scale based on cross-sectional topography according to claim 2, characterized in that: The determination of the flat beach river channel morphology adjustment mode based on the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width includes: Determine the radar chart between the relative adjustment intensity of the flat beach cross-sectional area, the relative adjustment intensity of the flat beach water depth, and the relative adjustment intensity of the flat beach river width; Coordinates are converted to the vertices in the radar map to respectively determine the relative adjustment intensity coordinate point of the flat beach water depth, the relative adjustment intensity coordinate point of the flat beach river width, and the vector between the center point of the radar map and the coordinate origin; Based on the relative adjustment intensity coordinate points of the flat beach water depth, the relative adjustment intensity coordinate points of the flat beach river width, and the vector between the center point of the radar chart and the coordinate origin, the flat beach river channel morphology adjustment mode is determined, and the flat beach river channel morphology adjustment mode includes at least one of the downcutting type adjustment, the widening type adjustment, and the comprehensive type adjustment.

4. The method for predicting river channel evolution scale based on cross-sectional topography according to claim 3, characterized in that: The determining of the flat beach river channel morphology adjustment mode based on the relative adjustment intensity coordinate point of the flat beach water depth, the relative adjustment intensity coordinate point of the flat beach river width, and the vector between the center point of the radar chart and the coordinate origin includes: Determine the first cross product value between the vector of the relative adjustment intensity coordinate point of the flat water depth and the coordinate origin, and the vector of the center point of the radar chart and the coordinate origin; Determine the second cross product value between the vector of the relative adjustment intensity coordinate point of Pingtan River width and the coordinate origin, and the vector of the center point of the radar chart and the coordinate origin; Determine a first ratio between the modulus of the vector of the relative adjustment intensity coordinate point of the flat beach water depth and the coordinate origin, and the modulus of the vector of the relative adjustment intensity coordinate point of the flat beach river width and the coordinate origin; The flat river channel morphology adjustment mode is determined based on the positive and negative values ​​of the first cross product value and the second cross product value and the size of the first ratio.

5. The method for predicting river channel evolution scale based on cross-sectional topography according to claim 1, characterized in that: Determining the quantitative relationship between the attenuation coefficient of the flat river channel morphological parameter at the river section scale and the adjustment equilibrium time of the flat river channel morphological parameter at the river section scale for the target river section includes: Determining, for each section in the target river reach, a second calculation formula for the flat-shoal river channel morphological parameters, wherein the second calculation formula includes a cross-sectional-scale flat-shoal river channel morphological parameter equilibrium value and a cross-sectional-scale flat-shoal river channel morphological parameter attenuation coefficient, both of which are calibrated based on the cross-sectional-scale flat-shoal river channel morphological parameters of the corresponding section in the target river reach and the operating life of the dam; For each section in the target river section, determining the time at which the adjustment amount of the flat river channel morphological parameter of the corresponding section in the target river section at the cross-sectional scale reaches the target adjustment amount, and using this time as the initial proposed adjustment equilibrium time, wherein the target adjustment amount is determined based on the difference between the equilibrium value of the flat river channel morphological parameter of the corresponding section in the target river section at the cross-sectional scale and the initial value of the flat river channel morphological parameter of the corresponding section in the target river section at the cross-sectional scale; The quantitative relationship is determined based on the attenuation coefficient of the flat river channel morphological parameters at the cross-sectional scale of all sections in the target river section and the initial adjustment equilibrium time.

6. The method for predicting river channel evolution scale based on cross-sectional topography according to claim 5, characterized in that: The quantitative relationship is determined based on the attenuation coefficient of the flat river channel morphological parameters of all sections in the target river section at the cross-sectional scale and the initial adjustment equilibrium time, including: Generate a scatter plot of the relationship between the attenuation coefficient of the flat river channel morphological parameter and the initial adjustment equilibrium time at the cross-sectional scale for all sections in the target river section; Determining upper and lower envelopes in the scatter plot; The data of the upper and lower envelopes are averaged and then curve fitting is performed to obtain the quantitative relationship.

7. The method for predicting river channel evolution scale based on cross-sectional topography according to claim 1, characterized in that: The first calculation formula is determined by the following formula: in, is the value of the flat river channel morphological parameters at (x, t), x is the distance from the dam to the corresponding section in the target river section, t is the operating life of the dam, Y e is the corresponding equilibrium value of the flat river channel morphological parameters, Φ(x, t) is the flat river channel morphological boundary function, Φ e is the equilibrium value of the flat river channel morphological boundary function, e is a natural constant, α is the attenuation coefficient of the corresponding flat river channel morphological parameters, v is the spatial propagation rate of the change of the corresponding flat river channel morphological parameters, and t0 is the initial value of the operating life of the dam.

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

9. A computer-readable storage medium, characterized in that A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps in the method for predicting river channel evolution scale based on cross-sectional topography according to any one of claims 1 to 7.

Images