GIS-based one-dimensional and two-dimensional hydrodynamic model lateral coupling association method and system

Abstract

The invention provides a GIS-based one-dimensional and two-dimensional hydrodynamic model lateral coupling association method and system, and relates to the technical field of hydrodynamic models.The method comprises the steps that a river channel water area surface, a river channel center line and a two-dimensional land area boundary vector file are collected; generating a two-dimensional partition boundary vector file based on the river water area surface, the river center line and the two-dimensional land boundary vector file, and generating a two-dimensional grid file consistent with the range of the two-dimensional partition boundary vector file through a grid generator; vectorizing the two-dimensional grid file to obtain a grid node file; extracting two-dimensional boundary line segmentation points; extracting a two-dimensional grid, a grid node, a two-dimensional boundary and a boundary grid node of each two-dimensional partition based on the two-dimensional boundary line segmentation point and the grid node file, and establishing a coupling relationship between the grid node and the boundary in each two-dimensional partition; and integrating the coupling relationship between the grid nodes and the boundaries in each two-dimensional partition to form a global coupling relationship.

Description

Technical Field

[0001] The present invention relates to the technical field of hydrodynamic models, and in particular to a GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method and system. Background Art

[0002] In hydrodynamic simulation and flood management, accurately simulating the dynamic exchange of water bodies across different spatial scales is crucial for accurate flood forecasting and prevention. Coupling one-dimensional and two-dimensional hydrodynamic models is a key approach for simulating complex water systems, such as river flooding, urban waterlogging, and dam failures. By effectively coupling one-dimensional river models with two-dimensional surface models, efficient characterization of the hydrodynamic interactions between linear water bodies (e.g., river channels) and planar water bodies (e.g., depressions and floodplains) within a watershed can be achieved.

[0003] Currently, mainstream one-dimensional and two-dimensional hydrodynamic coupling simulation platforms, such as the MIKE series, generally use lateral links to connect one-dimensional river channels with two-dimensional surface grids. In the MIKE series in particular, modelers often need to manually delineate the connection locations between the one-dimensional and two-dimensional models, defining "link structures" or "boundary conditions" to complete the construction of data interfaces and boundary information. Furthermore, current mainstream methods mostly rely on data such as topographic maps, river cross-sections, and remote sensing layers, requiring manual interpretation and demarcation of link locations to achieve the transmission of hydrodynamic information between river networks and surface areas.

[0004] However, existing methods rely on extensive manual work, which is not only inefficient but also prone to errors. This is especially complex and time-consuming in large-scale watershed models. Furthermore, manually defined coupling locations struggle to accurately reflect topographic variations, particularly in areas with sudden changes in topography, such as culverts and sluice gates. Furthermore, rigid connection mechanisms can easily lead to water imbalances and momentum transfer errors in highly non-uniform flow fields. Summary of the Invention

[0005] In view of the above shortcomings of the existing technology, the purpose of the embodiments of the present invention is to provide a GIS-based lateral coupling and association method for one- and two-dimensional hydrodynamic models. This method can address the technical issues that existing methods rely on extensive manual operations, which are not only inefficient but also prone to errors. This is particularly complex and time-consuming in large-scale watershed models. Furthermore, manually defined coupling locations fail to accurately reflect topographic changes, especially in areas of sudden topographic changes such as culverts and sluice gates. Furthermore, the rigid connection mechanism can easily cause water imbalances and momentum transfer errors in strongly non-uniform flow fields.

[0006] In a first aspect of an embodiment of the present invention, a GIS-based one- and two-dimensional hydrodynamic model lateral coupling correlation method is proposed, comprising:

[0007] S1: Collect river water area vector files, river centerline vector files and 2D land boundary vector files;

[0008] S2: Based on the river water area vector file, the river centerline vector file, and the two-dimensional land boundary vector file, a two-dimensional partition boundary vector file is generated, and a two-dimensional grid file is generated by a grid generator with a range consistent with the two-dimensional partition boundary vector file;

[0009] S3: vectorizing the two-dimensional grid file to obtain a grid node file;

[0010] S4: extracting two-dimensional boundary line segmentation points from the river water area vector file and the river centerline vector file;

[0011] S5: extracting the two-dimensional grid, grid nodes, two-dimensional boundary and boundary grid nodes of each of the two-dimensional partitions based on the two-dimensional boundary line segmentation points and the grid node file, and establishing a coupling relationship between the grid nodes and the boundary in each of the two-dimensional partitions;

[0012] S6: Integrate the coupling relationship between the grid nodes and the boundaries in each of the two-dimensional partitions to form a global coupling relationship.

[0013] A second aspect of an embodiment of the present invention provides a GIS-based one- and two-dimensional hydrodynamic model lateral coupling association system, comprising: a processor and a memory;

[0014] The memory stores programs or instructions that can be run on the processor, and when the programs or instructions are executed by the processor, the steps of the GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method as described in the first aspect are implemented.

[0015] According to a third aspect of an embodiment of the present invention, a readable storage medium is proposed, on which a program or instruction is stored. When the program or instruction is executed by a processor, the steps of the GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method as described in the first aspect are implemented.

[0016] The beneficial effects brought about by the technical solution provided by the embodiment of the present invention include at least:

[0017] In an embodiment of the present invention, by automatically processing the river water surface vector file, the river centerline vector file and the two-dimensional land boundary vector file, the system can automatically convert the two-dimensional grid and grid node file, eliminating the need for manual intervention, significantly improving data processing efficiency, and avoiding the inefficiency and errors caused by manual operation. At the same time, by automatically extracting the two-dimensional boundary line segmentation points in the river water surface and the river centerline, it can accurately reflect the terrain mutation area (such as culverts, gates, etc.), ensure the accuracy of the coupling position, and flexibly adapt to terrain changes. This method effectively avoids the water imbalance and momentum error that may be caused by the rigid connection mechanism in a strong non-uniform flow field, ensuring the efficiency and accuracy of the coupling process. BRIEF DESCRIPTION OF THE DRAWINGS

[0018] The accompanying drawings are only for the purpose of illustrating specific embodiments and are not to be considered as limiting the present invention. Throughout the drawings, the same reference symbols represent the same components. Obviously, the drawings described below are only some embodiments of the present invention. It is clear that those skilled in the art can derive other drawings based on these drawings without inventive effort.

[0019] Figure 1 1. It is a flow chart of a GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method provided by an embodiment of the present invention;

[0020] Figure 2 1. It is a flow chart of another GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method provided by an embodiment of the present invention;

[0021] Figure 3 The vectorized two-dimensional grid and grid node graph provided by the embodiment of the present invention;

[0022] Figure 4 This is a coupling association relationship output file format diagram provided by an embodiment of the present invention;

[0023] Figure 5 This is a comparison diagram of a dike and riverbank before and after treatment provided by an embodiment of the present invention;

[0024] Figure 6 It is a structural diagram of a lateral coupling association system of a one- and two-dimensional hydrodynamic model based on GIS provided in an embodiment of the present invention. DETAILED DESCRIPTION

[0025] In order to enable those skilled in the art to better understand the technical solutions in the embodiments of the present invention, the technical solutions of the present invention will be clearly and completely described below with reference to the accompanying drawings. Obviously, the embodiments described are part of the embodiments of the present invention, rather than all of the embodiments. It should be understood that these descriptions are merely exemplary and are not intended to limit the scope of the present invention. Based on the embodiments of the present invention, all other embodiments obtained by those of ordinary skill in the art without creative work should fall within the scope of protection of the present invention.

[0026] The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method provided by the embodiment of the present invention is described in detail below with reference to the accompanying drawings through specific embodiments and application scenarios.

[0027] Reference Manual Figure 1 , shows a flow chart of a GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method provided in an embodiment of the present invention.

[0028] Reference Manual Figure 2 A flow chart of another GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method provided by an embodiment of the present invention is shown.

[0029] Reference Manual Figure 3 A vectorized two-dimensional grid and a grid node graph provided by an embodiment of the present invention are shown.

[0030] Reference Manual Figure 4 A diagram showing a coupling association relationship output file format provided by an embodiment of the present invention is shown.

[0031] Reference Manual Figure 5 A comparison diagram of a levee and river bank before and after treatment provided by an embodiment of the present invention is shown.

[0032] The embodiment of the present invention provides a GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method, which may include the following steps:

[0033] S1: Collect river water area vector files, river centerline vector files, and two-dimensional partitioned land area vector files.

[0034] A river water surface vector file is a vector data file used to represent a river water surface (i.e., the water boundary area of ​​a river) in a GIS (Geographic Information System). This file typically contains information about the geometry, boundary location, and spatial characteristics of the river water area.

[0035] A river centerline vector file is a geographic data file used to represent the centerline or main flow of a river. In a GIS (Geographic Information System), river centerline vector files are typically stored as line segments, representing the main channel of a river or the direction of water flow.

[0036] The 2D Land Boundary Vector File is used to represent the land boundary data of the study area without surrounding rivers in a 2D hydrodynamic model. In a Geographic Information System (GIS), 2D land boundary vector files are typically stored as line segments, representing the land boundary of the study area without surrounding rivers.

[0037] In an embodiment of the present invention, the river water surface vector file, the river centerline vector file, and the two-dimensional land boundary vector file provide different types of spatial data required for the hydrodynamic model, so that the model can comprehensively describe the river and its surrounding area. By integrating these three different types of data, the direction, speed, flow rate, and water level changes of the water flow can be accurately simulated, thereby obtaining more reliable simulation results.

[0038] S2: Generate a two-dimensional partition boundary vector file based on the river water area vector file, the river centerline vector file and the two-dimensional land boundary vector file, and generate a two-dimensional grid file consistent with the range of the two-dimensional partition boundary vector file through the grid generator.

[0039] The Grid Generator is a software tool used to create computational grids in hydrodynamic models. It is primarily used to divide the study area into multiple small computational cells (grids). These grids can be used to perform numerical calculations in the model to simulate physical phenomena such as water flow, pollution diffusion, and heat transfer.

[0040] In a possible implementation, S2 specifically includes:

[0041] S201: Generate a partition surface file based on the river centerline vector file and the two-dimensional land boundary vector file.

[0042] Partition surface files are commonly used in hydrology, environmental simulation, geographic information systems (GIS), and hydrodynamic modeling to describe regional divisions or partition boundaries. Their primary function is to divide the model area into multiple subregions or "partitions," each representing an independent computational unit.

[0043] S202: Eliminate the river water area surface vector file from the partition surface file to obtain a two-dimensional partition boundary vector file.

[0044] S203: Based on the two-dimensional partition boundary vector file, generate a two-dimensional grid file consistent with the range of the two-dimensional partition boundary vector file through SMS grid generation software.

[0045] Among them, the Surface Water Modeling System (SMS) is a professional software widely used in the fields of hydrology and hydrodynamics, particularly for two-dimensional hydrodynamic modeling and mesh generation. SMS helps users build and manage complex hydrological and hydrodynamic models, especially those involving surface water bodies such as rivers, lakes, and oceans. By integrating multiple mesh generation methods, computational techniques, and post-processing tools, it provides precise solutions for problems such as water flow simulation, flood prediction, and pollutant dispersion.

[0046] In this embodiment of the present invention, by combining river water surface vector files, river centerline vector files, and two-dimensional land boundary vector files with SMS grid generation software, hydrological and hydrodynamic modeling can be efficiently and accurately performed, particularly for flood simulation, watershed management, and pollutant dispersion. This approach improves simulation accuracy, optimizes computational efficiency, and enhances flexibility and adaptability, providing an important scientific basis for disaster warning, environmental management, and decision support.

[0047] S3: Vectorize the two-dimensional grid file to obtain a grid node file.

[0048] The grid node file is a crucial data file in hydrological and hydrodynamic models. It contains detailed data describing each node in the computational grid. The grid node file stores the spatial location (coordinates), node number, and associated physical properties (such as elevation and flow velocity) of each grid node. These nodes are the fundamental units for numerical calculations.

[0049] In a possible implementation, S3 specifically includes:

[0050] S301: Reading node information and triangular mesh information from a two-dimensional mesh file, wherein the triangular mesh information includes the mesh name, and the node information includes the node name, node x-coordinate, node y-coordinate, node elevation, and node serial number.

[0051] S302: Arrange the node information and triangular mesh information into an Excel table.

[0052] S303: Use GIS software to convert the node coordinate information in the Excel table into a point vector file.

[0053] S304: De-duplication processing is performed on nodes with duplicate coordinates in the point vector file.

[0054] S305: Saving the point vector file after deduplication processing as a grid node file.

[0055] In this embodiment of the present invention, vectorizing a two-dimensional grid file and converting it into a grid node file helps improve data accuracy, efficiency, and compatibility. Furthermore, by precisely locating nodes, removing duplicate data, and standardizing storage formats, the accuracy and efficiency of hydrodynamic model calculations can be improved, and subsequent analysis, visualization, and data sharing can be streamlined.

[0056] S4: Extract two-dimensional boundary line segmentation points from the river water area vector file and the river centerline vector file.

[0057] 2D boundary line segmentation points (2D Boundary Line Segmentation Points) play an important role in hydrological and hydrodynamic models, particularly in simulating water flow and boundary conditions in rivers, lakes, reservoirs, and other areas. They are key elements in mesh generation, model calculations, and boundary condition definition, describing the geometry and location of boundaries within the model area.

[0058] In a possible implementation, S4 specifically includes:

[0059] S401: Convert the river water area vector file into the river water area shoreline, and remove the intersecting parts of the river water area shoreline to obtain the river embankment line.

[0060] S402: extracting the starting point and the end point of the river embankment line, and deleting the duplicated position points to obtain the first type of two-dimensional boundary line segmentation points.

[0061] S403: Calculate the spatial intersection of the river centerline vector file and the river water area vector file, and extract the intersection point feature.

[0062] S404: Delete duplicate location points in the intersection point elements to obtain the second type of two-dimensional boundary line segmentation points.

[0063] S405: Integrate the first type of two-dimensional boundary line segmentation points and the second type of two-dimensional boundary line segmentation points to form two-dimensional boundary line segmentation points.

[0064] In this embodiment, by extracting and integrating two-dimensional boundary segmentation points from river water area vector files and river centerline vector files, the accuracy, computational efficiency, and flexibility of hydrological and hydrodynamic models can be significantly improved. This eliminates duplicate data, precisely delineates boundaries, and defines boundary conditions, providing a solid foundation for subsequent grid generation, model calculations, disaster prediction, and decision support.

[0065] S5: Based on the 2D boundary line segmentation points and mesh node files, extract the 2D meshes, mesh nodes, 2D boundaries, and boundary mesh nodes of each 2D partition, and establish the coupling relationship between the mesh nodes and boundaries in each 2D partition.

[0066] In a possible implementation, S5 specifically includes:

[0067] S501: Extracting two-dimensional partition range vector data of the current two-dimensional partition.

[0068] S502: Based on the two-dimensional partition range, extract the river embankment line of the current two-dimensional partition, and eliminate the river embankment line within the 50m buffer range where the starting and ending points of the river embankment line are located to obtain the main river embankment line.

[0069] S503: Count the minimum side lengths of the two-dimensional grids in the current two-dimensional partition, generate a buffer zone inward with a buffer distance of 0.5 times the minimum side length, and extract two-dimensional grid cells associated with the current two-dimensional partition based on the generated buffer zone.

[0070] S504: Extracting grid nodes corresponding to the current two-dimensional partition according to the two-dimensional grid cells.

[0071] S505: Combining the two-dimensional grid cells, the two-dimensional boundary line segmentation points, and the main levee line of the river channel of the current two-dimensional partition, identifying and extracting the two-dimensional boundary of the current two-dimensional partition.

[0072] In a possible implementation, S505 specifically includes:

[0073] S5051: Fusing the two-dimensional grid cells of the current two-dimensional partition to generate a two-dimensional grid range surface, and extracting the corresponding two-dimensional grid range line based on the two-dimensional grid range surface.

[0074] S5052: The average grid side length calculated from the two-dimensional grid file is set as the tolerance value for boundary recognition, and the tolerance value is used as the search radius. The grid node with the smallest distance between the two-dimensional boundary line segmentation point and the grid node within the tolerance range is taken as the feature grid node.

[0075] S5053: Segment the two-dimensional grid range line according to the characteristic grid nodes to obtain multiple two-dimensional boundaries within the current two-dimensional partition.

[0076] S5054: Establish a one-to-one correspondence between each 2D boundary and the main embankment line of the river. If the direction of the 2D boundary is inconsistent with the direction of the main embankment line of the river, adjust the direction of the 2D boundary to make it consistent with the direction of the embankment line.

[0077] S5055: Label each 2D boundary. If a 2D boundary corresponds to a main river dike line, the current 2D boundary is marked as a lateral coupling boundary, and its number is consistent with the sequence number of the main river dike line. Otherwise, the current 2D boundary is marked as a land boundary and numbered using the LDBJ+sequence number rule.

[0078] In this embodiment of the present invention, the model dynamically adapts to varying terrain or structural changes by adjusting when the two-dimensional boundary misaligns with the river levee line. This flexibility enables the model to automatically update and adapt to changing circumstances, such as levee repairs or changes in river morphology. Furthermore, by calculating the average edge length of the grid and setting a tolerance for boundary recognition, the model automatically optimizes the division of grid cells, ensuring accurate boundaries without wasting computing resources.

[0079] S506: According to the many-to-many spatial relationship, the associated grid nodes of each two-dimensional boundary are extracted, and the associated grid nodes are associated with the boundary codes.

[0080] S507: Calculate the mileage value of each two-dimensional boundary grid node on the boundary line.

[0081] S508: Sort the boundary grid nodes in sequence according to the boundary name and mileage value, and complete the construction of the grid node and boundary coupling relationship of the current two-dimensional partition.

[0082] In this embodiment of the present invention, by extracting boundary grid nodes and associating them with the mileage values ​​along the boundary line, the model accurately defines the interaction between water flow and the boundary. Furthermore, associating grid nodes with boundary codes and sorting them by mileage ensures that boundary conditions (such as water level and flow velocity) are correctly applied to each node, avoiding calculation errors caused by inaccurate application of boundary conditions.

[0083] S509: Repeat S501 to S508 until the coupling relationship between the grid nodes and boundaries of all partitions is established.

[0084] In this embodiment, the accuracy, computational efficiency, and flexibility of hydrological and hydrodynamic models are significantly improved by precisely extracting and coupling 2D boundary segmentation points, grid nodes, 2D boundaries, and major river dike lines. This process ensures the accurate application of boundary conditions, optimizes the computational workflow, and provides strong support for large-scale, multi-region simulations.

[0085] S6: Integrate the coupling relationship between the grid nodes and the boundaries in each two-dimensional partition to form a global coupling relationship. In one possible implementation, S6 specifically includes:

[0086] S601: Merge the two-dimensional boundaries of each two-dimensional partition, and renumber the land boundaries that do not correspond to the main river embankment line.

[0087] S602: Merge the boundary grid nodes after sorting in each two-dimensional partition, and adjust the boundary number information of the corresponding nodes according to the unified land boundary number, complete the consistency matching of the whole domain boundary, and form a whole domain coupling relationship.

[0088] In this embodiment, by merging and renumbering the boundaries of each two-dimensional partition, the boundaries of different regions are ensured to be consistent across the entire model, avoiding simulation errors caused by inconsistent boundary definitions or non-standard numbering. Furthermore, by standardizing the numbering and adjusting the boundary nodes, the model's data structure becomes more standardized, facilitating subsequent data management, querying, and visualization.

[0089] In a possible implementation manner, after S6, the following steps are further included:

[0090] S7: Write the global coupling relationship into a text file, and output the text file as an inp format file and a dat format file.

[0091] The inp file format is a commonly used input data file format for hydrological and hydrodynamic models. In hydrodynamic simulation software, inp files are typically used to store all required model input data, such as model parameters, boundary conditions, initial conditions, and mesh data. They serve as the foundation for model operation and contain detailed information about model settings and environmental conditions.

[0092] DAT files are commonly used as data output files in hydrological and hydrodynamic models. After a model is run, the results are typically saved as DAT files, which store model calculations such as flow velocity, flood levels, pollutant concentrations, and sediment movement. These files contain temporal or spatially distributed data from the model simulation.

[0093] In a possible implementation, S7 specifically includes:

[0094] S701: Writing node information including node name, node x-coordinate, node y-coordinate, node elevation and node serial number in the two-dimensional grid file into a text file line by line.

[0095] S702: Writing the mesh names in the triangular mesh information into a text file line by line.

[0096] S703: Write the two-dimensional boundary and associated node information into a text file.

[0097] S704: Output the text file after writing the node information, grid name, two-dimensional boundary and associated node information into an inp format file and a dat format file.

[0098] In an embodiment of the present invention, node information, grid information, and boundary information are written separately and line by line into a text file, which can clearly distinguish between various types of data and facilitate subsequent data processing and analysis. Users can easily extract specific information (such as node location, grid type, boundary conditions, etc.) from the file and perform related operations. At the same time, by outputting standard dat files, it is easy to simulate and compare multiple scenarios, such as scenarios with different precipitation amounts and different river channel designs, and analyze the changes in water level, flow rate, etc. in each scenario.

[0099] The beneficial effects brought about by the technical solution provided by the embodiment of the present invention include at least:

[0100] In an embodiment of the present invention, by automatically processing the river water surface vector file, the river centerline vector file and the two-dimensional land boundary vector file, the system can automatically convert the two-dimensional grid and grid node file, eliminating the need for manual intervention, significantly improving data processing efficiency, and avoiding the inefficiency and errors caused by manual operation. At the same time, by automatically extracting the two-dimensional boundary line segmentation points in the river water surface and the river centerline, it can accurately reflect the terrain mutation area (such as culverts, gates, etc.), ensure the accuracy of the coupling position, and flexibly adapt to terrain changes. This method effectively avoids the water imbalance and momentum error that may be caused by the rigid connection mechanism in a strong non-uniform flow field, ensuring the efficiency and accuracy of the coupling process.

[0101] Reference Manual Figure 6 , shows a structural schematic diagram of a one- and two-dimensional hydrodynamic model lateral coupling association system based on GIS provided in an embodiment of the present invention.

[0102] The embodiment of the present invention provides a GIS-based one- and two-dimensional hydrodynamic model lateral coupling association system 20, comprising: a processor 201 and a memory 202;

[0103] The memory 202 stores programs or instructions that can be run on the processor 201. When the programs or instructions are executed by the processor 201, the steps of the above-mentioned GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method are implemented, and the same technical effect can be achieved. To avoid repetition, the present invention will not be repeated.

[0104] It should be understood that the processor 201 in the embodiment of the present invention may be a central processing unit (CPU), or may be other general-purpose processors, digital signal processors (DSP), application-specific integrated circuits (ASIC), field programmable gate arrays (FPGA), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. A general-purpose processor may be a microprocessor or any conventional processor, etc.

[0105] It should also be understood that the memory 202 in the embodiment of the present invention can be a volatile memory or a non-volatile memory, or can include both volatile and non-volatile memories. Among them, the non-volatile memory can be a read-only memory (ROM), a programmable read-only memory (PROM), an erasable programmable read-only memory (EPROM), an electrically erasable programmable read-only memory (EEPROM), or a flash memory. The volatile memory can be a random access memory (RAM), which is used as an external cache. By way of example and not limitation, many forms of random access memory (RAM) are available, such as static RAM (SRAM), dynamic random access memory (DRAM), synchronous DRAM (SDRAM), double data rate synchronous dynamic random access memory (DDR SDRAM), enhanced synchronous dynamic random access memory (ESDRAM), synchronous link DRAM (SLDRAM), and direct rambus RAM (DR RAM).

[0106] The above embodiments can be implemented in whole or in part through software, hardware (such as circuits), firmware, or any other combination. When implemented using software, the above embodiments can be implemented in whole or in part in the form of a computer program product. The computer program product includes one or more computer instructions or computer programs. When the computer instructions or computer program are loaded or executed on a computer, the processes or functions described in accordance with the embodiments of the present invention are generated in whole or in part. The computer can be a general-purpose computer, a special-purpose computer, a computer network, or other programmable device. The computer instructions can be stored in a computer-readable storage medium or transmitted from one computer-readable storage medium to another. For example, the computer instructions can be transmitted from one website, computer, server, or data center to another website, computer, server, or data center via a wired method (such as infrared, wireless, microwave, etc.). The computer-readable storage medium can be any available medium that can be accessed by a computer or a data storage device such as a server or data center that contains a collection of one or more available media. The available medium can be a magnetic medium (such as a floppy disk, hard disk, or magnetic tape), an optical medium (such as a DVD), or a semiconductor medium. The semiconductor medium can be a solid-state drive.

[0107] It should be understood that in various embodiments of the present invention, the size of the serial numbers of the above-mentioned processes does not mean the order of execution. The execution order of each process should be determined by its function and internal logic, and should not constitute any limitation on the implementation process of the embodiments of the present invention.

[0108] Those skilled in the art will appreciate that the units and algorithm steps of each example described in conjunction with the embodiments disclosed herein can be implemented in electronic hardware, or a combination of computer software and electronic hardware. Whether these functions are performed in hardware or software depends on the specific application and design constraints of the technical solution. Professionals and technicians can use different methods to implement the described functions for each specific application, but such implementation should not be considered beyond the scope of the present invention.

[0109] Those skilled in the art will clearly understand that, for the convenience and brevity of description, the specific working processes of the above-described equipment, devices and units can refer to the corresponding processes in the aforementioned method embodiments and will not be repeated here.

[0110] In the several embodiments provided by the present invention, it should be understood that the disclosed devices, apparatuses and methods can be implemented in other ways. For example, the device embodiments described above are merely illustrative. For example, the division of the units is merely a logical function division. In actual implementation, there may be other division methods, such as multiple units or components can be combined or integrated into another device, or some features can be ignored or not executed. Another point is that the mutual coupling or direct coupling or communication connection shown or discussed can be through some interface, indirect coupling or communication connection of the device or unit, which can be electrical, mechanical or other forms.

[0111] The units described as separate components may or may not be physically separate, and the components shown as units may or may not be physical units, that is, they may be located in one place or distributed across multiple network units. Some or all of these units may be selected to achieve the purpose of this embodiment according to actual needs.

[0112] In addition, each functional unit in each embodiment of the present invention may be integrated into one processing unit, or each unit may exist physically separately, or two or more units may be integrated into one unit.

[0113] If the functions are implemented in the form of software functional units and sold or used as independent products, they can be stored in a computer-readable storage medium. Based on this understanding, the technical solution of the present invention, or the part that contributes to the prior art, or the part of the technical solution, can be embodied in the form of a software product. The computer software product is stored in a storage medium and includes several instructions for enabling a computer device (which can be a personal computer, server, or network device, etc.) to execute all or part of the steps of the method described in each embodiment of the present invention. The aforementioned storage medium includes various media that can store program codes, such as a USB flash drive, a mobile hard disk, a read-only memory (ROM), a random access memory (RAM), a magnetic disk, or an optical disk.

[0114] An embodiment of the present invention provides a readable storage medium including: a program or instruction is stored on the readable storage medium, and when the program or instruction is executed by a processor, the steps of the above-mentioned GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method are implemented, and the same technical effect can be achieved. To avoid repetition, the present invention will not be described in detail.

[0115] Finally, it should be noted that the above embodiments are merely illustrative of the technical solutions of the embodiments of the present invention, and are not intended to limit them. Although the present invention has been described in detail with reference to the aforementioned embodiments, those skilled in the art should understand that they may still modify the technical solutions described in the aforementioned embodiments, or replace some of the technical features therein with equivalents; and such modifications or replacements do not deviate from the essence of the corresponding technical solutions from the spirit and scope of the technical solutions of the embodiments of the present invention. Any changes or replacements that can be easily conceived by those skilled in the art within the technical scope disclosed in the present invention should be covered by the scope of protection of the present invention.

Claims

1. A GIS-based one- and two-dimensional hydrodynamic model lateral coupling correlation method, characterized by: include: S1: Collect river water area vector files, river centerline vector files and 2D land boundary vector files; S2: Based on the river water area vector file, the river centerline vector file, and the two-dimensional land boundary vector file, a two-dimensional partition boundary vector file is generated, and a two-dimensional grid file is generated by a grid generator with a range consistent with the two-dimensional partition boundary vector file; S3: vectorizing the two-dimensional grid file to obtain a grid node file; S4: extracting two-dimensional boundary line segmentation points from the river water area vector file and the river centerline vector file; S5: extracting the two-dimensional grid, grid nodes, two-dimensional boundary and boundary grid nodes of each of the two-dimensional partitions based on the two-dimensional boundary line segmentation points and the grid node file, and establishing a coupling relationship between the grid nodes and the boundary in each of the two-dimensional partitions; S6: Integrate the coupling relationship between the grid nodes and the boundaries in each of the two-dimensional partitions to form a global coupling relationship.

2. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 1, characterized in that: The S2 specifically includes: S201: Generate a partition surface file based on the river centerline vector file and the two-dimensional land boundary vector file; S202: removing the river water area surface vector file from the partition surface file to obtain the two-dimensional partition boundary vector file; S203: Based on the two-dimensional partition boundary vector file, generate a two-dimensional grid file consistent with the range of the two-dimensional partition boundary vector file through SMS grid generation software.

3. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 1, characterized in that: The S3 specifically includes: S301: Reading node information and triangular mesh information in the two-dimensional mesh file; wherein the triangular mesh information includes a mesh name, and the node information includes a node name, a node x-coordinate, a node y-coordinate, a node elevation, and a node serial number; S302: Arranging the node information and the triangular mesh information into an Excel table; S303: Using GIS software to convert the node coordinate information in the Excel table into a point vector file; S304: performing deduplication processing on nodes with duplicate coordinates in the point vector file; S305: Saving the point vector file after deduplication processing as the grid node file.

4. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 1, characterized in that: The S4 specifically includes: S401: converting the river water area vector file into a river water area shoreline, and removing the intersecting parts of the river water area shoreline to obtain a river embankment line; S402: extracting the starting point and the end point of the river embankment line, and deleting duplicate position points to obtain the first type of two-dimensional boundary line segmentation points; S403: Calculating the spatial intersection of the river centerline vector file and the river water area vector file, and extracting the intersection point element; S404: Deleting duplicate points in the intersection point elements to obtain a second type of two-dimensional boundary line segmentation points; S405: Integrate the first type of two-dimensional boundary line segmentation points and the second type of two-dimensional boundary line segmentation points to form the two-dimensional boundary line segmentation points.

5. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 1, characterized in that: The S5 specifically includes: S501: extracting the two-dimensional partition range vector data of the current two-dimensional partition; S502: extracting the river embankment lines of the current two-dimensional partition based on the two-dimensional partition range, and removing the river embankment lines within the 50m buffer range where the starting and ending points of the river embankment lines are located, to obtain the main river embankment line; S503: Counting the minimum side length of the two-dimensional grid in the current two-dimensional partition, generating a buffer zone inward with a buffer distance of 0.5 times the minimum side length, and extracting the two-dimensional grid cells associated with the current two-dimensional partition based on the generated buffer zone; S504: Extracting a grid node corresponding to the current two-dimensional partition according to the two-dimensional grid unit; S505: identifying and extracting the two-dimensional boundary of the current two-dimensional partition by combining the two-dimensional grid cells of the current two-dimensional partition, the two-dimensional boundary line segmentation points, and the main levee line of the river channel; S506: extracting associated grid nodes of each two-dimensional boundary according to a many-to-many spatial relationship; and associating the associated grid nodes with boundary codes; S507: Calculating the mileage value of each of the two-dimensional boundary grid nodes on the boundary line; S508: Sort the boundary grid nodes in order according to the boundary name and the mileage value, and complete the construction of the coupling relationship between the grid nodes and the boundary of the current two-dimensional partition; S509: Repeat S501 to S508 until the coupling relationship between the grid nodes and boundaries of all partitions is established.

6. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 5, characterized in that: The S505 specifically includes: S5051: Fusing the two-dimensional grid cells of the current two-dimensional partition to generate a two-dimensional grid range surface, and extracting the corresponding two-dimensional grid range line based on the two-dimensional grid range surface; S5052: setting the average grid side length calculated from the two-dimensional grid file as a tolerance value for boundary recognition, and using the tolerance value as a search radius, taking the grid node with the smallest distance between the two-dimensional boundary line segmentation point and the grid node within the tolerance range as a feature grid node; S5053: Segmenting the two-dimensional grid range line according to the characteristic grid nodes to obtain multiple two-dimensional boundaries within the current two-dimensional partition; S5054: Establishing a one-to-one correspondence between each of the two-dimensional boundaries and the main embankment line of the river channel. If the directions of the two-dimensional boundaries and the main embankment line of the river channel are inconsistent, adjusting the direction of the two-dimensional boundaries to be consistent with the direction of the embankment line; S5055: Label each of the two-dimensional boundaries. If there is a corresponding main river embankment line for the two-dimensional boundary, the current two-dimensional boundary is recorded as a lateral coupling boundary, and the number is consistent with the serial number in the main river embankment line number; otherwise, the current two-dimensional boundary is recorded as a land boundary, and the land boundary is numbered according to the numbering rule of LDBJ+serial number.

7. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 6, characterized in that: The S6 specifically includes: S601: merging the two-dimensional boundaries of the two-dimensional partitions, and renumbering the land boundaries that do not correspond to the main embankment line of the river; S602: Merge the boundary grid nodes after sorting the two-dimensional partitions, and adjust the boundary number information of the corresponding nodes according to the unified land boundary number, complete the consistency matching of the global boundary, and form the global coupling relationship.

8. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 1, characterized in that: After S6, the method further includes: S7: writing the global coupling relationship into a text file, and outputting the text file into an inp format file and a dat format file.

9. The GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method according to claim 8, characterized in that: The S7 specifically includes: S701: writing node information including node name, node x-coordinate, node y-coordinate, node elevation and node serial number in the two-dimensional grid file into the text file line by line; S702: Writing the mesh names in the triangular mesh information into the text file line by line; S703: Writing the two-dimensional boundary and associated node information into the text file; S704: Outputting the text file after writing the node information, the grid name, the two-dimensional boundary and the associated node information into an inp format file and a dat format file.

10. A GIS-based one- and two-dimensional hydrodynamic model lateral coupling correlation system, characterized by: include: processor and memory; The memory stores programs or instructions that can be run on the processor, and when the programs or instructions are executed by the processor, the steps of the GIS-based one- and two-dimensional hydrodynamic model lateral coupling association method as described in any one of claims 1 to 9 are implemented.

Images