CN106897529B - Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation - Google Patents

Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation Download PDF

Info

Publication number
CN106897529B
CN106897529B CN201710120454.6A CN201710120454A CN106897529B CN 106897529 B CN106897529 B CN 106897529B CN 201710120454 A CN201710120454 A CN 201710120454A CN 106897529 B CN106897529 B CN 106897529B
Authority
CN
China
Prior art keywords
sub
basin
upstream
confluence
pipe network
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN201710120454.6A
Other languages
Chinese (zh)
Other versions
CN106897529A (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.)
China Institute of Water Resources and Hydropower Research
Original Assignee
China Institute of Water Resources and Hydropower Research
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 China Institute of Water Resources and Hydropower Research filed Critical China Institute of Water Resources and Hydropower Research
Priority to CN201710120454.6A priority Critical patent/CN106897529B/en
Publication of CN106897529A publication Critical patent/CN106897529A/en
Application granted granted Critical
Publication of CN106897529B publication Critical patent/CN106897529B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/20Design optimisation, verification or simulation
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06FELECTRIC DIGITAL DATA PROCESSING
    • G06F30/00Computer-aided design [CAD]
    • G06F30/10Geometric CAD
    • G06F30/18Network design, e.g. design based on topological or interconnect aspects of utility systems, piping, heating ventilation air conditioning [HVAC] or cabling
    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A20/00Water conservation; Efficient water supply; Efficient water use
    • Y02A20/152Water filtration

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Theoretical Computer Science (AREA)
  • General Physics & Mathematics (AREA)
  • Geometry (AREA)
  • Computer Hardware Design (AREA)
  • General Engineering & Computer Science (AREA)
  • Evolutionary Computation (AREA)
  • Computational Mathematics (AREA)
  • Pure & Applied Mathematics (AREA)
  • Mathematical Optimization (AREA)
  • Mathematical Analysis (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • Management, Administration, Business Operations System, And Electronic Commerce (AREA)

Abstract

The invention discloses a sewage pipe network topological relation-based sub-basin confluence arithmetic sequence calculation method, which comprises the following steps: dividing sub-basins according to a natural water system, and distributing unique codes for the sub-basins; constructing an upstream and downstream confluence topological relation table of a sub-basin where the natural water system is located according to the natural water system confluence path and the sub-basin codes; constructing an upstream and downstream confluence topological relation table of the sub-basin where the sewage pipe network is located according to the confluence path of the sewage pipe network and the serial number of the sub-basin; according to the two topological relation tables, circularly detecting the natural water system and the upstream sub-watersheds of the sewage pipe network one by one, and sequentially distributing continuous confluence calculation serial numbers to the upstream sub-watersheds according to the sequence of successfully finishing confluence calculation by the upstream sub-watersheds until all the sub-watersheds have confluence calculation serial numbers; and determining the confluence operation sequence of each sub-basin according to the confluence operation serial number. The invention solves the problem of repeated calculation or jump calculation of the calculation unit in the model convergence calculation process, and improves the operation efficiency of the convergence calculation.

Description

Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation
Technical Field
The invention relates to the field of hydrological models, in particular to a sub-basin confluence arithmetic sequence calculation method based on a sewage pipe network topological relation.
Background
Along with the development of the industrial and urbanization process, the urban pipe network coverage rate is larger and larger, the pipe network converging process is more and more complex, the distributed hydrological model and the pipe network model are commonly used for simulating urban waterlogging and urban water environment, how to reasonably carry out converging operation directly influences the simulation result and the model efficiency of water quantity and water quality. The conventional hydrological model is mainly based on digital elevation information to divide sub-watersheds and codes the sub-watersheds step by step according to the sequence from upstream to downstream, so that when river convergence calculation is performed, calculation is performed according to the number of the sub-watersheds from small to large, and the upstream sub-watersheds can be guaranteed to finish convergence calculation when the downstream sub-watersheds perform convergence calculation.
However, because the urban pipe network is complex in distribution, the sub-basin codes of the sewage path of the pipe network are not arranged from small to large, so that the pipe network convergence and the slope river convergence cannot be calculated simultaneously, or repeated calculation exists, and the risk of result error is caused. Therefore, whether the coding of pipe network convergence on the related sub-basin reasonably and directly influences the calculation efficiency and result of the water quantity and the water quality of the basin.
Disclosure of Invention
The invention aims to solve the technical problem that the existing hydrological model carries out confluence calculation according to the upstream and downstream relations of the natural water system, so that the confluence process of sewage pipes of each sub-basin has repeated calculation or jump calculation, thereby causing the wrong calculation result.
In order to solve the technical problems, the technical scheme adopted by the invention is to provide a sub-basin confluence operation sequence calculation method based on a sewage pipe network topological relation, which comprises the following steps:
step S10, dividing sub-domains according to the natural water system, and distributing a unique code for each sub-domain;
s20, constructing an upstream and downstream confluence topological relation table of the sub-basin where the natural water system is located according to the natural water system confluence path and the sub-basin codes;
s30, constructing an upstream and downstream confluence topological relation table of the sub basin where the sewage pipe network is located according to the sewage pipe network confluence path and the sub basin number;
step S40, circularly detecting the natural water system and the sewage pipe network upstream sub-basin one by one according to the upstream and downstream confluence topological relation table of the sub-basins where the natural water system and the sewage pipe network are located, and sequentially distributing continuous confluence calculation serial numbers for the sub-basins according to the sequence that the upstream sub-basins can smoothly complete confluence calculation until all the sub-basins have confluence calculation serial numbers;
and step S50, determining the convergence calculation sequence of each sub-basin in the distributed water quality model according to the convergence calculation sequence number.
In the above-mentioned method, the first step of the method,
the upstream and downstream confluence topological relation table of the sub-basin where the natural water system is located is constructed according to the upstream and downstream topological relation information of the natural water system relative to each sub-basin, which is obtained according to the topological relation of the confluence path of the natural water system;
the upstream and downstream confluence topological relation table of the sub-basin where the sewage pipe network is located is constructed according to the upstream and downstream topological relation information of the sewage pipe network relative to each sub-basin, which is obtained according to the topological relation of the sewage pipe network in the basin;
the topological relation table of upstream and downstream confluence of the sub-basin where the natural water system or the sewage pipe network is located comprises:
and (3) sub-stream domain coding: unique coding of all sub-domains in the system;
upstream "m" encoding: the coding value of the upstream sub-basin is 0, the coding value of the upstream sub-basin is no upstream sub-basin, and the coding value of the upstream sub-basin is greater than 0, so that the coding value of the upstream sub-basin corresponding to the sub-basin is represented;
number of upstream sub-watersheds: the number of the upstream sub-watersheds is equal to that of the upstream sub-watersheds with the coding value of m not 0 in the upstream and downstream confluence topological relation table of the sub-watersheds where the natural water system or the sewage pipe network is located.
In the above method, in step S40, the criterion that the upstream sub-basin can successfully complete the convergence calculation is:
the sub-basin has no upstream sub-basin;
alternatively, all the upstream sub-basins of the corresponding natural water system and sewage pipe network already have the confluence operation numbers.
In the above method, the step S40 includes the steps of:
step S41, setting the convergence calculation sequence number array runId (1: n) of all sub-domains in the entire convergence calculation system to 0, where the current convergence calculation sequence number nowId is 0 and n is the number of sub-domains in the entire convergence calculation system;
step S42, setting the traversal sequence number S of the current sub-basin to be 0, and starting to traverse the sub-basins one by one according to the upstream and downstream confluence topological relation table of the sub-basins where the natural water system and the sewage pipe network are located;
step S43, judging whether S is n, if yes, executing step S42; otherwise, go to step S44;
step S44, let S be S +1, and start detecting the next sub-basin;
step S45, determine whether the current convergence operation number runId (S) of the sub-basin is 0, if yes, execute step S46; otherwise, go to step S43;
step S46, determining whether the current sub-basin has a sub-basin upstream of the natural water system or an upstream sub-basin of the sewage pipe network, and an upstream sub-basin with a confluence algorithm serial number of 0, if yes, executing step S43; otherwise, go to step S47;
step S47, assigning the nowId to nowId +1, and assigning the nowId value to the merge operation serial number runid (S) of the current sub-domain;
step S48, judging whether the nowId is equal to the number n of the sub-watersheds, if so, executing step S49; otherwise, go to step S43;
step S49, the routine ends.
In the above method, the step S50 includes the steps of:
constructing a convergence operation sequence array runSub (1: n) and initializing to 0, wherein runSub (i) represents a sub-basin code of the ith operation, and a current operation serial number r is set to be 1;
step S52, setting the traversal sequence number S of the current sub-basin to 0;
step S53, let S be S +1, and determine whether the traversal sequence number S is greater than the number n of sub-domains, if yes, execute step S52; otherwise, go to step S54;
step S54, determine whether the current convergence operation sequence number runId (S) of the sub-basin is equal to the previous operation sequence number r, if yes, execute step S55; otherwise, go to step S53;
step S55, assigning the current sub-stream field coding to runsub (r), and making r ═ r + 1;
step S56, judging whether the current operation sequence number r is larger than the number n of the sub-watersheds, if so, executing step S57, otherwise, executing step S54;
step S57, the routine ends.
According to the invention, each sub-basin is uniquely coded, and each sub-basin is traversed according to the upstream and downstream confluence topological relation table of the sub-basin where the natural water system and the sewage pipe network are located, on the premise of ensuring that the natural water system confluence process and the sewage pipe network confluence process of the upstream sub-basin of each sub-basin are completed, confluence calculation is carried out on the sub-basins, so that the natural runoff and the sewage inflow amount obtained by calculating the upstream sub-basin can be directly used, the problem that the confluence calculation repeated calculation or jump calculation is carried out on the existing hydrological model according to the upstream and downstream relation of the natural water system is solved, and the method has the following specific beneficial effects:
(1) the invention can simultaneously meet the requirements of natural water system river confluence and sewage pipe network confluence algorithm of the sub-basin;
(2) the invention fully considers the topological relation of the sewage pipe network in the sub flow domain, and improves the operation efficiency;
(3) according to the topological relation of different pipe sections, the repeated calculation of sewage confluence is avoided.
Drawings
FIG. 1 is a flow chart of a method for calculating a sequence of basin convergence calculation based on a topological relation of a sewage pipe network according to the present invention;
FIG. 2 is a schematic view of an embodiment of the sub-basin natural water system and sewer network distribution of the present invention;
FIG. 3 is a flowchart illustrating step S40 according to the present invention;
fig. 4 is a flowchart illustrating step S50 in the present invention.
Detailed Description
The invention is described in detail below with reference to the figures and specific examples.
As shown in fig. 1, the method for calculating the sub-basin confluence algorithm sequence based on the topological relation of the sewage pipe network provided by the invention comprises the following steps:
step S10, the sub-domains are divided according to the natural water system and encoded so that each sub-domain has a unique code.
In step S10, the coding scheme of the sub-stream field is continuously increased from 1, and there is no repetition coding between the sub-stream fields. As shown in fig. 2, the whole example basin is divided into 14 sub-basins according to the natural water system, and each sub-basin is assigned with a unique number (numbers 1 to 14 in fig. 2), and it should be noted here that when the sub-basins are coded, each sub-basin code is not repeated, no matter how the sub-basin codes are carried out, different sub-basin coding systems are obtained, and the subsequent calculation of the method of the present invention is not affected.
And S20, constructing an upstream and downstream confluence topological relation table of the sub-basin where the natural water system is located according to the natural water system confluence path and the sub-basin codes, and obtaining the number information of the direct upstream sub-basin of each sub-basin relative to the natural water system.
In step S20, first, according to the topological relation of the natural water converging path (dashed line in fig. 2), the upstream topological relation information of each sub-basin with respect to the natural water system is obtained, as shown in table 1 (corresponding to fig. 2), where an upstream "m" code indicates an upstream sub-basin code merging into the current sub-basin in the natural water system, m indicates the several upstream sub-basins of the current sub-basin, a sub-basin code of "0" indicates that there is no upstream sub-basin, and a value greater than "0" indicates an encoding value of the upstream sub-basin corresponding to the sub-basin; the number of the upstream sub-stream domains is equal to the number of all upstream m coding values which are not 0 and are used for statistical calculation without influencing the calculation sequence.
And step S30, constructing an upstream and downstream confluence topological relation table of the sub-basin where the sewage pipe network is located according to the sewage pipe network confluence path and the sub-basin number, and obtaining the number information of the direct upstream sub-basin of each sub-basin relative to the sewage pipe network.
In step S30, first, according to the topological relationship of the sewage pipe network (the thick solid line in fig. 2) in the drainage basin, the upstream topological relationship information of the sewage pipe network relative to each sub-drainage basin is obtained, as shown in table 1, it should be noted that in the embodiment shown in table 1, the sewage pipe network system has only 2 upstream at most, and the natural water system has 3 upstream at most, which is related to a specific drainage basin, and the maximum number of upstream does not affect the subsequent calculation.
Table 1: and a topological relation table of upstream and downstream confluence of the sub-basin where the natural water system and the sewage pipe network are located.
Figure BDA0001236853000000061
And S40, circularly detecting the natural water system and the sewage pipe network upstream sub-basin one by one according to the upstream and downstream confluence topological relation table of the sub-basins where the natural water system and the sewage pipe network are located, and sequentially allocating continuous confluence calculation serial numbers to the sub-basins according to the sequence in which the upstream sub-basins can smoothly complete confluence calculation until all the sub-basins have confluence calculation serial numbers.
In the present invention, the fact that the upstream sub-basin can smoothly complete the convergence calculation means that the sub-basin has no upstream sub-basin; or all the upstream sub-watersheds of the corresponding natural water system and the sewage pipe network have the confluence operation serial number, namely when one sub-watersheds carries out confluence operation, the natural water system confluence process and the sewage pipe network confluence process of the upstream sub-watersheds are required to be completed, so that when the confluence operation is carried out on the current sub-watersheds, the natural runoff and the sewage inflow amount obtained by calculation of the upstream sub-watersheds can be directly used; and when the sub-basin has no upstream sub-basin, the upstream natural runoff and the sewage inflow amount are both 0, and the current sub-basin confluence algorithm is not influenced. If the sub-basin having no confluence operation number exists in the upstream sub-basin (the upstream sub-basin of the natural water system or the sewage pipe network), and the sub-basin cannot complete the confluence operation, the next sub-basin is directly skipped to detect, and when the next cycle is carried out, the upstream sub-basin is detected again.
In the present invention, as shown in fig. 3, step S40 includes the steps of:
step S41, setting the convergence calculation sequence number array runId (1: n) of all sub-domains in the whole convergence calculation system to 0, which indicates that the convergence calculation sequence number calculation is not completed in all sub-domains; the current convergence operation serial number nowId is set to 0, which means that no convergence operation serial number is assigned to any sub-domain, and n is the number of sub-domains in the entire convergence operation system.
And step S42, setting the traversal sequence number S of the current sub-basin to be 0, and starting to traverse the sub-basins one by one according to the upstream and downstream confluence topological relation table of the sub-basins where the natural water system and the sewage pipe network are located.
Step S43, judging whether S is n, if yes, executing step S42; otherwise, step S44 is executed.
In step S44, let S be S +1, and the next sub-basin starts to be detected.
Step S45, determine whether the current convergence operation number runId (S) of the sub-basin is 0, if yes, execute step S46; otherwise, step S43 is executed.
Step S46, determining whether the current sub-basin has a sub-basin upstream of the natural water system or an upstream sub-basin of the sewage pipe network, and an upstream sub-basin with a confluence algorithm serial number of 0, if yes, executing step S43; otherwise, step S47 is executed.
Step S47 represents adding one sub-domain in which the convergence operation can be smoothly completed, and assigning the value of nowId to the convergence operation number runid (S) of the current sub-domain.
Step S48, judging whether the nowId is equal to the number n of the sub-watersheds, if so, executing step S49; otherwise, step S43 is executed.
Step S49, the routine ends.
The calculation process of the confluence operation number will be described in detail below with reference to fig. 2 and table 1 as an example.
A first loop traversal is performed on table 1:
for the sub-basin 1, since the sub-basin 2 is the upstream sewage pipe network sub-basin and the current convergence operation serial number of the sub-basin 2 is 0, the convergence operation serial number of the sub-basin 1 cannot be set; the sub-basin traversal sequence number s is increased progressively to judge the sub-basin 2, and as for the sub-basin 2, no upstream sub-basin exists, the convergence operation sequence number of the sub-basin 2 can be set to be 1; the traversal sequence number of the incremental sub-basin is judged, and as for the sub-basins 3 to 6, because all the sub-basins exist in the upstream sub-basin and the corresponding upstream sub-basin has the condition that the convergence calculation sequence number is equal to 0, the convergence calculation sequence number cannot be assigned; for the sub-basin 7, since there is no upstream sub-basin, the number of the bus operation can be set, and since the sub-basin 2 is set with the number of the bus operation 1, the number of the bus operation of the sub-basin 7 is set to 2 in sequence; similarly, the convergence operation serial number of the sub-basin 12 is set to be 3; the convergence operation number of the sub-basin 13 is 4; the sub-basin 14 has an upstream sub-basin of the sewage pipe network, which is the sub-basin 3, and at this time, the sub-basin 3 has no converging operation serial number, so that the converging operation serial number cannot be assigned, and the first cycle traversal is finished.
A second loop traversal is performed on table 1:
for the sub-basin 1, since the confluence operation number of the sub-basin 2 is 1, the sub-basin 1 meets the standard that the operation can be completed smoothly, and the confluence operation number of the sub-basin 1 is set to be 5; similarly, the convergence operation serial number of the sub-basin 3 is set to 6; for the sub-basin 4, the upstream sub-basin 14 does not complete the setting of the serial number of the convergence operation, so that the assignment cannot be made; similarly, the convergence operation number of the sub-basin 5 is set to 7, the sub-basins 6, 8, 9, 10 cannot be set, and the convergence operation numbers of the sub-basins 11 and 14 are set to 8 and 9, respectively; and ending the second cycle traversal.
Note that in the second round of traversal, the sub-domains 2, 7, 12, and 13 do not need to be determined again because the bus operation numbers are already set, and therefore they are skipped directly in the second round of traversal.
And continuing to perform the third, fourth and fifth loop traversal of table 1 until all 14 sub-domains have the merge operation serial number, at which point the process ends.
And step S50, determining the convergence calculation sequence of each sub-basin in the distributed water quality model according to the convergence calculation sequence number.
In step S40, the calculated convergence operation sequence number is the model operation sequence, but in the actual model calculation, if the convergence operation sequence number is searched from 1 to 14 according to the sub-domain sequence numbers for operation, it is not efficient to go through continuously and search and calculate the sub-domains one by one according to the convergence operation sequence number, so the convergence operation sequence number sequence needs to be processed, and the actual calculation only needs to be performed for 14 (number of sub-domains) times;
in the present invention, step S50 includes the steps of:
step S51, constructing a convergence arithmetic sequence array runSub (1: n), and initializing to 0, wherein runSub (i) represents the sub-basin code of the ith operation; setting a current operation serial number r to be 1;
step S52, setting the traversal sequence number S of the current sub-basin to 0;
step S53, let S be S +1, and determine whether the traversal sequence number S is greater than the number n of sub-domains, if yes, execute step S52; otherwise, go to step S54;
step S54, determine whether the current convergence operation sequence number runId (S) of the sub-basin is equal to the previous operation sequence number r, if yes, execute step S55; otherwise, go to step S53;
step S55, assigning the current sub-stream field coding to runsub (r), and making r ═ r + 1;
step S56, judging whether the current operation sequence number r is larger than the number n of the sub-watersheds, if so, executing step S57, otherwise, executing step S54;
step S57, the routine ends.
The calculation process of the confluence operation number will be described in detail below with reference to fig. 2 and table 1 as an example.
Increasing in sequence according to r ═ 1 to n, and performing first cycle traversal on table 1: when r is 1, runSub (1) is 2; when r is 2, runSub (2) is 5; sequentially obtaining runSub (3) ═ 12 and runSub (4) ═ 13;
a second loop traversal is performed on table 1: obtaining runSub (5) ═ 1, runSub (6) ═ 3, runSub (7) ═ 5, runSub (8) ═ 11, runSub (9) ═ 14;
a third loop traversal is performed on table 1: obtaining runSub (10) ═ 4 and runSub (11) ═ 10;
a fourth loop traversal is performed on table 1: to obtain runSub (12) ═ 9;
a fifth round-robin traversal is performed on table 1: obtaining runSub (13) ═ 6 and runSub (14) ═ 8;
in the actual model calculation process, the sub-basin of the i-th calculation, runsub (i), is obtained in sequence from r 1 to 14.
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 (3)

1. A sub-basin confluence operation sequence calculation method based on a sewage pipe network topological relation is characterized by comprising the following steps:
step S10, dividing sub-domains according to the natural water system, and distributing a unique code for each sub-domain;
s20, constructing an upstream and downstream confluence topological relation table of the sub-basin where the natural water system is located according to the natural water system confluence path and the sub-basin codes;
s30, constructing an upstream and downstream confluence topological relation table of the sub basin where the sewage pipe network is located according to the sewage pipe network confluence path and the sub basin number;
step S40, circularly detecting the natural water system and the sewage pipe network upstream sub-basin one by one according to the upstream and downstream confluence topological relation table of the sub-basins where the natural water system and the sewage pipe network are located, and sequentially distributing continuous confluence calculation serial numbers for the sub-basins according to the sequence that the upstream sub-basins can smoothly complete confluence calculation until all the sub-basins have confluence calculation serial numbers;
step S50, determining the convergence calculation sequence of each sub-basin in the distributed water quality model according to the convergence calculation sequence number;
the upstream and downstream confluence topological relation table of the sub-basin where the natural water system is located is constructed according to the upstream and downstream topological relation information of the natural water system relative to each sub-basin, which is obtained according to the topological relation of the confluence path of the natural water system;
the upstream and downstream confluence topological relation table of the sub-basin where the sewage pipe network is located is constructed according to the upstream and downstream topological relation information of the sewage pipe network relative to each sub-basin, which is obtained according to the topological relation of the sewage pipe network in the basin;
the topological relation table of upstream and downstream confluence of the sub-basin where the natural water system or the sewage pipe network is located comprises:
and (3) sub-stream domain coding: unique coding of all sub-domains in the system;
upstream "m" encoding: the coding value of the upstream sub-basin is 0, the coding value of the upstream sub-basin is no upstream sub-basin, and the coding value of the upstream sub-basin is greater than 0, so that the coding value of the upstream sub-basin corresponding to the sub-basin is represented;
number of upstream sub-watersheds: the number of the upstream sub-basin is equal to that of the upstream sub-basins with the coding value of m not 0 in the upstream and downstream confluence topological relation table of the sub-basins where the natural water system or the sewage pipe network is located;
in step S40, the criterion that the upstream sub-basin can successfully complete the convergence calculation is:
the sub-basin has no upstream sub-basin;
alternatively, all the upstream sub-basins of the corresponding natural water system and sewage pipe network already have the confluence operation numbers.
2. The method of claim 1, wherein the step S40 includes the steps of:
step S41, setting the convergence calculation sequence number array runId (1: n) of all sub-domains in the entire convergence calculation system to 0, where the current convergence calculation sequence number nowId is 0 and n is the number of sub-domains in the entire convergence calculation system;
step S42, setting the traversal sequence number S of the current sub-basin to be 0, and starting to traverse the sub-basins one by one according to the upstream and downstream confluence topological relation table of the sub-basins where the natural water system and the sewage pipe network are located;
step S43, judging whether S is n, if yes, executing step S42; otherwise, go to step S44;
step S44, let S be S +1, and start detecting the next sub-basin;
step S45, determine whether the current convergence operation number runId (S) of the sub-basin is 0, if yes, execute step S46; otherwise, go to step S43;
step S46, determining whether the current sub-basin has a sub-basin upstream of the natural water system or an upstream sub-basin of the sewage pipe network, and an upstream sub-basin with a confluence algorithm serial number of 0, if yes, executing step S43; otherwise, go to step S47;
step S47, assigning the nowId to nowId +1, and assigning the nowId value to the merge operation serial number runid (S) of the current sub-domain;
step S48, judging whether the nowId is equal to the number n of the sub-watersheds, if so, executing step S49; otherwise, go to step S43;
step S49, the routine ends.
3. The method of claim 1, wherein the step S50 includes the steps of:
step S51, constructing a convergence operation sequence array runSub (1: n), and initializing to 0, where runSub (i) represents the sub-domain code of the ith operation, and sets the current operation serial number r to 1;
step S52, setting the traversal sequence number S of the current sub-basin to 0;
step S53, let S be S +1, and determine whether the traversal sequence number S is greater than the number n of sub-domains, if yes, execute step S52; otherwise, go to step S54;
step S54, determine whether the current convergence operation sequence number runId (S) of the sub-basin is equal to the previous operation sequence number r, if yes, execute step S55; otherwise, go to step S53;
step S55, assigning the current sub-stream field coding to runsub (r), and making r ═ r + 1;
step S56, judging whether the current operation sequence number r is larger than the number n of the sub-watersheds, if so, executing step S57, otherwise, executing step S54;
step S57, the routine ends.
CN201710120454.6A 2017-03-02 2017-03-02 Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation Active CN106897529B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN201710120454.6A CN106897529B (en) 2017-03-02 2017-03-02 Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN201710120454.6A CN106897529B (en) 2017-03-02 2017-03-02 Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation

Publications (2)

Publication Number Publication Date
CN106897529A CN106897529A (en) 2017-06-27
CN106897529B true CN106897529B (en) 2020-03-31

Family

ID=59185419

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201710120454.6A Active CN106897529B (en) 2017-03-02 2017-03-02 Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation

Country Status (1)

Country Link
CN (1) CN106897529B (en)

Families Citing this family (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN110766792B (en) * 2018-07-25 2021-03-12 中国水利水电科学研究院 Grid calculation order coding method based on ArcGIS hydrological analysis tool
CN110717231B (en) * 2019-09-11 2020-09-29 中国水利水电科学研究院 Sub-basin confluence simulation method based on slope channel river channel three-level structure
CN112612784A (en) * 2020-11-09 2021-04-06 中国再保险(集团)股份有限公司 Automatic river basin computing unit dividing method and device and computer equipment
CN113239502B (en) * 2021-04-23 2022-08-09 湖南中铭科技有限公司 Artificial intelligence image processing-based urban sewage pipe network simulation construction method
CN113128009B (en) * 2021-04-27 2021-12-07 中国水利水电科学研究院 Sub-basin unit division method considering mountain area plain landform difference
CN114882031B (en) * 2022-07-11 2022-09-06 江苏瑞立环保工程股份有限公司 Sewage treatment method and system based on activated sludge process
CN115577506B (en) * 2022-09-22 2023-08-25 长江生态环保集团有限公司 Method for pre-diagnosing and accurately checking problems of sewage pipe network

Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012127155A1 (en) * 2011-03-21 2012-09-27 Universite De Toulouse Ii Le Mirail Visuo-tactile atlas
CN103677826A (en) * 2013-12-09 2014-03-26 河海大学 Sub-basin dividing and information extracting method based on IDL and Mapinfo
CN104142812A (en) * 2014-07-30 2014-11-12 中国水利水电科学研究院 Distributed hydrological model parallel operation method

Patent Citations (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
WO2012127155A1 (en) * 2011-03-21 2012-09-27 Universite De Toulouse Ii Le Mirail Visuo-tactile atlas
CN103677826A (en) * 2013-12-09 2014-03-26 河海大学 Sub-basin dividing and information extracting method based on IDL and Mapinfo
CN104142812A (en) * 2014-07-30 2014-11-12 中国水利水电科学研究院 Distributed hydrological model parallel operation method

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
分布式水文模型子流域编码方法对比分析;刘佳嘉,周祖昊,贾仰文,王浩;《河海大学学报》;20170131;第45卷(第1期);22-29页 *
基于DEM河网干支拓扑关系的子流域编码规则;刘佳嘉,周祖昊,贾仰文,王浩;《河海大学学报》;20130731;第41卷(第4期);288-293页 *

Also Published As

Publication number Publication date
CN106897529A (en) 2017-06-27

Similar Documents

Publication Publication Date Title
CN106897529B (en) Sub-basin confluence operation sequence calculation method based on sewage pipe network topological relation
CN105260607A (en) Serial connection and parallel connection coupling multi-model hydrological forecasting method
Vojinovic et al. Multi-objective rehabilitation of urban drainage systems under uncertainties
CN108614915B (en) Hydrological model free construction strategy method based on scene driving
CN103729550A (en) Multi-model integrated flood forecasting method based on propagation time clustering analysis
CN111985106A (en) Distributed hydrological model parameter calibration method based on multipoint parallel correction
CN108734340A (en) A kind of river flood forecasting procedure generally changed based on big vast type
CN108647778A (en) Dynamic prediction method for drainage flow of drainage port of urban rainwater system
CN105303612A (en) Digital river network extraction method based on TIN model
CN115907536B (en) Reservoir characteristic curve correction method based on interval spurious rate
CN110782526B (en) Different sub-basin river channel grid calculation order coding method
WO2018078674A1 (en) Simulation device, simulation method, and recording medium
CN111275253B (en) Runoff probabilistic prediction method and system integrating deep learning and error correction
CN114841417B (en) High-precision salt tide forecasting method and system and readable storage medium
CN111327706B (en) Remote meter reading data processing method and device, gas meter system and gas cloud platform
Kwon et al. Optimal layout and pipe sizing of urban drainage networks to improve robustness and rapidity
Sood et al. Economic assessment of centralized and decentralized sewerage network systems: a case study of Ludhiana, India
CN108765190A (en) A kind of network of waterways data representation method towards the multistage nesting of large-scale parallel and the suitable network of waterways
CN115511354A (en) Pre-drainage scheduling method, system, equipment and medium for cross-basin water network communication project
CN114036709A (en) Method for optimizing and improving node ponding condition based on rainwater pipe network topological structure
CN113128008B (en) Sub-basin dividing and coding method reflecting subsurface bedding surface type characteristics
CN117390125B (en) Drainage basin hydrologic forecasting method with intelligent adaptation of flow production mode
CN116663211A (en) Method for converting complex river basin river channel in-out flow topological structure into sequence structure
CN113536504B (en) Well seat type determining method, well seat type determining device, computer equipment and readable storage medium
CN110362876B (en) Complex river network water level field interpolation calculation method based on limited measuring points

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant