CN-121659828-B - Dyke-breaking scene dynamic simulation method based on two-dimensional coupling model

CN121659828BCN 121659828 BCN121659828 BCN 121659828BCN-121659828-B

Abstract

The invention provides a dynamic simulation method of a breakwater scene based on a two-dimensional coupling model, which comprises the steps of firstly obtaining related data in a research area, constructing a one-dimensional river network model and a two-dimensional hydrodynamic model, coupling, setting simulation parameters and scene change trigger zone bits, separating the breakwater parameters into independent structures capable of being dynamically updated, distributing shared memory and defining an atomic access rule, determining time steps in simulation according to CFL conditions, advancing hydrodynamic calculation, detecting whether a user submits new breakwater parameters after each time step is finished, updating the parameters and triggering the zone bits if the user submits the new breakwater parameters, checking the zone bits by a main program in a non-blocking mode, correcting the breakwater elevation simulation breakwater after triggering, and continuously calculating reset zone bits. The method provided by the invention realizes the dynamic injection and thermal update of the breach scene, does not need restarting calculation, improves the timeliness of emergency decision, saves the computing resource, strengthens the decision support capability, and combines the simulation efficiency and the simulation precision.

Inventors

WANG WEIQI
QI XIN
GUO DANDAN
BI WUXIA
LIN WENQING
LU MINGDA
WANG FAN
ZHANG DAWEI

Assignees

中国水利水电科学研究院

Dates

Publication Date: 20260508
Application Date: 20251127

Claims (8)

1. A dyke-breaking scene dynamic simulation method based on a two-dimensional coupling model is characterized by comprising the following steps: step 1, acquiring data Acquiring a digital elevation model, river network section data, river network bank elevation data, land utilization type data and flood forecast data in a research area; step 2, constructing a coupling model Firstly, constructing a two-dimensional hydrodynamic model of a research area based on the digital elevation model and land utilization type data in the step 1; then, based on a spatial position relation, coupling the two-dimensional hydrodynamic model and the one-dimensional river network model in a lateral connection mode to complete construction of a two-dimensional coupling model of a research area; step 3, setting simulation parameters Setting simulation calculation parameters including calculation start time, calculation end time and output time step by taking flood forecast data obtained in the step 1 as inflow boundary conditions of a two-dimensional coupling model, setting an integer variable in a shared memory space as a scene change trigger flag bit, and setting an initial value as 0; Step 4, creating a crumple parameter structure In the initialization stage of the simulation system, the dike breach parameters are separated from a static data structure of a core algorithm, an independent 'breach parameter structure body' which can be updated dynamically is designed, a shared memory area which can be read and written by an interactive control program is allocated for the structure body, and an atomic access rule is defined for read-write operation; Step 5, hydrodynamic calculation and result output Judging and calculating the time step according to the flow field state of the current time t Calculating a main program based on flow field data at a time t, performing hydrodynamic deduction with a time step dt to obtain a flow field result at a time t+dt, updating a calculated flow field state according to the calculated result, and judging whether to output the result according to the output time step in the step 3, wherein the flow field state and the output result comprise water depth, x-direction flow velocity and y-direction flow velocity data of each grid unit in a two-dimensional hydrodynamic model, and water level and flow data of each river section in a one-dimensional river network model; Step6, the parameters of the crumple are updated and the flag bit is triggered After each calculation time step is finished, detecting whether a user sets a new crumple parameter through a graphical user interface or a command line tool and submits the new crumple parameter in the calculation time step process, if the user submits the new crumple parameter, receiving the crumple parameter by an interaction control program, overwriting the crumple parameter to a crumple parameter structure body in the step 4 through an atomicity operation, setting a scene change trigger flag bit in the step 3 to be 1, and if the user does not submit the crumple parameter, keeping the scene change trigger flag bit to be 0; step 7, scene change flag checking Checking the scene change trigger flag bit set in the step 3 in a non-blocking mode, if the value of the trigger flag bit is detected to be 1, immediately suspending the current calculation by the simulation calculation main program, and executing the step 8; step 8, reading parameters of the crumple After the simulation calculation main program obtains the access right of the shared memory, reading all parameter values in the updated crumple parameter structure in the step 6; step 9, embankment elevation correction and breach simulation The method comprises the steps of recording the center point coordinates of a burst opening according to a burst opening parameter structure body, determining the position of a river dike where the burst opening occurs according to the spatial position relation, correcting the elevation of the river dike in a one-dimensional river network model according to the width of the burst opening and the elevation information of the bottom of the burst opening, and simulating the burst opening state; step 10, bit zone reset and calculation recovery Resetting the trigger flag bit of scene change to 0, releasing the access right of the shared memory, and continuing to execute the step 5 to calculate the next time step; step 11, simulation termination judgment Judging whether the current simulation time is greater than the calculated termination time set in the step 3, returning to the step 5 to continue the time pushing cycle if the current time is not greater than the termination time, and executing the simulation program termination step if the current time is greater than the termination time; step 12, simulation program termination And releasing all the memories and terminating the simulation program.
2. The dynamic simulation method of a dyke-breaking situation according to claim 1, wherein in step 1, the acquired data are required to be a uniform geographic coordinate system and an elevation reference.
3. The dyke-breaking scene dynamic simulation method according to claim 1 is characterized in that in the step 2, a finite volume method is adopted to solve a two-dimensional shallow water equation set in the construction of the two-dimensional hydrodynamic model, a finite element method is adopted to solve a one-dimensional Saint Vinan equation set in the construction of the one-dimensional river network model, and a weir flow formula is adopted to calculate the exchange water volume between a one-dimensional river channel and a two-dimensional earth surface in the construction of the two-dimensional coupling model.
4. The dynamic simulation method of a dyke-breaking situation according to claim 1, wherein in the step 4, the parameter structure body of the dyke-breaking includes coordinates of a central point of the dyke-breaking, width of the dyke-breaking and elevation of the bottom of the dyke-breaking, and the atomic access rule adopts a mutual exclusion mechanism to protect a shared memory, so that only one process or thread can access the shared memory area at any moment, and mutual exclusivity and consistency of data access are ensured.
5. A method for dynamic simulation of a dyke-breaking scenario according to claim 1, characterized in that in step 5, the time step is The calculation formula of (2) is as follows: in the formula, As a safety factor, the safety factor of the device, 、 The space step length of the grid in the x direction and the y direction is respectively shown, u and v are the flow velocity in the x direction and the y direction respectively, g is the gravity acceleration; Is the water depth.
6. The dynamic simulation method of a breach scene according to claim 1, wherein in step 7, the non-blocking mode check means that after each calculation time step is finished, the main program judges whether a new breach scene needs to be processed by reading the value of the scene change trigger flag bit, without interrupting the current calculation flow.
7. The dynamic simulation method of a dike break situation according to claim 1, wherein the dike elevation correction and the dike break simulation in the step 9 specifically comprise the steps of modifying the elevation value of the dike position of the corresponding river channel into the elevation of the bottom of the dike break, and establishing the hydraulic connection between the one-dimensional river channel and the two-dimensional earth surface within the width range of the dike break.
8. The dynamic simulation method of a dike-break situation according to any one of claims 1 to 7, wherein the method supports dynamically introducing different dike-break situations multiple times in the simulation process, and repeatedly executing steps 6 to 10 each time a new dike-break situation is introduced.

Description

Dyke-breaking scene dynamic simulation method based on two-dimensional coupling model Technical Field The invention belongs to the technical field of flood simulation and flood control and disaster reduction, and particularly relates to a dyke-breaking scene dynamic simulation method based on a two-dimensional coupling model. Background Flood disasters are one of the most common and serious natural disasters worldwide, and particularly flood caused by embankment burst often causes huge casualties and economic losses. In order to effectively cope with flood risks, a flood evolution simulation technology becomes an important technical support means for flood control and disaster reduction work. At present, flood simulation is mainly realized by means of a hydrodynamic model, including a one-dimensional river network model and a two-dimensional hydrodynamic model. The one-dimensional river network model can efficiently simulate a river course flood evolution process by solving a Saint Vignan equation set, but is difficult to describe a complex flooding process of flood in a flooding area, and the two-dimensional hydrodynamic model can finely describe a flooding area flooding range and water depth distribution by solving a two-dimensional shallow water equation set, but has lower calculation efficiency. In order to achieve both calculation efficiency and simulation accuracy, researchers develop a two-dimensional coupling hydrodynamic model, generalize a river channel into a one-dimensional model, generalize a flooding area into a two-dimensional model, and realize full-river-basin flood-evolution simulation through coupling calculation. However, existing flood evolution simulation methods have significant limitations in handling burst scenarios such as dike breach. In the conventional simulation flow, key parameters such as time, spatial position, geometric form and the like of the crumple are usually used as initial conditions or static parameters, and must be preset before simulation calculation starts. In conventional models, once the simulation program is initiated, these crumple parameters are cured inside the model and are difficult to dynamically change during the calculation process. This "static preset" mode causes a number of inconveniences: Firstly, the flexibility is poor, and the timeliness is low. In a real emergency decision, a decision maker needs to quickly evaluate the difference of different breach schemes. The conventional method requires that the current simulation must be interrupted, the input file modified and then the calculation restarted from the beginning for each possible scenario, and the whole process is time-consuming and labor-consuming, and cannot meet the high requirement of emergency decision on timeliness. And secondly, the calculation resource is wasted. Flood evolution simulation, especially two-dimensional simulation for large areas, is huge in calculation amount and takes up to several hours or even days. The traditional method needs to restart calculation for each crumple scene, so that a large amount of completed calculation work is wasted, and serious calculation resource waste is caused. Thirdly, the decision support capability is limited. The assumption scene can not be dynamically introduced according to the flood real-time evolution situation in the simulation process, so that the potential of the simulation system as an interactive decision support tool is limited. For example, when the simulation shows that the pressure of a dike is great at a certain place, a decision maker wants to simulate the result of a break at the place immediately, and the prior art cannot realize the simulation in a 'response-on-time' mode. Therefore, how to overcome the defect of static preset of the breach scene in the traditional flood simulation method, and realize a flood simulation method which can dynamically, real-time and efficiently introduce and adjust the breach event in the simulation process, thereby providing stronger support for flood control emergency decision is a technical problem to be solved in the field. Disclosure of Invention Aiming at the problems that the existing flood simulation method is poor in flexibility, low in timeliness and incapable of dynamic interaction when the breach scene is processed, the invention aims to provide the dyke breach scene dynamic simulation method based on the two-dimensional coupling model, which can respond to a user instruction in real time and dynamically in the flood evolution simulation process, introduce or modify the flood evolution simulation of the dyke breach scene, realize 'thermal update' of simulation calculation, allow a user to dynamically inject new breach events when the simulation is carried out, and further remarkably improve the simulation analysis efficiency and the emergency decision support capability. The invention aims at realizing the following technical scheme: The invention provides a dyke-breaking scen