Flood wave continuous self-correction extraction method, system and medium

Abstract

The application relates to a flood wave continuous self-correction extraction method, a system and a medium, which comprises the following steps: forming a time sequence of flow / level fluctuation slope; respectively statistically analyzing period maximum value series of two stations on the section; judging whether m continuous flood peaks are independent; for the independent flood process, taking the flood peak position as a starting point, respectively searching forward and backward to divide the complete flood process; preliminarily determining the flood propagation time of the stations on the section; calculating the flood propagation characteristics information of the flood propagation time of the stations on the section, the rise difference coefficient and the like, and correcting the flood propagation time of the stations on the section according to the calculation result. The application can effectively realize the analysis of the flood propagation characteristics between the stations on the basin surface, and improve the understanding of the flood propagation law of the basin.

Description

Technical Field

[0001] This application relates to the field of autonomous extraction and analysis of flood waves, and more specifically, to a method, system and medium for continuous self-calibration extraction of flood waves. Background Technology

[0002] Climate change and severe human activities have significantly altered the flood evolution patterns in the Yangtze River Basin. Therefore, research on flood evolution patterns under these new conditions is crucial for improving our understanding of flood development and evolution in the Yangtze River Basin. However, traditional flood evolution analysis primarily employs graphical methods. This involves first plotting water level / discharge processes at upstream and downstream stations, then selecting relevant flood processes at these stations based on visual or empirical judgments, and finally analyzing flood evolution characteristic parameters. While accurate, this method is time-consuming and unsuitable for large-scale, inter-station flood process studies. Summary of the Invention

[0003] The purpose of this application is to provide a method, system and medium for continuous self-correction extraction of flood waves, so as to realize the automatic extraction and analysis of flood processes in hydrological sequences.

[0004] To achieve the above objectives, this application provides the following technical solution:

[0005] In a first aspect, embodiments of this application provide a flood wave continuity self-correction extraction method, comprising the following specific steps:

[0006] S1. Select any river section, identify the upstream and downstream hydrological / water level stations, use the station flood reporting data, arbitrarily select the flow / water level process data within a certain period of time, plot the flow / water level development process graph with time as the horizontal axis and flow / water level as the vertical axis, calculate the flow / water level slope between two adjacent time nodes, form a time series of flow / water level rise and fall slopes, select the time period maximum / minimum values ​​of the flow / water level series of upstream and downstream stations, and proceed to S2 for the next step of calculation;

[0007] S2. Statistically analyze the maximum value series of the two stations upstream and downstream of the cross section, and set the flood peak threshold for the two stations in three levels: small and medium flood, large flood, and flood exceeding the standard. Extract the sample that exceeds the flood peak threshold of a certain level from the maximum value of the flow / water level series as the flood peak flow, and transfer to S4 for the next step of calculation.

[0008] S3. If there are m consecutive flood peak flows / water levels, that is, there are no other non-flood peak maximum points between flood peaks, determine whether these m consecutive flood peaks are independent. Treat the non-independent consecutive flood peak flows / water levels as several flood peaks of an independent flood process, and proceed to S5 for the next step of calculation.

[0009] S4. For an independent flood process, starting from the flood peak location, search for the start and end times of the flood process forward and backward respectively to segment out the complete flood process, and then proceed to S6 for the next step of calculation.

[0010] S5. Preliminarily determine the flood propagation time of upstream and downstream stations on the cross section, and conduct similarity sliding analysis on the flood process of upstream and downstream stations on the cross section. If the flood events of the two stations on the cross section are similar, then the flood process of the downstream station is transmitted from the upstream flood.

[0011] S6. Calculate the flood propagation characteristics information such as the flood propagation time and the flood difference coefficient at upstream and downstream stations on the cross section, and correct the flood propagation time at upstream and downstream stations on the cross section in S5 based on the calculation results.

[0012] The implementation method of S1 is as follows: Select the flow / water level process of the upstream station of the cross section. / Flow / water level process at downstream stations of the cross section / Calculate the slope of flow / water level fluctuations between two adjacent time points to form a time series of slopes, including the flow / water level fluctuation slopes at upstream stations of the cross-section. The slope of flow / water level fluctuation at downstream stations of the cross section ,in, or , or The extreme values ​​of a time period can be determined using the following formula:

[0013] .

[0014] The implementation method of S2 is as follows: the maximum value of flow / water level at the upstream station of the cross section. / Maximum flow / water level at downstream stations of the cross-section / Based on the design flood data from upstream and downstream stations, determine the design flood discharge / water level corresponding to minor, major, and above-standard floods. / , / Statistical analysis was conducted to determine the flood type by identifying the points with the maximum flow / water level at upstream and downstream stations within a selected time period.

[0015]

[0016] or

[0017]

[0018] or .

[0019] The implementation of S3 is as follows: Based on the above analysis, all flood peaks within the study period can be selected. If m consecutive maximum points all belong to the flood category, it is necessary to determine whether these m consecutive flood peaks are independent. Assume the consecutive flood peaks are... The corresponding peak time is The time difference between the peaks of two adjacent flood peaks in a series of consecutive flood peaks is ,in If the following formula is satisfied, then m consecutive flood peaks are not independent and can be regarded as the same flood.

[0020]

[0021] In the formula, The base flow for site A; Refers to continuous flood peaks and The interval time; These are two consecutive consecutive flood peak flows; This refers to the minimum flow rate between two flood peaks; , This represents an empirical parameter.

[0022] The implementation of S4 is as follows: For an independent flood process, starting from the flood peak location, the start and end times of the flood process are searched forward and backward respectively to segment the complete flood process. The method for determining the start point of the flood process is as follows: if the flood peak is single-peak and the flood peak discharge is... The corresponding peak flood time is Then, search backward along the time axis for adjacent independent flood peaks, and count all the maximum points within the time interval between the occurrence of two independent flood peaks. ,in Represents the target flood peak Calculate the slope between two adjacent points with maximum flow rates. ,like If all are greater than 0, then The corresponding minimum point to the left of the maximum point The starting point of this flood, and the corresponding starting time is... Search forward along the time axis for adjacent independent flood peaks. ,in Represents the target flood peak Calculate the slope between two adjacent points with maximum flow rates. If the slope between two consecutive points of maximum flow is less than 0, then The corresponding minimum point to the right of the maximum value The receding point of this flood event, and the corresponding receding time is... ,

[0023]

[0024] In the formula, and These are empirical parameters. .

[0025] If the target flood peak is multi-peaked, starting from the first peak, search backwards along the time axis for adjacent independent flood peaks. If the searched independent flood peak is also a single peak, analyze the rising point between the two peaks using the above method. If the searched independent flood peak is multi-peaked, take the last peak among the multi-peaks as the object and analyze the rising point between the two peaks using the above method. If the target flood peak is multi-peaked, starting from the last peak, search forwards along the time axis for adjacent independent flood peaks. If the searched independent flood peak is a single peak, analyze the receding point between the two peaks using the above method. If the searched independent flood peak is multi-peaked, take the first peak among the multi-peaks as the object and analyze the receding point between the two peaks using the above method. Based on the time corresponding to the flood's rising and receding points, segment the complete flood process.

[0026] The same method was used to segment the flood flow / water level process at downstream stations of the cross section.

[0027] The implementation of S5 is as follows: For the flood flow / water level processes segmented at upstream and downstream stations on the cross-section, the peak occurrence time corresponding to the upstream flood process is... Flood peak ,magnanimity Based on the preliminary estimated upstream and downstream flood propagation time Analysis of downstream stations at the cross-section The flood process corresponding to an independent flood peak near a given time, where the flood peak occurs at a specific time. Flood peak ,magnanimity The following formula is used to analyze the similarity of two characteristic indicators, flood peak and flood volume, between two floods at upstream and downstream stations on the cross-section.

[0028]

[0029] In the formula, , Represents empirical factors. ; The closer the similarity is to 1, the higher the similarity of the peak values ​​and magnitudes of the two floods; conversely, the closer the similarity is to 1, the lower the similarity is to 1.

[0030] If the flood sequence lengths at upstream and downstream stations are the same, then the flood process at the upstream station can be represented as follows: The downstream flood process can be represented as The similarity between the two flood processes is calculated using the following formula:

[0031]

[0032] The smaller the similarity, the higher the similarity between the two floods, and vice versa.

[0033] If the flood sequence lengths at upstream and downstream stations are inconsistent, the flood process at the upstream station can be represented as follows: The downstream flood process can be represented as follows: The dynamic time adjustment algorithm was used to calculate the process similarity between the two floods.

[0034] After completing the flood similarity analysis of upstream and downstream stations on the cross section, the time axis of the downstream cross section is then moved forward. Continue to analyze the similarity of floods at upstream and downstream sections, and compare them. and The similarity between the downstream flood process and the upstream flood process at a given time section is used to determine the flood process transmitted from the upstream section to the downstream section.

[0035] The implementation method of S6 is as follows: calculate the flood propagation characteristic information such as the flood propagation time and the flood difference coefficient of the upstream and downstream stations, substitute the newly obtained flood propagation time into the time initially proposed in S5, take the average value of the two, and replace the flood propagation time initially proposed in S5, and carry out a new round of flood process similarity analysis.

[0036] Secondly, this application provides a flood wave continuous self-calibration extraction system. The system includes a memory and a processor. The memory includes a program for a flood wave continuous self-calibration extraction method. When the program for the flood wave continuous self-calibration extraction method is executed by the processor, it performs the following steps: S1. Select any river section, identify the upstream and downstream hydrological / water level stations, use the station flood reporting data, arbitrarily select the flow / water level process data within a certain period, plot the flow / water level development process graph with time as the horizontal axis and flow / water level as the vertical axis, calculate the flow / water level slope between two adjacent time nodes, form a time series of flow / water level rise and fall slopes, select the time period maximum / minimum values ​​of the flow / water level sequence of the upstream and downstream stations of the section, and proceed to S2 for the next calculation.

[0037] S2. Statistically analyze the maximum value series of the two stations upstream and downstream of the cross section, and set the flood peak threshold for the two stations in three levels: small and medium flood, large flood, and flood exceeding the standard. Extract the sample that exceeds the flood peak threshold of a certain level from the maximum value of the flow / water level series as the flood peak flow, and transfer to S4 for the next step of calculation.

[0038] S3. If there are m consecutive flood peak flows / water levels, that is, there are no other non-flood peak maximum points between flood peaks, determine whether these m consecutive flood peaks are independent. Treat the non-independent consecutive flood peak flows / water levels as several flood peaks of an independent flood process, and proceed to S5 for the next step of calculation.

[0039] S4. For an independent flood process, starting from the flood peak location, search for the start and end times of the flood process forward and backward respectively to segment out the complete flood process, and then proceed to S6 for the next step of calculation.

[0040] S5. Preliminarily determine the flood propagation time of upstream and downstream stations on the cross section, and conduct similarity sliding analysis on the flood process of upstream and downstream stations on the cross section. If the flood events of the two stations on the cross section are similar, then the flood process of the downstream station is transmitted from the upstream flood.

[0041] S6. Calculate the flood propagation characteristics information such as the flood propagation time and the flood difference coefficient at upstream and downstream stations on the cross section, and correct the flood propagation time at upstream and downstream stations on the cross section in S5 based on the calculation results.

[0042] Thirdly, embodiments of this application provide a computer-readable storage medium storing program code, which, when executed by a processor, implements the steps of the flood wave continuous self-correction extraction method as described above.

[0043] Compared with the prior art, the beneficial effects of the present invention are:

[0044] By developing an automated and standardized method for extracting flood waves (including both single-peak and multi-peak types) and employing self-calibration technology for key parameters, this invention enables simple, efficient, and accurate analysis of long-sequence hydrological processes at multiple stations. Compared to existing technologies, this invention is the first to propose a continuous self-calibrated flood wave extraction method, which can effectively analyze the flood propagation characteristics among numerous stations in a watershed, thereby enhancing the understanding of flood propagation patterns in the watershed. Attached Figure Description

[0045] To more clearly illustrate the technical solutions of the embodiments of this application, the accompanying drawings used in the embodiments of this application will be briefly introduced below. It should be understood that the following drawings only show some embodiments of this application and should not be regarded as a limitation of the scope. For those skilled in the art, other related drawings can be obtained based on these drawings without creative effort.

[0046] Figure 1 This is a flowchart of the method of the present invention.

[0047] Figure 2 The diagram shows the flow sequence of upstream station A and downstream station B at the cross-section.

[0048] Figure 3 The results are selected for the flow sequence at upstream station A and the flood peak at downstream station B.

[0049] Figure 4 The data are the flow sequences of upstream station A and downstream station B at the cross-section.

[0050] Figure 5 The starting point of the rise at upstream station A of the cross section is determined. Detailed Implementation

[0051] The technical solutions of the embodiments of this application will now be described with reference to the accompanying drawings. It should be noted that similar reference numerals and letters in the following drawings indicate similar items; therefore, once an item is defined in one drawing, it does not need to be further defined and explained in subsequent drawings.

[0052] The terms “comprising,” “including,” or any other variations thereof are intended to cover a non-exclusive inclusion, such that a process, method, article, or apparatus that comprises a list of elements includes not only those elements but also other elements not expressly listed, or elements inherent to such a process, method, article, or apparatus. Without further limitation, an element defined by the phrase “comprising one…” does not exclude the presence of other identical elements in the process, method, article, or apparatus that includes said element.

[0053] The terms “first,” “second,” etc., are used only to distinguish one entity or operation from another, and should not be construed as indicating or implying relative importance, nor as requiring or implying any such actual relationship or order between these entities or operations.

[0054] Please see Figure 1 A method for continuous self-correction extraction of flood waves includes the following steps:

[0055] S1, Calculation of upstream station flow / hydrological extreme values ​​at cross-section.

[0056] Taking the flow process as an example, the same applies below.

[0057] The flow sequences of upstream station A and downstream station B at the cross-section are as follows: Figure 2 As shown, the slope of flow / water level fluctuation between two adjacent time points is calculated to form a time series of slopes. The slope of flow fluctuation at the upstream station of the cross section is then transferred to S2 for the next step of calculation.

[0058]

[0059] S2, Selection of Peak Flow Rate

[0060] Assume the stratified flood thresholds for upstream stations A and B of the cross-section are as shown in the table below.

[0061] Site Small and medium floods Great Flood Floods exceeding standard levels A 5000 10000 20000 B 6000 12000 15

[0062] The peak flow selection for stations A and B is as follows: Figure 3 As shown

[0063] S3, Determination of Peak Flood Independence

[0064] like Figure 4 Taking the peak flow at upstream station A as an example, assume the four maximum values ​​are consecutive peaks, with values ​​of 10000 m³ / s, 11000 m³ / s, 8000 m³ / s, and 10000 m³ / s respectively, occurring at times of 3 hours, 12 hours, 36 hours, and 48 hours respectively. There are three minimum values ​​between the four peaks: 6000 m³ / s, 2000 m³ / s, and 2000 m³ / s respectively. for, Given a value of 0.5, determine whether the following formula is satisfied. From the calculation results, it can be seen that the flood peaks of 10000m3 / s and 11000m3 / s are not independent and are a double peak. The remaining flood peaks of 8000m3 / s and 10000m3 / s are independent and can be regarded as two independent flood peaks.

[0065]

[0066] S4, Flood Process Division

[0067] like Figure 5 The method is introduced using the flood discharge process segmentation at upstream station A as an example. Assume the first peak discharge of the double-peak flood is 10000 m³ / s, corresponding to a peak occurrence time of 8:30. The discharges at the first three maximum points are 3500 m³ / s, 3000 m³ / s, and 8000 m³ / s, corresponding to peak occurrence times of 6:30, 4:30, and 1:30. Calculate the slope between two adjacent discharge maximum points:

[0068]

[0069] Therefore, the starting point is the first minimum point to the left of the peak flow of 3000 m³ / s, assuming a corresponding flow rate of 2000 m³ / s, occurring at 4:00 AM. The first performance value point after the second peak of the double-peak pattern is the drawdown point, assuming a corresponding flow rate of 2000 m³ / s, occurring at 10:30 AM. ,satisfy:

[0070]

[0071] The same method was used to segment the flood flow / water level process at downstream stations of the cross section.

[0072] Similarity analysis of flood processes at upper and lower sections of S5

[0073] Assuming the peak flood discharge at the upstream section is 10000 m³ / s and the corresponding peak flood discharge at the downstream section is 10500 m³ / s, calculate the flood volume of the flood events at the upstream and downstream sections. , The following formula is used to analyze the similarity of two characteristic indicators, flood peak and flood volume, between two floods at upstream and downstream stations on the cross-section, assuming...

[0074]

[0075] In the formula, Indicates the upstream section flood process sequence, This represents the time series of flood events at the upstream section. Indicates the downstream section flood process sequence, This represents the time series of flood events at the downstream section; , Represents empirical factors. ; The closer it gets to 1.

[0076] Assumption , If the sequence lengths are consistent, then the upstream station flood events can be represented as follows: The downstream flood process can be represented as The similarity between the two flood processes is calculated using the following formula:

[0077]

[0078] In the formula The smaller the similarity, the higher the similarity between the two floods, and vice versa.

[0079] After completing the flood similarity analysis of upstream and downstream stations on the cross section, the time axis of the downstream cross section is then moved forward. We will continue to analyze the similarity of floods at upstream and downstream sections, and take the flood with the closest similarity as the flood process from the upstream section to the downstream section.

[0080] S6 parameter iterative update

[0081] Calculate flood propagation characteristics such as flood propagation time and flood difference coefficient at upstream and downstream stations. Substitute the newly obtained flood propagation time into the time initially proposed in S5, take the average of the two, and replace the flood propagation time initially proposed in S5. Conduct a new similarity analysis for the next round of flood events.

[0082] This application provides a flood wave continuous self-calibration extraction system. The system includes a memory and a processor. The memory includes a program for a flood wave continuous self-calibration extraction method. When the program for the flood wave continuous self-calibration extraction method is executed by the processor, it implements the steps of the flood wave continuous self-calibration extraction method as described above.

[0083] This application provides a computer-readable storage medium storing program code. When the program code is executed by a processor, it implements the steps of the flood wave continuous self-correction extraction method as described above.

[0084] Those skilled in the art will understand that embodiments of this application can be provided as methods, systems, or computer program products. Therefore, this application can take the form of a completely hardware embodiment, a completely software embodiment, or an embodiment combining software and hardware aspects. Furthermore, this application can take the form of a computer program product embodied on one or more computer-usable storage media (including but not limited to disk storage, CD-ROM, optical storage, etc.) containing computer-usable program code.

[0085] This application is described with reference to flowchart illustrations and / or block diagrams of methods, apparatus (systems), and computer program products according to embodiments of this application. It will be understood that each block of the flowchart illustrations and / or block diagrams, and combinations of blocks in the flowchart illustrations and / or block diagrams, can be implemented by computer program instructions. These computer program instructions can be provided to a processor of a general-purpose computer, special-purpose computer, embedded processor, or other programmable data processing apparatus to produce a machine, such that the instructions, which execute via the processor of the computer or other programmable data processing apparatus, generate instructions for implementing the flowchart... Figure 1 One or more processes and / or boxes Figure 1 A device that provides the functions specified in one or more boxes.

[0086] These computer program instructions may also be stored in a computer-readable storage medium that can direct a computer or other programmable data processing device to function in a particular manner, such that the instructions stored in the computer-readable storage medium produce an article of manufacture including instruction means, which are implemented in a process Figure 1 One or more processes and / or boxes Figure 1 The function specified in one or more boxes.

[0087] These computer program instructions may also be loaded onto a computer or other programmable data processing equipment to cause a series of operational steps to be performed on the computer or other programmable equipment to produce a computer-implemented process, thereby providing instructions that execute on the computer or other programmable equipment for implementing the process. Figure 1One or more processes and / or boxes Figure 1 The steps of the function specified in one or more boxes.

[0088] In a typical configuration, a computing device includes one or more processors (CPU), input / output interfaces, network interfaces, and memory.

[0089] Memory may include non-persistent memory in computer-readable media, such as random access memory (RAM) and / or non-volatile memory, such as read-only memory (ROM) or flash RAM. Memory is an example of computer-readable media.

[0090] Computer-readable media includes both permanent and non-permanent, removable and non-removable media that can store information using any method or technology. Information can be computer-readable instructions, data structures, modules of programs, or other data. Examples of computer storage media include, but are not limited to, phase-change memory (PRAM), static random access memory (SRAM), dynamic random access memory (DRAM), other types of random access memory (RAM), read-only memory (ROM), electrically erasable programmable read-only memory (EEPROM), flash memory or other memory technologies, CD-ROM, digital versatile optical disc (DVD) or other optical storage, magnetic tape, magnetic magnetic disk storage or other magnetic storage devices, or any other non-transferable medium that can be used to store information accessible by a computing device. As defined herein, computer-readable media does not include transient computer-readable media, such as modulated data signals and carrier waves.

[0091] The above description is merely an embodiment of this application and is not intended to limit the scope of protection of this application. Various modifications and variations can be made to this application by those skilled in the art. Any modifications, equivalent substitutions, improvements, etc., made within the spirit and principles of this application should be included within the scope of protection of this application.

Claims

1. A method for continuous self-correcting extraction of flood waves, characterized in that, The specific steps include the following: S1. Select any river section, identify the upstream and downstream hydrological / water level stations, use the station flood reporting data, arbitrarily select the flow / water level process data within a certain period of time, plot the flow / water level development process graph with time as the horizontal axis and flow / water level as the vertical axis, calculate the flow / water level slope between two adjacent time nodes, form a time series of flow / water level rise and fall slopes, select the time period maximum / minimum values ​​of the flow / water level series of upstream and downstream stations, and proceed to S2 for the next step of calculation; S2. Statistically analyze the maximum value series of the two stations upstream and downstream of the cross section, and set the flood peak threshold for the two stations in three levels: small and medium flood, large flood, and flood exceeding the standard. Extract the sample exceeding the flood peak threshold from the maximum value of the flow / water level series as the flood peak flow, and transfer to S4 for the next step of calculation. S3. If there are m consecutive flood peak flows / water levels, that is, there are no other non-flood peak maximum points between flood peaks, determine whether these m consecutive flood peaks are independent. Treat the non-independent consecutive flood peak flows / water levels as several flood peaks of an independent flood process, and proceed to S5 for the next step of calculation. S4. For an independent flood process, starting from the flood peak location, search for the start and end times of the flood process forward and backward respectively to segment out the complete flood process, and then proceed to S6 for the next step of calculation. S5. Preliminarily determine the flood propagation time of upstream and downstream stations on the cross section, and conduct similarity sliding analysis on the flood process of upstream and downstream stations on the cross section. If the flood events of the two stations on the cross section are similar, then the flood process of the downstream station is transmitted from the upstream flood. S6. Calculate the flood propagation characteristics information such as the flood propagation time and the flood difference coefficient at upstream and downstream stations on the cross section, and correct the flood propagation time at upstream and downstream stations on the cross section in S5 based on the calculation results.

2. The flood wave continuous self-correction extraction method according to claim 1, characterized in that, The implementation method of S1 is as follows: Select the flow / water level process of the upstream station of the cross section. / Flow / water level process at downstream stations of the cross section / Calculate the slope of flow / water level fluctuations between two adjacent time points to form a time series of slopes, including the flow / water level fluctuation slopes at upstream stations of the cross-section. The slope of flow / water level fluctuation at downstream stations of the cross section ,in, or , or The extreme values ​​of a time period can be determined using the following formula: 。 3. The flood wave continuous self-correction extraction method according to claim 1, characterized in that, The implementation method of S2 is as follows: the maximum value of flow / water level at the upstream station of the cross section. / Maximum flow / water level at downstream stations of the cross-section / Based on the design flood data from upstream and downstream stations, determine the design flood discharge / water level corresponding to minor, major, and above-standard floods. / , / Statistical analysis was conducted to determine the flood type by identifying the points with the maximum flow / water level at upstream and downstream stations within a selected time period. ， or , ， or .

4. The flood wave continuity self-correction extraction method according to claim 1, characterized in that, The implementation of S3 is as follows: Based on the above analysis, all flood peaks within the study period can be selected. If m consecutive maximum points all belong to the flood category, it is necessary to determine whether these m consecutive flood peaks are independent. Assume the consecutive flood peaks are... The corresponding peak time is The time difference between the peaks of two adjacent flood peaks in a series of consecutive flood peaks is ,in If the following formula is satisfied, then m consecutive flood peaks are not independent and can be regarded as the same flood. ， In the formula, The base flow for site A; Refers to continuous flood peaks and The interval time; These are two consecutive consecutive flood peak flows; This refers to the minimum flow rate between two flood peaks; , This represents an empirical parameter.

5. The flood wave continuous self-correction extraction method according to claim 1, characterized in that, The implementation of S4 is as follows: For an independent flood process, starting from the flood peak location, the start and end times of the flood process are searched forward and backward respectively to segment the complete flood process. The method for determining the start point of the flood process is as follows: if the flood peak is single-peak and the flood peak discharge is... The corresponding peak flood time is Then, search backward along the time axis for adjacent independent flood peaks, and count all the maximum points within the time interval between the occurrence of two independent flood peaks. ,in Represents the target flood peak Calculate the slope between two adjacent points with maximum flow rates. ,like If all are greater than 0, then The corresponding minimum point to the left of the maximum point The starting point of this flood, and the corresponding starting time is... Search forward along the time axis for adjacent independent flood peaks. ,in Represents the target flood peak Calculate the slope between two adjacent points with maximum flow rates. If the slope between two consecutive points of maximum flow is less than 0, then The corresponding minimum point to the right of the maximum value The receding point of this flood event, and the corresponding receding time is... , ， In the formula, and These are empirical parameters. , If the target flood peak is multi-peaked, starting from the first peak, search backwards along the time axis for adjacent independent flood peaks. If the searched independent flood peak is also a single peak, analyze the rising point between the two peaks using the above method. If the searched independent flood peak is multi-peaked, take the last peak among the multi-peaks as the object and analyze the rising point between the two peaks using the above method. If the target flood peak is multi-peaked, starting from the last peak, search forwards along the time axis for adjacent independent flood peaks. If the searched independent flood peak is a single peak, analyze the receding point between the two peaks using the above method. If the searched independent flood peak is multi-peaked, take the first peak among the multi-peaks as the object and analyze the receding point between the two peaks using the above method. Based on the time corresponding to the flood's rising and receding points, segment the complete flood process. The same method was used to segment the flood flow / water level process at downstream stations of the cross section.

6. The flood wave continuity self-correction extraction method according to claim 1, characterized in that, The implementation of S5 is as follows: For the flood flow / water level processes segmented at upstream and downstream stations on the cross-section, the peak occurrence time corresponding to the upstream flood process is... Flood peak ,magnanimity Based on the preliminary estimated upstream and downstream flood propagation time Analysis of downstream stations at the cross-section The flood process corresponding to an independent flood peak near a given time, where the flood peak occurs at a specific time. Flood peak ,magnanimity The following formula is used to analyze the similarity of two characteristic indicators, flood peak and flood volume, between two floods at upstream and downstream stations on the cross-section. ， In the formula, , Represents empirical factors. ; The closer the similarity is to 1, the higher the similarity of the peak values ​​and magnitudes of the two floods; conversely, the closer the similarity is to 1, the lower the similarity is to 1. If the flood sequence lengths at upstream and downstream stations are the same, then the flood process at the upstream station can be represented as follows: The downstream flood process can be represented as The similarity between the two flood processes is calculated using the following formula: ， The smaller the similarity, the higher the similarity between the two floods, and vice versa. If the flood sequence lengths at upstream and downstream stations are inconsistent, the flood process at the upstream station can be represented as follows: The downstream flood process can be represented as follows: The dynamic time adjustment algorithm was used to calculate the process similarity between the two floods. After completing the flood similarity analysis of upstream and downstream stations on the cross section, the time axis of the downstream cross section is then moved forward. Continue to analyze the similarity of floods at upstream and downstream sections, and compare them. and The similarity between the downstream flood process and the upstream flood process at a given time section is used to determine the flood process transmitted from the upstream section to the downstream section.

7. The flood wave continuity self-correction extraction method according to claim 1, characterized in that, The implementation method of S6 is as follows: calculate the flood propagation time and flood propagation characteristic information of upstream and downstream stations, substitute the newly obtained flood propagation time into the time initially proposed in S5, take the average value of the two, and replace the flood propagation time initially proposed in S5, and carry out a new round of flood process similarity analysis.

8. A flood wave continuous self-correcting extraction system, characterized in that, The system includes a memory and a processor. The memory includes a program for a continuous self-calibration extraction method for flood waves. When the program for the continuous self-calibration extraction method for flood waves is executed by the processor, the following steps are implemented: S1. Select any river section, identify the upstream and downstream hydrological / water level stations, use the station flood reporting data, arbitrarily select the flow / water level process data within a certain period of time, plot the flow / water level development process graph with time as the horizontal axis and flow / water level as the vertical axis, calculate the flow / water level slope between two adjacent time nodes, form a time series of flow / water level rise and fall slopes, select the time period maximum / minimum values ​​of the flow / water level series of upstream and downstream stations of the section, and proceed to S2 for the next step of calculation. S2. Statistically analyze the maximum value series of the two stations upstream and downstream of the cross section, and set the flood peak threshold for the two stations in three levels: small and medium flood, large flood, and flood exceeding the standard. Extract the sample exceeding the flood peak threshold from the maximum value of the flow / water level series as the flood peak flow, and transfer to S4 for the next step of calculation. S3. If there are m consecutive flood peak flows / water levels, that is, there are no other non-flood peak maximum points between flood peaks, determine whether these m consecutive flood peaks are independent. Treat the non-independent consecutive flood peak flows / water levels as several flood peaks of an independent flood process, and proceed to S5 for the next step of calculation. S4. For an independent flood process, starting from the flood peak location, search for the start and end times of the flood process forward and backward respectively to segment out the complete flood process, and then proceed to S6 for the next step of calculation. S5. Preliminarily determine the flood propagation time of upstream and downstream stations on the cross section, and conduct similarity sliding analysis on the flood process of upstream and downstream stations on the cross section. If the flood events of the two stations on the cross section are similar, then the flood process of the downstream station is transmitted from the upstream flood. S6. Calculate the flood propagation characteristics information such as the flood propagation time and the flood difference coefficient at upstream and downstream stations on the cross section, and correct the flood propagation time at upstream and downstream stations on the cross section in S5 based on the calculation results.

9. A computer-readable storage medium, characterized in that, The computer-readable storage medium stores program code that, when executed by a processor, implements the steps of the flood wave continuous self-correction extraction method as described in any one of claims 1 to 7.

Images