KR102328326B1 - Hydraulic hydrology analysis device applied with distributed model for advanced water quality prediction - Google Patents

Abstract

The present technology relates to a hydraulic hydrology analysis device applied with a distributed model for advanced water quality prediction. The hydraulic hydrology analysis device of the present technology includes: a grid setting unit that sets a grid resolution for a target basin; a hydrologic simulation performing unit that performs distributed hydrologic simulation by applying input data for the target basin corresponding to the set grid resolution to a hydrologic model; a hydraulic analysis performing unit that performs river hydraulic analysis by applying a hydraulic structure and a river cross-sectional shape for the target basin to a one-dimensional river hydraulic model; and a cooperation processing unit that allows the hydrological simulation performing unit to call a river hydraulic analysis step of the hydraulic analysis performing unit as a plurality of modularized functions so that an operation of the hydrologic simulation performing unit and an operation of the hydraulic analysis performing unit are cooperated with each other. The present technology can provide a hydraulic hydrology analysis device applied with a distributed model for advanced water quality prediction of basins or rivers.

Description

Hydraulic hydrology analysis device applied with distributed model for advanced water quality prediction}

The present invention relates to a hydraulic hydrological analysis apparatus to which a distributed model for advanced water quality prediction is applied, and more particularly, to a hydraulic hydrological analysis apparatus for advancing water quality prediction by linking distributed hydrological modeling and river hydraulic modeling.

The hydrologic model, which is a tool to analyze the water cycle in a watershed, can be divided into lumped, semi-distributed, and distributed models.

The centralized model views and interprets one watershed as one hydrologic unit. There is a limit to analyzing the cycle processes that are difficult to measure within the watershed.

The semi-distributed model overcomes the limitations of the centralized model by interpreting the subwatershed as a hydrological unit. However, there is a limitation that the circulation process within the subwatershed is also ignored.

The grid-based distribution model interprets individual grids as one hydrological unit. It can reflect the effects of various land uses, geology, and topography within the watershed.

Each of these three hydrologic models has their own strengths and weaknesses, but for integrated water management, a distributed hydrologic model that can physically grasp various hydrological cycles inside the watershed is suitable rather than focusing on the reaction at the exit of the basin or subcatchment. .

The distributed hydrology model can estimate the flow rate and filthy load at any point in the watershed, and because the process is physically traced, spatial distribution characteristics such as human activities or changes in natural conditions in the watershed are reflected in the calculation. There are advantages to doing it.

As a result of the four major rivers development project, 16 multi-function weirs were installed in the four major rivers, and the operation of the weirs significantly changed the river water quality compared to the past. In order to interpret the river water quality and ecosystem due to the influence of weirs, it is necessary to analyze the river hydraulics reflecting the operation of the weirs.

In addition, the issue of the effect of water quality improvement due to the opening of the weir is emerging as one of the major policy issues.

In the existing quasi-distributed and distributed hydrologic modeling, hydrological channel tracking functions such as kinematic wave, diffusion wave, and Muskingum have been utilized. A hydraulic tracing function is required. Therefore, in the present invention, it is intended to propose a method and apparatus for performing hydraulic channel tracking for a river section in a distributed hydrological model.

The inventor of the present invention has completed the present invention after long research and trial and error in order to solve these problems.

An embodiment of the present invention provides a hydraulic hydrological analysis apparatus to which a distribution model for advanced water quality prediction of watersheds and rivers is applied.

On the other hand, other objects not specified in the present invention will be additionally considered within the range that can be easily inferred from the following detailed description and effects thereof.

According to an embodiment of the present invention, there is provided a hydraulic hydrology analysis apparatus to which a distributed model for advanced water quality prediction is applied, comprising: a grid setting unit for setting a grid resolution for a target watershed; a sluice simulation performing unit for performing distributed sluice simulation by applying the input data for the target watershed corresponding to the set grid resolution to a sluice model; a hydraulic analysis performing unit for performing a river hydraulic analysis by applying a hydraulic structure and a river cross-sectional shape for the target watershed to a one-dimensional river hydraulic model; and a linkage processing unit that enables the hydrological simulation performing unit to call the river mathematical analysis step of the mathematical analysis performing unit as a plurality of modularized functions so that the operation of the hydrological simulation performing unit and the hydraulic analysis performing unit are linked with each other; may include

The linkage processing unit includes: a modularizing unit for modularizing the river hydraulic analysis step of the mathematical analysis performing unit into a function of several steps (K-RIVER functions); a conversion unit that converts the modularized K-RIVER functions into a library form so that they can be called by the sluice simulation execution unit; a calculation unit for classifying an outflow (slope) cell and a river cell in the sluice model, and calculating a flow rate (stream cell inflow) flowing into the river cell from the upstream cell; and an application unit that converts and applies the calculated river cell inflow into the lateral inflow of the river hydraulic model.

The link processing unit may perform simulation for the entire simulation time by repeatedly performing unit time simulation in the order of applying the conversion to the side inflow after calculating the river cell inflow.

At least one of the K-RIVER functions may be a function to determine the position of the river hydraulic model corresponding to the river cell in the sluice model, and to transmit the inflow of the river cell to the lateral inflow of the corresponding position.

At least one of the K-RIVER functions may be a function that determines a calculation time corresponding to the entire simulation time on the hydrological model and causes the river hydraulic analysis to be performed only at the corresponding calculation time.

The sluice simulation performing unit sets the initial soil moisture content of layers B and C among the soil conditions of the target watershed: layer A (surface flow), layer B (subsurface flow) and layer C (base runoff) to a fully saturated state. (a) step; (b) setting an observation-base runoff amount using observation data of a flow rate observatory located in the target watershed; (c) selecting a grid number corresponding to the location of the flow rate observatory; (d) performing runoff calculation (calculation-base runoff) of the entire inside of the target watershed per unit time under no rainfall condition; (e) comparing the calculated-basal runoff and the observed-basal runoff for the selected grid number; And, as a result of the comparison, if the calculation-base runoff is greater than the observed-base runoff, the initial soil function is automatically performed to improve the accuracy of the runoff simulation by returning to the step (d) and repeatedly calculating (f). It may include a correction function performing unit.

In the existing hydrological model for interpreting watershed runoff, hydrological channel tracing has been used to calculate runoff and track channel. For example, in rivers, hydraulic hydrologic analysis method and integrated device to which hydraulic channel tracking is applied can be provided.

1 is a diagram showing a detailed configuration of a hydrological analysis apparatus to which a distribution model for advanced water quality prediction is applied according to an embodiment of the present invention.
2 is a diagram illustrating a detailed configuration of a sluice simulator performing unit according to an embodiment of the present invention.
3 is a diagram showing a detailed configuration of a link processing unit according to an embodiment of the present invention.
4 is a view showing an example of topographical and geological data for hydrological analysis in the hydrological simulation performing unit.
5 is a view showing a rainfall observatory applied to hydrological analysis in the hydrological simulation performing unit.
6 is a view showing a water level observatory applied to the hydrological analysis in the hydrological simulation execution unit.
7 is a lattice-based distributed sluice model used by the sluice simulation performer according to an embodiment of the present invention, and is a conceptual diagram for explaining a K-DRUM.
8 is a view showing a flowchart of the operation of the initial soil function automatic correction function performing unit according to an embodiment of the present invention.
9 is a diagram illustrating a conceptual diagram of a groundwater simulation function performing unit performing a D-layer flow calculation according to an embodiment of the present invention.
10 is a view showing a process of reducing the calculation time by the MPI parallel code application unit according to an embodiment of the present invention.
11 is a view showing an example of the river bottom cross-section data applied to the river hydraulic analysis in the hydraulic analysis performing unit according to an embodiment of the present invention.
12 is a view showing a cross section and a schematic diagram of eight beams of Nakdonggang River applied to a river hydraulic analysis in a hydraulic analysis performing unit according to an embodiment of the present invention.
13 is a view showing an example of various hydraulic structures that can be considered in a K-River according to an embodiment of the present invention.
14 is a diagram illustrating a K-River governing equation according to an embodiment of the present invention.
15 is a diagram illustrating five-step functions modularized by a modularizer according to an embodiment of the present invention.
16 is a diagram illustrating a program and code flow diagram of a K-DRUM according to an embodiment of the present invention.
17 is a diagram illustrating a program and code flow diagram of a K-River according to an embodiment of the present invention.
18 is a diagram illustrating a K-River module call position by the link processing unit according to an embodiment of the present invention.
19 is a diagram illustrating a process in which a watershed runoff calculation, a river repair calculation, and a linked simulation are performed by repeatedly applying the calculation unit and the application unit to the unit time according to an embodiment of the present invention.
20 to 24 are diagrams showing the results of the distribution type hydrology model according to the hydraulic hydrologic analysis apparatus according to the embodiment of the present invention.
25 and 26 are diagrams showing the results of the simulation of the distribution type hydrology model and the river hydraulic model according to the hydraulic hydrologic analysis apparatus according to the embodiment of the present invention.
It is revealed that the accompanying drawings are exemplified by reference for understanding the technical idea of the present invention, and the scope of the present invention is not limited thereby.

In the following, the most preferred embodiment of the present invention is described. In the drawings, the thickness and the interval are expressed for convenience of description, and may be exaggerated compared to the actual physical thickness. In describing the present invention, well-known components that are not related to the gist of the present invention may be omitted. In adding reference numbers to the components of each drawing, it should be noted that only the same components are given the same number as possible even though they are indicated on different drawings.

1 is a diagram showing a detailed configuration of a hydrological analysis apparatus to which a distribution model for advanced water quality prediction is applied according to an embodiment of the present invention.

2 is a diagram illustrating a detailed configuration of a sluice simulator performing unit according to an embodiment of the present invention.

3 is a diagram illustrating a detailed configuration of a link processing unit according to an embodiment of the present invention.

First, referring to FIG. 1 , the hydro-hydraulic analysis apparatus 100 (hereinafter simply referred to as 'hydraulic-sluice analysis apparatus') to which the distributed model for advanced water quality prediction is applied includes a grid setting unit 110 and a hydrological simulation performing unit. 120 , a mathematical analysis performing unit 130 , and an associated processing unit 140 .

The grid setting unit 110 according to an embodiment of the present invention sets the grid resolution for the target watershed.

For example, when the watershed subject to hydrological analysis is the Nakdong River, it can be set from the minimum area and average area of administrative dongs of local governments belonging to the large Nakdong River area.

For example, if the minimum area of an administrative dong in Busan is 0.17 ㎢, it can be converted to a grid resolution to derive a minimum resolution of 412 m. If the average area of administrative dong in Busan is 3.37 km2, it is converted into grid resolution to derive an average resolution of 1.93 km. Accordingly, a predetermined value can be set as the grid resolution within the range of 412 m to 1.93 km.

If the minimum area of an administrative dong in Daegu Metropolitan City is 0.25㎢, it can be converted to a grid resolution to derive a minimum resolution of 500m. If the average area of administrative dong in Daegu Metropolitan City is 8.09 km2, it is converted into grid resolution to derive an average resolution of 2.84 km. Accordingly, a predetermined value can be set as the grid resolution within the range of 500m to 2.84km.

If the minimum area of administrative dong in Ulsan Metropolitan City is 0.5㎢, it can be converted to grid resolution to derive a minimum resolution of 707m. If the average area of administrative dong in Ulsan Metropolitan City is 18.18 km2, it is converted into grid resolution to derive an average resolution of 4.26 km. Accordingly, a predetermined value can be set as the grid resolution within the range of 707m to 4.26km.

As described above, the minimum resolution and average resolution can be derived for several administrative dongs belonging to the large Nakdong River region, and a predetermined value can be set as the grid resolution within these ranges. As an example, 1 km may be set as the grid resolution in consideration of model construction and calculation time. In this case, the large Nakdong River region can be composed of 39,928 grids. The grid resolution set in this way determines the basic grid size of the basic data.

Hereinafter, for convenience of explanation, an embodiment in which the target basin is the large Nakdong River region will be mainly examined, but the present invention is not limited thereto.

The sluice simulation performing unit 120 according to an embodiment of the present invention performs distributed sluice simulation by applying input data for a target watershed corresponding to the lattice resolution set by the lattice setting unit to the sluice model.

Input data may include topographical and geological data. The input data may be basically a grid-type digital elevation model (DEM).

For example, it may be spatial topography and geological data constructed by using a digital topographic map of the National Geographic Information Service (Ministry of Land, Infrastructure and Transport) to generate a 1km-resolution DEM and using a soil map of Soil Toram (Rural Development Administration). It is also possible to use the land cover map of the Environmental Geospatial Information Service (Ministry of Environment). 4 shows an example of topographical and geological data for hydrological analysis in the hydrological simulation performing unit. (a) is a DEM, (b) is a flow chart, (c) is a soil diagram, and (d) is a land use diagram, respectively.

For example, after collecting unit data from several rainfall observatories (eg, 212 locations) in the Nakdong River basin, it can be applied to a hydrological model as a grid-type rainfall through spatial distribution. 5 shows a rainfall observatory applied to hydrological analysis in the hydrological simulation execution unit. From this, rainfall observation and hydrological data are collected.

Also, for example, unit data can be collected from several water level observatories (74 locations) existing in the Nakdong River basin and applied to the hydrology model. 6 shows a water level observatory applied to hydrological analysis in the hydrological simulation execution unit. From this, water level observation and hydrological data are collected.

A K-water Distribution Runoff Model (K-DRUM) is applicable as a hydrology model according to an embodiment of the present invention. K-DRUM was developed as Fortran code by K-water. 7 is a lattice-based distributed sluice model used by the sluice simulation performing unit 120 according to an embodiment of the present invention, and is a conceptual diagram for explaining a K-DRUM.

As shown in Figure 7, K-DRUM calculates surface flow, subsurface flow, and base runoff by setting A layer, B layer, and C layer for each grid. . It is a model that allocates parameters such as land cover, soil depth, and roughness coefficient for each grid.

The vertical structure of the model consists of the surface layer, the shallow soil layers, A and B layers, and the groundwater layer, C layer. The tracing of the surface flow in the surface layer uses the universally applied kinematic wave, and the infiltration into the soil corresponding to the A and B layers uses the Green-Ampt formula, and the B layer And the flow of the groundwater layer C layer uses the linear storage method.

As shown in FIG. 2 , the hydrology simulation performing unit 120 according to the embodiment of the present invention performs an initial soil function automatic correction function performing unit 121, a rainfall spatial distribution interpolation function performing unit 122, and an evapotranspiration simulation function. It may include a unit 123 , a melting snow simulation function performing unit 124 , a groundwater simulation function performing unit 125 , and an MPI parallel code application unit 126 .

First, the initial soil function automatic correction function performing unit 121 improves the accuracy of the overall runoff simulation by reproducing the actual base runoff and applying the correct base runoff at the beginning of the rainfall runoff calculation.

To this end, the initial soil function automatic correction function performing unit may sequentially perform the following steps. It will be described with reference to FIG. 8 showing a flowchart of the operation of the initial soil function automatic correction function performing unit.

Step (a): Set the initial soil moisture content of layers B and C to full saturation among the soil conditions of layer A (surface flow), layer B (subsurface flow), and layer C (base runoff) of the target watershed.

Step (b): Observation-base runoff is set using the observation data of the flow station located in the target watershed.

(c) step: select the grid number corresponding to the location of the flow rate station

Step (d): Perform runoff calculation (calculation-base runoff) of the entire inside of the target watershed per unit time under no rainfall condition

Step (e): comparing the calculated-basal runoff and the observed-basal runoff for the selected grid number

Step (f): As a result of the comparison, if the calculated-basal runoff is greater than the observed-basal runoff, return to step (d) to repeat the calculation, and when the observed-basal runoff is the same, the initial soil moisture content is automatically corrected complete the course

The operation of the automatic initial soil function correction function performing unit uses the characteristic that the base runoff generated in each grid of K-DRUM decreases at a certain rate as the soil moisture content of the grid decreases in the fully saturated state. If uneven rainfall occurs throughout the watershed, the water content of the soil in the watershed shows an uneven distribution, but after about 3 to 5 days after the end of the rainfall, the decrease in the soil water content is due to the nature of the soil inside and the topographical and hydrological characteristics of the watershed. This is because it proceeds similarly to the decrease in base runoff.

The process of determining the initial conditions of K-DRUM provides convenience to the user through automated parameter setting rather than parameter correction based on the user's subjective judgment, and shortens the required time (trial and error method, etc.) Resolve any problems that may arise.

Next, the rainfall spatial distribution interpolation function performing unit 122 is a process of converting the point rainfall at the rainfall observatories scattered in the target watershed into the rainfall amount for each grid in order to convert the point rainfall amount to the input data of the watershed runoff model. Perform rainfall conversion.

The rainfall spatial distribution interpolation function execution unit converts the observation information of the rainfall observatory into spatial distribution grid rainfall to be utilized.

For this, one or more of Thiessen method, IDW method, and Kriging method can be selected to convert the spot rainfall.

Subsequently, the evapotranspiration simulation function performing unit 123 corrects the climate data of the observatory in the target watershed for each grid, calculates the potential evapotranspiration by the potential evapotranspiration calculation formula, and then applies the correction factor to calculate the actual evapotranspiration.

The evapotranspiration simulation function execution unit can utilize maximum temperature, minimum temperature, dew point temperature, sunshine time, and wind speed as climate factor data in this process. Here, it is possible to calculate the potential evapotranspiration for each grid by performing the correction by the daily temperature lapse rate and the correction by latitude (applied the FAO-56 Penman-Monteith calculation formula). Also, it is possible to calculate the actual evapotranspiration for each grid by applying a correction factor.

The melting snow simulation function performing unit 124 determines and calculates snowmelt if the atmospheric temperature is relatively low, and snowmelt if high, according to the reference temperature for classifying the melting snow for each grid. The process of estimating the melting snow is performed by calculating the atmospheric temperature for each grid by applying the temperature lapse rate accordingly.

The melting snow simulation function execution unit can utilize the maximum temperature, minimum temperature, dew point temperature, sunshine time, and wind speed as climate factor data in this process. Here, it is possible to calculate the amount of snow melt for each grid by performing correction by the daily temperature lapse rate (applied the Sugawara calculation formula).

The groundwater simulation function performing unit 125 sets the floor elevation of the D-floor simulating the groundwater outflow to be the same as sea level, and sets the depth of the D-layer according to the C-floor elevation to simulate the base runoff.

In K-DRUM, groundwater is simulated by the outflow of the D layer, and since the soil depth of the D layer is calculated to be constant, there is a limit to the simulation of the base runoff. Modify the depth to change according to the terrain shape. 9 is a conceptual diagram illustrating a groundwater simulation function performing unit performing a D-layer flow calculation. As shown in FIG. 9 , the bottom elevation of the D layer is set to be the same as the sea level, and the depth of the D layer is set according to the bottom elevation of the C layer.

The lateral runoff in the D-layer can be calculated according to the following equation.

Comparing the downstream grid and D-layer water level, if the upstream water level is high, the lateral runoff is applied to the downstream, and if the downstream water level is high, it means that there is no transverse runoff.

The MPI parallel code application unit 126 solves the computational problem, which is a disadvantage of the distributed model, through a parallelization processing technique.

To this end, the MPI parallel code application unit can perform task assignment and data distribution by using a method in which processes having a separate memory locally transmit, receive, and communicate data to share data. 10 is a view showing a process of shortening the calculation time by the MPI parallel code application unit. (a) shows region division for parallel calculation, and (b) shows a plurality of region division MPI algorithms. As shown in the figure, it can be divided into 7 regions, regions 2 to 7 are calculated independently of each other (that is, parallel calculation is applied), and region 1 receives information after the calculation of regions 2 to 7 is completed. It can be calculated by combining with MPI.

The hydraulic analysis performing unit 130 according to an embodiment of the present invention performs a river hydraulic analysis by applying a hydraulic structure and a river cross-sectional shape for a target watershed to a one-dimensional river hydraulic model.

For example, the river cross-section data constructed using cross-sectional information of about 700 points (for example, 710 pieces) at intervals of about 3 to 400 m for the main stream of the Nakdong River (Andong Dam ~ Nakdong Estuary Bank) can be applied. 11 shows an example of the riverbed cross-section data applied to the river hydraulic analysis in the hydraulic analysis performing unit.

Also, for example, the hydraulic structure data constructed from the specifications of 8 weirs (Sangju Weir, Nakdan Weir, Gumi Weir, Chilgok Weir, River Weir, Dalseong Weir, Hapcheon Weir, and Haman Weir) installed on the Nakdong River through the Four Major Rivers Maintenance Project are collected. can be applied. 12 shows a schematic diagram of the Nakdong River 8 beam cross-section and dimensions applied to the river hydraulic analysis in the hydraulic analysis performing unit.

As input data for the river hydraulic model, link data between the K-DRUM inflow point and the K-River river section number can be constructed. As an example, the relationship between the water level observatory and major points in the Nakdong River basin and the one-dimensional model can be constructed with the following form of K-DRUM inflow point - K-River river section number linkage data.

K-DRUM - KRiver connection point division tributary and water station One-dimensional river model point number note One Andong Joji Dam 699 2 semi-transition confluence 696 C12-C11 3 Joining Song Yacheon 679 C14-C13 4 Join Micheon 673 C16-C14 5 Confluence of Sinyeokcheon, Pungsancheon and Haacheon 637 △9+10+11 6 Gudam Watermark 620+407 Gudam Bridge 7 Confluence of Naeseongcheon Stream 579 Dalji-C20 8 Joining Yeonggang 566 C31-Several 9 private water table 559+082 Sangpung Bridge 10 Sangjubo 539+183 11 Joining Byung Seong-cheon 537 C38-C36 12 Join Wicheon 514 C43-C38 13 stepping stone 506+185 14 Nakdong Watermark 505+114 Nakdan Bridge 15 Joining Singokcheon 486 △33+34 16 Gumibo 470+195 17 Gamcheon Joining 468 Haepyeong-C49 18 Confluence of Bongamcheon, Daemangcheon, Seobuncheon, Innocheon, Suseongcheon, and Gumicheon 456 △37 19 Gumi water table 437+080 Gumi Bridge 20 Confluence of Igyecheon, Jangamcheon, and Dangpyeongcheon 431 △38 21 Chilgokbo 414+275 22 Waegwan water table 411+021 Nakdongganggu Railway Bridge 23 Baekcheon and Shincheon join 383 △39+40+41 24 Kang information 364+268 25 Confluence of Kumho River 363 C83-C62 26 Dalseongbo 322+295 27 Hyeonpung Water Table 314 28 convergence 270 C94+Δ64 29 Hapcheonbo 263+504 30 yellow river confluence 259 C104-C97 31 magic water table 239 32 Nam river confluence 202 C116-C107 33 vibration level table 191+110 34 Gaeseong Stream, etc. 182 △82, upstream of Nakdong Bridge 35 Hamanbo 174+445 36 Imhaejin water level table 163 37 Cheongdocheon Stream joins 150 △84 38 fisheries water table 138+020 Susan Bridge 39 Miryang River confluence 110 C130+△91 40 Samrangjin Watermark 107+195 Samrangjin Bridge 41 Confluence of Yangsancheon Stream 54 △92 42 Confluence of Daecheoncheon 42 C134-C132 43 Joining Deokcheoncheon 34 △94 44 Gupo water table 29+160 Gupyo Bridge 48 haggut 0

A K-River can be used as a one-dimensional river hydraulic model according to an embodiment of the present invention. It is a one-dimensional river hydraulic analysis model based on FVM developed by K-water with C++ code. In K-River, various hydraulic structures such as beams, weirs and sluice gates can be considered, and transverse structures such as beams including a plurality of weirs and sluice gates as shown in FIG. 13 can be considered. Also, in K-River, it is possible to consider the flow to analyze the flow in various culverts. As the K-River governing equation, the conservative Saint-Venant equation can be used. In addition, a well-balanced numerical technique can be applied to solve the discretization of FVM-based control volume and the numerical imbalance of flow/generation terms. 14 shows the K-River governing equation.

The linkage processing unit 140 according to an embodiment of the present invention performs the river hydraulic analysis step of the hydraulic simulation performing unit in the hydrological simulation performing unit so that the operation of the hydraulic simulation performing unit 110 and the hydraulic analysis performing unit 120 are linked with each other. can be called as multiple modularized functions.

To this end, the linkage processing unit 140 may include a modularizing unit 142 , a converting unit 144 , a calculating unit 146 , and an applying unit 148 as shown in FIG. 3 .

The modularization unit 142 modularizes the river mathematical analysis step of the mathematical analysis performing unit into a function of several steps (K-RIVER functions). For example, it can be modularized as a function of 5 steps. 15 is a diagram showing five steps of functions modularized by the modularization unit.

It corresponds to the process of modularizing K-River by grouping it by process.

This makes it possible to perform simulation as a single executable file by linking K-DRUM developed with Fortran code and K-River developed with C++ code through the conversion unit to be described later.

In detail, a program and code flow diagram of K-DRUM may be shown as shown in FIG. 16 . And, the program and code flow chart of K-River may be shown as shown in FIG. 17 . At this time, in order to link different codes, the calculation step of K-River can be modularized into 5 functions, that is, kriver_read_file, kriver_show_time_info, kriver_pre_proceed, kriver_proceed_hour, kriver_write_files_finished as shown in FIG. 15 .

The conversion unit 144 converts the modularized K-RIVER functions into a library form so that they can be called by the sluice simulation execution unit 120 .

It corresponds to the process of converting the modularized K-River into a library that can be called externally.

Specifically, the modularized K-River functions are built in the form of a library (libKRiver.a) and linked to be called inside K-DRUM. That is, the K-River module can be called from the K-DRUM code.

18 shows a K-River module call position by the link processing unit according to an embodiment of the present invention. (a) shows FIG. 16, and (b) visually shows the K-River module call position with respect to FIG.

Referring to FIG. 18 , the kriver_read_file function is called after input data of K-DRUM, and the kriver_show_time_info function and kriver_pre_proceed function are called after K-DRUM initial value correction. The kriver_proceed_hour function is included in the calculation loop of K-DRUM, and the kriver_set_lateral_flow function is called for the connection between K-DRUM and K-River before. Specifically, it is a function that identifies the cell for generating the inflow into the river hydraulic model on the K-DRUM, finds the location of the corresponding K-River, and transfers it to the lateral inflow. For example, the kriver_set_lateral_flow(river_i, station, Q) function is a function that transmits K-DRUM information (flow rate) to the K-River. By passing the calculation time to the K-River, only the required section (time) can perform river flow calculation. Meanwhile, the connection information may be input from an external file (connection_data.txt) file.

The calculation unit 146 divides the outflow (slope) cell and the river cell in the sluice model, and calculates the flow rate (stream cell inflow) flowing into the river cell from the upstream cell. In the river cell of the hydrologic model, it is assumed that no outflow to the downstream cell occurs.

The application unit 148 converts the calculated river cell inflow into the lateral inflow of the river hydraulic model and applies it.

19 shows that the watershed runoff by the calculation unit and the application unit according to the embodiment of the present invention is calculated (Step 1), the river repair is calculated (Step 2), and the steps of Step 1 and Step 2 are repeatedly performed for a unit time. It shows the process of performing linked simulation by applying (Step 3). (a) shows the block diagram, (b) shows each cell.

As shown in FIG. 19, the linkage simulation can be performed in the same way as above by separating outflow (slope) cells and river cells on the hydrologic model, and K-DRUM (Fortran code) and K-River (C++ code) It is possible to perform simulation with a single executable file by linking.

20 to 24 are diagrams showing the results of the distribution type hydrology model according to the hydraulic hydrologic analysis apparatus according to the embodiment of the present invention.

The simulation period of hydrological analysis using K-DRUM was calculated from 9:00 on 2018/08/22 to 23:00 on 2018/09/04. The analysis was performed by applying the observed rainfall to the 212 precipitation stations constructed previously.

The observed flow rates at nine dam inflow points (Namgang dam, Andong dam, Imha dam, Hapcheon dam, Miryang dam, Gunwi dam, Yeongcheon dam, Bohyeonsan dam, and Gimcheon Buhang dam) were compared with the calculation results of K-DRUM.

As shown in FIGS. 20 to 24, the simulation results show that the inflow and arrival times of Namgang Dam, Andong Dam, Hapcheon Dam, and Imha Dam with large watersheds are well simulated, but Bohyeonsan Dam, Yeongcheon Dam, and Gunwi with relatively small watersheds. It can be seen that all dams and the like tend to overestimate the amount of runoff. In general, when applying the distributed hydrology model, it is necessary to calibrate the parameters for each subwatershed. This can greatly affect the accuracy of the model results, and although the accuracy of the current model is not sufficient, it is judged that the accuracy of the distributed hydrology model can be improved through parameter optimization for each subwatershed in the future.

25 and 26 are diagrams showing the results of the simulation of the distribution type hydrology model and the river hydraulic model according to the hydraulic hydrologic analysis apparatus according to the embodiment of the present invention.

The simulation period was calculated for 10 days from 00:00 on 2018/08/25 to 00:00 on 2018/09/04. It was simulated by linking K-DRUM and KRiver, and KRiver is called as a modularized function inside K-DRUM. A total of 20 river tributaries were connected to KRiver, and the branches are shown in the table below.

Branch name River section number Branch name River section number half change 696 Kumho River 363 Song Yacheon 679 turnaround 271 Micheon 673 yellow steel 260 Naeseongcheon 580 Topyeongcheon 248 Yeonggang 566 Namgang 203 Byung Seong-cheon 536 Gwangnyeocheon 200 Wicheon 514 Cheongdocheon Stream 151 Gamcheon 468 Yangsancheon 117 agar 440 Miryang River 111 hundred thousand 380 Yangsancheon 55

The water level observed at eight water level stations (Gudam, Sabeol, Nakdong, Gumi, Hyeonpung, Jinjin, Samrangjin, and Gupo) was compared with the calculation results of the K-River. As shown in FIGS. 25 and 26, overall It seems that the water level rise and fall trends are well simulated, and the peak flow rate and arrival time are also well simulated. Considering that it is much more difficult to reproduce the water level during the flood season than during the normal season, it is judged that it is necessary to perform verification for the simulation period of about one year in the future. In addition, although simple operation rules are currently applied based on the level of management, it is expected that simulation accuracy will be improved if the simulation operation optimized for the target event is performed.

Embodiments of the present invention may be implemented in the form of program instructions that can be executed through various computer means and recorded in a computer-readable medium. The computer-readable medium may include program instructions, data files, data structures, etc. alone or in combination. The program instructions recorded on the medium may be specially designed and configured for the present invention, or may be known and available to those skilled in the art of computer software. Examples of the computer-readable recording medium include magnetic media such as hard disks, floppy disks and magnetic tapes, optical media such as CD-ROMs and DVDs, and magnetic such as floppy disks. - Examples of program instructions such as magneto-optical and ROM, RAM, and flash memory can be executed by a computer using an interpreter as well as machine code such as those created by a compiler. contains high-level language codes. The hardware device described above may be configured to operate as at least one software module to perform the operations of the embodiments of the present invention, and vice versa.

As described above, the present invention has been described with specific matters such as specific components and limited embodiments and drawings, but these are provided to help a more general understanding of the present invention, and the present invention is not limited to the above embodiments. , various modifications and variations are possible from these descriptions by those of ordinary skill in the art to which the present invention pertains. Therefore, the spirit of the present invention should not be limited to the described embodiments, and not only the claims to be described later, but also all those with equivalent or equivalent modifications to the claims will be said to belong to the scope of the spirit of the present invention. .

100: Hydraulic hydrological analysis device to which a distributed model for advanced water quality prediction is applied
110: grid setting unit
120: attendant of the sluice model
121: Initial soil function automatic correction function execution unit
122: rainfall spatial distribution interpolation function execution unit
123: evapotranspiration simulation function execution unit
124: fusion snow simulation function execution unit
125: groundwater simulation function execution unit
126: MPI parallel code application part
130: mathematical analysis execution unit
140: link processing unit
142: modularization part
144: conversion unit
146: calculator
148: application part

Claims (6)

As a hydro-hydraulic analysis device to which a distribution model is applied for advanced water quality prediction,
a grid setting unit for setting a grid resolution for a target watershed;
a sluice simulation performing unit for performing distributed sluice simulation by applying the input data for the target watershed corresponding to the set grid resolution to a sluice model;
a hydraulic analysis performing unit for performing a river hydraulic analysis by applying a hydraulic structure and a river cross-sectional shape for the target watershed to a one-dimensional river hydraulic model; and
A linkage processing unit that enables the hydrological simulation performing unit to call the river mathematical analysis step of the mathematical analysis performing unit as a plurality of modularized functions so that the operation of the hydrological simulation performing unit and the hydraulic analysis performing unit are linked with each other; but,
The link processing unit,
a modularization unit for modularizing the river hydraulic analysis step of the mathematical analysis performing unit into a function of several steps (K-RIVER functions);
a conversion unit that converts the modularized K-RIVER functions into a library form so that they can be called by the sluice simulation execution unit;
a calculation unit for classifying an outflow (slope) cell and a river cell in the sluice model, and calculating the flow rate (stream cell inflow) flowing into the river cell from the upstream cell; and
and an application unit that converts and applies the calculated river cell inflow into the lateral inflow of the river hydraulic model;
At least one of the K-RIVER functions,
Hydraulic hydrolog analysis apparatus, characterized in that it is a function to determine the position of the hydraulic model corresponding to the river cell in the hydrological model and to transmit the inflow of the river cell to the lateral inflow of the corresponding position. delete According to claim 1,
The link processing unit,
Hydrologistic analysis apparatus, characterized in that the simulation proceeds for the entire simulation time by repeatedly performing the unit time simulation in the order of applying the conversion to the side inflow after calculating the river cell inflow. delete As a hydro-hydraulic analysis device to which a distribution model is applied for advanced water quality prediction,
a grid setting unit for setting a grid resolution for a target watershed;
a sluice simulation performing unit for performing distributed sluice simulation by applying the input data for the target watershed corresponding to the set grid resolution to a sluice model;
a hydraulic analysis performing unit for performing a river hydraulic analysis by applying a hydraulic structure and a river cross-sectional shape for the target watershed to a one-dimensional river hydraulic model; and
A linkage processing unit that enables the hydrological simulation performing unit to call the river mathematical analysis step of the mathematical analysis performing unit as a plurality of modularized functions so that the operation of the hydrological simulation performing unit and the hydraulic analysis performing unit are linked with each other; but,
The link processing unit,
a modularization unit for modularizing the river hydraulic analysis step of the mathematical analysis performing unit into a function of several steps (K-RIVER functions);
a conversion unit that converts the modularized K-RIVER functions into a library form so that they can be called by the sluice simulation execution unit;
a calculation unit for classifying an outflow (slope) cell and a river cell in the sluice model, and calculating the flow rate (stream cell inflow) flowing into the river cell from the upstream cell; and
and an application unit that converts and applies the calculated river cell inflow into the lateral inflow of the river hydraulic model;
At least one of the K-RIVER functions,
Hydraulic hydrological analysis apparatus, characterized in that it is a function to determine the calculation time corresponding to the entire simulation time on the hydrologic model and to perform the river hydraulic analysis only in the corresponding calculation time. According to claim 1,
The execution unit of the sluice model,
(a) setting the initial soil moisture content of layer B and layer C to be fully saturated among the soil conditions of layer A (surface flow), layer B (subsurface flow) and layer C (base runoff) of the target watershed;
(b) setting an observation-base runoff amount using observation data of a flow rate observatory located in the target watershed;
(c) selecting a grid number corresponding to the location of the flow rate observatory;
(d) performing runoff calculation (calculation-base runoff) of the entire inside of the target watershed per unit time under no rainfall condition;
(e) comparing the calculated-basal runoff and the observed-basal runoff for the selected grid number; and
As a result of the comparison, if the calculated-base runoff is greater than the observed-base runoff, the initial soil function automatic correction to improve the accuracy of the runoff simulation by performing (f) returning to the step (d) and performing repeated calculations; Hydrographic analysis device comprising a; function performing unit.

Images