Automatic alluvial river pattern judgment method based on river facies coefficient and related equipment

Abstract

The invention discloses an alluvial river pattern automatic discrimination method based on a river facies coefficient and related equipment. The alluvial river pattern automatic judgment method based on the river facies coefficient comprises the steps that actually-measured flow achievement data of a target river reach in multiple historical time periods are obtained, wherein the actually-measured flow achievement data comprise actually-measured flow and actually-measured hydraulic geometric morphology parameters; a river facies relational expression corresponding to each historical time period is determined, the river facies relational expression comprises river facies parameters, and the river facies relational expression represents the incidence relation between the actually-measured flow and the actually-measured hydraulic geometric morphology parameters of the historical time periods on the basis of the river facies parameters; and determining the alluvial river pattern information of the target river reach based on the river facies parameters in the river facies relational expressions corresponding to the plurality of different historical time periods. According to the method, the alluvial river pattern can be distinguished more efficiently, quickly, simply and conveniently.

Description

Technical Field

[0001] The present application relates to the fields of river evolution simulation, river management and protection, and specifically to an automatic identification method for alluvial river types based on river phase coefficients and related equipment. Background Art

[0002] The river type of an alluvial river is a spatial form that adapts to external conditions and is formed in the long-term automatic adjustment process of an alluvial river. It is closely related to the dynamic characteristics of river flow, sediment transport characteristics, and the evolution of riverbed morphology. Accurately determining the river type and evolution trend is of great significance for predicting the evolution characteristics of the river channel and formulating appropriate regulation plans.

[0003] In related technologies, the identification of alluvial river types can be achieved by analyzing factors such as river phase sediment distribution, mainstream swing flow time series characteristics, the longitudinal transmission mechanism of upstream and downstream river water level and riverbed elevation changes, and the propagation characteristics of mainstream pendulum waves in the riverbed medium. However, the acquisition cost of river phase sediment distribution data is high and difficult, and the research time and space scale is large, while the analysis process of mainstream swing characteristics and longitudinal transmission mechanisms has problems such as complex operation and high basic data requirements. It can be seen that the current identification process of alluvial river types is complex and requires a high level of technical personnel. Summary of the invention

[0004] The embodiments of the present application provide an automated identification method for an alluvial river type based on a river phase coefficient and related equipment, aiming to make the identification of an alluvial river type more efficient, rapid and simple.

[0005] On the one hand, the present application provides an automatic identification method of alluvial river type based on river phase coefficient, the method comprising: Obtaining measured flow results data of a target river section in multiple historical time periods, wherein the measured flow results data includes measured flow and measured hydraulic geometry parameters; Determine a river phase relationship equation corresponding to each of the historical time periods, wherein the river phase relationship equation includes a river phase parameter, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameter; Based on the river phase parameters in the river phase relationship equations corresponding to the multiple different historical time periods, the alluvial river type information of the target river section is determined.

[0006] In some embodiments, the measured hydraulic geometry parameters include measured river width, measured water depth, and measured flow velocity, the river phase parameter in the river phase relationship equation representing the correlation between the measured flow rate and the measured river width in the historical time period is a river width parameter, the river phase parameter in the river phase relationship equation representing the correlation between the measured flow rate and the measured water depth in the historical time period is a water depth parameter, and the river phase parameter in the river phase relationship equation representing the correlation between the measured flow rate and the measured flow velocity in the historical time period is a flow velocity parameter; Based on the river phase parameters in the river phase relationship equations corresponding to the multiple different historical time periods, the alluvial river type information of the target river section is determined, including: Based on the river width parameters, water depth parameters and flow velocity parameters corresponding to the multiple different historical time periods, the alluvial river type information of the target river section is determined.

[0007] In some embodiments, the river width parameter includes a river width index, the water depth parameter includes a water depth index, and the flow velocity parameter includes a flow velocity index. The alluvial river type information of the target river section is determined based on the river width parameters, water depth parameters, and flow velocity parameters corresponding to a plurality of different historical time periods, including: Generate a triangular relationship diagram between the river width index, the water depth index, and the flow velocity index, wherein the triangular relationship diagram includes a plurality of position points, each of which corresponds to a river width index, a water depth index, and a flow velocity index corresponding to the historical time period; Based on the triangular relationship diagram, the alluvial river type information of the target river section is determined.

[0008] In some embodiments, determining the alluvial river type information of the target river section based on the triangular relationship diagram includes: In the triangular relationship diagram, determine the boundary lines of different alluvial river types; Based on the relative position relationship between the dividing line and a plurality of position points in the triangular relationship diagram, the alluvial river type information of the target river section is determined.

[0009] In some embodiments, the preset historical period includes a plurality of different river morphology adjustment cycles, the alluvial river type information includes the alluvial river type of the target river section in each of the river morphology adjustment cycles, and the alluvial river type information of the target river section is determined based on the relative position relationship between the dividing line and the plurality of position points in the triangular relationship diagram, including: For each of the river channel morphology adjustment cycles, determining, among the multiple position points in the triangular relationship diagram, position points corresponding to multiple historical time periods within the river channel morphology adjustment cycle; Based on the relative position relationship between the dividing line and the position points corresponding to multiple historical time periods within the river channel morphology adjustment cycle, the alluvial river type of the target river section in the river channel morphology adjustment cycle is determined.

[0010] In some embodiments, the river channel morphology adjustment period is determined by the following steps: Obtaining the runoff, sediment transport and cumulative scouring and silting volume per unit river length of the target river section in each historical time period in the preset historical period; Based on the runoff, the sediment transport and the cumulative scouring and silting amount per unit river length, the preset historical period is divided to obtain a plurality of different river channel morphology adjustment cycles.

[0011] In some embodiments, determining the river phase relationship corresponding to each of the historical time periods includes: Obtaining a preset power function relationship, wherein the power function relationship includes an undetermined exponent and an undetermined coefficient; The power function relationship is used to fit the correlation between the measured flow and the measured hydraulic geometric parameters of the historical time period to determine the value of the unknown index and the value of the unknown coefficient, and the river phase relationship corresponding to the historical time period is obtained, and the river phase parameter is at least one of the unknown index and the unknown coefficient.

[0012] In some embodiments, before determining the river phase relationship equation corresponding to each of the historical time periods, the method further includes: Delete the measured flow results data that meets the preset conditions, wherein the preset conditions include at least one of the following: The historical time period corresponding to the measured flow results data is within the preset dry season; The measured flow in the measured flow result data is less than a flow threshold, and the flow threshold is determined based on the flat beach flow in a historical time period corresponding to the measured flow.

[0013] On the other hand, an embodiment of the present application provides an automatic identification device for alluvial river type based on river phase coefficient, comprising: The first acquisition module is used to obtain the measured flow results data of the target river section in multiple historical time periods, wherein the measured flow results data includes the measured flow and the measured hydraulic geometry parameters; A first determination module is used to determine a river phase relationship equation corresponding to each of the historical time periods, wherein the river phase relationship equation includes a river phase parameter, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameter; The second determination module is used to determine the alluvial river type information of the target river section based on the river phase parameters in the river phase relationship equations corresponding to multiple different historical time periods.

[0014] On the other hand, the present application also provides a computer device, the computer device comprising: one or more processors; Memory; and One or more applications, wherein the one or more applications are stored in the memory and are configured to be executed by the processor to implement any step in the method for automatically distinguishing the alluvial river type based on river phase coefficient.

[0015] On the other hand, the present application also provides a computer-readable storage medium having a computer program stored thereon, wherein the computer program is loaded by a processor to execute the steps in any one of the methods for automatic identification of alluvial river types based on river phase coefficients.

[0016] On the other hand, an embodiment of the present application provides a computer program product, including a computer program or instructions, which are executed by a processor to implement the method for automatically distinguishing the alluvial river type based on the river phase coefficient as described in any of the above items.

[0017] The embodiments of the present application provide an automated method and related equipment for distinguishing the alluvial river type based on the river phase coefficient, by obtaining the measured flow and measured hydraulic geometry parameters of the target river section in multiple historical time periods, and then determining a river phase relationship equation that characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time periods. Based on the river phase parameters in the river phase relationship equation corresponding to multiple different historical time periods, the alluvial river type information of the target river section is determined. Compared with river type analysis schemes based on river phase sedimentation distribution, mainstream swing characteristics, longitudinal transfer mechanism, etc., the method can make the distinction of the alluvial river type more efficient, rapid, and simple, and has lower requirements on the level of technical personnel. BRIEF DESCRIPTION OF THE DRAWINGS

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

[0019] Figure 1 It is a schematic diagram of an embodiment of an automatic identification method of an alluvial river type based on a river phase coefficient provided in an embodiment of the present application; Figure 2It is a schematic diagram of the change process of the river width index, water depth index and flow velocity index over time corresponding to multiple different historical time periods provided in the embodiments of the present application; Figure 3 is a schematic diagram of a triangular relationship diagram including multiple position points in an embodiment of the present application; Figure 4 It is a schematic diagram of the change process of the runoff and sediment transport of the target river section in the preset historical period over time in the embodiment of the present application; Figure 5 It is a schematic diagram of the cumulative scouring and silting volume per unit river length of the target river section in the preset historical period over time in the embodiment of the present application; Figure 6 It is a schematic diagram of the structure of an embodiment of an automatic alluvial river type identification device based on river phase coefficient provided in an embodiment of the present application; Figure 7 It is a schematic diagram of the terminal structure of an embodiment of the computer device provided in the embodiments of the present application.

[0020] in, Figure 2 , 3 , 4 are color pictures to facilitate the distinction of different objects in the same picture. DETAILED DESCRIPTION

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

[0022] In the description of this application, the terms "first" and "second" are used for descriptive purposes only and should not be understood as indicating or implying relative importance or implicitly indicating the number of technical features indicated. Thus, a feature defined as "first" or "second" may explicitly or implicitly include one or more of the features. In the description of this application, "plurality" means two or more, unless otherwise clearly and specifically defined.

[0023] In this application, the term "in some embodiments" is used to mean "used as an example, illustration or description". Any embodiment described in this application as "in some embodiments" is not necessarily to be construed as being more preferred or advantageous than other embodiments. In order to enable any technician in the field to implement and use the present application, the following description is given. 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 examples, well-known structures and processes will not be 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 the present application.

[0024] 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, they are all corresponding data for processing by the computer device. The details will not be repeated here.

[0025] The embodiment of the present application provides an automatic identification method and related equipment for the river type of an alluvial river based on a river phase coefficient, which aims to obtain the relationship (y=αQ) between the measured hydraulic geometric parameters (measured river width B, measured water depth H, measured flow velocity V) and the measured flow Q based on the measured flow results data of the hydrological observation station of the target river section, and on the basis of analyzing the adjustment rules of the measured river width B, measured water depth H, and measured flow velocity V in different periods, using the cross-sectional river phase relationship in the form of a classic power function, and through nonlinear fitting. β ), by identifying the triangular relationship between the river width index β1, the water depth index β2 and the flow velocity index β3, the alluvial river type information of the target river section is determined, and the river type adjustment direction of the target river section is identified, so as to realize the efficient, fast and simple judgment of the alluvial river type and its evolution trend, which are described in detail below.

[0026] In one embodiment, referring to Figure 1 , the automatic identification methods of alluvial river types based on river phase coefficients include: 101. Obtain the measured flow results data of the target river section in multiple historical time periods, the measured flow results data including the measured flow and measured hydraulic geometry parameters.

[0027] In the embodiment of the present application, a hydrological observation station is set up in the target river section, and the hydrological observation station can obtain measured flow results data through actual observation. Each historical time period can be an observation year of the hydrological observation station, and the measured flow results data of multiple observation years form a long sequence of measured flow results data.

[0028] In some embodiments of the present application, the measured hydraulic geometry parameters include at least one of the measured river width B, the measured water depth H, and the measured flow velocity V.

[0029] In some embodiments of the present application, step 101 is illustrated. Specifically, step 101 may include: Step 1.1, collect and organize the measured flow results data of the representative hydrological observation stations of the target river section. For a specific hydrological observation station, set the observation year as yy and the number of observation points as nn, then the measured flow results data can be obtained, including the two-dimensional array of measured flow Q, measured river width B, measured water depth H, and measured flow velocity V:

[0030] Where Q is the measured flow rate; B is the measured river width; H is the measured water depth; V is the measured flow velocity; i is the observation year; j is the number of observation points of the measured flow data in each observation year. Where i=1, 2, …, yy; j=1, 2, …, nn.

[0031] Step 1.2, standardization of measured flow results data. Specifically, the measured flow results data of the same hydrological observation station in different observation years are stored in the same Excel file (forming a long sequence measured flow results table) in the format of observation year i, measured flow Q, measured river width B, measured water depth H and measured flow velocity V, while the measured flow results data of different hydrological observation stations are named after the hydrological observation station and stored in different Excel files.

[0032] 102. Determine a river phase relationship equation corresponding to each historical time period, wherein the river phase relationship equation includes river phase parameters, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameters.

[0033] In the embodiments of the present application, for each historical time period, it is necessary to determine the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period, and thus the corresponding river phase relationship equation can be determined by fitting the formula of the correlation relationship.

[0034] Among them, if the measured hydraulic geometric parameters include at least one of the measured river width B, the measured water depth H, and the measured flow velocity V, then the river phase parameter in the river phase relationship equation that characterizes the correlation between the measured flow rate Q and the measured river width B in the historical time period is the river width parameter, the river phase parameter in the river phase relationship equation that characterizes the correlation between the measured flow rate Q and the measured water depth H in the historical time period is the water depth parameter, and the river phase parameter in the river phase relationship equation that characterizes the correlation between the measured flow rate Q and the measured flow velocity V in the historical time period is the flow velocity parameter.

[0035] In some embodiments of the present application, the river phase relationship may be, for example, a formula form of a power function relationship, and the river phase parameter in the river phase relationship may be, for example, at least one of an undetermined exponent and an undetermined coefficient in a preset power function relationship. Specifically, step 102 may include: obtaining a preset power function relationship, the power function relationship including an undetermined exponent and an undetermined coefficient, for example, the power function relationship may be y=αQ β , y is the measured hydraulic geometry parameter, Q is the measured flow, α is the unknown coefficient, and β is the unknown exponent; the power function relationship is used to fit the correlation between the measured flow and the measured hydraulic geometry parameters in the historical time period to determine the value of the unknown coefficient α and the value of the unknown specified β, and the river phase relationship corresponding to the historical time period is obtained.

[0036] In some embodiments of the present application, before executing step 102, it may also include: deleting the measured flow result data that meets the preset conditions to filter the measured flow result data. The preset conditions include at least one of the following: the historical time period corresponding to the measured flow result data is within the preset dry season, and the preset dry season may be, for example, January to April and December; the measured flow in the measured flow result data is less than the flow threshold, and the flow threshold is determined based on the flat beach flow in the historical time period corresponding to the measured flow. For example, the flat beach flow can be multiplied by a preset ratio to obtain the flow threshold, and the preset ratio may be, for example, 15%. By filtering the measured flow result data under the preset conditions, the accuracy of distinguishing the river type of alluvial rivers can be further improved.

[0037] In some embodiments of the present application, step 102 and the filtering of the measured flow results data are illustrated. Specifically, considering that there are more data with smaller measured flow in the measured flow results data, which will affect the fitting result of the river phase relationship equation, before using the nonlinear fitting function to adopt the least squares fitting method to solve the river phase relationship equation, the measured flow results data from January to April and December and the measured flow results data with a measured flow less than 15% of the flat beach flow are eliminated. Based on the measured flow results data from May to November with a measured flow greater than 15% of the flat beach flow, the undetermined exponent and undetermined coefficient (i.e., B=α1Q β1 、H=α2Q β2 , V = α3Q β3The river width index β1, water depth index β2, flow velocity index β3, river width coefficient α1, water depth coefficient α2 and flow velocity coefficient α3 in the above-mentioned river width parameter include at least one of the river width index β1 and the river width coefficient α1, the water depth parameter include at least one of the water depth index β2 and the water depth coefficient α2, and the flow velocity parameter include at least one of the flow velocity index β3 and the flow velocity coefficient α3). The for loop statement in the programming language is used to calculate the values ​​of the unknown exponents and unknown coefficients in the power function relationship of the measured flow Q, the measured river width, the measured water depth and the measured flow velocity in different observation years.

[0038] In addition, using the river phase relationship corresponding to the specified historical time period, the value of the river phase parameter in the river phase relationship can be reversely calculated to obtain the specific flow rate (e.g. 3000m 3 / s) under specific hydraulic geometry parameters (such as specific river width, specific water depth and specific flow velocity).

[0039] 103. Based on the river phase parameters in the river phase relationship equations corresponding to multiple different historical time periods, the alluvial river type information of the target river section is determined.

[0040] In an embodiment of the present application, the alluvial river type information of the target river section can be determined based on the values ​​of the river phase parameters in the river phase relationship equation corresponding to multiple different historical time periods, that is, the alluvial river type information of the target river section is related to the river phase parameters.

[0041] In some embodiments of the present application, step 103 may include: determining the alluvial river type information of the target river section based on the river width parameters, water depth parameters and flow rate parameters corresponding to multiple different historical time periods, that is, the alluvial river type information of the target river section is related to the river width parameters, water depth parameters and flow rate parameters. Taking the example that the river width parameter includes the river width index β1, the water depth parameter includes the water depth index β2, and the flow rate parameter includes the flow rate index β3, a schematic diagram of the change of the river width index β1, the water depth index β2, and the flow rate index β3 corresponding to multiple different historical time periods over time is as follows: Figure 2 as shown in .

[0042] In some embodiments of the present application, in order to more conveniently determine the alluvial river type information of the target river section, a triangular relationship diagram may be used for analysis. Specifically, based on the river width parameters, water depth parameters, and flow velocity parameters corresponding to multiple different historical time periods, determining the alluvial river type information of the target river section may include: generating a triangular relationship diagram between the river width index, the water depth index, and the flow velocity index (e.g. Figure 3As shown in , the coordinate axes corresponding to the river width index, water depth index and flow velocity index are respectively one side of the triangle in the triangular relationship diagram, and the triangular relationship diagram includes multiple position points, each of which corresponds to the river width index, water depth index and flow velocity index corresponding to a historical time period, that is, the specific position of each position point is determined based on the values ​​of the river width index, water depth index and flow velocity index corresponding to the corresponding historical time period; based on the triangular relationship diagram, the alluvial river type information of the target river section is determined, thereby making the determination of the alluvial river type information simpler and more intuitive.

[0043] In some embodiments of the present application, determining the alluvial river type information of the target river section based on the triangular relationship diagram may include: determining the boundary line of different alluvial river types in the triangular relationship diagram, and the boundary line may be set based on experience. For example, when the alluvial river type includes wandering type and non-wandering type, the line in the triangular relationship diagram where the river width index is equal to the water depth index (i.e., β1=β2) may be used as the boundary line between the wandering type and the non-wandering type (e.g. Figure 3 Based on the relative position relationship between the dividing line and multiple points in the triangular relationship diagram, the alluvial river type information of the target river section can be determined. For example, Figure 3 The historical time periods corresponding to the position points in the area to the left of the dotted line in (β1>β2) are determined, and the alluvial river type of the target river section in these historical time periods is determined as the wandering type. For example, Figure 3 The historical time periods corresponding to the location points in the area to the right of the dotted line in (β1<β2) are used, and the alluvial river type of the target river section during these historical time periods is determined as non-wandering type.

[0044] In some embodiments of the present application, since the alluvial river type determined based on the position points corresponding to a single historical time period may not be accurate enough, and the alluvial river type is formed in the long-term adjustment process of the target river section, the position points corresponding to multiple historical time periods may be integrated to determine a more accurate alluvial river type. Specifically, all historical time periods are within a preset historical period (e.g., 1965 to 2015), so the preset historical period can be divided to obtain multiple different river morphology adjustment cycles in the preset historical period (e.g., 1965 to 1973, 1974 to 1980, 1981 to 1985, 1986 to 1999, and 2000 to 2015), that is, it is considered that the alluvial river type of the target river section in the same river morphology adjustment cycle is the same, so the alluvial river type information can include the alluvial river type of the target river section in each river morphology adjustment cycle. Accordingly, determining the alluvial river type information of the target river section based on the relative position relationship between the dividing line and the multiple position points in the triangular relationship diagram may include: for each river channel morphology adjustment period, determining the position points corresponding to the multiple historical time periods within the river channel morphology adjustment period from the multiple position points in the triangular relationship diagram; determining the alluvial river type of the target river section in the river channel morphology adjustment period based on the relative position relationship between the dividing line and the position points corresponding to the multiple historical time periods within the river channel morphology adjustment period, for example, Figure 3 In the figure, most of the position points corresponding to multiple historical time periods in the river morphology adjustment cycle from 1965 to 1973 (for example, more than 50%) are located in the area on the left side of the dotted line (β1>β2), so it can be determined that the alluvial river type of the target river section from 1965 to 1973 is a wandering type. For the four river morphology adjustment cycles from 1974 to 1980, 1981 to 1985, 1986 to 1999, and 2000 to 2015, most of the position points corresponding to multiple historical time periods in each river morphology adjustment cycle (for example, more than 50%) are located in the area on the right side of the dotted line (β1<β2), so it can be determined that the alluvial river type of the target river section in the four river morphology adjustment cycles from 1974 to 1980, 1981 to 1985, 1986 to 1999, and 2000 to 2015 is a non-wandering type.

[0045] In addition, the alluvial river type information of the target river section may also include the adjustment direction of the alluvial river type. For example, the adjustment direction of the alluvial river type may be determined based on the alluvial river type of the target river section in multiple river channel morphology adjustment cycles. The adjustment direction may be, for example, the alluvial river type of the target river section changing from a wandering type before 1973 to a non-wandering type after 1973.

[0046] In some embodiments of the present application, the division method of the preset historical period is described. Specifically, the river morphology adjustment period is determined by the following steps: obtaining the runoff, sediment transport and cumulative scouring and silting volume per unit river length of the target river section in each historical time period in the preset historical period, wherein the runoff and sediment transport of the target river section in each historical time period in the preset historical period are, for example, Figure 4 As shown in , the cumulative scouring and silting volume per unit river length of the target river section in each historical time period during the preset historical period is, for example, Figure 5 As shown in ; based on the runoff, sediment transport and cumulative scouring and silting per unit river length, the preset historical period is divided to obtain a plurality of different river morphology adjustment cycles. For example, based on the runoff, sediment transport and cumulative scouring and silting per unit river length, combined with the human activities of the target river section in each historical time period in the preset historical period, relevant technical personnel can identify the phased change characteristics of the water and sediment conditions and riverbed evolution of the target river section, and then divide the preset historical period into a plurality of different river morphology adjustment cycles (for example, the five river morphology adjustment cycles from 1965 to 1973, 1974 to 1980, 1981 to 1985, 1986 to 1999 and 2000 to 2015).

[0047] In the technical solution disclosed in this embodiment, by obtaining the measured flow and measured hydraulic geometry parameters of the target river section in multiple historical time periods, and then determining the river phase relationship equation that characterizes the correlation between the measured flow and the measured hydraulic geometry parameters in the historical time periods, based on the river phase parameters in the river phase relationship equation corresponding to multiple different historical time periods, the alluvial river type information of the target river section is determined. Compared with the river type analysis scheme based on river phase sedimentation distribution, mainstream swing characteristics, longitudinal transfer mechanism, etc., the identification of the alluvial river type can be made more efficient, rapid, and simple, and has lower requirements on the level of technical personnel.

[0048] Below, refer to Figure 1 , an example is given to illustrate the automatic identification method of alluvial river type based on river phase coefficient. Specifically, the evolution trend of the alluvial river type of the target river section under the change of water and sediment conditions is identified, including the following steps: Step S1, standardization processing of the long sequence measured flow result table to form a data set including measured flow Q, measured river width B, measured water depth H and measured flow velocity V.

[0049] Collect and organize the data of the long-sequence measured flow results table of representative hydrological observation stations in the target river section, set the observation time and observation point number, record the observation year yy=1965 to 2015, the number of observation points nn=142, and form a two-dimensional array including measured flow Q, measured river width B, measured water depth H and measured flow velocity V. The two-dimensional array can be, for example, Q(1965,142), B(1965,142), H(1965,142), V(1965,142). The measured flow results data of the same hydrological observation station in different observation years are stored in the same Excel file in the format of year, month, day, measured flow Q, measured river width B, measured water depth H and measured flow velocity V, as shown in Table 1 below: Table 1 Measured flow results Standardized data sets for each observation year

[0050] Step S2, river channel morphology adjustment period identification.

[0051] Collect and organize data on the long-term water and sediment conditions of the target river section (such as runoff and sediment transport), riverbed erosion and siltation (such as the cumulative erosion and siltation per unit river length), and human activities. Figure 4 As shown in the figure, the cumulative scouring and silting volume per unit river length of the target river section length sequence is as follows: Figure 5 As shown in the figure, combined with human activities such as reservoir construction in the target river section, the river channel evolution process of the target river section from 1965 to 2015 is divided into five different river channel morphology adjustment periods: 1965 to 1973, 1974 to 1980, 1981 to 1985, 1986 to 1999 and 2000 to 2015.

[0052] Step S3, fitting determines the values ​​of the exponents (such as the river width index β1, the water depth index β2, and the flow velocity index β3) and coefficients (such as the river width coefficient α1, the water depth coefficient α2, and the flow velocity coefficient α3) of the long sequence of river phase relationship equations, as well as the hydraulic geometric parameters under a specific flow rate.

[0053] Considering that there are too many data with smaller measured flow in the measured flow results data, which will affect the fitting results of the river phase relationship, the measured flow results data from January to April and December and the measured flow results data with a measured flow less than 15% of the flat beach flow are eliminated before solving the river phase relationship equation using the least squares fitting method using the nonlinear fitting function. Based on the measured flow results data from May to November with a measured flow greater than 15% of the flat beach flow, the undetermined exponent and undetermined coefficient in the power function relationship (i.e. B=α1Q β1 、H=α2Q β2 , V = α3Q β3The river width index β1, water depth index β2, flow velocity index β3, river width coefficient α1, water depth coefficient α2 and flow velocity coefficient α3, wherein the above-mentioned river width parameters include at least one of the river width index β1 and the river width coefficient α1, the above-mentioned water depth parameters include at least one of the water depth index β2 and the water depth coefficient α2, and the above-mentioned flow velocity parameters include at least one of the flow velocity index β3 and the flow velocity coefficient α3). The for loop statement in the programming language MATLAB is used to calculate the values ​​of the unknown exponents and unknown coefficients in the power function relationship of the measured flow rate Q, the measured river width, the measured water depth, and the measured flow velocity in different observation years. In addition, according to the river phase relationship corresponding to different observation years, the target river section can be calculated at a specific flow rate (for example, 3000m 3 / s) under specific hydraulic geometry parameters (such as specific river width, specific water depth and specific flow velocity).

[0054] Step S4, plotting the time variation of river phase parameters corresponding to different observation years.

[0055] Using the plot function in the programming language MATLAB, a schematic diagram of the river width index β1, water depth index β2, and flow velocity index β3 changing with time is shown, for example Figure 2 as shown in .

[0056] Step S5, determining the adjustment direction of the river channel type (ie, alluvial river type).

[0057] According to the river morphology adjustment cycle determined in step S2, the river width index β1, water depth index β2, and flow velocity index β3 are divided into five different stages: 1965 to 1973, 1974 to 1980, 1981 to 1985, 1986 to 1999, and 2000 to 2015.

[0058] The river width index β1, water depth index β2, and flow velocity index β3 of each observation year in the above five stages are plotted into the triangular relationship diagram between the river width index β1, water depth index β2, and flow velocity index β3, as shown in Figure 3 As shown. The dotted line β1=β2 is used as the characteristic line to distinguish the wandering type from the non-wandering type. The area on the left side of the dotted line (β1>β2) is the wandering type, and the area on the right side of the dotted line (β1<β2) is the non-wandering type. It can be seen that the alluvial river type of the target river section has basically changed from the wandering type before 1973 to the non-wandering type after 1973.

[0059] The embodiment of the present application provides an automated identification method for alluvial river types based on river phase coefficients. The method is based on the changing characteristics of the exponents in the classical river phase relationship equation in the form of a power function between the measured flow rate, measured river width, measured water depth and measured flow velocity in a long sequence of sections, and utilizes the triangular relationship between the exponents in each river phase relationship equation to explore the adjustment direction of the alluvial river type. The technology is based on the measured flow results data of the hydrological observation station, and has low difficulty in obtaining data, low cost, high accuracy, and low requirements on the level of technical personnel. The selected characteristic factors (such as river width index, water depth index, flow velocity index) take into account the lateral and vertical adjustment characteristics of the river channel and the longitudinal water flow dynamic characteristics, and can achieve efficient, rapid and simple identification of the alluvial river type and its evolution trend, which is of great significance for predicting the evolution characteristics of the river channel and formulating appropriate regulation plans.

[0060] In order to better implement the alluvial river type automatic discrimination method based on river phase coefficient in the embodiment of the present application, on the basis of the alluvial river type automatic discrimination method based on river phase coefficient, the embodiment of the present application also provides an alluvial river type automatic discrimination device based on river phase coefficient, such as Figure 6 As shown, the alluvial river type automatic identification device 600 based on river phase coefficient includes: The first acquisition module 601 is used to obtain the measured flow results data of the target river section in multiple historical time periods, and the measured flow results data includes the measured flow and the measured hydraulic geometry parameters; A first determination module 602 is used to determine a river phase relationship equation corresponding to each historical time period, wherein the river phase relationship equation includes river phase parameters, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameters; The second determination module 603 is used to determine the alluvial river type information of the target river section based on the river phase parameters in the river phase relationship equations corresponding to multiple different historical time periods.

[0061] The present application also provides a computer device that integrates any of the alluvial river type automatic identification devices based on river phase coefficients provided in the present application. Figure 7 As shown, it shows a schematic diagram of the structure of the computer device involved in the embodiment of the present application, specifically: 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 appreciate 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 the components differently. The processor 701 is the control center of the computer device. It uses various interfaces and lines to connect various parts of the entire computer device. By running or executing software programs and / or modules stored in the memory 702 and calling data stored in the memory 702, it executes various functions of the computer device and processes data, thereby monitoring the computer device as a whole. Optionally, the processor 701 may include one or more processing cores; preferably, the processor 701 may integrate an application processor and a modem processor, wherein the application processor mainly processes the operating system, user interface, and application programs, and the modem processor mainly processes wireless communications. It is understandable that the above-mentioned modem processor may not be integrated into the processor 701.

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

[0063] 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 through a power management system, so as to manage charging, discharging, power consumption and other functions through the power management system. The power supply 703 can also include one or more DC or AC power supplies, recharging systems, power failure detection circuits, power converters or inverters, power status indicators and other arbitrary components.

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

[0065] Although not shown, the computer device may further include a display unit, etc., which will not be described in detail herein. Specifically in this embodiment, the processor 701 in the computer device will load the executable files corresponding to the processes of one or more application programs into the memory 702 according to the following instructions, and the processor 701 will run the application programs stored in the memory 702, thereby realizing various functions, as follows: Obtain the measured flow results data of the target river section in multiple historical time periods, the measured flow results data including the measured flow and measured hydraulic geometry parameters; determine the river phase relationship equation corresponding to each historical time period, the river phase relationship equation includes river phase parameters, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameters; determine the alluvial river type information of the target river section based on the river phase parameters in the river phase relationship equation corresponding to multiple different historical time periods.

[0066] A person of ordinary skill in the art will appreciate that all or part of the steps in the various methods of the above embodiments may be completed 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.

[0067] To this end, the embodiment of the present application provides a computer-readable storage medium, which may include: a read-only memory (ROM), a random access memory (RAM), a disk or an optical disk, etc. A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps in any one of the methods for automatically distinguishing the river type of alluvial rivers based on river phase coefficients provided in the embodiment of the present application. For example, the computer program can be loaded by a processor to execute the following steps: Obtain the measured flow results data of the target river section in multiple historical time periods, the measured flow results data including the measured flow and measured hydraulic geometry parameters; determine the river phase relationship equation corresponding to each historical time period, the river phase relationship equation includes river phase parameters, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameters; determine the alluvial river type information of the target river section based on the river phase parameters in the river phase relationship equation corresponding to multiple different historical time periods.

[0068] The embodiment of the present application also provides a computer program product or a computer program, which includes a computer instruction stored in a computer-readable storage medium. A processor of an electronic device reads the computer instruction from the computer-readable storage medium and executes the computer instruction, so that the electronic device executes to implement the method for automatically distinguishing the river type of an alluvial river based on a river phase coefficient as described in any one of the above items, for example: Obtain the measured flow results data of the target river section in multiple historical time periods, the measured flow results data including the measured flow and measured hydraulic geometry parameters; determine the river phase relationship equation corresponding to each historical time period, the river phase relationship equation includes river phase parameters, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameters; determine the alluvial river type information of the target river section based on the river phase parameters in the river phase relationship equation corresponding to multiple different historical time periods.

[0069] In the above embodiments, the description of each embodiment has its own emphasis. 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.

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

[0071] The specific implementation of the above operations can be found in the previous embodiments, which will not be described in detail here.

[0072] The above is a detailed introduction to an automated identification method for alluvial river types based on river phase coefficients and related equipment provided in an embodiment of the present application. Specific examples are used herein to illustrate the principles and implementation methods of the present application. The description of the above embodiments is only used to help understand the method of the present application and its core idea. At the same time, for those skilled in the art, according to the idea 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. An automatic identification method of alluvial river type based on river phase coefficient, characterized in that: The method comprises: Obtaining measured flow results data of a target river section in multiple historical time periods, wherein the measured flow results data includes measured flow and measured hydraulic geometry parameters; Determine a river phase relationship equation corresponding to each of the historical time periods, wherein the river phase relationship equation includes a river phase parameter, and the river phase relationship equation characterizes the correlation between the measured flow and the measured hydraulic geometry parameters of the historical time period based on the river phase parameter; Based on the river phase parameters in the river phase relationship equations corresponding to the multiple different historical time periods, the alluvial river type information of the target river section is determined.

2. The method for automatically distinguishing the alluvial river type based on the river phase coefficient according to claim 1, characterized in that: The measured hydraulic geometry parameters include measured river width, measured water depth, and measured flow velocity. The river phase parameter in the river phase relationship equation that characterizes the correlation between the measured flow rate and the measured river width in the historical time period is a river width parameter. The river phase parameter in the river phase relationship equation that characterizes the correlation between the measured flow rate and the measured water depth in the historical time period is a water depth parameter. The river phase parameter in the river phase relationship equation that characterizes the correlation between the measured flow rate and the measured flow velocity in the historical time period is a flow velocity parameter. Based on the river phase parameters in the river phase relationship equations corresponding to the multiple different historical time periods, the alluvial river type information of the target river section is determined, including: Based on the river width parameters, water depth parameters and flow velocity parameters corresponding to the multiple different historical time periods, the alluvial river type information of the target river section is determined.

3. The method for automatically distinguishing the alluvial river type based on river phase coefficient according to claim 2, characterized in that: The river width parameter includes a river width index, the water depth parameter includes a water depth index, and the flow velocity parameter includes a flow velocity index. The alluvial river type information of the target river section is determined based on the river width parameters, water depth parameters, and flow velocity parameters corresponding to a plurality of different historical time periods, including: Generate a triangular relationship diagram between the river width index, the water depth index, and the flow velocity index, wherein the triangular relationship diagram includes a plurality of position points, each of which corresponds to a river width index, a water depth index, and a flow velocity index corresponding to the historical time period; Based on the triangular relationship diagram, the alluvial river type information of the target river section is determined.

4. The method for automatically distinguishing the alluvial river type based on river phase coefficient according to claim 3, characterized in that: The determining of the alluvial river type information of the target river section based on the triangular relationship diagram includes: In the triangular relationship diagram, determine the boundary lines of different alluvial river types; Based on the relative position relationship between the dividing line and a plurality of position points in the triangular relationship diagram, the alluvial river type information of the target river section is determined.

5. The method for automatically distinguishing the alluvial river type based on river phase coefficient according to claim 4, characterized in that: The preset historical period includes a plurality of different river morphology adjustment cycles, the alluvial river type information includes the alluvial river type of the target river section in each of the river morphology adjustment cycles, and the alluvial river type information of the target river section is determined based on the relative position relationship between the dividing line and a plurality of position points in the triangular relationship diagram, including: For each of the river channel morphology adjustment cycles, determining, among the multiple position points in the triangular relationship diagram, position points corresponding to multiple historical time periods within the river channel morphology adjustment cycle; Based on the relative position relationship between the dividing line and the position points corresponding to multiple historical time periods within the river channel morphology adjustment cycle, the alluvial river type of the target river section in the river channel morphology adjustment cycle is determined.

6. The method for automatically distinguishing the alluvial river type based on river phase coefficient according to claim 5, characterized in that: The river channel morphology adjustment period is determined by the following steps: Obtaining the runoff, sediment transport and cumulative scouring and silting volume per unit river length of the target river section in each historical time period in the preset historical period; Based on the runoff, the sediment transport and the cumulative scouring and silting amount per unit river length, the preset historical period is divided to obtain a plurality of different river channel morphology adjustment cycles.

7. The method for automatically distinguishing the alluvial river type based on river phase coefficient according to claim 1, characterized in that: Determining the river phase relationship corresponding to each of the historical time periods includes: Obtaining a preset power function relationship, wherein the power function relationship includes an undetermined exponent and an undetermined coefficient; The power function relationship is used to fit the correlation between the measured flow and the measured hydraulic geometric parameters of the historical time period to determine the value of the unknown index and the value of the unknown coefficient, and the river phase relationship corresponding to the historical time period is obtained, and the river phase parameter is at least one of the unknown index and the unknown coefficient.

8. The method for automatically distinguishing the alluvial river type based on river phase coefficient according to claim 1, characterized in that: Before determining the river phase relationship equation corresponding to each of the historical time periods, the method further includes: Delete the measured flow results data that meets the preset conditions, wherein the preset conditions include at least one of the following: The historical time period corresponding to the measured flow results data is within the preset dry season; The measured flow in the measured flow result data is less than a flow threshold, and the flow threshold is determined based on the flat beach flow in a historical time period corresponding to the measured flow.

9. 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 method for automatic identification of alluvial river types based on river phase coefficients as described in any one of claims 1 to 8.

10. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps in the method for automatically distinguishing the alluvial river type based on river phase coefficient as described in any one of claims 1 to 8.

Images