SWOT remote sensing data-based mountain torrent hydrological hydrodynamic model construction and parameter calibration method for data-deficient region

Abstract

The invention discloses an SWOT remote sensing data-based mountain torrent hydrological and hydrodynamic model construction and parameter calibration method for an area lacking data, and relates to the technical field of hydrological and hydrodynamic simulation and mountain torrent disaster early warning, and the method comprises the steps: calculating the observation water depth of each observation point in a target area according to SWOT satellite remote sensing data and soil static environment data; based on the cellular automaton, constructing an initial mountain torrent hydrological hydrodynamic model; constructing a microparameter learning framework based on the initial mountain torrent hydrological hydrodynamic model and a preset machine learning model; training the microparameter learning framework to optimize model parameters of a preset machine learning model, and determining corresponding optimal hydrological and hydrodynamic parameters; and determining a final mountain torrent hydrological hydrodynamic model according to the optimal hydrological hydrodynamic parameters and the initial mountain torrent hydrological hydrodynamic model. According to the method, efficient, accurate and automatic parameter calibration can be realized, so that the accuracy, timeliness and reliability of flood modeling are improved.

Description

Technical Field

[0001] This application relates to the technical fields of hydrological and hydrodynamic simulation and flash flood disaster warning, and particularly relates to a method for constructing a mountain flood hydrological and hydrodynamic model and calibrating parameters in data-deficient areas based on SWOT remote sensing data. Background Art

[0002] Flash flood disasters are natural disasters that seriously threaten the safety of lives and property in mountainous areas and along river areas, and are characterized by strong suddenness, great destructive power, and short duration. In the prediction and prevention of flash flood disasters, the construction of hydrological and hydrodynamic models is of crucial significance. However, the existing hydrological and hydrodynamic models are computationally complex and consume a large amount of resources, making it difficult to meet the rapid simulation requirements of large-scale regions. In the emergency situation of flash floods, traditional models often cannot provide real-time and accurate flash flood simulations. On the other hand, most of the existing hydrological model parameter calibration methods rely on traditional optimization algorithms, trial-and-error methods, or empirical formulas. These methods usually optimize model parameters based on single-point observation data and do not fully consider the spatial correlation of hydrological processes, resulting in low simulation accuracy for non-gauge areas. Moreover, traditional methods often rely on manual adjustment strategies, with a cumbersome optimization process and difficult to automate, and it is difficult to use large-scale satellite remote sensing data for global calibration.

[0003] In addition, obtaining hydrometeorological data also faces many challenges. Traditional hydrological stations are sparsely distributed and have limited observation coverage. Especially in mountainous and remote areas, due to complex terrain and environmental conditions, long-term continuous observation data is extremely lacking, it is difficult to obtain long-time series data for a single station, and the data missing problem is very serious. Traditional hydrological model parameter calibration methods usually rely on long-time series data of a single station and optimize the calibration by historical observation records of variables such as water level and flow velocity. However, in most non-gauge or data-deficient areas, since it is difficult to obtain continuous and stable long-time series data, traditional calibration methods are difficult to apply, resulting in a significant increase in the uncertainty of model parameters. In addition, even in some areas with observation data, the data records often have discontinuities, measurement errors, or insufficient time spans, further limiting the reliability of traditional calibration methods.

[0004] To make up for the above defects, data-deficient areas usually try to introduce parameter calibration methods based on machine learning. However, traditional optimization algorithms also require a large amount of computing resources, are prone to unstable calibration results, and it is difficult to ensure the rationality of parameters. The lack of these data directly limits the applicability and accuracy of traditional parameter calibration methods. Therefore, innovative technical means are urgently needed to address this challenge. Summary of the Invention

[0005] The objective of this application is to provide a method for constructing and parameter calibration of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data, which can achieve efficient, accurate, and automated parameter calibration, thereby improving the accuracy, timeliness, and reliability of flood modeling.

[0006] To achieve the above objective, this application provides the following solutions: This application provides a method for constructing and parameter calibration of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data, including: Obtain the soil static environment data, environmental driving sample data, and corresponding SWOT satellite remote sensing data of the target area; According to the SWOT satellite remote sensing data and the soil static environment data, calculate the observed water depth of each observation point in the target area; the soil static environment data, the environmental driving sample data, and the observed water depths of all corresponding observation points form a training sample, and multiple training samples form a training sample set; Based on the cellular automaton, construct an initial mountain flood hydrological and hydrodynamic model according to the soil static environment data; the initial mountain flood hydrological and hydrodynamic model includes hydrological and hydrodynamic parameters; Based on the initial mountain flood hydrological and hydrodynamic model and a preset machine learning model, construct a differentiable parameter learning framework; Use the training sample set to train the differentiable parameter learning framework to optimize the model parameters of the preset machine learning model and determine the corresponding optimal hydrological and hydrodynamic parameters; According to the optimal hydrological and hydrodynamic parameters and the initial mountain flood hydrological and hydrodynamic model, determine the final mountain flood hydrological and hydrodynamic model.

[0007] According to the specific embodiments provided in this application, the following technical effects are achieved: In this application, SWOT satellite remote sensing data is used for data processing. The high-resolution characteristics of SWOT satellite remote sensing data can be utilized to obtain accurate water depth and water surface change information in remote areas and mountainous areas where traditional observation means are difficult to cover, thus making up for the data missing problem. Based on cellular automata, an initial mountain flood hydrological and hydrodynamic model is constructed according to soil static environment data, and spatial distribution data such as water level and water depth within the entire basin can be obtained. These data can be highly matched with the SWOT satellite remote sensing observation data. The SWOT satellite remote sensing data obtains data on the basin surface, which is different from the time-series data at traditional observation stations. For this reason, this application constructs a differentiable parameter learning framework, combines machine learning algorithms, and uses the method of "exchanging space for time" to effectively calibrate the hydrological and hydrodynamic parameters in the initial mountain flood hydrological and hydrodynamic model with SWOT satellite remote sensing data. Moreover, through the differentiable parameter learning method, the repeated trial-and-error process in traditional model parameter calibration can be avoided, and the spatial similarity of the hydrological process is comprehensively considered, thereby improving the efficiency and accuracy of model parameter calibration.

[0008] In summary, this application uses cellular automata and a differentiable parameter learning framework to provide a new method for constructing a flood simulation model and automatic parameter calibration, which has higher efficiency, accuracy, and flexibility, improves the accuracy, timeliness, and reliability of flood modeling, and can be adapted to flood simulation in mountain flood areas lacking conventional monitoring data. BRIEF DESCRIPTION OF THE DRAWINGS

[0009] In order to more clearly illustrate the technical solutions in the embodiments of this application or in the prior art, the following will briefly introduce the drawings required to be used in the embodiments. Obviously, the drawings described below are only some embodiments of this application. For those of ordinary skill in the art, other drawings can be obtained based on these drawings without creative efforts.

[0010] Figure 1 It is a schematic flow chart of a method for constructing a mountain flood hydrological and hydrodynamic model and parameter calibration in data-deficient areas based on SWOT remote sensing data in an embodiment of this application.

[0011] Figure 2 It is a schematic diagram of grid runoff calculation based on infiltration capacity in an embodiment of this application.

[0012] Figure 3 It is a schematic diagram of the structure and parameter calibration of a differentiable parameter learning framework in an embodiment of this application.

[0013] Figure 4 It is a schematic diagram of the structure of a computer device provided in an embodiment of this application. DETAILED DESCRIPTION OF THE EMBODIMENTS

[0014] Next, in combination with the accompanying drawings in the embodiments of the present application, the technical solutions in the embodiments of the present application will be clearly and completely described. Obviously, the described embodiments are only a 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 those of ordinary skill in the art without creative efforts shall fall within the protection scope of the present application.

[0015] With the rapid development of remote sensing technology and artificial intelligence technology, it provides a new technical idea for flood modeling and model parameter calibration. These emerging technologies can not only solve the computational bottleneck in traditional methods, but also provide more and more accurate parameter calibration ideas. Model parameter calibration will no longer be limited by traditional manual parameter adjustment and data scarcity problems.

[0016] Based on this, the present application provides a technical solution by combining intelligent algorithms and high-resolution remote sensing data, realizing efficient, accurate, and automated parameter calibration, thereby greatly improving the accuracy, timeliness, and reliability of flood modeling. This will provide stronger technical support for flood prediction and disaster prevention and control.

[0017] To make the above objects, features, and advantages of the present application more obvious and understandable, the present application will be further described in detail below in conjunction with the accompanying drawings and specific embodiments.

[0018] In an exemplary embodiment, as Figure 1 shown, a method for constructing and parameter calibrating a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data is provided. This method can be executed by a computer device, specifically, it can be executed alone by a computer device such as a terminal or a server, or jointly executed by a terminal and a server. In the embodiments of the present application, the following steps 201 to step 206 are included.

[0019] Step 201, obtain the soil static environment data, environmental driving sample data, and corresponding SWOT satellite remote sensing data of the target area.

[0020] Among them, the soil static environment data includes DEM data, and may also include soil thickness data, land use data, and soil type data; the DEM data includes the actual terrain elevation of each observation point.

[0021] The SWOT satellite remote sensing data is obtained through the NASA official website. The corresponding SWOT remote sensing data includes features such as the water level height and range of large-scale water bodies. The obtained high-precision observation data can provide important support for the calibration of the subsequent mountain flood hydrological and hydrodynamic model. In an application example of the present application, the SWOT satellite remote sensing data includes the radar measurement elevation of each observation point (x, y).

[0022] Step 202: Calculate the observed water depth at each observation point in the target area according to the SWOT satellite remote sensing data and the soil static environment data; the soil static environment data, the environmental driving sample data, and the observed water depths of all corresponding observation points form a training sample, and multiple such training samples form a training sample set.

[0023] In an application example of the present application, the SWOT satellite can provide the water body elevation at each observation point in the basin area. Its basic principle is based on the radar-measured elevation plus a series of correction factors to eliminate the errors in radar measurement. Based on this, the calculation process of the observed water depth at each observation point in the target area includes: (21) For any observation point, calculate the observed water body elevation according to the radar-measured elevation, and the calculation formula is as follows: . Where is the radar-measured elevation; corrections is the correction data including earth curvature, atmospheric delay, orbital error, etc.; is the observed water body elevation of the SWOT satellite at time t.

[0024] (22) Calculate the difference between the observed water body elevation and the actual terrain elevation to obtain the observed water depth. That is, the inundation depth at each observation point is the difference between the water body elevation observed by the satellite and the actual terrain elevation of the observation point. The observed water depth at each observation point is: ; where is the actual terrain elevation of the observation point; is the inundation depth of the observation point, that is, the observed water depth.

[0025] After obtaining the observed water depths of all observation points, these points can be mapped to a regular grid through spatial interpolation, so as to obtain the spatial distribution of the entire water surface elevation, and then, in subsequent steps, compare with the predicted value output by the differentiable parameter learning framework to adjust the hydrological and hydrodynamic parameters in the initial mountain flood hydrological and hydrodynamic model to reduce the error.

[0026] Step 203: Based on Cellular Automaton (CA), construct an initial mountain flood hydrological and hydrodynamic model according to the soil static environment data; the initial mountain flood hydrological and hydrodynamic model includes hydrological and hydrodynamic parameters. Cellular Automaton, as the core tool for simulating the flood evolution process, divides the research area into two-dimensional grid cells (cells), and simulates the dynamic evolution process of floods in a highly flexible way. Each cell can independently carry geographical attributes (such as terrain elevation, water depth, flow velocity, etc.) and state change rules, so as to realize the integrated simulation of the entire hydrological and hydrodynamic process.

[0027] In an application example, step 203 includes the following (31)-(36).

[0028] (31) Use the HAND method to process the DEM data to obtain the terrain elevation data to be used. To optimize the hydrological and hydrodynamic simulation, the present application uses the HAND (Height Above Nearest Drainage) method to process the DEM data. The HAND method simplifies complex terrains by calculating the relative elevation of any point on the surface relative to the nearest neighboring river channel, can effectively identify the rainwater gathering areas, and provides the head information of the slope runoff process at the same time. Its mathematical definition is as follows: ; where H HAND is the HAND value of the current observation point, that is, the terrain elevation data to be used; Z C is the altitude of the current observation point; Z ND is the altitude of the river channel nearest to the current observation point. Here, "nearest neighbor" means the nearest along the runoff path, rather than the nearest in terms of straight-line distance.

[0029] The DEM data processed by the HAND method can accurately represent the water flow path and terrain features. Through the above processing, it not only helps to improve the flood inundation simulation and the accuracy of the basin runoff generation and concentration process, but also provides high-quality terrain data input for the cellular automaton model.

[0030] (32) Based on the terrain elevation data to be used, divide the target area into a cellular network; the cellular network includes multiple cells. Specifically, divide the target area into regular grids (cells) in meters, and assign a unique identifier or coordinates to each cell for dynamically recording the changes in its hydrological state.

[0031] (33) Define the state attributes corresponding to each cell; the state attributes include terrain elevation and soil characteristics; thus, key geographical attributes such as elevation are assigned to each cell. The state attributes corresponding to different target areas are different, and can also include initial conditions such as whether the terrain of the target area is steep, whether there is grass, and whether the soil is dry or wet.

[0032] (34) Use the Moore-type eight-neighborhood relationship to define the water flow direction and flow rules between the cells; specifically, use the flow rule driven by the height difference to simulate the water flowing from high to low.

[0033] (35) Set the initial conditions and boundary conditions of each cell; specifically, the initial conditions include the initial water depth, initial soil water content, outflow, etc., and the boundary conditions include drainage points, impermeable boundaries, etc.

[0034] The above settings provide an accurate spatial framework and physical basis for dynamically simulating the flood evolution process, that is, the water flow direction and flow rules in the cellular network, the state attributes of each cell, the initial conditions, and the boundary conditions constitute a two-dimensional cellular scenario; (36) Based on the two-dimensional cellular scenario, integrated runoff generation and concentration simulation are carried out to determine a set of water flow change functions; the two-dimensional cellular scenario and the set of water flow change functions constitute an initial mountain flood hydrological and hydrodynamic model; the set of water flow change functions is used to dynamically calculate the water flow direction, water flow velocity, water flow rate, and water depth of the cell; the hydrological and hydrodynamic parameters are determined by the set of water flow change functions.

[0035] Specifically, the cellular automaton adopts a grid-based spatial discretization method to divide the entire target area into interacting cells. Each cell can independently store and update the hydrological state and simulate the flood evolution process through local rules. Such discretization not only enables the runoff generated by rainfall (runoff generation process) to be precisely calculated in each cell but also spontaneously forms the concentration of water flow (concentration process) through the interaction between cells. Through integrated runoff generation and concentration calculation and analysis, the water flow direction, velocity, and water volume exchange are dynamically calculated to update the water level of each cell. Among them, the set of water flow change functions includes the cell soil infiltration capacity function, the inter-cell slope function, the water flow velocity and flow rate function, and the water balance function.

[0036] Rainfall is the main input of the surface hydrological cycle. A part of it will be lost, including infiltration, evaporation, vegetation interception, etc., and cannot be directly converted into surface runoff. The remaining part is called net rain, that is, the precipitation effectively participating in runoff generation. In this application, the infiltration capacity curve is used to calculate runoff generation to obtain the water depth of each cell, as Figure 2 shown, where the cell soil infiltration capacity function is: .

[0037] Among them, r is the soil infiltration rate, used to calculate the water depth of the cell; r f is the constant infiltration rate after soil saturation, that is, the stable infiltration rate; r0 is the infiltration rate at the initial rainfall, that is, the initial infiltration rate; β is the decay constant, indicating the rate at which the infiltration rate decreases with time; t is the time.

[0038] For the concentration part, by comparing the total water levels of the central cell and the neighboring cells, the path of water flow from the central cell to the low-potential cell can be determined. The inter-cell slope function is: .

[0039] Among them, is the slope from the central cell c to the neighboring cell n at time t, referring to the water flow direction; , are the water depths of the central cell c and the neighboring cell n at time t, respectively; , are the ground elevations of the central cell c and the neighboring cell n, respectively; is the horizontal distance from the central cell c to the neighboring cell n.

[0040] Manning's equation is an empirical formula used to calculate the flow velocity and discharge of water flow, which is widely applied in the fields of hydrology and hydraulics. By calculating the flow velocity using Manning's equation, the discharge can then be calculated. Therefore, the water flow velocity and discharge function is: .

[0041] .

[0042] Among them, is the water flow velocity from the central cell c to the neighboring cell n, is the Manning roughness coefficient; is the water flow discharge from the central cell c to the neighboring cell n; b is the cell width.

[0043] The water balance formula is the core formula in hydrology based on the law of conservation of mass. In the simulation of cellular automata, the water balance is reflected in that the water volume change in each cell is equal to the net change in the inflow water volume minus the outflow water volume, and the dynamic update of the water depth is calculated in combination with the area of the cell. The water balance function is: .

[0044] Among them, is the water depth of the i-th cell at time, is the water volume change duration, is the water depth of the i-th cell at time t; is the water flow discharge flowing into the central cell c; is the water flow discharge flowing out of the central cell c; A is the cell area.

[0045] Step 204: Based on the initial mountain flood hydrological and hydrodynamic model and the preset machine learning model, construct a differentiable parameter learning framework (dPL); among them, as a large-scale parameter learning method different from the traditional point-to-point parameter calibration method, the differentiable parameter learning framework will directly establish the global mapping relationship between the observed data and the model parameters, and at the same time use the model with physical mechanisms as the constraint conditions to achieve the effect of adaptive learning of the model parameters, as shown in Figure 3 .

[0046] The prediction of the initial mountain flood hydrological and hydrodynamic model based on cellular automata is related to external input conditions, static attributes, and parameters to be optimized. Therefore, the simulation values of cellular automata are dynamically updated based on the following formula: .

[0047] Among them, is the simulated water depth at the i-th cell at time t, which is the predicted value of the differentiable parameter learning framework; is the environmental driving sample data at the i-th cell at time t, which is generally an external driving variable related to time and location, such as meteorological forcing data (precipitation, temperature, etc.); is the static environmental data of the soil at the i-th cell, including static information related to spatial location such as terrain DEM, soil type, land cover, etc.; is the hydrological and hydrodynamic parameter of the i-th cell, which is an unobservable parameter specific to the location and needs to be determined separately through optimization. The hydrological and hydrodynamic parameters include the constant infiltration rate after soil saturation, the infiltration rate during initial rainfall, the decay constant, and the Manning roughness coefficient.

[0048] Compared with the traditional mountain flood hydrological and hydrodynamic model, in this application, the cellular automata uses local interaction rules to simulate the water flow evolution. Therefore, it is necessary to calibrate the parameters for each cell in detail to improve the simulation accuracy and stability. Among them, a large number of measured data required for calibration comes from the training sample set in step 202, and the parameters to be calibrated are hydrological and hydrodynamic parameters, as shown in Table 1 below.

[0049] Table 1

[0050] The predicted value of the differentiable parameter learning framework The relationship with the observed value can be expressed by a mapping function, and the function formula is as follows: .

[0051] Among them, the observed value is the observed water depth determined by SWOT satellite remote sensing data in the training sample set, is the mapping function, which is used to convert the predicted value of the model into the same unit or scale as the observed value; is the error between the predicted value and the observed value, and the closer it is to 0, the better.

[0052] Step 205: Use the training sample set to train the differentiable parameter learning framework to optimize the model parameters of the preset machine learning model and determine the corresponding optimal hydrological and hydrodynamic parameters. In an application example, step 205 includes the following steps (51)-(54).

[0053] (51) During a training iteration, the soil thickness data, the land use data, and the soil type data are input into the differentiable parameter learning framework, and the preset machine learning model generates predicted hydro - hydrodynamic parameters.

[0054] In a specific application, the preset machine learning model adopts an LSTM (Long Short - Term Memory) network. LSTM was initially developed for learning sequential data in the field of artificial intelligence and has now become a common choice for processing hydrological time - series data. Different from ordinary recurrent neural networks (such as RNN), LSTM has two states (cell state and hidden state) and three gating mechanisms (input gate, forget gate, output gate). By finding a set of parameters , such that the output of the physical model is as close as possible to the observed value . Therefore, by capturing the time - series characteristics through the LSTM network, mapping the dynamic input to the state update of the hydrological model, and back - propagating the error of the model output value to adjust the model parameters, automatic parameter calibration can be achieved.

[0055] (52) The predicted hydro - hydrodynamic parameters are input into the initial mountain flood hydro - hydrodynamic model in the differentiable parameter learning framework for parameter update, and then, driven by the environmental driving sample data, mountain flood hydro - dynamic simulation is carried out to obtain the simulated water depths at all observation points.

[0056] (53) When the preset training stop condition is not reached, based on the observed water depths and the simulated water depths at all observation points corresponding to the environmental driving sample data, a loss function is calculated, and the gradient is calculated using the loss function to optimize the model parameters of the preset machine learning model, and then the next iteration training is entered; different iteration trainings correspond to different training samples.

[0057] Among them, the loss function is used to measure the error between the model prediction value and the true observed value, and at the same time combines physical constraints to optimize the model performance to ensure that the model conforms to physical laws; the formula of the loss function is: .

[0058] Among them, L total is the value of the loss function; L data is the data - term loss value, which is used to quantify the error between the observed water depth and the simulated water depth; is the physical - constraint - term loss value, which is used to ensure mass conservation and momentum conservation; is the weight parameter, which is used to balance the priority of data fitting and physical consistency.

[0059] By adjusting the parameters , the model prediction is made to match the observed data most closely, thereby optimizing the model performance. Correspondingly, the functional formula for the data item loss value is: .

[0060] where t is the time, i is the index of the cell, is the simulated water depth at time t of the i-th cell, is the observed water depth, is the mapping function between the simulated water depth and the observed water depth; is the environmental driving sample data at time t of the i-th cell, is the static environmental data of the soil of the i-th cell, is the hydro-hydraulic parameter of the i-th cell; is the square of the Euclidean distance, measuring the difference between the predicted value and the observed value.

[0061] The functional formula for the loss value of the physical constraint term is: .

[0062] where q is the flow rate; y is the water depth; is the divergence of the flow rate, indicating the degree of change of the flow rate with space; is the rate of change of the water depth with time.

[0063] In practical applications, when performing gradient calculation, the weights of the neural network are optimized through gradient backpropagation, indirectly optimizing the parameters of the model. The gradient with respect to the parameters is expressed as: .

[0064] where is the gradient of the loss function L with respect to the parameter θ.

[0065] When updating the parameters, the gradient descent method is used to iteratively adjust the parameters: . Where is the learning rate, used to control the step size during each parameter update.

[0066] (54) When the preset training stop condition is reached, the predicted hydro-hydraulic parameters are marked as the optimal hydro-hydraulic parameters. That is, the differentiable parameter learning framework can be used to optimize the parameters. By repeating the training process and gradually reducing the value of the loss function until the objective function converges or reaches the preset threshold, the optimal hydro-hydraulic parameters are obtained.

[0067] In addition, the present application selects the Nash-Sutcliffe efficiency coefficient (NSE) as the main error measurement index. NSE measures the fitting degree of the model prediction value to the observed value, and the closer it is to 1, the better the model performance. The model performance is verified through the error index until the error converges. The calculation formula of NSE is as follows: .

[0068] where is the mean value of the observed values at all time steps t of the i-th cell.

[0069] Step 206: Determine the final mountain flood hydrodynamics model according to the optimal hydrographic and hydrodynamic parameters and the initial mountain flood hydrodynamics model.

[0070] Compared with the prior art, the present application has the following advantages: (1) The SWOT satellite is based on the interferometric synthetic aperture radar technology, breaking through the limitation of traditional satellites that only rely on point observations. It can provide large-scale and grid-based hydrographic data and accurately measure the height, width, and slope of floods in three dimensions. Compared with other satellites, SWOT has higher spatial and temporal resolutions and can obtain synchronous water level information over a wide area in a short time, forming a continuous planar data set. Therefore, the present application uses SWOT data to provide a new idea of "exchanging space for time". Although it cannot provide long-term observational data at a single point, through synchronous observations of a large area, relying on spatially distributed data to make up for the lack of time data.

[0071] (2) As a hydrographic and hydrodynamic simulation method, the cellular automaton has a grid-based computational unit in essence, which can be highly matched with grid-based remote sensing observation data (such as SWOT satellite data), enabling the model to provide a more refined hydrodynamic state simulation at different spatial positions and showing higher computational efficiency and flexibility when dealing with discrete dynamic hydrographic processes.

[0072] (3) In the process of parameter calibration of the cellular automaton, the core challenge of traditional methods lies in the need for complete grid-level data that matches it. However, the parameter calibration of traditional flood models is mainly based on single-point observational data such as water level and flow velocity at hydrological stations, and usually adopts overall or local adjustment methods, which are difficult to meet the refined calibration requirements of the CA model at the grid scale. In addition, due to the evolution mechanism of CA depending on the interaction between grids, simply relying on limited point observational data for interpolation or approximate optimization is difficult to accurately characterize the spatial heterogeneity of the flood process, resulting in limited calibration accuracy.

[0073] In response to this, the present application conducts simulations of the mountain flood hydrological and hydrodynamic processes based on the cellular automata model, uses the discrete characteristics of cellular automata to carry out flood process calculations, and dynamically updates the cell states to simulate the basin runoff generation and concentration processes under different hydrodynamic conditions. Among them, the differentiable parameter learning framework is mainly based on deep learning and gradient optimization mechanisms, and can optimize the parameters of all grid cells simultaneously on a global scale, thereby effectively improving the efficiency and accuracy of parameter calibration. Its core advantage lies in that the high-resolution gridded hydrological data provided by the SWOT satellite highly matches the grid cells of the CA model, enabling the model to directly calibrate using the observed data of each cell without relying on interpolation or approximate optimization based on limited point observations. In addition, the present application regards the model parameters as trainable variables, uses gradient descent optimization and combines with the backpropagation algorithm for dynamic adjustment, so that the model prediction results can match the SWOT observation data to the greatest extent, realizing more refined and efficient parameter optimization.

[0074] Considering the time-varying nature of the flood process, a long short-term memory network is further introduced to capture the spatio-temporal variation patterns of the parameters, ensuring that the calibrated parameters can adapt to the hydrodynamic characteristics under different time steps and different basin conditions. This method not only breaks through the limitation of traditional calibration methods relying on limited observation points, but also greatly improves the spatial resolution, time adaptability and computational efficiency of calibration, providing a more accurate and physically constrained solution for flood process simulation.

[0075] In summary, the present application is different from the traditional point-to-point parameter calibration method. It directly establishes a global mapping relationship between the observed data and the model parameters, and uses a model with a physical mechanism as a constraint condition to realize the adaptive adjustment of the parameters, thereby improving the accuracy and automation of the model. At the same time, it can also fully consider the similarity between geographical processes, that is, transfer the knowledge learned in one area to other areas, so as to effectively utilize large-scale observation data. The present application constructs a differentiable parameter learning framework, effectively making up for the computational bottleneck of the traditional model calibration method, realizing that the parameter calibration of the mountain flood model is no longer limited by the traditional manual parameter adjustment and data scarcity problems. At the same time, through the combination of intelligent algorithms and high-resolution remote sensing data, efficient, accurate and automated parameter calibration is achieved, thereby greatly improving the accuracy, timeliness and reliability of flood modeling, which will provide stronger technical support for flood prediction and disaster prevention and control.

[0076] In an exemplary embodiment, a computer device is provided. The computer device can be a server or a terminal, and its internal structure diagram can be as Figure 4As shown in the figure. The computer device includes a processor, a memory, an input / output interface (Input / Output, abbreviated as I / O), and a communication interface. Among them, the processor, the memory, and the input / output interface are connected through a system bus, and the communication interface is connected to the system bus through the input / output interface. Among them, the processor of the computer device is used to provide computing and control capabilities. The memory of the computer device includes a non-volatile storage medium and an internal memory. The non-volatile storage medium stores an operating system, a computer program, and a database. The internal memory provides an environment for the operation of the operating system and the computer program in the non-volatile storage medium. The input / output interface of the computer device is used to exchange information between the processor and external devices. The communication interface of the computer device is used to communicate with an external terminal through a network connection. When the computer program is executed by the processor, it realizes a method for constructing and parameter calibration of a mountain flood hydrological hydrodynamic model in data-deficient areas based on SWOT remote sensing data.

[0077] Those skilled in the art can understand that Figure 4 the structure shown in the figure is only a block diagram of some structures related to the solution of the present application, and does not constitute a limitation on the computer device to which the solution of the present application is applied. The specific computer device may include more or fewer components than those shown in the figure, or combine some components, or have different component arrangements.

[0078] In an exemplary embodiment, a computer device is provided, including a memory and a processor. A computer program is stored in the memory, and when the processor executes the computer program, the steps in the above method embodiments are implemented.

[0079] In an exemplary embodiment, a computer-readable storage medium is provided, storing a computer program, and when the computer program is executed by the processor, the steps in the above method embodiments are implemented.

[0080] In an exemplary embodiment, a computer program product is provided, including a computer program, and when the computer program is executed by the processor, the steps in the above method embodiments are implemented.

[0081] It should be noted that the user information (including but not limited to user device information, user personal information, etc.) and data (including but not limited to data for analysis, stored data, displayed data, etc.) involved in the present application are all information and data authorized by the user or fully authorized by all parties, and the collection, use, and processing of relevant data need to comply with relevant regulations.

[0082] Those of ordinary skill in the art can understand that all or part of the processes in the methods of the above embodiments can be completed by instructing relevant hardware through a computer program. The computer program can be stored in a non-volatile computer-readable storage medium. When the computer program is executed, it can include the processes of the embodiments of the above methods. Among them, any reference to a memory, database, or other medium used in the embodiments provided in the present application can include at least one of non-volatile and volatile memories. Non-volatile memories can include read-only memory (ROM), magnetic tapes, floppy disks, flash memories, optical memories, high-density embedded non-volatile memories, resistive random access memories (ReRAM), magnetoresistive random access memories (MRAM), ferroelectric random access memories (FRAM), phase change memories (PCM), graphene memories, etc. Volatile memories can include random access memory (RAM) or external cache memories, etc. By way of illustration and not limitation, RAM can be in various forms, such as static random access memory (SRAM) or dynamic random access memory (DRAM), etc.

[0083] The databases involved in the embodiments provided in the present application can include at least one of relational databases and non-relational databases. Non-relational databases can include distributed databases based on blockchain, etc., without limitation. The processors involved in the embodiments provided in the present application can be general-purpose processors, central processors, graphics processors, digital signal processors, programmable logics, data processing logics based on quantum computing, etc., without limitation.

[0084] The technical features of the above embodiments can be combined arbitrarily. For the sake of concise description, not all possible combinations of the technical features in the above embodiments are described. However, as long as there is no contradiction in the combination of these technical features, it should be considered to be within the scope described in this specification.

[0085] In this article, specific examples are used to elaborate on the principles and implementation manners of the present application. The description of the above embodiments is only used to help understand the method and its core idea of the present application; at the same time, for those of ordinary skill in the art, according to the idea of the present application, there will be changes in the specific implementation manners and application scopes. In summary, the content of this specification should not be construed as a limitation to the present application.

Claims

1. A method for constructing a mountain flood hydrological and hydrodynamic model and calibrating parameters in data-deficient areas based on SWOT remote sensing data, characterized in that, The method includes: Obtaining soil static environment data, environmental driving sample data, and corresponding SWOT satellite remote sensing data of the target area; Calculating the observed water depth of each observation point in the target area according to the SWOT satellite remote sensing data and the soil static environment data; the soil static environment data, the environmental driving sample data, and the observed water depths of all corresponding observation points form a training sample, and multiple such training samples form a training sample set; Based on the cellular automaton, constructing an initial mountain flood hydrological and hydrodynamic model according to the soil static environment data; the initial mountain flood hydrological and hydrodynamic model includes hydrological and hydrodynamic parameters; Based on the initial mountain flood hydrological and hydrodynamic model and a preset machine learning model, constructing a differentiable parameter learning framework; Using the training sample set to train the differentiable parameter learning framework to optimize the model parameters of the preset machine learning model and determine the corresponding optimal hydrological and hydrodynamic parameters; Determining the final mountain flood hydrological and hydrodynamic model according to the optimal hydrological and hydrodynamic parameters and the initial mountain flood hydrological and hydrodynamic model.

2. The method for constructing and parameter calibration of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data according to claim 1, characterized in that, The SWOT satellite remote sensing data includes the radar measurement elevation of each observation point; the soil static environment data includes DEM data, and the DEM data includes the actual terrain elevation of each observation point; The calculation process of the observed water depth of each observation point in the target area includes: For any observation point, calculating the observed water body elevation according to the radar measurement elevation; Calculating the difference between the observed water body elevation and the actual terrain elevation to obtain the observed water depth.

3. A method for constructing and calibrating parameters of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data according to claim 1, characterized in that, The soil static environment data includes DEM data; The construction process of the initial mountain flood hydrological and hydrodynamic model includes: Using the HAND method to process the DEM data to obtain the terrain elevation data to be used; Based on the terrain elevation data to be used, dividing the target area into a cellular network; the cellular network includes multiple cells; Defining the state attributes corresponding to each cell; the state attributes include terrain elevation and soil characteristics; Using the Moore-type eight-neighborhood relationship to define the water flow direction and flow rules between the cells; Setting the initial conditions and boundary conditions of each cell; the water flow direction and flow rules in the cellular network, the state attributes of each cell, the initial conditions, and the boundary conditions constitute a two-dimensional cellular scenario; Based on the two-dimensional cellular scenario, performing integrated runoff and confluence simulation to determine a set of water flow change functions; the two-dimensional cellular scenario and the set of water flow change functions constitute the initial mountain flood hydrological and hydrodynamic model; the set of water flow change functions is used to dynamically calculate the water flow direction, water flow velocity, water flow rate, and water depth of the cells; the hydrological and hydrodynamic parameters are determined by the set of water flow change functions.

4. A method for constructing and calibrating parameters of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data according to claim 3, characterized in that, The set of water flow change functions includes a cell soil infiltration capacity function, a cell-to-cell slope function, a water flow velocity and flow rate function, and a water balance function; The cell soil infiltration capacity function is: ； where r is the soil infiltration rate, which is used to calculate the water depth of the cell; r f is the constant infiltration rate after the soil is saturated; r0 is the infiltration rate at the initial rainfall; β is the decay constant; t is the time; The cell-to-cell slope function is: ； Among them, is the slope from the central cell c to the neighboring cell n at time t, referring to the water flow direction; and are the water depths of the central cell c and the neighboring cell n at time t, respectively; and are the ground elevations of the central cell c and the neighboring cell n, respectively; is the horizontal distance from the central cell c to the neighboring cell n; The water flow velocity and flow rate function is: ； ； Among them, is the water flow velocity from the central cell c to the neighboring cell n, is the Manning roughness coefficient; is the water flow rate from the central cell c to the neighboring cell n; b is the cell width. The water balance function is: ； Among them, is the water depth of the i-th cell at moment, is the water volume change duration, is the water depth of the i-th cell at time t; is the water flow rate flowing into the central cell c; is the water flow rate flowing out of the central cell c; A is the cell area.

5. A method for constructing and calibrating parameters of a mountain flood hydrological hydrodynamic model in data - scarce areas based on SWOT remote sensing data according to claim 4, characterized in that, The hydrological and hydrodynamic parameters include the constant infiltration rate after soil saturation, the infiltration rate at the initial rainfall, the decay constant, and the Manning roughness coefficient.

6. A method for constructing and calibrating parameters of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data according to claim 1, characterized in that, The differentiable parameter learning framework is expressed as: ； wherein, is the simulated water depth of the i-th cell at time t; is the environmental driving sample data of the i-th cell at time t; is the soil static environmental data of the i-th cell; is the hydrological and hydrodynamic parameter of the i-th cell.

7. A method for constructing and calibrating parameters of a mountain flood hydrological hydrodynamic model in data-deficient areas based on SWOT remote sensing data according to claim 1, characterized in that, The soil static environment data includes soil thickness data, land use data, and soil type data; Using the training sample set, train the differentiable parameter learning framework to optimize the model parameters of the preset machine learning model and determine the corresponding optimal hydrological and hydrodynamic parameters, including: In a training iteration process, input the soil thickness data, the land use data, and the soil type data into the differentiable parameter learning framework, and the preset machine learning model generates predicted hydrological and hydrodynamic parameters; Input the predicted hydrological and hydrodynamic parameters into the initial mountain flood hydrological and hydrodynamic model in the differentiable parameter learning framework for parameter update, and then, driven by the environmental driving sample data, perform mountain flood hydrological and hydrodynamic simulation to obtain the simulated water depths at all observation points; When the preset training stop condition is not reached, calculate the loss function based on the observed water depths at all observation points corresponding to the environmental driving sample data and the simulated water depths at all observation points, and use the loss function to calculate the gradient to optimize the model parameters of the preset machine learning model, and then enter the next iteration training; different iteration trainings correspond to different training samples; When the preset training stop condition is reached, mark the predicted hydrological and hydrodynamic parameters as the optimal hydrological and hydrodynamic parameters.

8. A method for constructing and calibrating parameters of a mountain flood hydrological hydrodynamic model in data - scarce areas based on SWOT remote sensing data according to claim 7, characterized in that, The formula of the loss function is: ； Among them, L total is the value of the loss function; L data is the data item loss value, which is used to quantify the error between the observed water depth and the simulated water depth; is the physical constraint item loss value, which is used to ensure mass conservation and momentum conservation; is the weight parameter.

9. A method for constructing and calibrating parameters of a mountain flood hydrological and hydrodynamic model in data-deficient areas based on SWOT remote sensing data according to claim 8, characterized in that, The function formula of the data item loss value is: ； where \(t\) is the time, \(i\) is the index of the cell, is the simulated water depth of the \(i\)-th cell at time \(t\), is the observed water depth, is the mapping function between the simulated water depth and the observed water depth; is the environmental driving sample data of the \(i\)-th cell at time \(t\), is the soil static environmental data of the \(i\)-th cell, is the hydro-hydraulic parameter of the \(i\)-th cell; is the square of the Euclidean distance.

10. A method for constructing and calibrating parameters of a mountain flood hydrological and hydrodynamic model in data - scarce areas based on SWOT remote sensing data according to claim 9, characterized in that, The function formula of the physical constraint item loss value is: ； where q is the flow rate; y is the water depth; is the divergence of the flow rate; is the rate of change of the water depth with time.

Images