CN107145678B - Method for calibrating roughness of river channel plane two-dimensional water flow model - Google Patents
Method for calibrating roughness of river channel plane two-dimensional water flow model Download PDFInfo
- Publication number
- CN107145678B CN107145678B CN201710365078.7A CN201710365078A CN107145678B CN 107145678 B CN107145678 B CN 107145678B CN 201710365078 A CN201710365078 A CN 201710365078A CN 107145678 B CN107145678 B CN 107145678B
- Authority
- CN
- China
- Prior art keywords
- roughness
- water
- water level
- value
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
Images
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F30/00—Computer-aided design [CAD]
- G06F30/20—Design optimisation, verification or simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Measuring Volume Flow (AREA)
Abstract
The invention discloses a method for calibrating the roughness of a two-dimensional water flow model of a river channel plane, which comprises the steps of firstly defining roughness zones according to the characteristics of the river channel plane and assigning initial values to the roughness in each zone; selecting the current flow rate of which the roughness is to be calibrated and the corresponding river reach outlet water level; carrying out hydrodynamic calculation by using a river channel plane two-dimensional water flow mathematical model; determining the error between the model calculated value and the measured value of the water level of each monitoring point; judging whether the water level error is small enough; determining a water passing area according to the water depth, and estimating the water level error of each grid point in the water passing area by the water level error interpolation of each monitoring point; adjusting the roughness values of the grid points of the water passing area by a bisection method according to the water level errors of the grid points in the water passing area; further judging whether all the waiting rate fixed flow rates are calibrated; if yes, the roughness rating is completed, and the result of the roughness rating is stored. The method can reduce the man-made interference in the iterative trial calculation process and overcome the defects of unclear flow, harsh conditions, difficult application and the like in the prior art.
Description
Technical Field
The invention relates to the technical field of river dynamics research, in particular to a method for calibrating the roughness of a two-dimensional river channel plane water flow model.
Background
Roughness is an important parameter of a two-dimensional water flow mathematical model (hereinafter referred to as a "two-dimensional water flow model") on a river channel plane, the two-dimensional water flow model can be used for calculating factors such as flow velocity, water level and water depth of a natural river, and the hydraulic factors can be obtained to provide hydraulic conditions for solving the problems of water conservancy, ecology, environment and the like. The two-dimensional water flow model is realized on a computer, is called as river channel plane two-dimensional water flow mathematical model software, the calculation process is essentially a process of solving a discrete plane two-dimensional shallow water equation, besides initial conditions and boundary conditions are required to be given, relevant parameters are also required to be given, wherein the roughness is an important parameter of the two-dimensional water flow model, the river channel roughness is not only related to the roughness of a bed surface, but also related to the river flow plane shape, sediment forming accumulation body and water flow characteristics, and is a comprehensive hydraulic friction coefficient, and the parameter directly determines the rationality of a model calculation result.
In the prior art, three methods are mainly used for calibrating the roughness of a two-dimensional water flow model of a river channel plane: the method comprises the steps of iterative trial algorithm, empirical distribution method and theoretical calculation method. However, the three methods have the defects of unclear flow, harsh conditions, difficult application and the like, no experience or theoretical formula can directly calculate the roughness distribution of the planar two-dimensional bed surface at present, and the traditional iterative trial calculation process mainly depends on manual intervention, has no certain rule and strong subjectivity, and roughness values obtained by different personnel in a ratio are different even if the same model is used for processing the same river reach.
Disclosure of Invention
The invention aims to provide a method for calibrating the roughness of a two-dimensional water flow model of a river channel plane, which can reduce artificial interference in an iterative trial calculation process and overcome the defects of unclear flow, harsh conditions, difficult application and the like in the prior art.
The embodiment of the invention provides a method for calibrating the roughness of a two-dimensional water flow model of a river channel plane, which comprises the following steps:
step 1: defining roughness partitions according to the characteristics of the river channel plane, and assigning initial values to the roughness in each partition by using a roughness initial value function;
step 2: selecting the current flow of which the roughness is to be calibrated and the corresponding river reach outlet water level, and using the current flow as a trial calculation input condition of a river course plane two-dimensional water flow mathematical model;
and step 3: carrying out hydrodynamic calculation by using a river channel plane two-dimensional water flow mathematical model according to the trial calculation input conditions to obtain model calculation values of the water levels of all monitoring points;
and 4, step 4: determining the error between the model calculated value and the measured value of the water level of each monitoring point;
and 5: judging whether the water level error of each monitoring point in the step 4 is small enough; if yes, executing step 8; if not, executing the step 6;
step 6: determining a water passing area according to the water depth, and estimating the water level error of each grid point in the water passing area by the water level error interpolation of each monitoring point;
and 7: determining the adjustment direction of the roughness according to the size and the sign of the water level error of each grid point in the water passing area, adjusting the roughness values of the grid points in the water passing area by a bisection method, and returning to the step 3 for processing;
and 8: further judging whether all the waiting rate fixed flow rates are calibrated; if yes, executing step 9; if not, returning to the step 2, and carrying out roughness rating of the next flow level;
and step 9: the roughness rating was completed and the results of the roughness rating were saved.
According to the technical scheme provided by the invention, the method can reduce the artificial interference in the iterative trial calculation process, overcomes the defects of unclear flow, harsh conditions, difficult application and the like in the prior art, and improves the roughness calibration efficiency while ensuring the reasonable roughness value of the two-dimensional water flow model in the river plane.
Drawings
In order to more clearly illustrate the technical solutions of the embodiments of the present invention, the drawings needed to be used in the description of the embodiments are briefly introduced below, and it is obvious that the drawings in the following description are only some embodiments of the present invention, and it is obvious for those skilled in the art to obtain other drawings based on the drawings without creative efforts.
Fig. 1 is a schematic flow chart of a method for calibrating roughness of a two-dimensional river plane water flow model according to an embodiment of the present invention;
fig. 2 is a schematic diagram of water levels at two monitoring points in an example of the present invention.
Detailed Description
The technical solutions in the embodiments of the present invention are clearly and completely described below with reference to the drawings in the embodiments of the present invention, and it is obvious that the described embodiments are only a part of the embodiments of the present invention, and not all embodiments. All other embodiments, which can be derived by a person skilled in the art from the embodiments of the present invention without making any creative effort, shall fall within the protection scope of the present invention.
The embodiment of the present invention will be further described in detail with reference to the accompanying drawings, and as shown in fig. 1, a schematic flow chart of a method for rating the roughness of a two-dimensional river channel plane water flow model provided by the embodiment of the present invention is shown, where the method includes:
step 1: defining roughness partitions according to the characteristics of the river channel plane, and assigning initial values to the roughness in each partition by using a roughness initial value function;
in the step, the defining of the roughness partition according to the characteristics of the river channel plane is to define the roughness partition according to the bed sand composition and vegetation type factors of the river channel plane to be treated, and then to assign initial values to the roughness in each partition according to a roughness initial value function.
For natural river channels, the bed surface characteristics of river channels, islands, beaches and cardiac beaches are often different, so that roughness partitions are divided according to the bed surface characteristics, river reach satellite pictures can be referred to, and if site survey data exists, the roughness partitions can be accurately performed.
After the roughness division is carried out on the researched river reach, the initial roughness values n of different areas can be determined according to the existing hydraulic achievement0The value ranges of (A) are shown in Table 1 belowThe following steps:
TABLE 1 Natural river roughness n0Value of
After the rough rate partition is carried out, for each partition with different rough rate value intervals, the process of assigning initial values to the rough rates in each partition by using a rough rate initial value function is as follows:
the specific form of the roughness initial value function established is as follows:
in the above formula, n0(x, y) is the initial roughness value at the midpoint (x, y) of the plane; zx,yIs the bed surface elevation at point (x, y), ZmaxAnd ZminThe maximum bed surface elevation value and the minimum bed surface elevation value of the current roughness partition are respectively; n ismaxAnd nminThe maximum and minimum acceptable roughness values for the current roughness partition, respectively.
The essence of the above formula is that the elevation of the bed surface is taken as a basis, linear interpolation is carried out on the roughness of each coordinate point in a roughness partition, and the roughness values are considered to be smaller at the lower part of the elevation of the bed surface in the same roughness partition.
Step 2: selecting the current flow of which the roughness is to be calibrated and the corresponding river reach outlet water level, and using the current flow as a trial calculation input condition of a river course plane two-dimensional water flow mathematical model;
in this step, a typical flow rate is selected as the current flow rate of the roughness to be measured, and the river reach outlet water level corresponding to the flow rate is determined, and the current flow rate and the water level are known measured values, which are used as trial input conditions of the mathematical model of the two-dimensional water flow on the river course plane.
For the same overflowing area in a river reach, when the overflowing flow (water depth) is different, the represented roughness values are also different, so that the roughness values of different incoming water flows need to be considered in the calculation of a two-dimensional river plane water flow mathematical model. After typical flows with different sizes are selected as grading flows, the embodiment of the invention respectively carries out roughness rating on each flow grade, and the grading roughness corresponding to the grading flows can be used for interpolation calculation of roughness values of other flows.
And step 3: carrying out hydrodynamic calculation by using a river channel plane two-dimensional water flow mathematical model according to the trial calculation input conditions to obtain model calculation values of the water levels of all monitoring points;
and 4, step 4: determining the error between the model calculated value and the measured value of the water level of each monitoring point;
in the step, the specific process is as follows: after the model calculation value of the water level of each monitoring point is obtained, the model calculation value is compared with the actually measured water level of each monitoring point, and an error value between the model calculation value and the actually measured water level is calculated.
And 5: judging whether the water level error of each monitoring point in the step 4 is small enough; if yes, executing step 8; if not, executing the step 6;
step 6: determining a water passing area according to the water depth, and estimating the water level error of each grid point in the water passing area by the water level error interpolation of each monitoring point;
in the step, only the roughness of the water passing area influences the water level of the monitoring point, so that the water passing area under the current flow condition can be determined according to the water depth value output in the trial calculation process of the two-dimensional water flow model, and thus, the roughness adjusting area on the river channel plane can be accurately determined.
The water level error of each grid point is the basis of roughness adjustment, and because only the monitoring point has the calculated error value of the water level, and the water level errors of other non-monitoring points in the water passing area need to be interpolated and estimated by using a two-dimensional interpolation method, the embodiment of the invention adopts an inverse distance weighting method to estimate the water level error of each grid point in the water passing area, and the specific calculation formula is as follows:
in the above formula, e0Is an estimate of the water level error at non-monitored points in the water passing region, eiThe water level error calculation value of the ith monitoring point of the water passing area is obtained, and I is the total number of the monitoring points in the water passing area; d0,iIs the distance between the currently processed non-monitoring point and the ith monitoring point; k is an index whose value is such that e0The estimation error of (2) is small, and k can be 1 in the practical application process.
And 7: determining the roughness adjusting direction according to the size and the sign of the water level error of each grid point in the water passing area, adjusting the roughness value of the grid point in the water passing area by a bisection method, and returning to the step 3 for processing;
in this step, after obtaining the water level error of each grid point in the water passing region, the roughness value of each grid point in the water passing region may be calibrated by using a dichotomy roughness adjustment formula, where the dichotomy roughness adjustment formula is as follows:
in the above formula, nx,yRepresents the roughness value at the current point (x, y), and k +1 in the superscript represent the number of iterations; n isi,minRepresents the minimum advisable roughness value of the ith roughness partition to which the current point (x, y) belongs; n isi,maxRepresenting the maximum admissible roughness value of the i-th roughness partition to which the current point (x, y) belongs;
when the calculated value of the water level is higher, the point roughness is reduced; when the calculated water level value is lower, the point roughness is increased.
And after the roughness values of all the points in the water passing area are adjusted, repeating the steps 3-7 until the water level calculation errors of all the monitoring points meet the precision requirement under the current flow rate condition.
And 8: further judging whether all the waiting rate fixed flow rates are calibrated; if yes, executing step 9; if not, returning to the step 2, and carrying out roughness rating of the next flow level;
and step 9: the roughness rating was completed and the results of the roughness rating were saved.
The following detailed description will be made of the application of the above calibration method by specific examples:
firstly, a plane two-dimensional water flow mathematical model of a bending branch river reach is established, two monitoring points exist in the river reach, and then water levels and corresponding flow rates of different scenes in the river reach are calculated, as shown in the following tables 2 and 3:
TABLE 2 calculation of design flood levels in river reach
TABLE 3 calculation of river reach design flood flow
After the calculation conditions are complete, the roughness of each level of flow in table 3 is calibrated by the method described in the embodiment of the invention. The water levels at the two monitoring points under each flow are calculated by a well-calibrated model, as shown in fig. 2, which is a schematic diagram of the water levels at the two monitoring points in the example of the present invention, as can be seen from fig. 2: after the roughness is rated by the method of the embodiment of the invention, the calculated values of the water levels of the monitoring points 1 and 2 are basically consistent with the measured values, which shows that the method for rating the roughness of the two-dimensional water flow model of the river channel plane is feasible.
In summary, the method of the embodiment of the present invention has the following advantages:
(1) the method has clear process for rating the roughness, loose conditions and easy implementation, can realize the roughness rating by a machine, has less artificial interference, and greatly reduces the subjectivity of manual trial calculation;
(2) according to the method, the roughness of the plane river channel is divided into regions according to the characteristics of the bed surface, the roughness of each region is selected from an initial value and is rated in a roughness reasonable value range specified by hydraulics, and the reasonability of the rated roughness can be comprehensively ensured;
(3) after each round of trial calculation, the water passing area is determined according to the water depth information, the adjustment space of the roughness can be accurately determined, the calculated amount is greatly reduced, meanwhile, a dichotomy roughness adjustment formula is provided based on the monotonic relation between the roughness and the water level, and the convergence speed is higher.
The above description is only for the preferred embodiment of the present invention, but the scope of the present invention is not limited thereto, and any changes or substitutions that can be easily conceived by those skilled in the art within the technical scope of the present invention are included in the scope of the present invention. Therefore, the protection scope of the present invention shall be subject to the protection scope of the claims.
Claims (6)
1. A method for calibrating roughness of a two-dimensional water flow model of a river channel plane is characterized by comprising the following steps:
step 1: defining roughness partitions according to the characteristics of the river channel plane, and assigning initial values to the roughness in each partition by using a roughness initial value function;
step 2: selecting the current flow of which the roughness is to be calibrated and the corresponding river reach outlet water level, and using the current flow as a trial calculation input condition of a river course plane two-dimensional water flow mathematical model;
and step 3: carrying out hydrodynamic calculation by using a river channel plane two-dimensional water flow mathematical model according to the trial calculation input conditions to obtain model calculation values of the water levels of all monitoring points;
and 4, step 4: determining the error between the model calculated value and the measured value of the water level of each monitoring point;
and 5: judging whether the water level error of each monitoring point in the step 4 is smaller than a set threshold value or not; if yes, executing step 8; if not, executing the step 6;
step 6: determining a water passing area according to the water depth, and estimating the water level error of each grid point in the water passing area by the water level error interpolation of each monitoring point;
and 7: determining the roughness adjusting direction according to the size and the sign of the water level error of each grid point in the water passing area, adjusting the roughness value of the grid point in the water passing area by a bisection method, and returning to the step 3 for processing;
and 8: further judging whether all the waiting rate fixed flow rates are calibrated; if yes, executing step 9; if not, returning to the step 2, and carrying out roughness rating of the next flow level;
and step 9: the roughness rating was completed and the results of the roughness rating were saved.
2. A method for calibrating roughness of a two-dimensional river channel plane water flow model according to claim 1, wherein in the step 1, the defining of the roughness partition according to the characteristics of the river channel plane specifically comprises:
and defining a roughness zone according to the bed sand composition and vegetation type factors of the river channel plane to be treated.
3. A method for calibrating roughness of a two-dimensional river channel plane water flow model according to claim 1, wherein in the step 1, the initial roughness value function is used to assign initial roughness values to the partitions, and specifically comprises:
the specific form of the roughness initial value function established is as follows:
in the above formula, n0(x, y) is the initial roughness value at the midpoint (x, y) of the plane; z is a radical ofx,yIs the bed surface elevation at point (x, y), zmaxAnd zminThe maximum bed surface elevation value and the minimum bed surface elevation value of the current roughness partition are respectively; n ismaxAnd nminThe maximum and minimum acceptable roughness values for the current roughness partition, respectively.
4. A method for calibrating roughness of a two-dimensional water flow model of a river channel plane according to claim 1, wherein in the step 4, the determining an error between the calculated value and the measured value of the model of the water level at each monitoring point specifically comprises:
after the model calculation value of the water level of each monitoring point is obtained, the model calculation value is compared with the actually measured water level of each monitoring point, and an error value between the model calculation value and the actually measured water level is calculated.
5. A method for calibrating roughness of a two-dimensional water flow model in a river course plane according to claim 1, wherein in the step 6, the estimating of the water level error of each grid point in the water passing area by interpolation of the water level error of each monitoring point specifically comprises:
estimating the water level error of each grid point in the water passing area by adopting an inverse distance weight method, wherein the specific calculation formula is as follows:
in the above formula, e0Is an estimate of the water level error at non-monitored points in the water passing region, eiThe calculated value of the water level error of the ith monitoring point in the water passing area is obtained, and l is the total number of the monitoring points in the water passing area; d0,iIs the distance between the currently processed non-monitoring point and the ith monitoring point; k is an index whose value is such that e0The estimation error of (2) takes a small value.
6. A method for calibrating roughness of a two-dimensional river channel plane water flow model according to claim 1, wherein in the step 7, the adjusting roughness values of grid points of the water passing area by dichotomy specifically comprises:
the dichotomy roughness adjustment formula is as follows:
in the above formula, nx,yRepresents the roughness value at the current point (x, y), and k +1 in the superscript represent the number of iterations; n isi,minRepresents the minimum advisable roughness value of the ith roughness partition to which the current point (x, y) belongs; n isi,maxRepresenting the maximum admissible roughness value of the i-th roughness partition to which the current point (x, y) belongs;
when the calculated value of the water level is higher, the point roughness is reduced; when the calculated water level value is lower, the point roughness is increased.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710365078.7A CN107145678B (en) | 2017-05-22 | 2017-05-22 | Method for calibrating roughness of river channel plane two-dimensional water flow model |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710365078.7A CN107145678B (en) | 2017-05-22 | 2017-05-22 | Method for calibrating roughness of river channel plane two-dimensional water flow model |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107145678A CN107145678A (en) | 2017-09-08 |
CN107145678B true CN107145678B (en) | 2020-10-23 |
Family
ID=59778393
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710365078.7A Active CN107145678B (en) | 2017-05-22 | 2017-05-22 | Method for calibrating roughness of river channel plane two-dimensional water flow model |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107145678B (en) |
Families Citing this family (11)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN107451383B (en) * | 2017-09-29 | 2020-10-23 | 中国水利水电科学研究院 | Calibration method for initial bed sand gradation of plane two-dimensional water sand mathematical model |
CN108021762B (en) * | 2017-12-14 | 2021-07-16 | 黑龙江省水利水电勘测设计研究院 | Method and device for determining width of channel bottom |
CN108345738A (en) * | 2018-02-06 | 2018-07-31 | 广州地理研究所 | A kind of self rating method of middle Storm flood of small basins confluence Runoff Model parameter |
CN109145396B (en) * | 2018-07-27 | 2020-10-02 | 同济大学 | River course roughness zoning calibration method based on vegetation distribution |
CN109033725B (en) * | 2018-09-14 | 2020-05-05 | 中国水利水电科学研究院 | Estimation method for large-area bed surface shear stress of fixed bed river model test |
CN109840368B (en) * | 2019-01-17 | 2023-03-24 | 广东省水利水电科学研究院 | Irrigation area channel water level flow relation calibration method |
CN109979172A (en) * | 2019-04-09 | 2019-07-05 | 南京信息工程大学 | A kind of dynamic mountain torrents Critical Rainfall forecasting procedure based on Xinanjiang model |
CN110851765B (en) * | 2019-11-01 | 2022-02-18 | 华北电力大学 | Method for acquiring river course roughness in gradual flow state |
CN111753417B (en) * | 2020-06-18 | 2022-06-17 | 清华大学 | Parameter sensitivity analysis method based on uncertainty separation |
CN111859820A (en) * | 2020-06-23 | 2020-10-30 | 河海大学 | Roughness partition calibration method and system based on posterior distribution |
CN114818320B (en) * | 2022-04-25 | 2022-11-08 | 珠江水利委员会珠江水利科学研究院 | Giant river network hydrodynamic force mathematical model calibration method, system and storage medium |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090011735A (en) * | 2007-07-27 | 2009-02-02 | 인하대학교 산학협력단 | Developing analysis system of river plane change using air-photos and a gis and method of the same |
CN103590363A (en) * | 2013-11-26 | 2014-02-19 | 武汉大学 | Control node locating protecting method of undercutting channel segment water level lowering amplitude |
CN104750972A (en) * | 2015-01-30 | 2015-07-01 | 长江水利委员会水文局 | Complex water network area wedge-shaped water volume computing method |
CN106021854A (en) * | 2016-05-06 | 2016-10-12 | 西安理工大学 | A natural river course cross section stage-discharge relation determining method |
CN106372375A (en) * | 2016-10-26 | 2017-02-01 | 武汉大学 | Method for predicting diversion ratio change of distributary channels of two-stage bifurcated riverway in the downstream of reservoir |
-
2017
- 2017-05-22 CN CN201710365078.7A patent/CN107145678B/en active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
KR20090011735A (en) * | 2007-07-27 | 2009-02-02 | 인하대학교 산학협력단 | Developing analysis system of river plane change using air-photos and a gis and method of the same |
CN103590363A (en) * | 2013-11-26 | 2014-02-19 | 武汉大学 | Control node locating protecting method of undercutting channel segment water level lowering amplitude |
CN104750972A (en) * | 2015-01-30 | 2015-07-01 | 长江水利委员会水文局 | Complex water network area wedge-shaped water volume computing method |
CN106021854A (en) * | 2016-05-06 | 2016-10-12 | 西安理工大学 | A natural river course cross section stage-discharge relation determining method |
CN106372375A (en) * | 2016-10-26 | 2017-02-01 | 武汉大学 | Method for predicting diversion ratio change of distributary channels of two-stage bifurcated riverway in the downstream of reservoir |
Non-Patent Citations (2)
Title |
---|
复杂地形河道的二维数值模拟及问题处理;沈锐燕 等;《河北工程大学学报( 自然科学版)》;20130630;第30卷(第2期);62-64 * |
河道糙率率定方法存在问题与对策措施;滕凯 等;《西北水电》;20170131(第1期);1-4,12 * |
Also Published As
Publication number | Publication date |
---|---|
CN107145678A (en) | 2017-09-08 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107145678B (en) | Method for calibrating roughness of river channel plane two-dimensional water flow model | |
Tang et al. | Parameter identification for modeling river network using a genetic algorithm | |
CN107401140B (en) | The determination method of tributary lower reaches section channel forming disahcge | |
CN114329950B (en) | Dynamic generalization-based slope type submerged dike wave hydrodynamic influence numerical simulation method | |
CN106021854A (en) | A natural river course cross section stage-discharge relation determining method | |
CN109271672B (en) | River channel water surface line calculation method under interaction of river-lake-pump station | |
CN105069287A (en) | Prediction method of forest growth of overhead transmission line passageway | |
CN103886141A (en) | Method for computing fluvial cross section information in automatic-batch mode | |
CN111462450A (en) | Mountain torrent early warning method considering rainfall spatial heterogeneity | |
CN109902267B (en) | River channel safety discharge amount calculation method influenced by downstream lake jacking | |
CN111104746B (en) | River flood beach elevation determination method based on wavelet analysis | |
CN106446547A (en) | Generalized second-category beta distribution based flood frequency analysis method | |
CN113486556B (en) | Improved efficient automatic history fitting method for oil and gas reservoir | |
CN106320255A (en) | Overbank flow floodplain and channel excess flow computing method | |
CN103295203A (en) | Silicon wafer surface morphology reconstruction method based on B-spline | |
CN105277974A (en) | Stratum data interpolation method | |
CN103336907A (en) | Method for quickly calculating static storage capacity of reservoir based on DSI technology | |
Fortunato et al. | Automatic improvement of unstructured grids for coastal simulations | |
KR101273932B1 (en) | Method for interpolation using inverse distance weighting | |
CN109325254A (en) | Drain calculation method, system and the device of channel sand holding ability | |
CN103942409A (en) | Dredging channel back-silting amount predicting method based on timing sequence analysis-Markov chain method | |
CN111079065B (en) | Nested self-stabilization confluence calculation method for space-time dynamic optimization | |
CN109063306B (en) | Soil infiltration capacity space dispersion method of gridding Hebei model | |
CN113869350A (en) | Ocean current prediction method and system based on spatial feature difference | |
CN108345738A (en) | A kind of self rating method of middle Storm flood of small basins confluence Runoff Model parameter |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |