Method and device for calculating surge amplitude of bifurcated channel

Abstract

The invention discloses a method and a device for calculating the surge amplitude of a bifurcated channel. The method comprises the following steps: acquiring a surge wave of a surge before shunting at a bifurcated opening of the bifurcated channel; the method comprises the following steps: performing simplified fitting on a surge wave in front of a bifurcated opening to form an initial surge waveform; decomposing and dispersing the initial surge waveform into independent surge waves; on the basis of the initial surge waveform and the initial height of the independent surge wave, the amplitude of the independent surge wave at the bifurcated opening position at the moment t1 before collision reflection occurs in the shunting process is calculated; according to the movement and reflection path of each independent swell wave in the bifurcated mouth, the position of the independent swell wave at the moment t2 after collision reflection of the independent swell wave occurs is obtained; overlapping the amplitudes of the independent surge waves at the same position at the t2 moment to obtain the overlapped amplitude of the independent surge waves; and determining the superposed wave amplitude meeting a preset condition as the wave amplitude of the surge in the bifurcated channel. The invention provides a bifurcated channel amplitude calculation method.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of surge wave amplitude, and particularly relates to a method and device for calculating surge wave amplitude of a bifurcated river. BACKGROUND

[0002] The surge formula calculation system is the most rapid calculation method for evaluating surge disasters. Existing research shows that when the surge generated by river embankment landslide and the like propagates in a single river, the propagation process can be divided into a ring propagation process and a parallel propagation process. The ring propagation area can be defined as a ring area with the landslide water entry point as the center and the river width as the radius, and the parallel propagation area is the remaining area affected by the surge, which has a larger area and range.

[0003] During the propagation process of the surge, the bifurcation of the river is encountered, and the bifurcation causes the surge to be divided, thereby changing the propagation direction of the surge and the surge wave amplitude.

[0004] Existing surge calculation methods are mostly for a single river, and the surge wave is simplified as a constant flow, without considering the influence of surge division on the calculation of the surge wave amplitude, so that the wave shape change before and after the surge division cannot be accurately calculated, and the surge wave amplitude at the bifurcation position cannot be calculated. SUMMARY

[0005] The present application provides a method and device for calculating the surge wave amplitude of a bifurcated river to solve the problem that the existing technology cannot calculate the surge wave amplitude at the bifurcation position.

[0006] The first aspect of the embodiment of the present application provides a method for calculating the surge wave amplitude of a bifurcated river, comprising: obtaining the surge wave before the bifurcation of the bifurcated river; simplifying and fitting the surge wave before the bifurcation to form an initial surge wave shape; discretely decomposing the initial surge wave shape in a direction perpendicular to the original propagation direction of the surge wave into a plurality of independent surge waves with the same wavelength and different initial positions; based on the initial surge wave shape and the initial height of the independent surge wave, calculating the wave amplitude of the independent surge wave at the bifurcation position at a time t1 before the collision and reflection of the independent surge wave during the division process; calculating the motion and reflection path of each independent surge wave in the bifurcation, and obtaining the position of the independent surge wave at a time t2 after the collision and reflection, wherein the time t2 is located after the time t1 in time sequence; superimposing the wave amplitudes of the independent surge waves at the same position at the time t2 to obtain the superimposed wave amplitude of the independent surge wave, wherein the wave amplitude of the independent surge wave at the time t2 is determined based on the wave amplitude of the independent surge wave at the time t1. The amplitude of the superimposed wave meeting the preset condition is determined as the amplitude of the surge in the branch channel.

[0007] In some optional embodiments of the first aspect of the embodiments of the present application, the method further includes: obtaining position information of the surge; determining the initial amplitude of the surge at the position of the surge according to the instability mode of the landslide and a corresponding calculation formula; determining the surge wave before the surge is divided at the branch mouth of the branch channel according to the initial amplitude of the surge, the distance between the position of the surge and the branch mouth, and a generalized model of the branch channel.

[0008] In some optional embodiments of the first aspect of the embodiments of the present application, the method further includes: simplifying and fitting the surge wave before the branch mouth according to a wave equation to form an initial surge wave shape, wherein the wave equation includes a Fourier series equation and a cosine wave equation.

[0009] In some optional embodiments of the first aspect of the embodiments of the present application, the method further includes: determining the amplitude of the independent surge wave at the position of the branch mouth at the time t1 before the collision and reflection of the independent surge wave during the division process according to the initial height of the independent surge wave and the initial height of the initial surge wave shape. determining the amplitude of the independent surge wave at the position of the branch mouth at the time t1 before the collision and reflection of the independent surge wave during the division process c i i = 1, 2, 3,...), wherein, the initial height of each independent surge wave, x l the wave range of the surge in the division direction during the division process, S 0 is the riverbed slope, and g is the acceleration of gravity.

[0010] In some optional embodiments of the first aspect of the embodiments of the present application, the method further includes: defining the propagation direction of the surge in the original propagation channel as the y direction, defining the propagation direction of the surge in the branch channel as the x direction, and establishing a calculation grid at the position of the branch mouth of the branch channel according to the width of the original propagation channel and the width of the branch channel; calculating the movement and reflection path of each independent surge wave in the branch mouth of the branch channel based on the mirror reflection principle according to the calculation grid. ​​

[0011] With reference to the first aspect, in some possible embodiments, the determining the superimposed wave amplitude that meets the preset condition as the wave amplitude of the surge in the braided river channel comprises: determining the maximum superimposed wave amplitude at the position of the center line in the braided river channel as the wave amplitude of the surge in the braided river channel.

[0012] With reference to the first aspect, in some possible embodiments, the method further comprises: determining the influence range and damage degree of the surge on the river channel according to the determined wave amplitude of the surge in the braided river channel.

[0013] The second aspect of the present application provides a device for calculating the wave amplitude of a surge in a braided river channel, comprising: a first obtaining module configured to obtain the wave of the surge before the surge is branched at a branching point of the braided river channel; a forming module configured to form an initial surge wave shape by simplifying and fitting the wave of the surge before the surge is branched at the branching point of the braided river channel; a decomposing and discretizing module configured to decompose and discretize the initial surge wave shape into a preset number of independent surge waves with the same wavelength and different initial positions in a direction perpendicular to the original propagation direction of the surge; a first calculating module configured to calculate the wave amplitude of each independent surge wave at the branching point at a time t1 before the independent surge waves collide and reflect during the branching process based on the initial surge wave shape and the initial height of the independent surge waves; a second calculating module configured to calculate the movement and reflection path of each independent surge wave in the branching point, and obtain the position of each independent surge wave at a time t2 after the independent surge waves collide and reflect, wherein the time t2 is later than the time t1 in time sequence; a second obtaining module configured to superimpose the wave amplitudes of the independent surge waves at the same position at the time t2 to obtain the superimposed wave amplitude of the independent surge waves, wherein the wave amplitude of the independent surge waves at the time t2 is determined based on the wave amplitude of the independent surge waves at the time t1. a determining module configured to determine the superimposed wave amplitude that meets the preset condition as the wave amplitude of the surge in the braided river channel.

[0014] The third aspect of the present application further provides an electronic device, comprising a processor, a memory, and a program stored in the memory and executable on the processor, and the program, when executed by the processor, implements the steps of the method for calculating the wave amplitude of a surge in a braided river channel provided by the embodiments of the present application.

[0015] The fourth aspect of the present application further provides a computer readable storage medium, and the computer readable storage medium stores a computer program, and the computer program, when executed by a processor, implements the steps of the method for calculating the wave amplitude of a surge in a braided river channel provided by the embodiments of the present application.

[0016] The technical scheme provided by the embodiments of the present application has at least the following beneficial effects: In the embodiments of the present application, the wave height of the surge in the bifurcation river is calculated by obtaining the wave of the surge before the bifurcation of the bifurcation river; the initial wave shape of the surge is formed by simplifying and fitting the wave of the surge before the bifurcation; the initial wave shape of the surge is decomposed and dispersed into a plurality of independent surge waves with the same wavelength and different initial positions in the direction perpendicular to the original propagation direction of the surge; the wave amplitude of the independent surge wave at the bifurcation position at time t1 before the collision and reflection of the independent surge wave during the bifurcation process is calculated based on the initial wave shape of the surge and the initial height of the independent surge wave; the movement and reflection path of each independent surge wave in the bifurcation are calculated, and the position of the independent surge wave at time t2 after the collision and reflection is obtained, wherein t2 is in time sequence after t1; the wave amplitudes of the independent surge waves at the same position at time t2 are superimposed to obtain the superimposed wave amplitude of the independent surge wave, wherein the wave amplitude of the independent surge wave at time t2 is determined based on the wave amplitude of the independent surge wave at time t1; and the superimposed wave amplitude meeting the preset condition is determined as the wave amplitude of the surge in the bifurcation river. Thus, a bifurcation river surge wave amplitude calculation method is provided, which considers the wave shape characteristics of the surge wave and can accurately solve the wave shape and amplitude of the surge after the bifurcation in the bifurcation river, thereby providing technical guarantee for more accurate evaluation of surge disasters. BRIEF DESCRIPTION OF DRAWINGS

[0017] Figure 1 Flowchart of the bifurcation river surge wave amplitude calculation method provided by the embodiments of the present application Figure One ; Figure 2 Schematic diagram of the fitting and dispersion processing of the wave shape of the surge wave when the surge wave propagates to the bifurcation position in the bifurcation river surge wave amplitude calculation method provided by the embodiments of the present application Figure 3 Schematic diagram of the surge wave at time t1 after entering the bifurcation in the bifurcation river surge wave amplitude calculation method provided by the embodiments of the present application Figure 4 Schematic diagram of the wave amplitude calculation results of the surge wave at different times during the movement process of the bifurcation in the bifurcation river surge wave amplitude calculation method provided by the embodiments of the present application Figure 5 Schematic diagram of the actual application embodiment of the bifurcation river surge wave amplitude calculation method provided by the embodiments of the present application Figure 6 Flowchart of the bifurcation river surge wave amplitude calculation method provided by the embodiments of the present application Figure Two . DETAILED DESCRIPTION

[0018] In order to make the purposes, technical solutions and advantages of the embodiments of the present application clearer, the technical solutions of the embodiments of the present application will be described clearly and completely below with reference to the drawings of the embodiments of the present application. Obviously, the described embodiments are part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by a person of ordinary skill in the art without creative work fall within the protection scope of the present application.

[0019] The terms "first", "second", and the like in the specification and claims of the present application are used to distinguish similar objects, and are not used to describe a specific order or sequence. It should be understood that the data used in this way can be interchanged under appropriate circumstances, so that the embodiments of the present application can be implemented in an order other than those illustrated or described herein, and the objects distinguished by "first", "second", etc. are generally a category and do not limit the number of objects, for example, the first object can be one or more. In addition, "and / or" in the specification and claims indicates at least one of the connected objects, and the character " / ", generally indicates that the front and rear associated objects are in an "or" relationship.

[0020] In addition, the technical features involved in different embodiments of the present application described below can be combined with each other as long as there is no conflict between them.

[0021] The branched river channel surge wave amplitude calculation method provided by the embodiments of the present application will be described in detail below with reference to the drawings, through specific embodiments and their application scenarios.

[0022] As shown in the figure, the branched river channel surge wave amplitude calculation method can specifically include: Figure 1 Step 101, obtaining the surge wave before the surge wave is divided at the branch mouth of the branched river channel. Step 101, obtaining the surge wave before the surge wave is divided at the branch mouth of the branched river channel.

[0023] The surge wave is usually generated by the collapse of the river bank, the falling of the object into the water, and the like. In order to accurately obtain the amplitude information of the surge wave, the technical solution of the present application selects the corresponding calculation formula such as the river channel generalization model, the wave making formula, the propagation wave formula and the wave climbing height formula according to the collapse mode of the landslide such as falling, dumping, sliding and the like, and calculates and determines the initial amplitude of the surge wave.

[0024] For example, for the landslide surge wave, the following calculation formula is used to calculate and determine the maximum wave amplitude of the surge wave:

[0025] The specific parameters include: landslide speed v, gravity acceleration g, landslide water entry angle a , original propagation river width B0, landslide thickness s, landslide water entry length l sw , water depth hw , landslide into water width w w The equalization parameter is used to calculate the maximum surge wave amplitude, and the maximum surge wave amplitude is determined as the initial wave amplitude of the surge wave at the disaster surge wave occurrence position.

[0026] After obtaining the initial wave amplitude of the surge wave, the maximum wave amplitude value of the surge wave generated by the landslide body instability into water is calculated according to the selected wave making formula and the propagation wave formula a max , that is, the initial wave amplitude of the surge wave at the surge wave occurrence position, and the wave amplitude of the surge wave propagating to the bifurcation a c0max .

[0027] Taking the landslide surge wave propagation calculation formula as an example, the wave amplitude of the surge wave before the bifurcation is calculated through the distance l of the bifurcation from the landslide position and the water depth hw a c0max .

[0028] It can be seen that in the technical scheme of the present application, the surge wave wave before the bifurcation of the bifurcated river channel can include: obtaining the surge wave occurrence position information; According to the landslide body instability mode, the corresponding calculation formula is selected to calculate and determine the initial wave amplitude of the surge wave at the surge wave occurrence position; According to the initial wave amplitude of the surge wave, the distance between the surge wave occurrence position and the bifurcation, and the bifurcated river channel generalization model, the wave of the surge wave before the bifurcation of the bifurcated river channel is determined.

[0029] Step 102, the initial surge wave shape is formed by simplifying and fitting the surge wave before the bifurcation.

[0030] In the embodiment of the present application, the wave shape characteristics of the surge wave are considered, so that the wave shape and amplitude of the surge wave after the bifurcation of the bifurcated river channel can be accurately solved, and the technical defects of the prior art that lack research on the wave shape of the surge wave are overcome. At the same time, the surge wave calculation requirements of the surge wave bifurcation calculation method in the present patent are met, which provides important technical support for the surge wave prediction and disaster prevention and reduction of complex river channels such as reservoir areas.

[0031] Specifically, the surge wave before the bifurcation can be simplified and fitted by a Fourier series or a cosine wave equation to form an initial surge wave shape.

[0032] In a specific embodiment, the wavelength of the surge wave needs to be determined during the fitting process l which is much higher than the bifurcated river channel width B.

[0033] In another specific embodiment, the maximum wave amplitude of the simplified surge wave can be equal to the wave amplitude of the surge wave before the bifurcation a c0max .

[0034] Step 103: Decompose the initial surge waveform into a preset number of independent surge waves with the same wavelength but different initial positions in a direction perpendicular to the original propagation direction of the surge.

[0035] Specifically, the simplified initial surge waveform before the split is decomposed and discretized into 100 waveforms with the same wavelength perpendicular to the original propagation direction of the surge. l / 100, and independent surge waves with different initial positions, each with an initial height of . c0 (i) In a specific embodiment, the maximum amplitude of all independent surge waves can be: a c0max .

[0036] The schematic diagram of the wave waveform fitting and discretization process when the wave propagates to the bifurcation point in the wave amplitude calculation method of the bifurcation channel provided in this application embodiment can be seen as follows: Figure 2 As shown.

[0037] Step 104: Based on the initial surge waveform and the initial height of the independent surge wave, calculate the amplitude of the independent surge wave at the bifurcation point at time t1 before the collision reflection occurs during the diversion process.

[0038] In this step, the specific details can be determined according to... Calculate the amplitude of the independent surge wave at the bifurcation point at time t1 before collision reflection during the flow splitting process. c ( i )( i = 1,2,3...), where, The initial height of each individual surge wave, x l This refers to the extent of the surge's reach in the direction of the diversion process. S 0 represents the riverbed slope, and g represents the acceleration due to gravity.

[0039] The schematic diagram of the surge amplitude calculation method for bifurcation channels provided in this application at time t1 after the surge enters the bifurcation is shown in the following embodiment. Figure 3 As shown.

[0040] Step 105: Calculate and obtain the motion and reflection path of each independent surge wave within the branching point, and obtain the position (x, y) of the independent surge wave at time t2 after the collision and reflection.

[0041] The original propagation direction of the surge is defined as the y direction, and the direction in which the deflection occurs to generate the flow splitting is defined as the x direction. A calculation grid of the location of the branch mouth of the branch channel is established according to the grid size, which is the original propagation channel width B0 / 100 multiplied by the branch channel width B / 100. Based on the mirror reflection principle, the movement and reflection path of each column of independent surge waves in the branch mouth are calculated and obtained. The above t2 time is located in time sequence after the t1 time.

[0042] In other embodiments, the grid size can also be established according to other specific parameters.

[0043] The position y0 of the independent surge wave at the t2 time after entering the branch mouth can be calculated by y0=v0t0, where v0 is the wave speed of the surge before the flow splitting, t0 is the time of the surge wave in the branch mouth, and t2 is the time of the surge wave in the branch mouth. o is the time of the surge wave in the branch mouth.

[0044] The calculation results of the wave amplitude of the surge at different times in the process of the movement of the surge in the branch mouth of the branch channel provided by the embodiments of the present application can be as shown in Figure 4

[0045] Step 106, superimposing the wave amplitudes of the independent surge waves at the same position at the t2 time to obtain the superimposed wave amplitude of the independent surge waves, wherein the wave amplitude of the independent surge wave at the t2 time is determined based on the wave amplitude of the independent surge wave at the t1 time.

[0046] According to the position information in the calculation grid, the wave amplitudes of the independent surge waves at the same position at the t2 time are superimposed to obtain the superimposed wave amplitude of the independent surge waves.

[0047] The wave amplitude of the independent surge wave at the t2 time can be specifically determined based on the wave amplitude of the independent surge wave at the t1 time before the collision and reflection of the surge.

[0048] For example, the wave amplitude of the independent surge wave at the t2 time can be directly determined as the wave amplitude of the independent surge wave at the t1 time before the collision and reflection of the surge; or, the wave amplitude of the independent surge wave at the t2 time can also be determined based on the wave amplitude of the independent surge wave at the t1 time before the collision and reflection of the surge according to mature and feasible ways such as corresponding attenuation standards and actual use experience. For this, the present application does not make specific limitations.

[0049] Step 107, determining the superimposed wave amplitude that meets the preset condition as the wave amplitude of the surge in the branch channel.

[0050] The maximum surge wave amplitude at the center line position in the branch channel can be taken as the maximum surge wave amplitude in the subsequent propagation process a c1 , and the maximum surge wave amplitudes in all branch channels and the original propagation channel at the time when the flow splitting is completed are calculated. Wherein, the subscript 1 represents the branch channel number. ​

[0051] Alternatively, based on references such as the different degrees of damage corresponding to different surge amplitudes, surges with amplitudes higher than a preset height value can be defined as the surge amplitude in the bifurcation channel.

[0052] The wave amplitude calculation method for swells in bifurcation channels provided in the above embodiments of the present invention can accurately solve the waveform and amplitude of swells after they split in bifurcation channels by combining the waveform characteristics of swells, thus providing a guarantee for more accurate evaluation of swell disasters.

[0053] The wave amplitude calculation method for bifurcation channels provided in this invention can achieve significant technical effects in practical applications, for example, as... Figure 5 As shown (wave amplitude distribution of the surge generated by the Gongjiafang landslide at the confluence of the Yangtze River and the Daning River), based on the technical solution of this application, the surge of the Gongjiafang landslide that occurred in Wushan was calculated. The calculated surge amplitude at the confluence of the Daning River and the Yangtze River in Wushan County was 0.36 m, which is only 10% different from the actual survey result of 0.4 m.

[0054] In another specific embodiment, such as Figure 6 As shown, the technical solution of this application also includes: Step 108: Determine the extent of the impact of the surge on the river channel and the degree of damage based on the determined surge amplitude in the bifurcation channel.

[0055] Specifically, after the surge is diverted, assuming that the attenuation of the surge during its subsequent propagation still follows the propagation attenuation law within a single river channel, a c1 replace a max Substitute return Calculations are performed to obtain the surge amplitude during subsequent propagation, where... l This is the distance from the calculated point along the river to the landslide.

[0056] Obtain the wave amplitude distribution across the entire river channel affected by the swell, and calculate the swell run-up R at all locations along the bank using the swell run-up formula. For example, based on the wave amplitude 'a' at the corresponding location in the river channel c ( l ), bank slope inclination angle β Calculate the wave height along the swell path.

[0057] After determining the maximum surge wave amplitude, the embodiment of the present application can solve the propagation attenuation of the surge wave based on the real situation of the river channel through a corresponding mode, calculate the falling dangerous rock or landslide by using a surge wave attenuation prediction model such as a side-limited or non-side-limited propagation area, and divide the intensity of the surge wave in the region, for example, the surge wave height A>2m, 2m>A>1.5m, 1.5m>4>1m, 1m>A>0.5m, and A≤0.5m are defined as "extremely high surge wave area", "very high surge wave area", "high surge wave area", "medium surge wave area", and "low surge wave area" respectively, and determine the influence range and damage degree of the surge wave on the river channel based on the above surge wave area.

[0058] The method for calculating the surge wave amplitude of the branched river channel provided by the embodiment can provide important technical support for surge wave prediction and disaster prevention and reduction in complex river channels such as reservoir areas.

[0059] The embodiment of the present application also provides a device for calculating the surge wave amplitude of a branched river channel, comprising: A first acquisition module is configured to acquire the surge wave amplitude before the surge wave is divided at the branch mouth of the branched river channel. A forming module is configured to form an initial surge wave shape by simplifying and fitting the surge wave amplitude before the surge wave is divided at the branch mouth of the branched river channel. A decomposition and discretization module is configured to decompose and discretize the initial surge wave shape into a plurality of independent surge waves with the same wavelength and different initial positions in a direction perpendicular to the original propagation direction of the surge wave. A first calculation module is configured to calculate the wave amplitude of each independent surge wave at the branch mouth at a time t1 before the collision and reflection of the independent surge wave during the division process based on the initial surge wave shape and the initial height of the independent surge wave. A second calculation module is configured to calculate the movement and reflection path of each independent surge wave in the branch mouth, and acquire the position (x, y) of the independent surge wave at a time t2 after the collision and reflection of the independent surge wave, wherein the time t2 is located after the time t1 in time sequence. A second acquisition module is configured to superimpose the wave amplitudes of the independent surge waves at the same position at the time t2 to obtain the superimposed wave amplitude of the independent surge wave, wherein the wave amplitude of the independent surge wave at the time t2 is determined based on the wave amplitude of the independent surge wave at the time t1. A determination module is configured to determine the superimposed wave amplitude that meets a preset condition as the wave amplitude of the surge wave in the branched river channel.

[0060] The device for calculating the surge wave amplitude of the branched river channel provided by the embodiment of the present application has the same technical process and technical effects as the method for calculating the surge wave amplitude of the branched river channel provided by the embodiment of the present application, and thus repeated description is omitted here.

[0061] The embodiment of the present application further provides an electronic device, comprising a processor, a storage, a computer program stored in the storage and executable on the processor, which, when executed by the processor, implements each process of the embodiment of the bifurcation river channel wave amplitude calculation method and achieves the same technical effects. To avoid repetition, details are not described herein.

[0062] The embodiment of the present application further provides a computer readable storage medium, which stores a computer program, and the computer program, when executed by a processor, implements each process of the embodiment of the bifurcation river channel wave amplitude calculation method and achieves the same technical effects. To avoid repetition, details are not described herein. The computer readable storage medium is, for example, a Read-Only Memory (ROM), a Random Access Memory (RAM), a magnetic disk or an optical disk.

[0063] It should be noted that, in this document, the terms "comprising", "including", or any other variant thereof are intended to cover non-exclusive inclusions, such that processes, methods, articles, or devices that include a series of elements not only include those elements, but also include other elements not explicitly listed, or inherent to such processes, methods, articles, or devices. Without more limitations, the element defined by the phrase "including a" does not exclude the presence of another identical element in the process, method, article, or device that includes the element.

[0064] From the above description of the embodiments, those skilled in the art can clearly understand that the above embodiment methods can be implemented by means of software plus a necessary general hardware platform, and of course, can also be implemented by hardware, but in many cases, the former is a better embodiment. Based on such understanding, the technical solutions of the present application can be embodied in the form of a software product, which is stored in a storage medium (such as a ROM / RAM, a magnetic disk, or an optical disk) and includes a number of instructions for causing a terminal (which can be a mobile phone, a computer, a server, an air conditioner, or a network device) to execute the methods described in the various embodiments of the present application.

[0065] Obviously, the above embodiments are merely examples for clear illustration, and are not limitations on the embodiments. For those skilled in the art, other different forms of changes or modifications can be made on the basis of the above description. Here, it is not necessary or possible to exhaust all the embodiments. The obvious changes or modifications derived therefrom are still within the scope of protection of the present application.

Claims

1. A method for calculating the wave amplitude of a braided river, characterized by, The method comprises the following steps: obtaining a surge wave before the surge wave is divided at a bifurcation of a bifurcated river channel; forming an initial surge wave shape by simplifying and fitting the surge wave before the surge wave is divided at the bifurcation; discretely dividing the initial surge wave shape in a direction perpendicular to the original propagation direction of the surge wave into a preset number of independent surge waves with the same wavelength and different initial positions; calculating the amplitude of the independent surge waves at the bifurcation at a time t1 before the independent surge waves collide and reflect during the division process based on the initial surge wave shape and the initial height of the independent surge waves; calculating the movement and reflection path of each independent surge wave in the bifurcation to obtain the position of the independent surge wave at a time t2 after the independent surge wave collides and reflects, wherein the time t2 is later than the time t1 in time sequence; superimposing the amplitudes of the independent surge waves at the same position at the time t2 to obtain a superimposed amplitude of the independent surge waves, wherein the amplitude of the independent surge wave at the time t2 is determined based on the amplitude of the independent surge wave at the time t1; determining the superimposed amplitude that meets a preset condition as the amplitude of the surge wave in the bifurcated river channel.

2. The method of claim 1, wherein, The method for obtaining the surge wave before the surge wave is divided at the bifurcation of the bifurcated river channel comprises the following steps: obtaining position information of the surge wave; selecting a corresponding calculation formula to calculate the initial amplitude of the surge wave at the position of the surge wave according to a collapse and slide instability mode; determining the surge wave before the surge wave is divided at the bifurcation of the bifurcated river channel based on the initial amplitude of the surge wave, the distance between the position of the surge wave and the bifurcation, and a bifurcated river channel generalization model.

3. The method of claim 1, wherein, The method for forming the initial surge wave shape by simplifying and fitting the surge wave before the surge wave is divided at the bifurcation comprises the following steps: simplifying and fitting the surge wave before the surge wave is divided at the bifurcation based on a wave shape equation to form the initial surge wave shape, wherein the wave shape equation comprises a Fourier series equation and a cosine wave equation.

4. The method of claim 1, wherein, The method for calculating the amplitude of the independent surge waves at the bifurcation at the time t1 before the independent surge waves collide and reflect during the division process based on the initial surge wave shape and the initial height of the independent surge waves comprises the following steps: according to Calculate the amplitude of the independent surge wave at the bifurcation point at time t1 before collision reflection during the flow splitting process. c ( i )( i = 1,2,3...), where, The initial height of each individual surge wave, x l This refers to the extent of the surge's reach in the direction of the diversion process. S 0 represents the riverbed slope, and g represents the acceleration due to gravity.

5. The method of claim 1, wherein, The method for calculating the movement and reflection path of each independent surge wave in the bifurcation comprises the following steps: defining the propagation direction of the surge wave in the original propagation river channel as the y direction, defining the propagation direction of the surge wave in the bifurcated river channel as the x direction, and establishing a calculation grid at the position of the bifurcation of the bifurcated river channel based on the width of the original propagation river channel and the width of the bifurcated river channel; calculating the movement and reflection path of each independent surge wave in the bifurcation based on the mirror reflection principle according to the calculation grid.

6. The method of claim 1, wherein, The method for determining the superimposed amplitude that meets the preset condition as the amplitude of the surge wave in the bifurcated river channel comprises the following steps: determining the maximum superimposed amplitude at the center line position in the bifurcated river channel as the amplitude of the surge wave in the bifurcated river channel.

7. The method according to claim 1 or 6, wherein The method further comprises the following steps: determining the influence range and damage degree of the surge wave on the river channel based on the determined amplitude of the surge wave in the bifurcated river channel.

8. A device for calculating the surge amplitude of a braided river channel, characterized in that, The method comprises the following steps: a first obtaining module for obtaining a surge wave before the surge wave is divided at a bifurcation of a bifurcated river channel; a forming module for forming an initial surge wave shape by simplifying and fitting the surge wave before the surge wave is divided at the bifurcation; The decomposition discrete module is configured to decompose the initial surge wave into preset independent surge waves with the same wavelength and different initial positions in a direction perpendicular to the original propagation direction of the surge wave; The first calculation module is configured to calculate the amplitude of the independent surge wave at the bifurcation position at a time t1 before the collision reflection based on the initial surge wave and the initial height of the independent surge wave; The second calculation module is configured to calculate the movement and reflection path of each independent surge wave in the bifurcation position, and obtain the position of the independent surge wave at a time t2 after the collision reflection, wherein the time t2 is located after the time t1 in time sequence; The second acquisition module is configured to superimpose the amplitudes of the independent surge waves at the same position at the time t2 to obtain the superimposed amplitude of the independent surge wave, wherein the amplitude of the independent surge wave at the time t2 is determined based on the amplitude of the independent surge wave at the time t1; The determination module is configured to determine the superimposed amplitude that meets the preset condition as the amplitude of the surge in the bifurcation river channel.

9. An electronic device, comprising: The processor, the memory, and the program stored in the memory and executable on the processor are provided, and the program is executed by the processor to implement the steps of the bifurcation river channel surge amplitude calculation method according to any one of claims 1 to 7. The computer program is stored in the computer readable storage medium, and the computer program is executed by the processor to implement the steps of the bifurcation river channel surge amplitude calculation method according to any one of claims 1 to 7.

10. A computer-readable storage medium, characterized in that, ​

Images