An automated data processing method, device and medium for flat river channels
Through automated data processing methods, historical observation data are used to determine the beach lip position and the beach water level of the flat beach river trough, solving the problem of low efficiency and accuracy of flat beach water level recognition in the existing technology, and achieving more efficient and accurate flat beach water level recognition.
Patent Information
- Application Number
- CN202411818397.5
- Authority / Receiving Office
- CN · China
- Patent Type
- Patents(China)
- Current Assignee / Owner
- Filing Date
- 2024-12-11
- Publication Date
- 2025-03-14
- Estimated Expiration
- 2044-12-11
AI Technical Summary
In the prior art, the identification efficiency and accuracy of flat beach water levels are low, making it difficult to meet the needs of actual production and life, especially when multi-sectional processing.
An automated data processing method is provided, by obtaining the starting point distance and riverbed elevation of multiple observation points of the target section at multiple historical observations, determining the lowest point elevation, and then determining the beach lip position and flat beach water level, and determining the actual flat beach water level in multiple historical observations.
The recognition efficiency and accuracy of flat beach water level are improved, making the identified actual flat beach water level more accurate, and the recognition efficiency is higher than that of manual judgment.
Smart Images

Figure CN119290113B_ABST
Abstract
Description
Technical Field
[0001] The present application relates to the technical field of flat river channels, and in particular to an automated data processing method, device and medium for flat river channels. Background Art
[0002] A flat river channel refers to a characteristic river channel morphology under a flat water level. It is closely related to the flow capacity of the river channel and the stability of the riverbed. It is an important basis for the design of stable river channel geometry, river management and restoration. However, affected by human or natural factors such as climate change and reservoir construction, the water and sediment conditions entering the river channel are often changing, causing drastic adjustments in the flat river channel morphology, affecting the ecological health of the river and the safety of the project. Therefore, it is of great theoretical and practical value to accurately and quickly predict the morphological parameters of the flat river channel.
[0003] Accurate and rapid prediction of flat river channel morphological parameters involves accurate identification of flat water level. In related technologies, the identification of flat water level is generally based on human experience to judge the river section. Although this method can give full play to the accumulated practical experience of practitioners, it is highly subjective and has low work efficiency, which is difficult to meet the needs of actual production and life (actual work often involves dozens or even hundreds of sections). It can be seen that the current identification efficiency and accuracy of flat water level are low. Summary of the invention
[0004] The embodiments of the present application provide an automated data processing method, device and medium for a flat river channel, aiming to improve the efficiency and accuracy of identifying the flat water level.
[0005] On the one hand, the present application provides an automated data processing method for a flat river channel, the method comprising:
[0006] Acquire historical observation data of a target section in a target river section, wherein the historical observation data includes starting point distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times;
[0007] For each of the historical observation times, determining the lowest point elevation among the riverbed elevations of the plurality of observation points at the historical observation times;
[0008] Determine the beach lip position of the target section at the time of the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the time of the historical observation times;
[0009] Determine the flat beach water level of the target section at the historical observation times based on the riverbed elevation at the beach lip position of the target section at the historical observation times;
[0010] The actual flat water level of the target section is determined according to the flat water level of the target section at multiple historical observation times.
[0011] In some embodiments, the determining the beach lip position of the target section at the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times includes:
[0012] The observation point where the lowest point elevation is located is used as the deep point of the target section at the historical observation times;
[0013] According to the assumed calculated water level and the riverbed elevations of the plurality of observation points at the historical observation times, the observation point sequence is divided to obtain a plurality of subsequences, wherein the observation point sequence includes the plurality of observation points at the historical observation times sorted according to the distances from the starting points;
[0014] Determine a target subsequence where the deep point is located among the plurality of subsequences;
[0015] Based on the riverbed elevations of the plurality of observation points in the target subsequence, the beach lip position of the target section at the time of the historical observation times is determined.
[0016] In some embodiments, determining the beach lip position of the target section at the historical observation times based on the riverbed elevations of the plurality of observation points in the target subsequence includes:
[0017] Determine a maximum starting point distance and a minimum starting point distance among the starting point distances of a plurality of observation points in the target subsequence;
[0018] Detecting whether there is a target observation point that satisfies a first preset condition or a second preset condition, wherein the first preset condition includes that the target observation point does not exist in the target subsequence, the target observation point is adjacent to the observation point where the maximum starting point distance is located in the observation point sequence, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the maximum starting point distance is located satisfy a first preset size relationship, and the second preset condition includes that the target observation point does not exist in the target subsequence, the target observation point is adjacent to the observation point where the minimum starting point distance is located in the observation point sequence, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the minimum starting point distance is located satisfy a second preset size relationship;
[0019] If the target observation point exists, adding the target observation point to the target subsequence to obtain an expanded target subsequence;
[0020] Based on the multiple observation points in the expanded target subsequence, the beach lip position of the target section at the historical observation times is determined.
[0021] In some embodiments, determining the beach lip position of the target section at the historical observation times based on the multiple observation points in the expanded target subsequence includes:
[0022] Determine, among the multiple observation points of the expanded target subsequence, a first observation point with the largest starting point distance and a second observation point with the smallest starting point distance;
[0023] Determine the beach lip position of the target section at the historical observation times, wherein the beach lip position includes the starting point distance and riverbed elevation of the first observation point, and the starting point distance and riverbed elevation of the second observation point.
[0024] In some embodiments, determining the actual flat water level of the target section according to the flat water level of the target section at a plurality of the historical observation times includes:
[0025] Generate a cross-sectional view of the target section based on the flat beach water level of the target section at multiple historical observation times and the beach lip position of the target section at multiple historical observation times;
[0026] In the cross-sectional diagram, the actual flat water level of the target section is determined.
[0027] In some embodiments, generating a cross-sectional view of the target section based on the flat beach water level of the target section at multiple historical observation times and the beach lip position of the target section at multiple historical observation times includes:
[0028] For each of the historical observation times, determining a plurality of third observation points among the plurality of observation points of the expanded target subsequence, wherein the third observation point is adjacent to the first observation point or the second observation point, and a starting point distance of each of the third observation points is smaller than a starting point distance of the first observation point and larger than a starting point distance of the second observation point;
[0029] Obtaining the difference between the flat water level of the target section at the historical observation times and the riverbed elevation of each of the third observation points;
[0030] If there is a fourth observation point whose difference is less than the preset difference among the plurality of third observation points, the first observation point and / or the second observation point are re-determined with the goal of eliminating the fourth observation point, so as to determine the coordinate range of the flat river channel;
[0031] Based on the flat water level of the target section at multiple historical observation times and the flat river channel coordinate range of the target section at multiple historical observation times, a cross-sectional view of the target section is generated.
[0032] In some embodiments, after determining the flat beach water level of the target section at the historical observation times based on the beach lip position of the target section at the historical observation times, the method further includes:
[0033] Determining the flat river channel morphological parameters of the target section at the time of the historical observation times based on the flat river channel coordinate range of the target section at the time of the historical observation times;
[0034] Obtaining a preset morphological parameter calculation formula, wherein the morphological parameter calculation formula includes the flat river channel morphological parameters at adjacent moments and a preset first undetermined parameter;
[0035] Using the flat river channel morphological parameters of the target section at a plurality of the historical observation times and the morphological parameter calculation formula, the value of the first undetermined parameter is determined to update the morphological parameter calculation formula;
[0036] The updated morphological parameter calculation formula is used to determine the predicted flat river channel morphological parameters of the target section at a future time.
[0037] In some embodiments, the morphological parameter calculation formula also includes a morphological parameter adjustment strength attenuation coefficient, and the morphological parameter adjustment strength attenuation coefficient represents the influence of the river water and sediment factors at the historical observation times on the flat beach river channel morphological parameter adjustment process. The flat beach river channel morphological parameters of the target section at multiple historical observation times and the morphological parameter calculation formula are used to determine the value of the first undetermined parameter, including:
[0038] Obtaining a preset morphological parameter adjustment intensity attenuation coefficient calculation formula, wherein the morphological parameter adjustment intensity attenuation coefficient calculation formula includes the morphological parameter adjustment intensity attenuation coefficient, the river flow data and the river sediment content data of the target section at a single historical observation number, and a preset second parameter to be determined;
[0039] The values of the first and second parameters to be determined are determined by using the river flow data and river sediment content data of the target section at multiple historical observation times, the calculation formula for the intensity attenuation coefficient of the morphological parameter adjustment, the flat river channel morphological parameters of the target section at multiple historical observation times, and the morphological parameter calculation formula.
[0040] On the other hand, an embodiment of the present application provides an automated data processing device for a flat river channel, comprising:
[0041] A first acquisition module is used to acquire historical observation data of a target section in a target river section, wherein the historical observation data includes starting distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times;
[0042] A first determination module is used to determine the lowest point elevation among the riverbed elevations of multiple observation points at the historical observation times for each of the historical observation times;
[0043] The second determination module is used to determine the beach lip position of the target section at the time of the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the time of the historical observation times;
[0044] A third determination module is used to determine the flat beach water level of the target section at the historical observation times based on the riverbed elevation of the beach lip position of the target section at the historical observation times;
[0045] The fourth determination module is used to determine the actual flat water level of the target section according to the flat water level of the target section at multiple historical observation times.
[0046] On the other hand, the present application also provides a computer device, the computer device comprising:
[0047] one or more processors;
[0048] Memory; and
[0049] 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 automated data processing for flat river channels.
[0050] 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 of any one of the methods for automated data processing for flat river channels.
[0051] 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 automated data processing method for flat river channels as described in any of the above items.
[0052] The automated data processing method, device and medium for flat beach river channels provided in the embodiments of the present application determine the lowest point elevation by acquiring the starting point distance and riverbed elevation of multiple observation points of the target section at multiple historical observation times, and determine the beach lip position of the target section at the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevation of multiple observation points at the historical observation times, and then determine the flat beach water level based on the riverbed elevation at the beach lip position, and determine the actual flat beach water level by combining the flat beach water levels at multiple historical observation times, so that the identified actual flat beach water level can be more accurate and the identification efficiency is higher than that of manual judgment. BRIEF DESCRIPTION OF THE DRAWINGS
[0053] 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.
[0054] Figure 1 It is a schematic diagram of an embodiment of an automated data processing method for a flat river channel provided in an embodiment of the present application;
[0055] Figure 2 It is a schematic diagram of an embodiment of a typical alluvial river cross-section beach channel morphology provided in the embodiments of the present application;
[0056] Figure 3a It is a schematic diagram of an embodiment of the flat beach area of the target section at different historical observation times in the embodiment of the present application;
[0057] Figure 3b It is a schematic diagram of an embodiment of the Pingtan River width of the target section at different historical observation times in the embodiment of the present application;
[0058] Figure 3c It is a schematic diagram of an embodiment of the flat water depth of the target section at different historical observation times in the embodiment of the present application;
[0059] Figure 3d It is a schematic diagram of an embodiment of the river phase coefficient of the flat river channel of the target section at different historical observation times in the embodiment of the present application;
[0060] Figure 4 It is a schematic diagram of an embodiment of the prediction result of the river phase coefficient of the flat river channel in the flat river channel morphological parameters of the embodiment of the present application;
[0061] Figure 5 It is a schematic diagram of the structure of an embodiment of an automated data processing device for a flat river channel provided in an embodiment of the present application;
[0062] Figure 6 It is a schematic diagram of the terminal structure of an embodiment of the computer device provided in the embodiments of the present application. DETAILED DESCRIPTION
[0063] 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.
[0064] 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.
[0065] 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.
[0066] 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.
[0067] In the related technologies, the efficient simulation of the evolution law of the morphological parameters of the flat beach channel involves two problems: one is the accurate identification of the flat beach water level, and the other is the construction of the calculation method of the long-term series evolution trend of the flat beach channel.
[0068] Regarding the first question, there are two main methods for determining the flat water level in current research. The first is to judge the flat water level of the river section based on human experience; the second is to achieve this by identifying the mutation points of the relationship curve between characteristic hydraulic parameters such as river channel storage, water flow sediment entrainment, hydraulic radius and water level. Taking the Yellow River, a typical alluvial river in my country, as an example, its river section is wide and shallow, its branches are intertwined, and its artificially built grid dikes are intertwined. It has complex and changeable characteristics, and it is very difficult to identify the flat water level. Although the first method mentioned above can give full play to the accumulated practical experience of practitioners, it is highly subjective and has low work efficiency, and it is difficult to meet the needs of actual production and life (actual work often involves dozens or even hundreds of sections); in the second method, the selected single characteristic hydraulic parameter is not representative enough to fully reflect the flow and stability characteristics of the flat river channel.
[0069] Regarding the second question, for the simulation of the long-term evolution trend of flat river channels, in addition to the mathematical models based on the water-sediment continuity equation and the motion equation (although it can better simulate the water flow, sediment movement and riverbed scouring and siltation process, it is often limited to small areas and short durations due to the difficulties in closing the control equations and the large amount of calculation), the commonly used description methods currently include nonlinear equation method, hyperbolic equation method, sliding average and delayed response model method. Although fast calculation is achieved, there are still defects such as considering a single influencing factor and focusing more on the evolution results while ignoring the process.
[0070] In response to the above problems, the embodiments of the present application provide an automated data processing method, device and medium for flat beach channels, which is intended to be based on the physical meaning of flat beach channels, and on the basis of summarizing and refining the beach channel morphological characteristics, compare the similarities and differences between the cross-sectional morphologies of different measurements, and use multiple cross-sectional morphologies for comprehensive judgment to achieve rapid and accurate identification of the beach lip position and flat beach water level. In addition, based on this, considering the influence of factors such as incoming water and incoming sand on the adjustment potential and adjustment rate of flat beach channels, a more complete calculation formula for the evolution law of the cross-sectional morphology of flat beach channels is proposed, so as to achieve efficient prediction of the long-term series evolution trend of flat beach channels under changing water and sand conditions (that is, to determine the predicted flat beach channel morphological parameters of the target section at a future moment). The automated data processing method, device and medium for flat beach channels provided in the embodiments of the present application can be widely used for rapid and accurate judgment of riverbed evolution trends in multiple fields such as riverbed evolution prediction, river channel management and protection, which are described in detail below.
[0071] In one embodiment, the automated data processing method for a flat river channel includes:
[0072] 101. Obtain historical observation data of a target section in a target river section, wherein the historical observation data includes starting point distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times.
[0073] In the embodiment of the present application, the target river section is a designated river section, the target section is a designated section in the target river section, and the historical observation data is the observation data actually measured for the target section.
[0074] In an embodiment of the present application, the historical observation data is standardized data. In the historical observation data, the format of the starting point distance and the riverbed elevation can be a three-dimensional array, for example, it can be: X(i, j, k), Y(i, j, k), where X is the measured distance from the observation point to the starting point of the section, that is, the starting point distance, Y is the measured riverbed elevation of the observation point, i is the serial number of the target section, j is the serial number of the historical observation times, and k is the serial number of the observation point. The starting point distance and riverbed elevation of multiple observation points of the same section at multiple historical observation times are stored in the same file in the format of the above three-dimensional array, and the data of different sections are named after the section name and stored in different files for reading.
[0075] 102. For each historical observation number, the lowest point elevation is determined from the riverbed elevations of multiple observation points during the historical observation number.
[0076] In the embodiment of the present application, the lowest point elevation is the elevation of the deep point, and the deep point refers to the deepest observation point in the target section. For each historical observation number, among the riverbed elevations of multiple observation points at the historical observation number, the riverbed elevation with the smallest value is determined, that is, the lowest point elevation.
[0077] 103. Using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times, the beach lip position of the target section at the historical observation times is determined.
[0078] In the embodiment of the present application, for each historical observation number, the corresponding beach lip position is determined using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times. Among them, the observation point where the lowest point elevation is located is the deep point, and the beach lip position refers to the edge of the river beach adjacent to the water flow.
[0079] In some embodiments of the present application, step 103 may include: taking the observation point where the lowest point elevation is located as the deep throbbing point of the target section at the historical observation times; dividing the observation point sequence according to the assumed calculated water level and the riverbed elevations of multiple observation points at the historical observation times to obtain multiple subsequences, the observation point sequence including multiple observation points sorted according to the starting point distance at the historical observation times, for example, the riverbed elevations of multiple observation points in the observation point sequence may be compared with the assumed calculated water level, thereby screening out multiple subsequences, the riverbed elevation of each observation point in each subsequence being less than the assumed calculated water level; determining the target subsequence where the deep throbbing point is located among the multiple subsequences, the target subsequence where the deep throbbing point is located is the subsequence that is more likely to contain the beach lip position; determining the beach lip position of the target section at the historical observation times based on the riverbed elevations of multiple observation points in the target subsequence.
[0080] The observation point sequence can be obtained through the previous data sorting. For example, when sorting the data, the multiple observation points at the historical observation times can be sorted from small to large according to the starting distances of the multiple observation points at the historical observation times to obtain the observation point sequence. It can be seen that the order of the multiple observation points in the observation point sequence is the same as their order along the horizontal direction in the target section.
[0081] In a further embodiment, the observation point sequence is divided according to the assumed calculated water level ans_Z and the riverbed elevations of multiple observation points at the historical observation times to obtain multiple subsequences, which may include: determining the initial value H_mean of the average water depth of the target section at the historical observation times based on the riverbed elevations of multiple observation points of the target section at the historical observation times; taking multiple values at a specified step size (e.g., 0.05 meters) within the target water level range to obtain multiple assumed calculated water levels ans_Z, wherein the target water level range may be 0.8 times the initial value of the average water depth of the section above the deep fen point ~ the maximum value of the riverbed elevation, for example, [the elevation of the deep fen point + 0.8 H_mean, Y_max]; using each assumed calculated water level ans_Z, the observation point sequence is divided respectively to obtain multiple subsequences under each calculated water level ans_Z. Accordingly, determining the target subsequence where the deep point is located in the multiple subsequences may include: determining the target subsequence where the deep point is located under each calculated water level ans_Z in the multiple subsequences under each calculated water level.
[0082] In a further embodiment, due to the rule that the left edge point of the beach lip is higher than its adjacent right observation point, and the right edge point of the beach lip is higher than its adjacent left observation point, the target subsequence where the deep thalassemia point is located at each calculated water level can be expanded respectively. Specifically, based on the riverbed elevations of multiple observation points in the target subsequence, determining the beach lip position of the target section at the time of historical observation times can include: for each target subsequence where the deep thalassemia point is located at the calculated water level, determining the maximum starting point distance and the minimum starting point distance among the starting point distances of multiple observation points in the target subsequence; detecting whether there is a target observation point that meets the first preset condition or the second preset condition; if there is a target observation point, adding the target observation point to the target subsequence until there is no target observation point that meets the first preset condition or the second preset condition, and obtaining an expanded target subsequence; based on the multiple observation points in the expanded target subsequence, determining the beach lip position of the target section at the time of historical observation times.
[0083] Among them, the first preset condition includes that the target observation point does not exist in the target subsequence, the target observation point is adjacent to the observation point with the maximum starting point distance in the observation point sequence, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point with the maximum starting point distance satisfy the first preset size relationship. The first preset size relationship can, for example, be that the riverbed elevation of the target observation point is greater than the riverbed elevation of the observation point with the maximum starting point distance, that is, it does not conform to the rule that the right edge point of the beach lip is higher than its adjacent left observation point. Therefore, it means that the observation point with the maximum starting point distance is not the actual right edge point of the beach lip, and therefore the target observation point can be added to the target subsequence. Similarly, the second preset condition includes that the target observation point does not exist in the target subsequence, the target observation point is adjacent to the observation point where the minimum starting point distance is located in the observation point sequence, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the minimum starting point distance is located satisfy the second preset size relationship. The second preset size relationship may be, for example, that the riverbed elevation of the target observation point is greater than the riverbed elevation of the observation point where the minimum starting point distance is located, that is, it does not conform to the rule that the left edge point of the beach lip is higher than its adjacent right observation point, thus indicating that the observation point where the minimum starting point distance is located is not the left edge point of the actual beach lip, and therefore the target observation point can be added to the target subsequence. After the target observation point is added to the target subsequence, the starting point distances of multiple observation points in the target subsequence still satisfy the order of arrangement from small to large.
[0084] In a further embodiment, determining the beach lip position of the target section at the historical observation times based on multiple observation points in the expanded target subsequence may include: for each expanded target subsequence under the calculated water level, determining the first observation point node_right with the largest starting point distance and the second observation point node_left with the smallest starting point distance among the multiple observation points in the expanded target subsequence; determining the beach lip position of the target section at the historical observation times, wherein the beach lip position includes the starting point distance and riverbed elevation of the first observation point and the starting point distance and riverbed elevation of the second observation point in at least one expanded target subsequence.
[0085] In a further embodiment, determining the beach lip position of the target section at the historical observation times may include: for each expanded target subsequence under the calculated water level, obtaining the minimum value among the calculated water level, the riverbed elevation Y (node_right) of the first observation point, and the riverbed elevation Y (node_left) of the second observation point as the updated calculated water level; obtaining the cross-section river width and cross-section water depth under the updated calculated water level; for each expanded target subsequence under the calculated water level, obtaining the riverbed elevation Y (node_right) of the first observation point, The minimum value of the riverbed elevation Y (node_left) of the second observation point, and then obtain the absolute value of the difference between the minimum value and the updated calculated water level; among multiple calculated water levels, determine the calculated water level when the absolute value is the minimum; if the cross-sectional river width under the calculated water level when the absolute value is the minimum is within the preset cross-sectional river width range, and the cross-sectional water depth under the calculated water level when the absolute value is the minimum is within the preset cross-sectional water depth range, it is determined that the beach lip position corresponding to the calculated water level when the absolute value is the minimum is more likely to be the beach lip position of the target section at the historical observation times. At this time, the beach lip position of the target section at the historical observation times may include: the starting point distance and riverbed elevation of the first observation point in the expanded target subsequence under the calculated water level when the absolute value is the minimum, and the starting point distance and riverbed elevation of the second observation point, so that the determined beach lip position is more accurate.
[0086] Among them, obtaining the cross-sectional river width and cross-sectional water depth under the updated calculated water level may include: calculating the difference between the starting point distance of the first observation point and the starting point distance of the second observation point for the expanded target subsequence under the updated calculated water level, and using it as the cross-sectional river width; calculating the cross-sectional area by the trapezoidal area summation method between adjacent observation points based on the starting point distances and riverbed elevations of adjacent observation points; and using the ratio of the cross-sectional area to the cross-sectional river width as the cross-sectional water depth.
[0087] 104. Based on the riverbed elevation of the beach lip position of the target section at the historical observation times, determine the flat beach water level of the target section at the historical observation times.
[0088] In an embodiment of the present application, at the shoal lip position of the target section at the historical observation times, the smaller value of the riverbed elevation at the first observation point and the riverbed elevation at the second observation point can be directly used as the flat beach water level of the target section at the historical observation times.
[0089] 105. Determine the actual flat water level of the target section based on the flat water level of the target section at multiple historical observation times.
[0090] In some embodiments of the present application, after obtaining the flat water level of the target section at multiple historical observation times, the flat water level of the target section at multiple historical observation times can also be tested and optimized based on the rule that in the year when the flood occurred, the siltation of the beach will cause the flat water level to rise to a large extent, while in the year when the flood did not occur, the flat water level will not rise to a large extent. Specifically, it can include: determining whether the target section has a flood in the time period of each historical observation number; for each historical observation number, if the target section has a flood in the time period of the historical observation number, it indicates that the flat water level of the target section at the historical observation number should be increased relative to the flat water level of the target section at the previous historical observation number. If an increase is detected, it is determined that the flat water level of the target section at the historical observation number is correct; if no increase is detected, it is determined that the flat water level of the target section at the historical observation number is incorrect, so the flat water level of the target section at the historical observation number can be adjusted to make the flat water level of the target section at the historical observation number conform to the rule of increase. In addition, after adjusting the flat beach water level, the flat beach river channel coordinate range of the target section at the time of the historical observation times can also be updated.
[0091] In some embodiments of the present application, step 105 may also include: generating a cross-sectional diagram of the target section based on the flat beach water level of the target section at multiple historical observation times and the beach lip position of the target section at multiple historical observation times. For example, the flat beach water level of the target section at multiple historical observation times, the beach lip position of the target section at multiple historical observation times, and the observation points between the beach lip positions may all be drawn into the same preset cross-sectional diagram to obtain a cross-sectional diagram of the target section; in the cross-sectional diagram, the actual flat beach water level of the target section is determined. For example, the cross-sectional diagram may be displayed, and relevant technical personnel may determine a more accurate actual flat beach water level of the target section by viewing the data drawn in the cross-sectional diagram. Compared with direct manual judgment of the river section, this method is simpler, more efficient and more accurate.
[0092] In some embodiments of the present application, taking into account the complexity of the cross-sectional morphology of alluvial rivers, there may be areas with smaller water depths on the left and right sides of the beach lip position determined in step 103, and the flat beach river channel is the main flood discharge channel of the river channel. The smaller water depth has little effect on the flood discharge capacity, but has a greater impact on the width of the flat beach river. Therefore, based on the beach lip position determined in step 103, the smaller water depths on the left and right sides of the flat beach river channel can be eliminated, thereby optimizing the flat beach river channel.
[0093] Specifically, based on the flat beach water level of the target section at multiple historical observation times and the beach lip position of the target section at multiple historical observation times, generating a cross-sectional diagram of the target section can include: for each historical observation time, among the multiple observation points of the expanded target subsequence (under the calculated water level at the above minimum value), determining multiple third observation points, wherein the third observation point is adjacent to the first observation point or the second observation point, and the starting point distance of each third observation point is smaller than the starting point distance of the first observation point and larger than the starting point distance of the second observation point, that is, the third observation point is located between the first observation point and the second observation point, and is adjacent to the first observation point or the second observation point; obtaining the difference between the flat beach water level of the target section at the historical observation times and the riverbed elevation of each third observation point; if If there is a fourth observation point whose difference is less than the preset difference (e.g., 1 meter) among the multiple third observation points, it indicates that the fourth observation point is a part with a smaller water depth. Therefore, with the goal of eliminating the current fourth observation point, the first observation point and / or the second observation point are re-determined (e.g., if the fourth observation point is located on the right side of the deep point and adjacent to the first observation point, the first observation point can be moved to the left by at least one observation point, and / or if the fourth observation point is located on the left side of the deep point and adjacent to the second observation point, the second observation point is moved to the right by at least one observation point), so as to determine the corresponding flat beach channel coordinate range (the flat beach channel coordinate range includes at least the re-determined first observation point and / or second observation point), and the flat beach channel coordinate range is updated to obtain an optimized and more accurate flat beach channel. Then, based on the flat beach water level of the target section at multiple historical observation times, the flat beach channel coordinate range of the target section at multiple historical observation times, and the observation points within the flat beach channel coordinate range, a cross-sectional view of the target section is generated.
[0094] In the technical solution disclosed in this embodiment, the lowest point elevation is determined by obtaining the starting point distance and riverbed elevation of multiple observation points of the target section at multiple historical observation times, and the beach lip position of the target section at the historical observation times is determined by using the observation point where the lowest point elevation is located and the riverbed elevation of multiple observation points at the historical observation times. The flat beach water level is then determined based on the riverbed elevation at the beach lip position, and the actual flat beach water level is determined by combining the flat beach water levels at multiple historical observation times, so that the identified actual flat beach water level can be more accurate and the identification efficiency is higher than that of manual judgment.
[0095] In some embodiments of the present application, after step 104, it may also include: determining the flat beach channel morphological parameters of the target section at the historical observation times based on the flat beach channel coordinate range of the target section at the historical observation times, the flat beach channel morphological parameters may include, for example, flat beach area, flat beach river width, flat beach water depth, flat beach river channel river phase coefficient, etc.; obtaining a preset morphological parameter calculation formula, the morphological parameter calculation formula includes the flat beach channel morphological parameters at adjacent moments and a preset first pending parameter; using the flat beach channel morphological parameters and morphological parameter calculation formula of the target section at multiple historical observation times, determining the value of the first pending parameter to update the morphological parameter calculation formula; using the updated morphological parameter calculation formula to determine the predicted flat beach channel morphological parameters of the target section at future times, thereby achieving efficient prediction of the long-term series evolution trend of the flat beach channel.
[0096] Among them, based on the flat beach channel coordinate range of the target section at the historical observation times, determining the flat beach channel morphological parameters of the target section at the historical observation times may include: obtaining multiple fifth observation points within the flat beach channel coordinate range from multiple observation points at the historical observation times, and determining the flat beach area A based on the trapezoidal area summation method between adjacent fifth observation points; determining the flat beach river width B based on the difference between the starting point distance of the (redetermined) first observation point within the flat beach channel coordinate range and the starting point distance of the (redetermined) second observation point within the flat beach channel coordinate range; determining the flat beach water depth H based on the ratio of the flat beach area to the flat beach river width, for example, H=A / B; determining the flat beach channel river phase coefficient ξ based on the flat beach river width and the flat beach water depth, for example, ξ=B 0.5 / H. The river phase coefficient ξ of the flat river channel can directly reflect the narrowness and depth of the flat river channel cross-section, and is closely related to the flow capacity of the target section. It is selected here as the key variable for the evolution of the flat river channel morphology.
[0097] In a further embodiment, the morphological parameter calculation formula also includes a morphological parameter adjustment intensity attenuation coefficient, which characterizes the influence of the river water and sediment factors at the historical observation times on the adjustment process of the flat river channel morphological parameters. Therefore, using the flat river channel morphological parameters and the morphological parameter calculation formula of the target section at multiple historical observation times to determine the value of the first parameter to be determined can include: obtaining a preset morphological parameter adjustment intensity attenuation coefficient calculation formula, which includes the morphological parameter adjustment intensity attenuation coefficient, the river flow data and river sediment content data of the target section at a single historical observation time, and the preset second parameter to be determined; using the river flow data and river sediment content data of the target section at multiple historical observation times, the morphological parameter adjustment intensity attenuation coefficient calculation formula, the flat river channel morphological parameters of the target section at multiple historical observation times, and the morphological parameter calculation formula to determine the value of the first parameter to be determined and the value of the second parameter to be determined.
[0098] The preset morphological parameter calculation formula may be, for example:
[0099] λ(i)=(1-e αdt )λe(i-1)+e αdt λ(i-1), formula (1)
[0100] Among them, λ(i) and λ(i-1) are the flat river channel morphological parameters of the i-th historical observation number and the i-1-th historical observation number sorted in chronological order, that is, the flat river channel morphological parameters at adjacent moments, the preset first undetermined parameter includes α (characterizing the attenuation coefficient of non-equilibrium evolution), and dt is the calculation time step. In addition, the preset morphological parameter calculation formula can also include the morphological parameter adjustment intensity attenuation coefficient.
[0101] The preset morphological parameter adjustment intensity attenuation coefficient calculation formula may be, for example:
[0102] λe(i)=k1S(i) b / Q(i) a , formula (2)
[0103] Among them, Q(i) is the river channel (average) flow data of the target section at the i-th historical observation number, S(i) is the river channel (average) sediment content data of the target section at the i-th historical observation number, k1, a, b are the preset second undetermined parameters, which are used to reflect the influence of basin climate characteristics, vegetation conditions, sediment material composition, etc. on the river channel morphology. It can be seen that formula (2) reflects the corresponding relationship between the morphological parameter adjustment intensity attenuation coefficient and the characteristic value of river water and sediment.
[0104] The non-equilibrium evolution attenuation coefficient α in the preset first undetermined parameter can be calculated by the following formula:
[0105] , formula (3)
[0106] in, is the average value of the river flow data of the target section at multiple historical observation times, is the average value of the river sediment content data of the target section at multiple historical observation times, and k2, c, and d are the preset third unknown parameters.
[0107] , formula (4)
[0108] , Formula (5)
[0109] By using the river flow data and river sediment content data of the target section at multiple historical observation times, the flat river channel morphological parameters of the target section at multiple historical observation times, and the above formulas (1), (2), (3), (4), and (5), the unknown parameters k1, k2, a, b, c, and d are calibrated, and the morphological parameter calculation formula can be updated.
[0110] It can be seen that the prediction method for the long-term series evolution trend of the flat river channel provided in the embodiment of the present application comprehensively considers the influence of incoming water, incoming sand, etc. on the adjustment potential and adjustment rate of the flat river channel, overcomes the defects of insufficient consideration of factors in previous studies and the emphasis on evolution results while ignoring the process, and better realizes the efficient prediction of the long-term series evolution trend of the cross-sectional morphological parameters of the flat river channel under changing water and sand conditions, and serves the rapid judgment of river channel evolution trends in multiple fields such as riverbed evolution prediction, river management and protection.
[0111] Below, the above embodiments are specifically described by taking the large cross section measured in Gaocun, a wide and shallow section of the lower Yellow River, as an example. Figure 1 , which may specifically include the following steps:
[0112] Step 1: Standardization processing of observation data from 1986 to 2015 from Gaocun, a large section of the wide and shallow lower reaches of the Yellow River.
[0113] The long-term flood information of the target section was collected, and the year of flood was determined to be 1996. The historical observation data of the target section measured before and after the flood season were sorted out, and the starting point distance and riverbed elevation of multiple observation points at different historical observation times of the same target section were stored in the same file in the format of starting point distance X and elevation Y, as shown in Table 1 below.
[0114] Table 1 Standardized data set of historical observation times of measured target sections (taking 4 measurements from 1986 to 1987 as an example)
[0115]
[0116] In Table 1, "198620" in the second row represents before the flood season of 1986, and "198660" represents after the flood season of 1986; the first column is the serial number of the observation point; the second and third columns are the starting point distance and elevation before the flood season of 1986, respectively; the fourth and fifth columns are the starting point distance and elevation after the flood season of 1986, respectively; the sixth and seventh columns are the starting point distance and elevation before the flood season of 1987, respectively; the eighth and ninth columns are the starting point distance and elevation after the flood season of 1987, respectively. Figure 2 This is an example of a cross-section diagram of the lower Yellow River at Huayuankou. Figure 2 In the figure, the horizontal axis is the distance from the starting point and the vertical axis is the riverbed elevation.
[0117] Step 2: Identify the beach lip position of the target section.
[0118] 2.1. According to the channel morphology characteristics of the target section, the approximate range of the river width and water depth of the flat beach of the target section, as well as the initial value of the average water depth of the section, are determined based on experience. Then, for each historical observation number, the riverbed elevation and position of the deep point are first found. The following is an example of the post-flood data of Gaocun section in 2015:
[0119] The approximate ranges of Pingtan River width and Pingtan water depth are B=[200,1500] and H=[1.5,10] respectively. The initial value of the average water depth of the section is H_mean=4m, and the riverbed elevation at the deep point is h_thaw=min(Y)=55.99;
[0120] At each different calculated water level, the area where the riverbed elevation is less than the calculated water level is divided to obtain multiple subsequences. Among them, the value range of the calculated water level is [the elevation of the deep well point + 0.8 H_mean, Y_max], and the step length between different calculated water levels is 0.05 m.
[0121] Among the multiple subsequences, determine the target subsequence where the deep point is located. Then, based on the rule that the left edge point of the lip is higher than its adjacent right observation point, and the right edge point of the lip is higher than its adjacent left observation point, traverse multiple observation points outside the target subsequence, and add the target observation points that meet the first preset condition and / or the second preset condition to the target subsequence until there is no target observation point that meets the first preset condition or the second preset condition, and obtain the expanded target subsequence.
[0122] Step 3: Determine the flat water level and calculate the flat river channel morphological parameters (including flat area A, flat river width B, flat water depth H, and flat river channel river phase coefficient ξ).
[0123] For each expanded target subsequence under the calculated water level, the minimum value among the calculated water level, the riverbed elevation of the first observation point, and the riverbed elevation of the second observation point is obtained as the updated calculated water level. Then, for each expanded target subsequence under the calculated water level, the absolute value of the difference between the minimum value among the riverbed elevation of the first observation point and the riverbed elevation of the second observation point and the updated calculated water level is obtained, and the results shown in Table 2 are obtained as follows:
[0124] Table 2
[0125]
[0126] It can be seen from Table 2 that the minimum absolute value is 0.020, and the corresponding calculated water level is 62.74 meters. Since the cross-section river width at the calculated water level of 62.74 is within the preset cross-section river width range, and the cross-section water depth at the calculated water level of 62.74 is within the preset cross-section water depth range, the beach lip position of the target section at the time of the historical observation can be further determined.
[0127] After obtaining the flat beach water level of the target section at multiple historical observation times, the flat beach water level of the target section at multiple historical observation times can be tested and optimized based on the rule that in the years when the flood occurs, the siltation of the beach will cause the flat beach water level to rise to a large extent, while in the years when the flood does not occur, the flat beach water level will not rise to a large extent. After optimizing the flat beach water level of the target section at multiple historical observation times, the actual flat beach water level of the target section can be determined based on the flat beach water level of the target section at multiple historical observation times.
[0128] Taking into account the complexity of the cross-sectional morphology of alluvial rivers, there may be areas with smaller water depths on the left and right sides of the beach lip position determined above. As the flat beach river channel is the main flood discharge channel of the river channel, the smaller water depth has little effect on the flood discharge capacity, but has a greater impact on the width of the flat beach river. Therefore, based on the beach lip position, the areas with smaller water depths on the left and right sides of the flat beach river channel can be eliminated to optimize the flat beach river channel.
[0129] Based on the coordinate range of the flat beach channel, the corresponding flat beach channel morphological parameters can be determined. The flat beach channel morphological parameters include flat beach area, flat beach river width, flat beach water depth, and flat beach channel river phase coefficient. Figure 3a As shown in Figure 1, the flat beach area of the target section changes with the year; Figure 3b As shown in the figure, the relationship between the width of the Pingtan River in the target section and the years changes; Figure 3c As shown in the figure, the relationship between the flat water depth of the target section and the year changes; Figure 3d As shown, the relationship between the flat river channel river phase coefficient of the target section and the years.
[0130] Step 4: Based on the calculation results in step 3, an efficient prediction is made on the long-term evolution trend of the flat river channel.
[0131] Among them, the average value of the river flow data of the target section at multiple historical observation times can be, for example, 752.76m 3 / s, the average value of the river sediment content data of the target section at multiple historical observation times can be, for example, 12.94kg / m 3 .
[0132] Based on the above formulas (1), (2), (3), (4), and (5), with the minimum variance between the calculated value and the measured value as the judgment principle, the least squares method is used to calibrate the unknown parameters k1, k2, a, b, c, and d, and the results are k1=764.30, k2=-60.40, a=0.72, b=0.20, c=-0.29, and d=-0.72.
[0133] By using the parameter values obtained by calibration and combining the above formulas (1), (2), (3), (4), and (5), the predicted flat river channel morphological parameters of the target section at the future time can be determined. Taking the prediction of the flat river channel phase coefficient as an example, the predicted flat river channel phase coefficient of the target section at the future time is Figure 4 The river phase coefficient of the flat river channel in the middle prediction area is shown in Figure 2. Figure 4 It is a color picture to distinguish the measured and calculated values of the beach-channel river phase coefficient.
[0134] In order to better implement the automated data processing method for a flat river channel in the embodiment of the present application, on the basis of the automated data processing method for a flat river channel, an automated data processing device for a flat river channel is also provided in the embodiment of the present application, such as Figure 5 As shown, the automatic data processing device 500 for flat river channels includes:
[0135] The first acquisition module 501 is used to acquire historical observation data of a target section in a target river section, wherein the historical observation data includes starting distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times;
[0136] The first determination module 502 is used to determine the lowest point elevation among the riverbed elevations of multiple observation points at the historical observation times for each historical observation time;
[0137] The second determination module 503 is used to determine the beach lip position of the target section at the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times;
[0138] The third determination module 504 is used to determine the flat beach water level of the target section at the historical observation times based on the riverbed elevation of the beach lip position of the target section at the historical observation times;
[0139] The fourth determination module 505 is used to determine the actual flat water level of the target section according to the flat water level of the target section at multiple historical observation times.
[0140] The present application also provides a computer device that integrates any of the automatic data processing devices for flat river channels provided in the present application. Figure 6 As shown, it shows a schematic diagram of the structure of the computer device involved in the embodiment of the present application, specifically:
[0141] The computer device may include components such as a processor 601 with one or more processing cores, a memory 602 with one or more computer-readable storage media, a power supply 603, and an input unit 604. Those skilled in the art will appreciate that Figure 6 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.
[0142] The processor 601 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 602 and calling data stored in the memory 602, it executes various functions of the computer device and processes data, thereby monitoring the computer device as a whole. Optionally, the processor 601 may include one or more processing cores; preferably, the processor 601 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 601.
[0143] The memory 602 can be used to store software programs and modules. The processor 601 executes various functional applications and data processing by running the software programs and modules stored in the memory 602. The memory 602 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 602 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 602 may also include a memory controller to provide the processor 601 with access to the memory 602.
[0144] The computer device also includes a power supply 603 for supplying power to various components. Preferably, the power supply 603 can be logically connected to the processor 601 through a power management system, so as to manage charging, discharging, and power consumption through the power management system. The power supply 603 can also include any components such as one or more DC or AC power supplies, recharging systems, power failure detection circuits, power converters or inverters, and power status indicators.
[0145] The computer device may further include an input unit 604, 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.
[0146] 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 601 in the computer device will load the executable files corresponding to the processes of one or more application programs into the memory 602 according to the following instructions, and the processor 601 will run the application programs stored in the memory 602, thereby realizing various functions, as follows:
[0147] Obtain historical observation data of a target section in a target river section, wherein the historical observation data includes starting point distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times; for each historical observation time, determine the lowest point elevation among the riverbed elevations of multiple observation points at the historical observation times; determine the beach lip position of the target section at the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times; determine the flat beach water level of the target section at the historical observation times based on the riverbed elevation of the beach lip position of the target section at the historical observation times; determine the actual flat beach water level of the target section based on the flat beach water levels of the target section at multiple historical observation times.
[0148] 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.
[0149] 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 automated data processing methods for flat river channels provided in the embodiment of the present application. For example, the computer program can be loaded by a processor to execute the following steps:
[0150] Obtain historical observation data of a target section in a target river section, wherein the historical observation data includes starting point distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times; for each historical observation time, determine the lowest point elevation among the riverbed elevations of multiple observation points at the historical observation times; determine the beach lip position of the target section at the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times; determine the flat beach water level of the target section at the historical observation times based on the riverbed elevation of the beach lip position of the target section at the historical observation times; determine the actual flat beach water level of the target section based on the flat beach water levels of the target section at multiple historical observation times.
[0151] 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 automatic data processing method for a flat river channel as described in any one of the above items, for example:
[0152] Obtain historical observation data of a target section in a target river section, wherein the historical observation data includes starting point distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times; for each historical observation time, determine the lowest point elevation among the riverbed elevations of multiple observation points at the historical observation times; determine the beach lip position of the target section at the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times; determine the flat beach water level of the target section at the historical observation times based on the riverbed elevation of the beach lip position of the target section at the historical observation times; determine the actual flat beach water level of the target section based on the flat beach water levels of the target section at multiple historical observation times.
[0153] 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.
[0154] 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.
[0155] The specific implementation of the above operations can be found in the previous embodiments, which will not be described in detail here.
[0156] The above is a detailed introduction to an automated data processing method, device and medium for flat river channels provided in the embodiments 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 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. An automated data processing method for flat river channels, characterized in that: The method comprises: Acquire historical observation data of a target section in a target river section, wherein the historical observation data includes starting point distances and riverbed elevations of multiple observation points of the target section at multiple historical observation times; For each of the historical observation times, determining the lowest point elevation among the riverbed elevations of the plurality of observation points at the historical observation times; Determine the beach lip position of the target section at the time of the historical observation times by using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the time of the historical observation times; Based on the riverbed elevation of the beach lip position of the target section at the historical observation times, determine the flat beach water level of the target section at the historical observation times, wherein the beach lip position of the target section at the historical observation times includes the starting distance and riverbed elevation of the first observation point, the starting distance and riverbed elevation of the second observation point, the starting distance of the first observation point is greater than the starting distance of the second observation point, and the smaller value of the riverbed elevation of the first observation point and the riverbed elevation of the second observation point is used as the flat beach water level of the target section at the historical observation times; Determining the actual flat water level of the target section according to the flat water level of the target section at a plurality of the historical observation times; Wherein, determining the actual flat beach water level of the target section according to the flat beach water level of the target section at multiple historical observation times comprises: generating a cross-sectional view of the target section based on the flat beach water level of the target section at multiple historical observation times and the beach lip position of the target section at multiple historical observation times; determining the actual flat beach water level of the target section in the cross-sectional view; The method of using the observation point where the lowest point elevation is located and the riverbed elevations of multiple observation points at the historical observation times to determine the beach lip position of the target section at the historical observation times includes: taking the observation point where the lowest point elevation is located as the deep throat point of the target section at the historical observation times; dividing the observation point sequence according to the assumed calculated water level and the riverbed elevations of multiple observation points at the historical observation times to obtain multiple subsequences, the observation point sequence including multiple observation points sorted according to the starting point distance at the historical observation times; determining the target subsequence where the deep throat point is located among the multiple subsequences; and determining the beach lip position of the target section at the historical observation times based on the riverbed elevations of multiple observation points in the target subsequence.
2. The automated data processing method for flat river channels according to claim 1, characterized in that: The determining, based on the riverbed elevations of the plurality of observation points in the target subsequence, the beach lip position of the target section at the time of the historical observation times comprises: Determine a maximum starting point distance and a minimum starting point distance among the starting point distances of a plurality of observation points in the target subsequence; Detecting whether there is a target observation point that satisfies a first preset condition or a second preset condition, wherein the first preset condition includes that the target observation point does not exist in the target subsequence, the target observation point is adjacent to the observation point where the maximum starting point distance is located in the observation point sequence, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the maximum starting point distance is located satisfy a first preset size relationship, and the second preset condition includes that the target observation point does not exist in the target subsequence, the target observation point is adjacent to the observation point where the minimum starting point distance is located in the observation point sequence, and the riverbed elevation of the target observation point and the riverbed elevation of the observation point where the minimum starting point distance is located satisfy a second preset size relationship; If the target observation point exists, adding the target observation point to the target subsequence to obtain an expanded target subsequence; Based on the multiple observation points in the expanded target subsequence, the beach lip position of the target section at the historical observation times is determined.
3. The automated data processing method for flat river channels according to claim 2, characterized in that: The determining, based on the multiple observation points in the expanded target subsequence, the beach lip position of the target section at the historical observation times comprises: Determine, among the multiple observation points of the expanded target subsequence, a first observation point with the largest starting point distance and a second observation point with the smallest starting point distance; Determine the beach lip position of the target section at the historical observation times, wherein the beach lip position includes the starting point distance and riverbed elevation of the first observation point, and the starting point distance and riverbed elevation of the second observation point.
4. The automated data processing method for flat river channels according to claim 3, characterized in that: The generating of the cross-sectional view of the target section based on the flat beach water level of the target section at multiple historical observation times and the beach lip position of the target section at multiple historical observation times comprises: For each of the historical observation times, determining a plurality of third observation points among the plurality of observation points of the expanded target subsequence, wherein the third observation point is adjacent to the first observation point or the second observation point, and a starting point distance of each of the third observation points is smaller than a starting point distance of the first observation point and larger than a starting point distance of the second observation point; Obtaining the difference between the flat water level of the target section at the historical observation times and the riverbed elevation of each of the third observation points; If there is a fourth observation point whose difference is less than the preset difference among the plurality of third observation points, the first observation point and / or the second observation point are re-determined with the goal of eliminating the fourth observation point, so as to determine the coordinate range of the flat river channel; Based on the flat water level of the target section at multiple historical observation times and the flat river channel coordinate range of the target section at multiple historical observation times, a cross-sectional view of the target section is generated.
5. The automated data processing method for flat river channels according to any one of claims 1 to 4, characterized in that: After determining the flat beach water level of the target section at the historical observation times based on the beach lip position of the target section at the historical observation times, the method further includes: Determine the flat river channel morphological parameters of the target section at the time of the historical observation times based on the flat river channel coordinate range of the target section at the time of the historical observation times; Obtaining a preset morphological parameter calculation formula, wherein the morphological parameter calculation formula includes the flat river channel morphological parameters at adjacent moments and a preset first undetermined parameter; Using the flat river channel morphological parameters of the target section at a plurality of the historical observation times and the morphological parameter calculation formula, the value of the first undetermined parameter is determined to update the morphological parameter calculation formula; The updated morphological parameter calculation formula is used to determine the predicted flat river channel morphological parameters of the target section at a future time.
6. The automated data processing method for flat river channels according to claim 5, characterized in that: The morphological parameter calculation formula also includes a morphological parameter adjustment strength attenuation coefficient, which represents the influence of the river water and sediment factors at the historical observation times on the flat beach river channel morphological parameter adjustment process. The flat beach river channel morphological parameters of the target section at multiple historical observation times and the morphological parameter calculation formula are used to determine the value of the first undetermined parameter, including: Obtaining a preset morphological parameter adjustment intensity attenuation coefficient calculation formula, wherein the morphological parameter adjustment intensity attenuation coefficient calculation formula includes the morphological parameter adjustment intensity attenuation coefficient, the river flow data and the river sediment content data of the target section at a single historical observation number, and a preset second parameter to be determined; The values of the first and second parameters to be determined are determined by using the river flow data and river sediment content data of the target section at multiple historical observation times, the calculation formula for the intensity attenuation coefficient of the morphological parameter adjustment, the flat river channel morphological parameters of the target section at multiple historical observation times, and the morphological parameter calculation formula.
7. A computer device, characterized in that: The computer device comprises: one or more processors; Memory; and one or more applications, wherein the one or more applications are stored in the memory and configured to be executed by the processor to implement the steps in the automated data processing method for flat river channels as described in any one of claims 1 to 6.
8. A computer-readable storage medium, characterized in that: A computer program is stored thereon, and the computer program is loaded by a processor to execute the steps in the automated data processing method for flat river channels as claimed in any one of claims 1 to 6.
Citation Information
Patent Citations
Rapid calculation method for river channel flat beach water level under strong scouring condition
CN116356752A