CN116167242A - Hydrologic function communication assessment method suitable for multi-gate dam plain urban river network - Google Patents
Hydrologic function communication assessment method suitable for multi-gate dam plain urban river network Download PDFInfo
- Publication number
- CN116167242A CN116167242A CN202310311253.XA CN202310311253A CN116167242A CN 116167242 A CN116167242 A CN 116167242A CN 202310311253 A CN202310311253 A CN 202310311253A CN 116167242 A CN116167242 A CN 116167242A
- Authority
- CN
- China
- Prior art keywords
- hydrologic
- flow
- gate dam
- water
- index
- 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.)
- Pending
Links
- 238000004891 communication Methods 0.000 title claims abstract description 75
- 238000000034 method Methods 0.000 title claims abstract description 51
- 238000011156 evaluation Methods 0.000 claims abstract description 13
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims description 108
- 238000004590 computer program Methods 0.000 claims description 10
- 238000003860 storage Methods 0.000 claims description 9
- 238000005315 distribution function Methods 0.000 claims description 5
- 238000004422 calculation algorithm Methods 0.000 claims description 4
- 230000011218 segmentation Effects 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 9
- 238000011160 research Methods 0.000 description 7
- 238000009826 distribution Methods 0.000 description 5
- 239000002131 composite material Substances 0.000 description 4
- 230000001186 cumulative effect Effects 0.000 description 4
- 239000004576 sand Substances 0.000 description 4
- 238000013508 migration Methods 0.000 description 3
- 230000005012 migration Effects 0.000 description 3
- 238000004458 analytical method Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000013210 evaluation model Methods 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 238000011144 upstream manufacturing Methods 0.000 description 2
- 238000001276 Kolmogorov–Smirnov test Methods 0.000 description 1
- 238000009825 accumulation Methods 0.000 description 1
- 230000031018 biological processes and functions Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004140 cleaning Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000001419 dependent effect Effects 0.000 description 1
- 230000006866 deterioration Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000007613 environmental effect Effects 0.000 description 1
- 239000003673 groundwater Substances 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000011002 quantification Methods 0.000 description 1
- 230000001105 regulatory effect Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 238000011064 split stream procedure Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 239000000126 substance Substances 0.000 description 1
- 238000010998 test method Methods 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 238000012876 topography Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
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
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/10—Services
- G06Q50/26—Government or public services
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2111/00—Details relating to CAD techniques
- G06F2111/08—Probabilistic or stochastic CAD
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F2113/00—Details relating to the application field
- G06F2113/08—Fluids
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A10/00—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE at coastal zones; at river basins
- Y02A10/40—Controlling or monitoring, e.g. of flood or hurricane; Forecasting, e.g. risk assessment or mapping
Landscapes
- Business, Economics & Management (AREA)
- Engineering & Computer Science (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- Physics & Mathematics (AREA)
- Economics (AREA)
- General Physics & Mathematics (AREA)
- Tourism & Hospitality (AREA)
- Development Economics (AREA)
- Educational Administration (AREA)
- Strategic Management (AREA)
- General Business, Economics & Management (AREA)
- Marketing (AREA)
- Entrepreneurship & Innovation (AREA)
- Game Theory and Decision Science (AREA)
- Computer Hardware Design (AREA)
- Evolutionary Computation (AREA)
- Operations Research (AREA)
- Quality & Reliability (AREA)
- Geometry (AREA)
- General Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Primary Health Care (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
Aiming at the limitation of the prior art, the invention provides a hydrologic function communication evaluation method suitable for a multi-gate dam plain urban river network; according to the method, through acquiring the overflow area and the acquired flow data, the actual flow and the ecological base flow of the river basin are obtained, after the actual flow and the ecological base flow are converted into hydrologic flow indexes, the acquired hydrologic continuous indexes are combined, and hydrologic communication comprehensive indexes serving as hydrologic function communication evaluation indexes are output; dynamic changes of hydrologic function communication can be represented in real time, and more scientific and reasonable method support is provided for a gate dam dispatching scheme.
Description
Technical Field
The invention relates to the technical field of hydrologic communication assessment, in particular to a hydrologic function communication assessment method suitable for a multi-gate dam plain city river network.
Background
Connectivity is a key to research hydrology, topography, environment, ecology, flood, and has important influence on connection between water systems, water and sand migration, energy flow, biological migration and the like. Hydrologic communication research is receiving more and more attention because it is closely related to physical, chemical and biological processes of materials, energy and biology transfer and river systems.
From the meaning of hydrologic communication, hydrologic communication is divided into structural communication and functional communication: structural communication (static communication) refers to the spatial pattern of the river system and the degree of physical connection between rivers and lakes; functional communication (dynamic communication) refers to the degree to which different water systems are interconnected between river systems by some dynamic process. Wherein, the functional communication more embodies the essence of hydrologic communication. The current hydrologic function communication state is mainly characterized by using the gate dam communication degree, although the hydrologic function communication level of the water system can be primarily analyzed by using the gate dam communication degree, for example, in the chinese invention application of publication No. 2017.02.22: in the river network connectivity assessment and gate dam optimization method based on comprehensive resistance weights and graph theory, the thought of the comprehensive resistance weights and graph theory is combined, a river network gate dam system is generalized into a river network directed graph model which takes a gate dam as a vertex and takes a river section as a side, the comprehensive resistance weights covering hydrology, hydrodynamic force, water quality and hydraulic engineering items are calculated, and the comprehensive connectivity of the river network is obtained by combining a shortest path algorithm. However, the functional connectivity assessment model based on graph theory and hydrologic connectivity functions cannot obtain dynamic changes of river network hydrologic connectivity, such as connectivity changes caused by gate dam scheduling. In addition, the prior art has certain limitations because hydrologic data acquisition difficulty, hydrologic nonlinear variation, complexity of hydrologic scheduling rules and the like limit the quantification of hydrologic functional connectivity.
Disclosure of Invention
Aiming at the limitation of the prior art, the invention provides a hydrologic function communication assessment method applicable to a multi-gate dam plain city river network, which adopts the following technical scheme:
a hydrologic function communication assessment method suitable for a multi-gate dam plain city river network comprises the following steps:
s1, acquiring data about water level, flow and flow velocity acquired and recorded at a preset main gate dam point;
s2, river channel data of the main gate dam point are obtained, and the area of the main gate dam point on the water cross section, namely the overflow area, is obtained according to the river channel data;
s3, obtaining the actual flow and the ecological base flow of the main gate dam point on the water cross section according to the overflow area and the data about the flow and the flow speed obtained by collecting and recording;
s4, acquiring a hydrological flow index according to the actual flow of the main gate dam point on the water cross section and the ecological base flow;
s5, acquiring a hydrologic continuous index according to the data about the water level obtained by the acquisition record;
and S6, obtaining the hydrologic communication comprehensive index serving as the hydrologic function communication evaluation index according to the hydrologic flow index and the hydrologic continuous index.
Compared with the prior art that the hydrologic function communication state is characterized by mainly utilizing the gate dam communication degree, although the hydrologic function communication level of a water system can be primarily analyzed by utilizing the gate dam communication degree, the technical problem that the dynamic change of the hydrologic communication of a river network cannot be obtained by a functional communication evaluation model based on graph theory and hydrologic communication functions is solved, the actual flow and the ecological base flow of a river basin are obtained by obtaining the overflow area and the collected flow data, and after the actual flow and the ecological base flow are converted into hydrologic flow indexes, the collected hydrologic continuous indexes are combined, and the hydrologic communication comprehensive indexes serving as hydrologic function communication evaluation indexes are output; dynamic changes of hydrologic function communication can be represented in real time, and more scientific and reasonable method support is provided for a gate dam dispatching scheme.
As a preferred embodiment, the flow-through area is obtained by:
determining the starting point distances of the water edges of the left and right banks of the water section through a river channel midpoint distance algorithm, and extracting the water surface widths between the adjacent sounding plumb lines and the plumb lines according to the distance range of the starting points of the water edges; dividing the water cross section into a plurality of parts by adopting an average segmentation method and taking the sounding vertical lines as boundaries, wherein the distance between the adjacent vertical lines is a part width, and multiplying the average value of the water depths of the adjacent vertical lines to obtain the area of each part of the water cross section; and calculating the sum of the areas of all the parts to obtain the flow-through area.
As a preferred embodiment, the hydrological flow index Φ is obtained by the following formula:
wherein Q is 0 For the actual flow rate of the main gate dam point on the water cross section, Q e And eta is the dividing/converging ratio for the ecological base flow of the main gate dam point on the water cross section.
Further, the ecological base flow of the main gate dam point on the water cross section is obtained by the following steps:
based on the month average natural flow of a plurality of years, a P-III type curve is utilized, and the minimum month average natural flow under a preset frequency is selected as the ecological base flow.
Further, for the two branched river channels in the split state, in the case of diversion split in the split state, the score/confluence ratio η is obtained by:
Q 1 ＝Q 0 ；
Q 2 ＝(1-η)Q 0 ；
wherein Q is 0 For dry flow, Q 1 ，Q 2 Respectively, the tributary flows.
Further, for the two branched river channels in the split state, in the case that the split state is natural split, the score/confluence ratio η is obtained by:
wherein A is L1 And A L2 For the area of the tributary, θ 1 And theta 2 Is the angle between the substream and the main stream.
As a preferred embodiment, the hydrographic continuous index is obtained by:
wherein F (x) is a probability distribution function, x represents the water level at time/day/month, and the unit is m; pi is hydrologic continuous index, pi is more than or equal to 0 and less than or equal to 1.
As a preferred embodiment, the hydrologic communication comprehensive index C is obtained as follows:
C＝Φ×Π；
wherein Φ is a hydrologic flow index, and pi is a hydrologic continuous index.
The invention also includes the following:
a storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the hydrologic function connectivity assessment method as described above for a multi-gate dam plain urban river network.
A computer device comprising a storage medium, a processor and a computer program stored in the storage medium and executable by the processor, which when executed by the processor performs the steps of the hydrologic function connectivity assessment method as described above for a multi-gate dam plain urban river network.
Drawings
Fig. 1 is a flow chart of a hydrologic function communication evaluation method suitable for a multi-gate dam plain city river network provided in embodiment 1 of the present invention.
Detailed Description
The drawings are for illustrative purposes only and are not to be construed as limiting the present patent;
it should be understood that the described embodiments are merely some, but not all, of the embodiments of the present application. All other embodiments, which can be made by one of ordinary skill in the art without undue burden from the embodiments of the present application, are within the scope of the embodiments of the present application.
The terminology used in the embodiments of the application is for the purpose of describing particular embodiments only and is not intended to be limiting of the embodiments of the application. As used in this application and the appended claims, the singular forms "a," "an," and "the" are intended to include the plural forms as well, unless the context clearly indicates otherwise. It should also be understood that the term "and/or" as used herein refers to and encompasses any or all possible combinations of one or more of the associated listed items.
When the following description refers to the accompanying drawings, the same numbers in different drawings refer to the same or similar elements, unless otherwise indicated. The implementations described in the following exemplary examples are not representative of all implementations consistent with the present application. Rather, they are merely examples of apparatus and methods consistent with some aspects of the present application as detailed in the accompanying claims. In the description of this application, it should be understood that the terms "first," "second," "third," and the like are used merely to distinguish between similar objects and are not necessarily used to describe a particular order or sequence, nor should they be construed to indicate or imply relative importance. The specific meaning of the terms in this application will be understood by those of ordinary skill in the art as the case may be.
Furthermore, in the description of the present application, unless otherwise indicated, "a plurality" means two or more. "and/or", describes an association relationship of an association object, and indicates that there may be three relationships, for example, a and/or B, and may indicate: a exists alone, A and B exist together, and B exists alone. The character "/" generally indicates that the context-dependent object is an "or" relationship. The invention is further illustrated in the following figures and examples.
In order to solve the limitations of the prior art, the present embodiment provides a technical solution, and the technical solution of the present invention is further described below with reference to the drawings and the embodiments.
Example 1
Referring to fig. 1, the hydrologic function communication assessment method suitable for the multi-gate dam plain city river network comprises the following steps:
s1, acquiring data about water level, flow and flow velocity acquired and recorded at a preset main gate dam point;
s2, river channel data of the main gate dam point are obtained, and the area of the main gate dam point on the water cross section, namely the overflow area, is obtained according to the river channel data;
s3, obtaining the actual flow and the ecological base flow of the main gate dam point on the water cross section according to the overflow area and the data about the flow and the flow speed obtained by collecting and recording;
s4, acquiring a hydrological flow index according to the actual flow of the main gate dam point on the water cross section and the ecological base flow;
s5, acquiring a hydrologic continuous index according to the data about the water level obtained by the acquisition record;
and S6, obtaining the hydrologic communication comprehensive index serving as the hydrologic function communication evaluation index according to the hydrologic flow index and the hydrologic continuous index.
Compared with the prior art that the hydrologic function communication state is characterized by mainly utilizing the gate dam communication degree, although the hydrologic function communication level of a water system can be primarily analyzed by utilizing the gate dam communication degree, the technical problem that the dynamic change of the hydrologic communication of a river network cannot be obtained by a functional communication evaluation model based on graph theory and hydrologic communication functions is solved, the actual flow and the ecological base flow of a river basin are obtained by obtaining the overflow area and the collected flow data, and after the actual flow and the ecological base flow are converted into hydrologic flow indexes, the collected hydrologic continuous indexes are combined, and the hydrologic communication comprehensive indexes serving as hydrologic function communication evaluation indexes are output; dynamic changes of hydrologic function communication can be represented in real time, and more scientific and reasonable method support is provided for a gate dam dispatching scheme.
Specifically, the proposal and establishment of hydrologic function communication indexes are key to connectivity assessment, and the selection principle of the indexes is to show effectiveness, availability, objectivity, independence and the like. From the connotation and mechanism analysis of hydrologic communication, the fluidity and continuity of water system are the basic conditions for the existence of river, and are the most important elements of hydrologic communication. Therefore, the flow rate and the passing rate of the artificial structures are mainly considered in establishing the hydrologic function communication index of the water system. Wherein, the flow rate reflects the fluidity of the water flow, and the artificial structure passing rate reflects the continuity of the water flow. However, in real life, the state of water flow is rarely calibrated by the flow rate, and the flow rate is generally used for expressing the overflow capacity of a water system. Meanwhile, the flow data is easier to acquire, and can be classified into different types according to the statistical characteristics. If the flow is divided into annual runoff, lunar runoff, daily runoff, timely runoff and the like according to the time period; the flow rate is classified into a runoff rate, a flood peak flow rate, a minimum flow rate, and the like according to the flow rate characteristics (see table 1). In addition, the connectivity level of the water conservancy junction only characterizes the density and the connection degree of the gate dams between the river channels, and the real-time state of water flow passing through the water conservancy junction cannot be reflected. Therefore, in order to quantify hydrologic functional connectivity, a water flow passing rate evaluation method needs to be established starting from the scheduling rules of the water conservancy junction.
TABLE 1 flow index
The water system overflow capacity refers to the flow rate of a river channel under a certain water level, and reflects the flowing state of the river channel. The calculation formula is as follows:
Q c ＝A L U；
in which Q c M is the overcurrent capacity 3 /s；A L For the area of overcurrent, m 2 The method comprises the steps of carrying out a first treatment on the surface of the U is flow rate, m/s.
In general, flood control flow, flat beach flow, ecological flow and the like can effectively reflect the characteristic flow of a river channel.
Flood control flow: the flood control flow mainly refers to the maximum flood peak flow which can safely pass through the river channel. For the river channel, the maximum overflow flow cannot exceed the flood control flow. The calculation method comprises a model method, a statistical method, a formula method and the like.
Flat beach flow rate: the flat beach flow is the flow when the river channel water level reaches the beach lip, and the flow reflects the overflow capacity and sand conveying capacity of the main channel of the river channel. The calculation method comprises a water level flow method and a Manning formula method.
Ecological flow rate: in order to protect the biodiversity, prevent river channel from breaking and shrinking, maintain the self-cleaning capability of water body, ensure reasonable groundwater level, exert the water-sand migration function of river, the water system needs a certain base flow, namely the ecological flow of river. The ecological flow is mainly influenced by factors such as river hydrology, an ecological system and the like, such as upstream inflow, rainfall, biological species and the like. The main calculation method of the ecological flow comprises a hydrology method, a hydraulics method, an ecological simulation method, an ecological hydraulic radius method and the like.
The water flow section refers to the area surrounded by the water surface line and the river bottom line at a certain research moment, and the overflow area has the same meaning as the water flow section, and can be represented by the water flow section.
The traditional water cross-section area measurement mode is complex, and the water cross-section area is measured and calculated by mainly using a computer program algorithm. Taking the calculation of the cross-sectional area of water from the middle point distance of a river channel as an example:
the middle point distance of the river channel refers to the starting point distance corresponding to the lowest point of the measured elevation in each water section of the river channel.
When the method for analyzing the center distance of the river channel is adopted, the left bank water edge and the right bank water edge of the river channel section can be rapidly and accurately identified by the computer program, so that the water cross section range is locked, the corresponding depth perpendicular under each starting point distance is extracted, and the water cross section area of the river channel is calculated. The program algorithm adopts a method specified by technical standards of hydrologic industry, when in calculation, the depth-measuring vertical line is taken as a boundary, a graph formed by the water surface and the river bottom between adjacent depth-measuring vertical lines, the water surface between the vertical lines is regarded as a trapezoid, the water surface width between the adjacent vertical lines is taken as the height of the trapezoid, and the area of each part is calculated respectively. The partial areas of the two banks are calculated according to the triangular area, and the sum of the partial areas is the cross section area of the river channel.
Thus, as a preferred embodiment, the flow-through area is obtained by:
determining the starting point distances of the water edges of the left and right banks of the water section through a river channel midpoint distance algorithm, and extracting the water surface widths between the adjacent sounding plumb lines and the plumb lines according to the distance range of the starting points of the water edges; dividing the water cross section into a plurality of parts by adopting an average segmentation method and taking the sounding vertical lines as boundaries, wherein the distance between the adjacent vertical lines is a part width, and multiplying the average value of the water depths of the adjacent vertical lines to obtain the area of each part of the water cross section; and calculating the sum of the areas of all the parts to obtain the flow-through area.
Specifically, the fluidity of a river is an important measure of hydrologic communication, and the fluidity cannot be estimated by considering the upstream inflow water flow or the overflow capacity of a river channel alone. Therefore, the relative relationship between the actual flow rate of the section and the river channel flow capacity is the important point of consideration, namely, the comprehensive consideration is taken into consideration in terms of coordination and collocation between the two, so as to reflect the fluidity of the water flow. In summary, the present invention defines the ratio between the actual flow rate and the flow capacity of the section as the hydrologic flow index, and introduces the split/confluence ratio η while considering the influence of the split/confluence.
Thus, as a preferred embodiment, the hydrological flow index Φ is obtained by the following formula:
wherein Q is 0 For the actual flow rate of the main gate dam point on the water cross section, Q e And eta is the dividing/converging ratio for the ecological base flow of the main gate dam point on the water cross section.
Further, the ecological base flow of the main gate dam point on the water cross section is obtained by the following steps:
based on the month average natural flow of several years (data of 30 years or more are selected in the embodiment, the data can be obtained from relevant units of hydrologic bureau and the like, or can be monitored by a flow rate water level measuring instrument), the P-III type curve is utilized, and the minimum month average natural flow at a preset frequency (calculation frequency of P=90% is selected in the embodiment) is selected as the ecological base flow.
Specifically, the P-III type curve, also called the Pearson III (P-III) type curve, is a common method in hydrologic frequency estimation, and the related formula is as follows:
in the above formula, x is the average flow in months, Γ (α) is the gamma function of α; alpha, beta, a 0 Shape dimension and position parameter, alpha, of pearson III type distribution respectively>0，β>0; p is x is greater than or equal to x p Is a cumulative frequency of (a) for a plurality of frequency bands.
α、β、a 0 And (3) with
There is a relation among (flow average), cv (variation coefficient) and Cs (deviation coefficient), wherein Cv and Cs can be obtained by querying a pearson m-type frequency curve coefficient table.
In particular, the tributaries are an important type of water communication, and when the tributaries exist in the main stream, the over-current capacity of the river channel is significantly affected. Thus, the split/split ratio is an important parameter for flow index calculation. In the present embodiment, the shunt state is divided into two cases: diversion and natural diversion.
More specifically, for the diversion state, the diversion amount of the main flow can be controlled by utilizing the dispatching of the diversion gate, and the diversion flow is obtained through the flow monitoring after the gate, so that the diversion ratio is calculated.
Therefore, further, for the two-branched river channel in which the split state occurs, in the case where the split state is the diversion split, the score/confluence ratio η is obtained by:
Q 1 ＝ηQ 0 ；
Q 2 ＝(1-η)Q 0 ；
wherein Q is 0 For dry flow, Q 1 ，Q 2 Respectively, the tributary flows.
More specifically, the split flow of each branch is uncontrollable in the natural split state, and it is difficult to determine the split ratio of each branch. The existing calculation formula of the split ratio mainly comprises a water level difference equivalence method, a momentum balance method, a sand content method and the like. However, since the plain city river is severely channelized, the tributary and the dry stream roughness are basically consistent, the present embodiment derives the split/confluence ratio of the natural split stream by:
Q 1 U 1 sinθ 1 ＝Q 2 U 2 sinθ 2 ；
by means of Q c ＝A L U, converting the above formula into
Let Q 1 ＝ηQ 0 ，Q 2 ＝(1-η)Q 0 Substituting the above formula to obtain
After further conversion from the above formula, for the two branch river channel in the split state, the score/confluence ratio η is obtained by:
wherein A is L1 And A L2 For the area of the tributary, θ 1 And theta 2 Is the angle between the substream and the main stream.
From this, it is clear that the split/confluence ratio in the natural split state has a close relationship with the water passing area of the branch, the branch angle, and the like.
When the branch flows are in a confluent state, the branch/confluent ratio generally obtains the flow rates of the branch flows and the main flow according to actual measurement data, model prediction and the like.
It should be noted that the value of the split/confluence ratio should be determined according to the actual situation, and if the cross section has no split or confluence state, the corresponding split/confluence ratio η=0, the hydroflow index Φ becomes
At this time, when Q 0 Less than Q e When the flow rate of the water flow is lower than the minimum flow rate for maintaining the ecological system, the hydrodynamic condition is weaker, the normal ecological requirement cannot be maintained, and the hydrologic connectivity is extremely poor.
In particular, for rivers, the continuity of water flow is the basis for the communication of hydrologic functions. The plain city river network area is generally located in a densely populated and more developed city area, is more easily interfered by human activities, the water system is controlled by a large number of gate dams to prevent flood and environmental deterioration, and the construction of the large number of gate dams and the complex scheduling rules bring uncertain factors for the continuity of water flow. In order to analyze the influence of the gate dam on hydrologic function communication, scholars have made a great deal of research. For example, the water head difference between adjacent gate dam points is used for representing the magnitude of hydrologic functional communication of a research area; analyzing longitudinal functional connectivity of the river by using tree-like river connectivity index (DCI); and analyzing the hydrologic function communication state with 50% gate dam opening probability in a mode of combining hydrologic communication functions and graph theory. Because the gate dams in plain river network areas are mainly limited by complex scheduling rules, such as water level of a specific section, rainfall and the like. Although the above-mentioned research method can analyze the communication degree of the gate dams to a certain extent, it is not suitable for plain river network areas with complex scheduling rules. In summary, the floodgate is regulated by specific scheduling rules to ensure that the water level does not exceed the local warning water level. Currently, the warning water level is a key that affects the opening or closing of the floodgate, so the occurrence probability of the warning water level represents the cumulative probability of the floodgate opening. Therefore, under the sluice scheduling rule of the research area, the embodiment defines the safety interval accumulation probability of the water level as the sluice passing probability, namely the hydrologic continuous index pi.
As a preferred embodiment, the hydrographic continuity index is obtained by:
wherein F (x) is a probability distribution function, x represents the water level at time/day/month, and the unit is m; pi is hydrologic continuous index, pi is more than or equal to 0 and less than or equal to 1.
Specifically, the present embodiment introduces a probability distribution function F (x) for quantifying the continuous index pi, and as an alternative embodiment, F (x) examines the optimal probability distribution function of the water level sequence by the Kolmogorov-Smirnov (K-S) method. The Kolmogorov-Smirnov test is one of the most commonly used normal test methods, which mainly calculates the distance between an empirical distribution and a theoretical distribution, and takes the largest distance (difference) as a test statistic. The theoretical cumulative frequency distribution of the classified data is first obtained, compared with the observed empirical cumulative frequency distribution, their maximum deviation value is determined, and then it is checked at a given level of significance whether such deviation value is occasional.
The analysis process of this example shows that the main implications of connectivity are fluidity and continuity of water flow. Thus, the present embodiment defines a hydrologic connectivity composite index as the product of a hydrologic flow index and a hydrologic continuity index, and creates a composite index formula, characterizing the hydrologic functional connectivity level with the connectivity composite index, the hydrologic connectivity composite index C being obtained as follows:
C＝Φ×Π；
wherein Φ is a hydrologic flow index, and pi is a hydrologic continuous index.
The traditional hydrologic function communication evaluation method mainly comprises a graph theory method, a hydrologic communication function method, a water level difference method and the like, and the method cannot accurately represent the real-time state and the actual scheduling situation of water flow. The hydrologic communication comprehensive index formula established by the invention is improved aiming at the defects, so that the model is suitable for hydrologic function communication evaluation of complex scheduling rules and more accords with the actual situation of hydrologic communication.
Example 2
A storage medium having stored thereon a computer program which, when executed by a processor, implements the steps of the hydrologic function connectivity assessment method applicable to multi-gate dam plain urban river networks as described in embodiment 1.
Example 3
A computer device comprising a storage medium, a processor and a computer program stored in the storage medium and executable by the processor, which when executed by the processor implements the steps of the hydrologic function connectivity assessment method applicable to multi-gate dam plain urban river network as described in embodiment 1.
It is to be understood that the above examples of the present invention are provided by way of illustration only and not by way of limitation of the embodiments of the present invention. Other variations or modifications of the above teachings will be apparent to those of ordinary skill in the art. It is not necessary here nor is it exhaustive of all embodiments. Any modification, equivalent replacement, improvement, etc. which come within the spirit and principles of the invention are desired to be protected by the following claims.
Claims (10)
1. The hydrologic function communication assessment method suitable for the multi-gate dam plain urban river network is characterized by comprising the following steps of:
s1, acquiring data about water level, flow and flow velocity acquired and recorded at a preset main gate dam point;
s2, river channel data of the main gate dam point are obtained, and the area of the main gate dam point on the water cross section, namely the overflow area, is obtained according to the river channel data;
s3, obtaining the actual flow and the ecological base flow of the main gate dam point on the water cross section according to the overflow area and the data about the flow and the flow speed obtained by collecting and recording;
s4, acquiring a hydrological flow index according to the actual flow of the main gate dam point on the water cross section and the ecological base flow;
s5, acquiring a hydrologic continuous index according to the data about the water level obtained by the acquisition record;
and S6, obtaining the hydrologic communication comprehensive index serving as the hydrologic function communication evaluation index according to the hydrologic flow index and the hydrologic continuous index.
2. The hydrologic function connectivity assessment method applicable to multi-gate dam plain city river network according to claim 1, wherein the overflow area is obtained by:
determining the starting point distances of the water edges of the left and right banks of the water section through a river channel midpoint distance algorithm, and extracting the water surface widths between the adjacent sounding plumb lines and the plumb lines according to the distance range of the starting points of the water edges; dividing the water cross section into a plurality of parts by adopting an average segmentation method and taking the sounding vertical lines as boundaries, wherein the distance between the adjacent vertical lines is a part width, and multiplying the average value of the water depths of the adjacent vertical lines to obtain the area of each part of the water cross section; and calculating the sum of the areas of all the parts to obtain the flow-through area.
3. The method for evaluating hydrologic functional connectivity of a multi-gate dam plain urban river network according to claim 1, wherein the hydrologic flow index Φ is obtained by the following formula:
wherein Q is 0 For the actual flow rate of the main gate dam point on the water cross section, Q e And eta is the dividing/converging ratio for the ecological base flow of the main gate dam point on the water cross section.
4. The hydrologic function communication assessment method applicable to multi-gate dam plain city river network according to claim 3, wherein the ecological base flow of the main gate dam point at the water cross section is obtained by the following method:
based on the month average natural flow of a plurality of years, a P-III type curve is utilized, and the minimum month average natural flow under a preset frequency is selected as the ecological base flow.
5. The hydrologic function communication evaluation method applicable to multi-gate dam plain urban river network according to claim 3, wherein for the two branch river channels in the split state, the score/confluence ratio η is obtained by:
Q 1 ＝ηQ 0 ；
Q 2 ＝(1-η)Q 0 ；
wherein Q is 0 For dry flow, Q 1 ，Q 2 Respectively, the tributary flows.
6. The hydrologic function communication evaluation method applicable to multi-gate dam plain city river network according to claim 3, wherein for the two branch river channels in the split state, the score/confluence ratio η is obtained by:
wherein A is L1 And A L2 For the area of the tributary, θ 1 And theta 2 Is the angle between the substream and the main stream.
7. The method for evaluating hydrologic function connectivity of a multi-gate dam plain urban river network according to claim 1, wherein the hydrologic continuity index is obtained by:
wherein F (x) is a probability distribution function, x represents the water level at time/day/month, and the unit is m; pi is hydrologic continuous index, pi is more than or equal to 0 and less than or equal to 1.
8. The hydrologic function connectivity assessment method applicable to multi-gate dam plain city river network according to claim 1, wherein the hydrologic connectivity comprehensive index C is obtained as follows:
C＝Φ×Π；
wherein Φ is a hydrologic flow index, and pi is a hydrologic continuous index.
9. A storage medium having a computer program stored thereon, characterized by: the computer program, when executed by a processor, implements the steps of the hydrologic function connectivity assessment method according to any one of claims 1 to 8, applicable to multi-gate dam plain urban river networks.
10. A computer device, characterized by: comprising a storage medium, a processor, a computer program stored in the storage medium and executable by the processor, the computer program when executed by the processor implementing the steps of the hydrologic function connectivity assessment method according to any one of claims 1 to 8, applicable to multi-gate dam plain urban river network.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310311253.XA CN116167242A (en) | 2023-03-28 | 2023-03-28 | Hydrologic function communication assessment method suitable for multi-gate dam plain urban river network |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310311253.XA CN116167242A (en) | 2023-03-28 | 2023-03-28 | Hydrologic function communication assessment method suitable for multi-gate dam plain urban river network |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN116167242A true CN116167242A (en) | 2023-05-26 |
Family
ID=86422095
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310311253.XA Pending CN116167242A (en) | 2023-03-28 | 2023-03-28 | Hydrologic function communication assessment method suitable for multi-gate dam plain urban river network |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116167242A (en) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117057272A (en) * | 2023-08-15 | 2023-11-14 | 长江水利委员会长江科学院 | Method, system and medium for calculating hydrologic connectivity of gated lake and river |
-
2023
- 2023-03-28 CN CN202310311253.XA patent/CN116167242A/en active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117057272A (en) * | 2023-08-15 | 2023-11-14 | 长江水利委员会长江科学院 | Method, system and medium for calculating hydrologic connectivity of gated lake and river |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN110009002B (en) | Mountain torrent rapid early warning and forecasting method based on multi-dimensional rainfall characteristic space optimal decision | |
| CN116304161B (en) | Hydrologic connectivity calculation method, system and medium for natural connected river-lake | |
| CN116167242A (en) | Hydrologic function communication assessment method suitable for multi-gate dam plain urban river network | |
| Chen et al. | Real-time discharge measurement in tidal streams by an index velocity | |
| Ensign et al. | A method for using shoreline morphology to predict suspended sediment concentration in tidal creeks | |
| Kamel | Application of a hydrodynamic MIKE 11 model for the Euphrates River in Iraq | |
| CN110847112B (en) | River flood discharge early warning method based on hydraulics simulation | |
| Bina et al. | Experimental study of discharge coefficient and trapping ratio in mesh-panel bottom rack for sediment and non-sediment flow and supercritical approaching conditions | |
| Azarang et al. | Simulating the erosion and sedimentation of Karun Alluvial River in the Region of Ahvaz (southwest of Iran) | |
| Westergard et al. | Regionalized equations for bankfull-discharge and channel characteristics of streams in New York State—Hydrologic Region 5 in central New York | |
| Snell et al. | Application of the meta‐channel concept: Construction of the meta‐channel hydraulic geometry for a natural catchment | |
| Zorndt et al. | The influence of hydrodynamic boundary conditions on characteristics, migration, and associated sand transport of sand dunes in a tidal environment: A long-term study of the Elbe Estuary | |
| Negm et al. | Impacts of GERD on the accumulated sediment in Lake Nubia using machine learning and GIS techniques | |
| CN115511291A (en) | Karst depression water passing capacity analysis method based on inland inundation water depth monitoring | |
| Demissie et al. | Hydrology, hydraulics, and sediment transport, Kankakee and Iroquois Rivers | |
| CN113887053A (en) | Municipal drainage data quality assessment method and system for pipe network water flow calculation | |
| Navratil et al. | Test of three methods to detect the overbank flow from water level time‐series analysis | |
| CN118133726B (en) | River channel dead water level decline inhibition node identification and classification method | |
| Nama et al. | Review and state of the art for the hydro-morphological modeling of transboundary rivers, Tigris River as a case study | |
| Milišić et al. | Estimation of channel and flood plain roughness using HEC-RAS model: A case study of the Veseočica River, Bosnia and Herzegovina | |
| CN117057272B (en) | Method, system and medium for calculating hydrologic connectivity of gated lake and river | |
| Raileanu et al. | An Evaluation of the Dynamics of Some Meteorological and Hydrological Processes along the Lower Danube. Sustainability 2023, 15, 6087 | |
| Smith | Investigating Ice-Impacted Fluvial Processes in a Small Riffle-Pool Sequence in Southern Ontario | |
| Cilli | Evaluation of bedload yield in two small sand-bed rivers of the Ravenna province (Italy) | |
| Shrestha et al. | Stone Quarries and Hydrological Regime of Roshi with Emphasis on Sedimentation |
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 |