CN113221366A - Method and system for calculating dynamic sediment transport water amount of river channel - Google Patents

Method and system for calculating dynamic sediment transport water amount of river channel Download PDF

Info

Publication number
CN113221366A
CN113221366A CN202110549830.XA CN202110549830A CN113221366A CN 113221366 A CN113221366 A CN 113221366A CN 202110549830 A CN202110549830 A CN 202110549830A CN 113221366 A CN113221366 A CN 113221366A
Authority
CN
China
Prior art keywords
river channel
river
sand
data
acquisition
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.)
Granted
Application number
CN202110549830.XA
Other languages
Chinese (zh)
Other versions
CN113221366B (en
Inventor
张金良
刘继祥
鲁俊
谢亚光
仝亮
付健
罗秋实
崔振华
宋伟华
慕红伟
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Yellow River Engineering Consulting Co Ltd
Original Assignee
Yellow River Engineering Consulting Co Ltd
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Yellow River Engineering Consulting Co Ltd filed Critical Yellow River Engineering Consulting Co Ltd
Priority to CN202110549830.XA priority Critical patent/CN113221366B/en
Publication of CN113221366A publication Critical patent/CN113221366A/en
Application granted granted Critical
Publication of CN113221366B publication Critical patent/CN113221366B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • G—PHYSICS
    • G06—COMPUTING; CALCULATING OR COUNTING
    • G06F—ELECTRIC DIGITAL DATA PROCESSING
    • G06F2113/00—Details relating to the application field
    • G06F2113/08—Fluids

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)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention provides a method and a system for calculating the dynamic sand transportation water amount of a river channel, wherein the method comprises the following steps: collecting basic data, and calculating the corresponding river channel section average water depth under the condition of actually measured flow; fitting a power function relation of the average water depth and the flow of the river channel, combining a preset formula, deducing an expression of the instantaneous sand conveying capacity of the river channel, defining a non-uniform coefficient in the process of introducing the flow, deducing a calculation method of the dynamic sand conveying capacity of the river channel, taking measured data of the year that the erosion and deposition of the river channel is close to the erosion and deposition balance, and calibrating parameters in the calculation method of the dynamic sand conveying capacity of the river channel; establishing a relation between the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount to obtain a calculation method of the river channel dynamic sand conveying water amount; and analyzing the change rule of the non-uniform coefficient phi and the erosion and deposition parameter lambda in the river course flow process by counting the basic data to obtain the calculation method of the dynamic sediment transport water volume of the river course. The method has the advantages of simple steps, reliable results, convenient calculation and easy operation.

Description

Method and system for calculating dynamic sediment transport water amount of river channel
Technical Field
The invention relates to the technical field of engineering sediment, in particular to a method and a system for calculating dynamic sediment transport water quantity of a river channel.
Background
The total length of the yellow river dry flow is 5464km, which is divided into upstream, midstream and downstream, the river estuary town from river source to Tokton county of inner Mongolia is the yellow river upstream, the river estuary town to Henan Zhengzhou and Huayu are the yellow river midstream, and the yellow river downstream from below the Tayu to Haishu. The upstream of the yellow river is the main water source area of the yellow river, the midstream of the yellow river is the main sand producing area of the yellow river, and the downstream of the yellow river is the key river segment for flood control of the yellow river.
The research on the sand transporting water amount of the river channel at the downstream of the yellow river has important significance on the treatment of the downstream of the yellow river, and is the key for constructing a high-efficiency sand transporting channel at the downstream of the yellow river and ensuring that the river channel at the downstream of the yellow river does not silt. The characteristics of different river sections of the river course at the downstream of the yellow river have obvious differences, the sand transportation water volumes of the river sections with different characteristics have larger differences, and the influence of the morphological characteristics of the cross section of the river course on the sand transportation water volume of the river course is lack of research in the past achievements.
In order to solve the problems, the invention deeply analyzes the relation between the average water depth and the flow of the river channel according to the combination of the actual measurement data and the theory, and establishes a method and a system for calculating the dynamic sand transportation water volume of the river channel. The method can be used for universally calculating the dynamic sand conveying capacity, has simple and clear steps, reliable results, simple and convenient calculation and easy operation, and is a simple and convenient method which is easy to master and use by basic-level science and technology workers.
Disclosure of Invention
The invention provides a method and a system for calculating the dynamic sand transportation water amount of a river channel,
the invention provides a method for calculating dynamic sand transporting water amount of a river channel, which comprises the following steps:
step 1: collecting basic data related to the river channel;
step 2: calculating the corresponding river channel section average water depth under the condition of actually measured flow according to the basic data;
and step 3: fitting a first relational expression of the average water depth of the river channel section and the river channel flow, and combining a preset formula to obtain an instantaneous river channel sand conveying capacity expression;
defining an uneven coefficient phi and a scouring parameter lambda in the process of introducing flow, and acquiring a dynamic sediment transport capacity calculation method of the river channel based on the expression of the instantaneous sediment transport capacity of the river channel;
meanwhile, acquiring measured data related to river channel erosion and deposition balance, and determining an uneven coefficient phi and an erosion and deposition parameter lambda of a required parameter flow process in a river channel dynamic sand transportation capacity calculation method based on the measured data;
and 4, step 4: establishing an incidence relation between the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount;
and 5: and determining the change rule of the non-uniform coefficient phi and the erosion parameter lambda in the flow process based on the basic data, obtaining a calculation method of the dynamic sediment transport water volume of the river channel based on the change rule, the association relation and the determined required parameters, obtaining the dynamic sediment transport water volume of the river channel according to the calculation method of the dynamic sediment transport water volume of the river channel, and outputting and displaying the dynamic sediment transport water volume.
In one possible implementation, collecting the basic data related to the river includes:
collecting actually measured water and sand process data of a river channel hydrological measurement station and actually measured topographic data of river channel hydrological sections at the downstream of the yellow river in the flood season and the non-flood season in the previous year;
wherein, the data in the actually measured water and sand process data comprises: the daily average flow rate of the river channel, the daily average sand conveying rate of the river channel and the median particle size of the river channel sediment;
wherein, after collecting the basic data related to the river channel, the method further comprises the following steps:
and actually measuring topographic data based on the hydrological section of the river channel, and calculating the river channel silt flushing amount of each river section in the river channel in different time periods by adopting a section method.
In a possible implementation manner, the step of calculating the average water depth of the river section under the measured flow condition according to the basic data includes:
acquiring actual measurement flood factors based on the basic data;
acquiring an actually measured section area and a river width related to the actually measured flow condition;
and calculating the corresponding river channel section average water depth based on the actually measured flood factors, the actually measured section area and the river width.
In a possible implementation manner, the specific steps of step 3 include:
fitting a first relational expression of the average water depth of the river cross section and the river flow, wherein the first relational expression is as follows: h alpha Qβ；
Wherein h represents the average water depth of the river cross section, and the unit is m; q represents the river flow in m3/s；
Combining a preset formula, wherein the preset formula is as follows, and acquiring an expression of the instantaneous sediment transport capacity of the river channel;
Qs＝QS；
Figure BDA0003074995070000031
Figure RE-GDA0003146384470000032
wherein Q issExpressing the sand conveying rate of the river channel, wherein the unit is kg/s; s is the sand content of river water flow, and the unit is kg/m3(ii) a U represents the average water flow velocity of the river cross section and is unit m/s; n represents the roughness coefficient of the river channel; j represents the hydraulic ratio drop of the river channel; g is gravity acceleration, and is 9.81m2S; omega represents mud of river channelThe average settling velocity of the sand is m/s; k is the water flow sand-carrying force coefficient, and m is the water flow sand-carrying force index;
defining an uneven coefficient phi and a scouring and silting parameter lambda in the process of introducing flow as follows, and obtaining a river channel dynamic sand conveying capacity calculation method based on the river channel instantaneous sand conveying capacity expression, wherein the formula is as follows:
Figure BDA0003074995070000033
Figure RE-GDA0003146384470000034
in the above formula:
Figure BDA0003074995070000035
is the average flow over a period of time, in m3/s；ΔtiIs the time interval in units of s; and N is the number of time segments. Meanwhile, actual measurement data related to the silt flushing balance of the silt channel is obtained, and required parameters K, m in the method for calculating the dynamic sediment transport capacity of the silt channel are determined based on the actual measurement data.
In a possible implementation manner, the correlation between the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount is established based on the following formula:
λ(Ws-Ws *)＝ΔWs
wherein, WsRepresenting the amount of sand coming from the river channel, and the unit is kg; ws *Expressing the dynamic sand conveying capacity in kg; Δ WsThe unit of the river channel silt is kg.
In one possible implementation manner, the method further includes: the method for calculating the dynamic sediment transport water volume of the river channel is calculated according to the following formula:
Figure BDA0003074995070000041
wherein, WwShows the dynamic sand transportation water quantity m of the river3。
In a possible implementation manner, in step 1, before collecting the basic data related to the river, the method further includes:
determining historical acquisition information of each equipment acquisition interface on acquisition equipment for acquiring basic data, wherein the historical acquisition information comprises: the interface acquisition instruction and the corresponding interface acquisition data;
determining the acquisition attribute of the acquisition interface of the corresponding acquisition equipment according to the historical acquisition information and factory setting information of the acquisition interface of the equipment;
meanwhile, according to the historical acquisition information, determining first interconnectivity among different equipment acquisition interfaces, transmitting an interface acquisition instruction corresponding to the historical acquisition information of the first interface to a second interface, and acquiring data by the second interface according to the corresponding interface acquisition instruction to obtain data to be compared;
acquiring the difference between the data to be compared and the corresponding interface acquisition data;
determining a second interconnectivity of the second interface and the first interface according to the acquisition attributes of the second interface and the acquired differences;
selecting a standby interface from the equipment acquisition interfaces according to the first interconnectivity and the second interconnectivity;
collecting basic data related to a river channel from a main acquisition interface in the process of collecting the basic data related to the river channel, and simultaneously monitoring the main acquisition interface in real time;
and judging the reliability of the data collected by the main collection interface according to the real-time monitoring result so as to determine whether to switch to the standby interface for data collection.
In a possible implementation manner, the determining, according to a real-time monitoring result, reliability of data collected by the main acquisition interface includes:
determining a first data set and a second data set acquired by the main acquisition interface;
acquiring a first sub-item of the first data set, and inputting the first sub-item into a comparison unit table for first calibration;
acquiring a second sub-item of the second data set, and inputting the second sub-item into a comparison unit table for second calibration;
comparing the same sub-items in the first calibration result and the second calibration result one by one, and determining difference information between the same sub-items;
meanwhile, determining the actual acquisition period of each second sub-item, and judging whether each actual acquisition period is consistent with the standard acquisition period of the corresponding first sub-item;
if the data are consistent, the main acquisition interface is judged to be reliable;
if the two are inconsistent, judging that the main acquisition interface is unreliable;
if partial inconsistency exists, extracting the sub-collected data of the inconsistent second sub-item, retrieving the item details related to the inconsistent second sub-item from a storage database, and extracting the abnormal variable of the inconsistent second sub-item based on the item details;
calculating a comprehensive abnormal value Z according to all the extracted abnormal variables and the following formula;
Figure BDA0003074995070000051
wherein Z1 represents the number of inconsistent second sub-items and has a value range of [1, Z11 ]](ii) a Z2 represents the number of abnormal variables in each inconsistent second sub-item and has a value in the range of [1, Z22 ]]；Xz1(fz2) The z2 abnormal variable f in the z1 th inconsistent second sub-itemz2Comparing the outliers corresponding to the standard variables; exp () represents an exponential function; zetaz1The abnormal adjustment coefficient of the second sub item inconsistent with the z1 th item is represented, and the value range is (0.6, 1)]；ζz2The abnormal adjustment coefficient of the z2 th abnormal variable is shown, and the value range is (0.4, 0.8)]；ιz1z2Z2 abnormal changes in the z1 th inconsistent second sub-itemThe weight value of the quantity is (0, 1)]；
And when the comprehensive abnormal value Z is smaller than a preset threshold value, judging that the main acquisition interface is reliable, otherwise, judging that the main acquisition interface is unreliable, and automatically switching to the standby interface for data collection.
The invention provides a river course dynamic sediment outflow calculation system, which comprises:
the collection module is used for collecting basic data related to the river channel;
the calculation module is used for calculating the corresponding river channel section average water depth under the condition of actually measured flow according to the basic data;
the first acquisition module is used for fitting a first relational expression of the average water depth of the river channel section and the river channel flow and acquiring an instantaneous river channel sand transportation capacity expression by combining a preset formula;
defining an uneven coefficient phi and a scouring parameter lambda in the process of introducing flow, and acquiring a dynamic sediment transport capacity calculation method of the river channel based on the expression of the instantaneous sediment transport capacity of the river channel;
meanwhile, acquiring measured data related to river channel erosion and deposition balance, and determining an uneven coefficient phi and an erosion and deposition parameter lambda of a required flow process in a river channel dynamic sand transportation capacity calculation method based on the measured data; the establishment module is used for establishing the correlation among the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount;
and the second acquisition module is used for determining the change rule of the non-uniform coefficient phi and the parameter lambda in the flow process based on the basic data, obtaining a calculation method of the dynamic sediment transport water volume of the river channel based on the change rule, the association relation and the determined required parameter, acquiring the dynamic sediment transport water volume of the river channel according to the calculation method of the dynamic sediment transport water volume of the river channel, and outputting and displaying the acquired dynamic sediment transport water volume.
Compared with the prior art, the invention has the following advantages and beneficial effects:
the invention provides a dynamic sediment transport water quantity expression of each river section at the downstream of the yellow river by collecting actually measured water and sand data of a representative hydrological station of the river channel and the topography change condition of the section of the downstream river channel, deeply analyzing the relation between the average water depth and the flow of the river channel and combining with a sediment transport theory. Compared with the prior art, the method has the advantages of simple and clear steps, reliable results, simple and convenient calculation and easy operation, and is a simple and convenient method which is easy to master and use by basic-level science and technology workers.
Additional features and advantages of the invention will be set forth in the description which follows, and in part will be obvious from the description, or may be learned by practice of the invention. The objectives and other advantages of the invention will be realized and attained by the structure particularly pointed out in the written description and claims hereof as well as the appended drawings.
The technical solution of the present invention is further described in detail by the accompanying drawings and embodiments.
Drawings
The accompanying drawings, which are included to provide a further understanding of the invention and are incorporated in and constitute a part of this specification, illustrate embodiments of the invention and together with the description serve to explain the principles of the invention and not to limit the invention. In the drawings:
FIG. 1 is a flow chart of a method for calculating the amount of dynamically transported sand water in a river according to the present invention;
FIG. 2 is a graph of the average water depth h of a high village and the flow Q;
FIG. 3 shows the dynamic sand-transporting capacity W of the high river reach of the inventions *And
Figure BDA0003074995070000071
a relationship diagram of (1);
FIG. 4 is a graph showing the relationship between the average water depth h and the flow Q of the moxa cone according to the present invention;
FIG. 5 shows the dynamic sand transportation capacity W of the high-Artemisia river reach of the inventions *And
Figure BDA0003074995070000072
a relationship diagram of (1);
FIG. 6 is a graph showing the relationship between the average water depth h and the flow Q of the Airy river reach according to the present invention;
FIG. 7 shows the dynamic sand transportation ability W of the Aili river reach of the present inventions *And
Figure BDA0003074995070000073
a relationship diagram of (1);
fig. 8 is a structural diagram of a system for calculating a dynamic sand transportation amount of a river according to an embodiment of the present invention.
Detailed Description
The preferred embodiments of the present invention will be described in conjunction with the accompanying drawings, and it will be understood that they are described herein for the purpose of illustration and explanation and not limitation.
The invention provides a method for calculating the dynamic sand transportation water amount of a river channel, which comprises the following steps of:
step 1: collecting basic data related to the river channel;
step 2: calculating the corresponding river channel section average water depth under the condition of actually measured flow according to the basic data;
and step 3: fitting a first relational expression of the average water depth of the river channel section and the river channel flow, and combining a preset formula to obtain an instantaneous river channel sand conveying capacity expression;
defining an uneven coefficient phi and a scouring parameter lambda in the process of introducing flow, and acquiring a dynamic sediment transport capacity calculation method of the river channel based on the expression of the instantaneous sediment transport capacity of the river channel;
meanwhile, acquiring measured data related to river channel erosion and deposition balance, and determining an uneven coefficient phi and an erosion and deposition parameter lambda of a required parameter flow process in a river channel dynamic sand transportation capacity calculation method based on the measured data;
and 4, step 4: establishing an incidence relation between the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount;
and 5: and determining the change rule of the non-uniform coefficient phi and the erosion parameter lambda in the flow process based on the basic data, obtaining a calculation method of the dynamic sediment transport water volume of the river channel based on the change rule, the association relation and the determined required parameters, obtaining the dynamic sediment transport water volume of the river channel according to the calculation method of the dynamic sediment transport water volume of the river channel, and outputting and displaying the dynamic sediment transport water volume.
The beneficial effects of the above technical scheme are: the invention provides a dynamic sediment transport water quantity expression of each river section at the downstream of the yellow river by collecting actually measured water and sand data of a representative hydrological station of the river channel and the topography change condition of the section of the downstream river channel, deeply analyzing the relation between the average water depth and the flow of the river channel and combining with a sediment transport theory. Compared with the prior art, the method has the advantages of simple and clear steps, reliable results, simple and convenient calculation and easy operation, and is a simple and convenient method which is easy to master and use by basic-level science and technology workers.
For the above 1-5 steps, see in particular the following examples:
calculating the dynamic sediment transport water volume of the river channel from the garden mouth of the downstream of the yellow river to the river reach of the high village (called the flower high river reach for short);
and collecting water and sand data of the hydrological survey station and terrain change of the section of river channel by taking the hydrological station in the high village as the representative station of the section of river channel. Calculating the corresponding flow Q (m) according to the actual flood element data3Calculating the corresponding average water depth h (m) of the representative section by the actually measured section area and the river width B (m) under the condition of/s), drawing a point into a graph by taking the flow Q as an abscissa and the average water depth h as an ordinate, and fitting a power function relation expression h ═ alpha Q of the representative section average water depth h and the flow QβReferring to fig. 2, according to the fitting result, α is 0.141 and β is 0.350, i.e., h is 0.141Q0.350。
Using the above relation formula and sand conveying rate formula QsQS, manning formula
Figure RE-GDA0003146384470000081
Substituting zuelai sand-holding force formula
Figure RE-GDA0003146384470000082
And arranging to obtain an expression of the instantaneous sand conveying capacity of the high river reach
Figure RE-GDA0003146384470000083
Figure RE-GDA0003146384470000091
The flow process non-uniformity coefficient to be defined is
Figure RE-GDA0003146384470000092
The sand conveying capacity of the floricome river reach in a period of time can be obtained after the sand conveying capacity is put into an expression of the instantaneous sand conveying capacity of the floricome river reach and is finished
Figure RE-GDA0003146384470000093
Figure RE-GDA0003146384470000094
According to the topographic change of the river section, the year-to-year silt flushing amount of the high river section (hydrology) is counted, and the year of the river section close to the silt flushing balance is selected (1972, 1974, 1975, 1978, 1979, 1980, 1984, 1985, 1993, 1995, 1997, 1999).
And (3) according to the measured data of the years, calibrating parameters K, m in the dynamic sand-transporting capacity expression of the high river reach, wherein K is 0.033 and m is 0.92, so as to obtain the dynamic sand-transporting capacity expression of the high river reach:
Figure RE-GDA0003146384470000095
Figure RE-GDA0003146384470000096
according to the statistical analysis of the measured data, the non-uniform coefficient phi of the flow process of the high river reach is 34607591 in the average value for many years, the parameter lambda is 0.95, and then the dynamic sand water delivery quantity W of the high river reach is obtainedWIs expressed as
Figure BDA0003074995070000097
Calculating the dynamic sand conveying capacity of a river channel of a high village-ai mountain river section (a high ai river section for short) at the downstream of a yellow river;
and collecting water and sand data of the hydrological survey station and terrain changes of the section of river channel by taking the moxa-mountain hydrological station as the representative station of the section of river channel. Calculating the corresponding flow Q (m) according to the actual flood element data3The actually measured section area and the river width B (m) under the condition of/s) calculate the average water depth h (m) corresponding to the representative sectionTaking the flow Q as an abscissa and the average water depth h as an ordinate, after points are drawn into a graph, fitting a power function relation h which represents the average water depth h of a section and the flow Q and is alpha QβReferring to fig. 4, according to the fitting result, α is 0.099, β is 0.470, that is, h is 0.099Q0.470。
Using the above relation formula and sand conveying rate formula QsQS, manning formula
Figure RE-GDA0003146384470000101
Substituting zuelai sand-holding force formula
Figure RE-GDA0003146384470000102
And arranging to obtain an expression of the instantaneous sand transporting capacity of the high-moxa river reach
Figure RE-GDA0003146384470000103
Figure RE-GDA0003146384470000104
A process non-uniformity coefficient to be defined as
Figure RE-GDA0003146384470000105
The sand conveying capacity of the high-moxa river reach within a period of time can be obtained after the sand conveying capacity is put into the expression of the instantaneous sand conveying capacity of the high-moxa river reach and is sorted
Figure RE-GDA0003146384470000106
Figure RE-GDA0003146384470000107
Wherein the flow process non-uniformity coefficient is
Figure RE-GDA0003146384470000108
Figure RE-GDA0003146384470000109
According to the topographic change of the section of the river channel, the year-to-year silt-flushing amount of the high-altitude ai river reach (hydrology), the year of the section of the river channel close to the silt-flushing balance is selected (1972, 1974, 1975, 1976, 1977, 1978, 1979, 1980, 1981, 1982, 1983, 1986, 1987, 1988, 1989, 1991, 1992, 1993, 1994, 1995, 1996, 1997, 1998, 1999, 2000, 2001, 2002, 2003, 2008, 2009, 2010, 2011, 2012), and the parameter K, m in the expression of the dynamic sand-transporting capacity of the high-altitude ai river reach is determined according to the measured data of the year, wherein K is 0.0222, and m is 0.92 (see figure 5), namely the expression of the dynamic sand-transporting capacity of the high-altitude ai river reach:
Figure RE-GDA00031463844700001010
according to the statistical analysis of the measured data, the non-uniform coefficient phi of the high-moxa flow process is 36544904 in the average value for many years, and the parameter lambda is 0.95, so that the dynamic sand conveying water quantity W of the high-flower river reachWIs expressed as
Figure BDA00030749950700001011
Calculating the dynamic sediment transport capacity of a river channel of a downstream Aishan-Rijin river reach (Airy river reach for short);
and collecting water and sand data of the hydrological measuring station and terrain change of the river channel section by taking the Lijin hydrological station as the representative station of the river channel section. Calculating the corresponding flow Q (m) according to the actual flood element data3Calculating the corresponding average water depth h (m) of the representative section by the actually measured section area and the river width B (m) under the condition of/s), drawing a point into a graph by taking the flow Q as an abscissa and the average water depth h as an ordinate, and fitting a power function relation expression h ═ alpha Q of the representative section average water depth h and the flow QβReferring to fig. 6, according to the fitting result, α is 0.135, β is 0.413, that is, h is 0.135Q0.413。
Using the above relation formula and sand conveying rate formula QsQS, manning formula
Figure RE-GDA0003146384470000111
Substituting zuelai sand-holding force formula
Figure RE-GDA0003146384470000112
And arranging to obtain an expression of the instantaneous sand transporting capacity of the high-moxa river reach
Figure RE-GDA0003146384470000113
Figure RE-GDA0003146384470000114
The flow process non-uniformity coefficient to be defined is
Figure RE-GDA0003146384470000115
The sand conveying capacity of the high-moxa river reach within a period of time can be obtained after the sand conveying capacity is put into the expression of the instantaneous sand conveying capacity of the high-moxa river reach and is sorted
Figure RE-GDA0003146384470000116
According to the topographic change of the river section, the year-to-year silt flushing amount of the Gaotai river reach (hydrology year) is counted, the year in which the river section of the Gaotai river reach is close to the silt flushing balance is selected (1967, 1972, 1973, 1974, 1976, 1979, 1980, 1983, 1984, 1985, 1989, 1991, 1993, 1994, 1995, 1996, 1998, 1999, 2000, 2001, 2004, 2005, 2006, 2007, 2008, 2009, 2014, 2015, 2016, 2017 and 2018), parameters K, m in the expression of the dynamic sand transportation capacity of the Gaotai river reach are determined according to the measured data of the year, wherein K is 0.0343, and m is 0.92 (see figure 7), namely the expression of the dynamic sand transportation capacity of the Gaotai river reach:
Figure RE-GDA0003146384470000117
according to the statistical analysis of the measured data, the non-uniform coefficient phi of the flow process is 40013477 in the average value of many years, the parameter lambda is 0.95, and then the dynamic sand conveying water quantity W of the high river reachWIs expressed as
Figure BDA0003074995070000118
Figure BDA0003074995070000121
Through the embodiments, the dynamic sand conveying quantity can be calculated conveniently and effectively.
In one embodiment, step 1, before collecting the basic data related to the river, further includes:
determining historical acquisition information of each equipment acquisition interface on acquisition equipment for acquiring basic data, wherein the historical acquisition information comprises: the interface acquisition instruction and the corresponding interface acquisition data;
determining the acquisition attribute of the acquisition interface of the corresponding acquisition equipment according to the historical acquisition information and factory setting information of the acquisition interface of the equipment;
meanwhile, according to the historical acquisition information, determining first interconnectivity among different equipment acquisition interfaces, transmitting an interface acquisition instruction corresponding to the historical acquisition information of the first interface to a second interface, and acquiring data by the second interface according to the corresponding interface acquisition instruction to obtain data to be compared;
acquiring the difference between the data to be compared and the corresponding interface acquisition data;
determining a second interconnectivity of the second interface and the first interface according to the acquisition attributes of the second interface and the acquired differences;
selecting a standby interface from the equipment acquisition interfaces according to the first interconnectivity and the second interconnectivity;
collecting basic data related to a river channel from a main acquisition interface in the process of collecting the basic data related to the river channel, and simultaneously monitoring the main acquisition interface in real time;
and judging the reliability of the data collected by the main collection interface according to the real-time monitoring result so as to determine whether to switch to the standby interface for data collection.
In the embodiment, since the collected data is an application basis of subsequent mathematical computation, in the embodiment, a series of judgments are performed on the device acquisition interface for collecting the data, and finally, the standby interface is screened out, so that the validity of data collection is ensured, and the possibility of data computation errors caused by incomplete data collection is avoided.
In this embodiment, the interface acquisition data is also associated with the underlying data.
In this embodiment, factory setting information, such as the type of data collected, data that can be collected, and the like, is set.
In this embodiment, the collection attributes relate to data type, data volume, data category, and the like.
In this embodiment, the first interconnectivity is determined according to a data type and a data type, and is used for the interconnectivity between every two interfaces on the acquisition device.
In this embodiment, the first interface is any one of all the device acquisition interfaces in the acquisition device, and the second interface excludes all the remaining interfaces of the first interface from all the corresponding device acquisition interfaces;
in this embodiment, the data to be compared is obtained by controlling the second interface to perform data acquisition according to the interface acquisition instruction of the first interface.
In this embodiment, the difference refers to difference information between the data to be compared and interface acquisition data acquired by the corresponding first interface according to the interface acquisition instruction.
In this embodiment, the second association refers to an association condition between the first interface and all the second interfaces.
The beneficial effects of the above technical scheme are: the standby interfaces are screened out by judging a series of judgment on the acquisition interfaces of the data collection equipment and determining the first relevance and the second relevance, so that the effectiveness of data collection is ensured, and the possibility of data calculation errors caused by incomplete data collection is avoided.
In one embodiment, determining the reliability of the data collected by the main collection interface according to the real-time monitoring result includes:
determining a first data set and a second data set acquired by the main acquisition interface;
acquiring a first sub-item of the first data set, and inputting the first sub-item into a comparison unit table for first calibration;
acquiring a second sub-item of the second data set, and inputting the second sub-item into a comparison unit table for second calibration;
comparing the same sub-items in the first calibration result and the second calibration result one by one, and determining difference information between the same sub-items;
meanwhile, determining the actual acquisition period of each second sub-item, and judging whether each actual acquisition period is consistent with the standard acquisition period of the corresponding first sub-item;
if the data are consistent, the main acquisition interface is judged to be reliable;
if the two are inconsistent, judging that the main acquisition interface is unreliable;
if partial inconsistency exists, extracting the sub-collected data of the inconsistent second sub-item based on the difference information, retrieving the item details related to the inconsistent second sub-item from a storage database, and extracting the abnormal variable of the inconsistent second sub-item based on the item details;
calculating a comprehensive abnormal value Z according to all the extracted abnormal variables and the following formula;
Figure BDA0003074995070000141
wherein Z1 represents the number of inconsistent second sub-items and has a value range of [1, Z11 ]](ii) a Z2 represents the number of abnormal variables in each inconsistent second sub-item and has a value in the range of [1, Z22 ]]；Xz1(fz2) The z2 abnormal variable f in the z1 th inconsistent second sub-itemz2Comparing the outliers corresponding to the standard variables; exp () represents an exponential function; zetaz1The abnormal adjustment coefficient of the second sub item inconsistent with the z1 th item is represented, and the value range is (0.6, 1)]；ζz2The abnormal adjustment coefficient of the z2 th abnormal variable is shown, and the value range is (0.4, 0.8)]；ιz1z2The weight value of the z2 th abnormal variable in the z1 th inconsistent second sub-item is represented, and the value range is (0),1]；
And when the comprehensive abnormal value Z is smaller than a preset threshold value, judging that the main acquisition interface is reliable, otherwise, judging that the main acquisition interface is unreliable, and automatically switching to the standby interface for data collection.
In this embodiment, the monitoring result is formed by the second data set, and the first data set is formed by collecting data when the primary collection interface is qualified.
In this embodiment, the comparison unit table, the first sub-item, and the second sub-item are all related to the acquired objects and parameter variables, for example, the objects in the sub-items are water and sand process data, the flood season in the past year, and the hydrological section of the yellow river downstream of the non-flood season, and the corresponding parameter variables are topographic data, the daily average flow rate of the river, the daily average sand rate of the river, and the median diameter of the sediment in the river.
In this embodiment, the sub-collected data of the inconsistent second sub-item is, for example, that some parameter variables or missing of the sub-item exist in the collection process, and the abnormal variable may also be regarded as missing of some parameters or data abnormality, etc.
In this embodiment, item details, such as an extra sub-item or a missing sub-item, may be considered item details.
The beneficial effects of the above technical scheme are: the method comprises the steps of inputting different data sets of a main acquisition interface into a comparison unit table for calibration, comparing the same sub-items one by one, and improving the reliability of comparison.
The invention provides a river channel dynamic sediment outflow calculation system, as shown in fig. 8, comprising:
the collection module is used for collecting basic data related to the river channel;
the calculation module is used for calculating the corresponding river channel section average water depth under the condition of actually measured flow according to the basic data;
the first acquisition module is used for fitting a first relational expression of the average water depth of the river channel section and the river channel flow and acquiring an instantaneous river channel sand transportation capacity expression by combining a preset formula;
defining an uneven coefficient phi and a scouring parameter lambda in the process of introducing flow, and acquiring a dynamic sediment transport capacity calculation method of the river channel based on the expression of the instantaneous sediment transport capacity of the river channel;
meanwhile, acquiring measured data related to river channel erosion and deposition balance, and determining an uneven coefficient phi and an erosion and deposition parameter lambda of a required flow process in a river channel dynamic sand transportation capacity calculation method based on the measured data;
the establishment module is used for establishing the correlation among the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount;
and the second acquisition module is used for determining the change rule of the non-uniform coefficient phi and the parameter lambda in the flow process based on the basic data, obtaining a calculation method of the dynamic sediment transport water volume of the river channel based on the change rule, the association relation and the determined required parameter, acquiring the dynamic sediment transport water volume of the river channel according to the calculation method of the dynamic sediment transport water volume of the river channel, and outputting and displaying the acquired dynamic sediment transport water volume.
The beneficial effects of the above technical scheme are: the invention provides a dynamic sediment transport water quantity expression of each river section at the downstream of the yellow river by collecting actually measured water and sand data of a representative hydrological station of the river channel and the topography change condition of the section of the downstream river channel, deeply analyzing the relation between the average water depth and the flow of the river channel and combining with a sediment transport theory. Compared with the prior art, the method has the advantages of simple and clear steps, reliable results, simple and convenient calculation and easy operation, and is a simple and convenient method which is easy to master and use by basic-level science and technology workers.
It will be apparent to those skilled in the art that various changes and modifications may be made in the present invention without departing from the spirit and scope of the invention. Thus, if such modifications and variations of the present invention fall within the scope of the claims of the present invention and their equivalents, the present invention is also intended to include such modifications and variations.

Claims (9)

1. A method for calculating the dynamic sand transporting water amount of a river channel is characterized by comprising the following steps:
step 1: collecting basic data related to the river channel;
step 2: calculating the corresponding river channel section average water depth under the condition of actually measured flow according to the basic data;
and step 3: fitting a first relational expression of the average water depth of the river channel section and the river channel flow, and combining a preset formula to obtain an instantaneous river channel sand conveying capacity expression;
defining an uneven coefficient phi and a scouring parameter lambda in the process of introducing flow, and acquiring a dynamic sediment transport capacity calculation method of the river channel based on the expression of the instantaneous sediment transport capacity of the river channel;
meanwhile, acquiring measured data related to river channel erosion and deposition balance, and determining an uneven coefficient phi and an erosion and deposition parameter lambda of a required parameter flow process in a river channel dynamic sand transportation capacity calculation method based on the measured data;
and 4, step 4: establishing an incidence relation between the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount;
and 5: and determining the change rule of the non-uniform coefficient phi and the erosion parameter lambda in the flow process based on the basic data, obtaining a calculation method of the dynamic sediment transport water volume of the river channel based on the change rule, the association relation and the determined required parameters, obtaining the dynamic sediment transport water volume of the river channel according to the calculation method of the dynamic sediment transport water volume of the river channel, and outputting and displaying the dynamic sediment transport water volume.
2. The computing method of claim 1, wherein collecting the channel-related basic data comprises:
collecting actually measured water and sand process data of a river channel hydrological measurement station and actually measured topographic data of river channel hydrological sections at the downstream of the yellow river in the flood season and the non-flood season in the previous year;
wherein, the data in the actually measured water and sand process data comprises: the daily average flow rate of the river channel, the daily average sand conveying rate of the river channel and the median particle size of the river channel sediment;
wherein, after collecting the basic data related to the river channel, the method further comprises the following steps:
and actually measuring topographic data based on the hydrological section of the river channel, and calculating the river channel silt flushing amount of each river section in the river channel in different time periods by adopting a section method.
3. The method of claim 1, wherein the step of calculating the corresponding average water depth of the river cross section under the measured flow rate condition according to the basic data comprises:
acquiring actual measurement flood factors based on the basic data;
acquiring an actually measured section area and a river width related to the actually measured flow condition;
and calculating the corresponding river channel section average water depth based on the actually measured flood factors, the actually measured section area and the river width.
4. The computing method according to claim 1, wherein the specific steps of step 3 include:
fitting a first relational expression of the average water depth of the river cross section and the river flow, wherein the first relational expression is as follows: h ═ α Qβ；
Wherein h represents the average water depth of the river cross section, and the unit is m; q represents the river flow in m3S; alpha and beta are parameters;
combining a preset formula, wherein the preset formula is as follows, and acquiring an expression of the instantaneous sediment transport capacity of the river channel; qs＝QS；
Figure RE-FDA0003146384460000021
Figure RE-FDA0003146384460000022
Wherein Q issExpressing the sand conveying rate of the river channel, wherein the unit is kg/s; s is the sand content of river water flow, and the unit is kg/m3(ii) a U represents the average water flow velocity of the river cross section and is unit m/s; n represents the roughness coefficient of the river channel; j represents the hydraulic ratio drop of the river channel; g is gravity acceleration, and is 9.81m2S; omega represents the average settling velocity of the silt in the riverway, and the unit is m/s; k is the water flow sand-carrying force coefficient, and m is the water flow sand-carrying force index;
defining an uneven coefficient phi and a scouring and silting parameter lambda in the process of introducing flow as follows, and obtaining a river channel dynamic sand conveying capacity calculation method based on the river channel instantaneous sand conveying capacity expression, wherein the formula is as follows:
Figure RE-FDA0003146384460000023
Figure RE-FDA0003146384460000031
in the above formula:
Figure RE-FDA0003146384460000032
is the average flow over a period of time, in m3/s；ΔtiIs the time interval in units of s; and N is the number of time segments. And meanwhile, acquiring measured data related to the erosion and deposition balance of the river channel, and determining a water flow sand-entrainment force coefficient K and a water flow sand-entrainment force index m in the calculation method of the dynamic sand transport capacity of the river channel based on the measured data.
5. The calculation method according to claim 1, wherein the correlation between the amount of sand coming from the river, the dynamic sand transport capacity and the amount of silt flowing from the river is established based on the following formula:
λ(Ws—Ws *)＝ΔWs
wherein, WsRepresenting the amount of sand coming from the river channel, and the unit is kg; ws *Expressing the dynamic sand conveying capacity in kg; Δ WsThe unit of the river channel silt is kg.
6. The computing method of claim 1, further comprising: the method for calculating the dynamic sediment transport water volume of the river channel is calculated according to the following formula:
Figure FDA0003074995060000033
wherein, WwExpresses the dynamic sand transportation water quantity of the river channel, and the unit is m3。
7. The computing method of claim 1, wherein step 1, before collecting the basic data related to the river, further comprises:
determining historical acquisition information of each equipment acquisition interface on acquisition equipment for acquiring basic data, wherein the historical acquisition information comprises: the interface acquisition instruction and the corresponding interface acquisition data;
determining the acquisition attribute of the acquisition interface of the corresponding acquisition equipment according to the historical acquisition information and factory setting information of the acquisition interface of the equipment;
meanwhile, according to the historical acquisition information, determining first interconnectivity among different equipment acquisition interfaces, transmitting an interface acquisition instruction corresponding to the historical acquisition information of the first interface to a second interface, and acquiring data by the second interface according to the corresponding interface acquisition instruction to obtain data to be compared;
acquiring the difference between the data to be compared and the corresponding interface acquisition data;
determining a second interconnectivity of the second interface and the first interface according to the acquisition attributes of the second interface and the acquired differences;
selecting a standby interface from the equipment acquisition interfaces according to the first interconnectivity and the second interconnectivity;
collecting basic data related to a river channel from a main acquisition interface in the process of collecting the basic data related to the river channel, and simultaneously monitoring the main acquisition interface in real time;
and judging the reliability of the data collected by the main collection interface according to the real-time monitoring result so as to determine whether to switch to the standby interface for data collection.
8. The computing method of claim 7, wherein determining the reliability of the data collected by the primary collection interface based on the real-time monitoring comprises:
determining a first data set and a second data set acquired by the main acquisition interface;
acquiring a first sub-item of the first data set, and inputting the first sub-item into a comparison unit table for first calibration;
acquiring a second sub-item of the second data set, and inputting the second sub-item into a comparison unit table for second calibration;
comparing the same sub-items in the first calibration result and the second calibration result one by one, and determining difference information between the same sub-items;
meanwhile, determining the actual acquisition period of each second sub-item, and judging whether each actual acquisition period is consistent with the standard acquisition period of the corresponding first sub-item;
if the data are consistent, the main acquisition interface is judged to be reliable;
if the two are inconsistent, judging that the main acquisition interface is unreliable;
if partial inconsistency exists, extracting the sub-collected data of the inconsistent second sub-item, retrieving the item details related to the inconsistent second sub-item from a storage database, and extracting the abnormal variable of the inconsistent second sub-item based on the item details;
calculating a comprehensive abnormal value Z according to all the extracted abnormal variables and the following formula;
Figure FDA0003074995060000051
wherein Z1 represents the number of inconsistent second sub-items and has a value range of [1, Z11 ]](ii) a z2 represents eachThe number of abnormal variables in the inconsistent second sub-item is in the range of [1, Z22 ]]；Xz1(fz2) The z2 abnormal variable f in the z1 th inconsistent second sub-itemz2Comparing the outliers corresponding to the standard variables; exp () represents an exponential function; zetaz1The abnormal adjustment coefficient of the second sub item inconsistent with the z1 th item is represented, and the value range is (0.6, 1)]；ζz2The abnormal adjustment coefficient of the z2 th abnormal variable is shown, and the value range is (0.4, 0.8)]；ιz1z2The weight value of the z2 th abnormal variable in the z1 th inconsistent second sub-item is represented, and the value range is (0, 1)]；
And when the comprehensive abnormal value Z is smaller than a preset threshold value, judging that the main acquisition interface is reliable, otherwise, judging that the main acquisition interface is unreliable, and automatically switching to the standby interface for data collection.
9. A riverway dynamic sediment transport water volume calculation system is characterized by comprising:
the collection module is used for collecting basic data related to the river channel;
the calculation module is used for calculating the corresponding river channel section average water depth under the condition of actually measured flow according to the basic data;
the first acquisition module is used for fitting a first relational expression of the average water depth of the river channel section and the river channel flow and acquiring an instantaneous river channel sand transportation capacity expression by combining a preset formula;
defining an uneven coefficient phi and a parameter lambda in an inflow process, and acquiring a river channel dynamic sand conveying capacity calculation method based on the river channel instantaneous sand conveying capacity expression;
meanwhile, acquiring measured data related to river channel erosion and deposition balance, and determining an uneven coefficient phi and an erosion and deposition parameter lambda of a required flow process in a river channel dynamic sand transportation capacity calculation method based on the measured data; the establishment module is used for establishing the correlation among the river channel sand inflow amount, the dynamic sand conveying capacity and the river channel silt flushing amount;
and the second acquisition module is used for determining the change rule of the non-uniform coefficient phi and the erosion parameter lambda in the flow process based on the basic data, obtaining a calculation method of the dynamic sediment transport water amount of the river channel based on the change rule, the association relation and the determined required parameters, acquiring the dynamic sediment transport water amount of the river channel according to the calculation method of the dynamic sediment transport water amount of the river channel, and outputting and displaying the acquired dynamic sediment transport water amount of the river channel.
CN202110549830.XA 2021-05-20 2021-05-20 Method and system for calculating dynamic sediment transport water amount of river channel Active CN113221366B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110549830.XA CN113221366B (en) 2021-05-20 2021-05-20 Method and system for calculating dynamic sediment transport water amount of river channel

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110549830.XA CN113221366B (en) 2021-05-20 2021-05-20 Method and system for calculating dynamic sediment transport water amount of river channel

Publications (2)

Publication Number Publication Date
CN113221366A true CN113221366A (en) 2021-08-06
CN113221366B CN113221366B (en) 2023-01-31

Family

ID=77093764

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110549830.XA Active CN113221366B (en) 2021-05-20 2021-05-20 Method and system for calculating dynamic sediment transport water amount of river channel

Country Status (1)

Country Link
CN (1) CN113221366B (en)

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113737710A (en) * 2021-08-12 2021-12-03 中国长江三峡集团有限公司 Estimation method of sediment transport rate of natural river bed load

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136636A (en) * 2013-03-15 2013-06-05 中广核工程有限公司 Nuclear power project progress measuring method and system
CN108629055A (en) * 2017-03-16 2018-10-09 中交上海航道勘察设计研究院有限公司 A kind of sandy cruiseway silt amount forecasting procedure based on saturation sediment transport principle
US20190354873A1 (en) * 2018-02-16 2019-11-21 Lucas Pescarmona Analysis system and hydrology management for basin rivers
AU2020101063A4 (en) * 2020-06-19 2020-07-23 Jilin Province Water Resource and Hydropower Consultative Company of P.R.CHINA A Calculation Method of River Ecological Flow
CN112149314A (en) * 2020-10-16 2020-12-29 黄河勘测规划设计研究院有限公司 Multi-sand reservoir capacity silt flushing simulation method based on virtual reservoir capacity correction
CN112287528A (en) * 2020-10-16 2021-01-29 黄河勘测规划设计研究院有限公司 Flood control high water level determination method for sandy river reservoir based on high beach trough
CN112487732A (en) * 2020-11-18 2021-03-12 黄河勘测规划设计研究院有限公司 Novel method and system for simulating erosion-deposition river channel
CN112595489A (en) * 2020-11-26 2021-04-02 黄河勘测规划设计研究院有限公司 Method and system for calculating dynamic sand transporting capacity of river channel

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103136636A (en) * 2013-03-15 2013-06-05 中广核工程有限公司 Nuclear power project progress measuring method and system
CN108629055A (en) * 2017-03-16 2018-10-09 中交上海航道勘察设计研究院有限公司 A kind of sandy cruiseway silt amount forecasting procedure based on saturation sediment transport principle
US20190354873A1 (en) * 2018-02-16 2019-11-21 Lucas Pescarmona Analysis system and hydrology management for basin rivers
AU2020101063A4 (en) * 2020-06-19 2020-07-23 Jilin Province Water Resource and Hydropower Consultative Company of P.R.CHINA A Calculation Method of River Ecological Flow
CN112149314A (en) * 2020-10-16 2020-12-29 黄河勘测规划设计研究院有限公司 Multi-sand reservoir capacity silt flushing simulation method based on virtual reservoir capacity correction
CN112287528A (en) * 2020-10-16 2021-01-29 黄河勘测规划设计研究院有限公司 Flood control high water level determination method for sandy river reservoir based on high beach trough
CN112487732A (en) * 2020-11-18 2021-03-12 黄河勘测规划设计研究院有限公司 Novel method and system for simulating erosion-deposition river channel
CN112595489A (en) * 2020-11-26 2021-04-02 黄河勘测规划设计研究院有限公司 Method and system for calculating dynamic sand transporting capacity of river channel

Non-Patent Citations (4)

* Cited by examiner, † Cited by third party
Title
SHIFENG HUANG等: "Dynamic monitoring of Yellow River estuary based on remote sensing in the recent ten years and analysis of correlation with flow and sediment condition", 《IEEE XPLORE》 *
吴保生等: "黄河内蒙古河段输沙量与淤积量计算方法", 《水科学进展》 *
杨丽丰等: "渭河下游输沙用水量研究", 《泥沙研究》 *
马睿等: "黄河流域典型沙质河段冲淤量预估方法及应用", 《水利学报》 *

Cited By (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113737710A (en) * 2021-08-12 2021-12-03 中国长江三峡集团有限公司 Estimation method of sediment transport rate of natural river bed load

Also Published As

Publication number Publication date
CN113221366B (en) 2023-01-31

Similar Documents

Publication Publication Date Title
CN105042339B (en) One kind is based on nondimensional leakage of finished oil pipeline amount estimating system and method
Di Stefano et al. Flow resistance equation for rills
CN109902111B (en) Earth and rockfill dam safety information mining method
CN108256177B (en) Parameter optimization method and system for river water sand model
Gabrielle et al. Analysis and field evaluation of the CERES models water balance component
CN112595489B (en) Method and system for calculating dynamic sand transporting capacity of river channel
CN113095694B (en) Rainfall sand transportation model construction method suitable for multiple landform type areas
CN110929359A (en) Pipe network siltation risk prediction modeling method based on PNN neural network and SWMM technology
CN105809249B (en) A kind of PM2.5 Concentration Testings based on amphineura network and forecasting system and method
CN111310976B (en) Quick water quality forecasting and early warning method and system based on distributed computing
CN113221366B (en) Method and system for calculating dynamic sediment transport water amount of river channel
CN114997548B (en) Water movement geological seepage real-time early warning method based on genetic neural network
CN105486358A (en) Gas-liquid two-phase flow parameter measuring method based on double-differential pressure of Venturi tube
CN110160435A (en) The moisture content increment of load response ratio Prediction Parameters that come down and method
WO2022257661A1 (en) Data processing method of dredging and transport system on long-distance pipeline transport site
CN101673096B (en) Soft-measuring method for density in concentration process of salvia miltiorrhiza injection production
CN109783934A (en) A kind of mean velocity in section fitting rating method based on H-ADCP
CN115587546A (en) Method for identifying flood level variation trend and cause of multiple river branches under dam
CN109299830A (en) A kind of hydrologic forecast system and method based on Soil Moisture Monitoring
CN106320257B (en) Method is determined based on the lake and reservoir channel storage curve of hydrological observation
CN113361130B (en) Construction method of high-sand-content flood post-river reach bed form adjustment prediction model
CN107507396A (en) A kind of method for early warning of the dangerous three-dimensional multi-point multi objective of rain-induced landslide
CN112818495A (en) Novel dynamic correction method for pipeline pressure drop measurement and calculation algorithm parameters
CN116541678B (en) Pressure monitoring method and device for gas station safety pipeline
CN107655593A (en) A kind of intelligent calorimeter is in Line synthesis error calibration method

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