CN111667168B - Method for diagnosing running state of drainage system based on liquid level monitoring - Google Patents

Method for diagnosing running state of drainage system based on liquid level monitoring Download PDF

Info

Publication number
CN111667168B
CN111667168B CN202010497087.3A CN202010497087A CN111667168B CN 111667168 B CN111667168 B CN 111667168B CN 202010497087 A CN202010497087 A CN 202010497087A CN 111667168 B CN111667168 B CN 111667168B
Authority
CN
China
Prior art keywords
pipeline
flow
drainage
liquid level
level monitoring
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
CN202010497087.3A
Other languages
Chinese (zh)
Other versions
CN111667168A (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.)
Tongji University
Original Assignee
Tongji University
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 Tongji University filed Critical Tongji University
Priority to CN202010497087.3A priority Critical patent/CN111667168B/en
Publication of CN111667168A publication Critical patent/CN111667168A/en
Application granted granted Critical
Publication of CN111667168B publication Critical patent/CN111667168B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Classifications

    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Administration; Management
    • G06Q10/06Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
    • G06Q10/063Operations research, analysis or management
    • G06Q10/0639Performance analysis of employees; Performance analysis of enterprise or organisation operations
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01MTESTING STATIC OR DYNAMIC BALANCE OF MACHINES OR STRUCTURES; TESTING OF STRUCTURES OR APPARATUS, NOT OTHERWISE PROVIDED FOR
    • G01M99/00Subject matter not provided for in other groups of this subclass
    • G01M99/005Testing of complete machines, e.g. washing-machines or mobile phones
    • GPHYSICS
    • G06COMPUTING; CALCULATING OR COUNTING
    • G06QINFORMATION 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/00Systems or methods specially adapted for specific business sectors, e.g. utilities or tourism
    • G06Q50/10Services
    • G06Q50/26Government or public services
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/10Detection; Monitoring
    • GPHYSICS
    • G16INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
    • G16YINFORMATION AND COMMUNICATION TECHNOLOGY SPECIALLY ADAPTED FOR THE INTERNET OF THINGS [IoT]
    • G16Y40/00IoT characterised by the purpose of the information processing
    • G16Y40/20Analytics; Diagnosis

Abstract

The invention discloses a method for diagnosing an operation state of a drainage system based on liquid level monitoring, and belongs to the technical field of drainage system operation state diagnosis. The method comprises the following steps: (1) Optimally arranging liquid level monitoring points, and dividing a drainage system into a plurality of drainage subareas by taking the liquid level monitoring points as boundaries; (2) Collecting and preprocessing collected drainage pipe network liquid level monitoring data; (3) Taking the preprocessed liquid level monitoring data as an input condition, and calculating a pipeline flow value through a Manning formula or solving a Saint Vigna equation set; (4) Checking the pipeline flow by using the portable flowmeter, and determining the fouling degree of the pipeline; (5) And calculating the sunny day average flow, the underground water infiltration amount and the domestic sewage mixing amount of each drainage partition based on the checked pipeline flow, and performing detailed investigation on the areas with serious underground water infiltration and serious domestic sewage mixing. The invention has higher diagnosis efficiency and wider application range, and is widely applicable to large-scale urban drainage systems in China.

Description

Method for diagnosing running state of drainage system based on liquid level monitoring
Technical Field
The invention relates to a method for diagnosing the running state of a drainage system, in particular to a method for diagnosing the running state of the drainage system based on liquid level monitoring, and belongs to the technical field of running state diagnosis of the drainage system.
Background
At present, most urban drainage pipe network architectures in China are basically formed, but a series of factors such as coexistence of new pipe networks and old pipe networks, design and construction dislocation, improper maintenance and management exist in the urban development process, so that a plurality of problems exist in the running process of the urban drainage pipe network system, such as: the domestic sewage is mixed and connected into a rainwater system to cause the rainwater pump station to be discharged from river in dry days; the rainwater is mixed and connected into a sewage pipe network, so that the water inflow of a sewage treatment plant in a rainy day is higher, and the water inflow concentration is lower; the damage of the pipeline leads to larger infiltration amount of groundwater, and reduces the drainage capacity of a drainage pipe network; the rainwater pipeline is seriously deposited, and the river course is black and odorous in rainy days caused by the pollution of initial rain and the overflow of a pipe network. How to efficiently diagnose the running state problem of the drainage pipe network system, ensure the healthy running of the drainage pipe network, and is the key for improving the service quality and efficiency of the urban drainage pipe network and reducing the water environment pollution of the river course.
There are a great deal of researches on a diagnosis method of the running state of a drainage pipe network at home and abroad. Most of the existing researches focus on rain and sewage mixed connection and groundwater infiltration, the main research methods mainly comprise four types of flow investigation method, tracer method, water quality characteristic factor method and pipeline detection technology, and the following are representative researches:
the literature ([ 1]: field R, pitt R, lalor M, et al, investments of dry-weather pollutant entries into storm-drainage systems, journal of Environmental engineering 1994, vol.120 (05): 1044-1066: [2] Almeida M C, brito R.S. System diagnostics using flow data: quantifying sources and opportunities for performance improvement.Portland:2002.) discloses the use of flow surveys to diagnose the operating condition of a drainage network, primarily by installing flow meters to monitor the flow of the pipes in pipes or important drainage mains where rain and sewage mixing may occur, constructing a drainage network monitoring system with pipe flow as a core, and using these flow data to quantify the rain and sewage mixing or groundwater infiltration. The method has the advantages of intuitively monitoring the flow change of the pipeline and quantitatively analyzing the non-rainwater inflow rate in the rainwater system, but has the following defects: (1) The flowmeter is expensive, is generally only installed near important nodes or discharge ports, and is difficult to popularize in a pipe network; (2) Because of the influence of pipeline uncertainty factors such as silt, the measurement accuracy of the flowmeter is difficult to guarantee. Although the flow survey method is practically applied in many cases, the diagnosis result is not ideal enough, and it is difficult to meet the target requirements of most drainage companies.
The literature (1]Gokhale S,Graham J A.A new development in locating leaks in sanitary sewers.Tunnelling and Underground Space Technology.2004,Vol.19 (01): 85-96. [2]Jewell C.A systematic methodology for the identification and remediation of illegal connections.Proceedings of the Water Environment Federation.2001,Vol.2001 (02): 669-683.) discloses the use of tracer methods for diagnosing the running state of a drainage pipe network, wherein the methods are to accurately judge the illegal inflow condition of a certain section of pipeline or the whole pipe network by adding tracer substances (such as dyes, isotopes and the like) on the upstream of the drainage pipe network and comparing the concentration differences of the upstream and downstream tracer substances. The method has the advantages that the illegal inflow rate in a certain pipeline or the whole pipeline network can be accurately diagnosed, but due to the fact that the drain pipeline network is complicated, the situation that trace substances are difficult to track possibly exists, the accuracy of calculation results is affected, meanwhile, the method has high requirements on operators, is complex to operate and high in detection cost, and the method is limited to be widely popularized.
The literature (1]Lilly L,Stack B P,Caraco D.Pollution loading from illicit sewage discharges in two mid-Atlantic subwatersheds and implications for nutrient and bacterial total maximum daily loads. Watershed Science bulletin.2012, vol.3 (01): 7-17. [2] Xu Zuxin, wang Lingling, yin Hailong, etc.) discloses a characteristic factor-based groundwater infiltration analysis method for a drainage pipe network, university of the same university (Nature Science edition), 2016, vol.44 (04): 593-599.), which is used for diagnosing the running state of the drainage pipe network by using a water quality characteristic factor method, wherein the method utilizes the concentration differences of water quality characteristic factors in different mixed sources such as domestic sewage, industrial wastewater, groundwater and the like, and a chemical mass balance equation is established by monitoring the concentration differences of the water quality characteristic factors at different points in the pipe network and combining the water quality characteristic factor concentrations in the mixed sources, so as to solve the proportion of each mixed source. The method has the advantages that the method is most widely studied in the literature, and the method is based on the solution of a chemical mass conservation equation, has clearer theoretical support, and can analyze the proportions of water from different sources in a system. But it is disadvantageous in that: (1) The difference of the spatial-temporal distribution of the water quality characteristic factors causes the unclosed analysis result and influences the accuracy of the analysis result; (2) The method can only analyze the proportion of each mixed source in the drainage system, and can determine the water quantity of each mixed source by combining a flow investigation method; (3) When there are mixed sources in the system that are not considered, the deviation of the calculation results is large.
The method adopts a pipeline detection technology to diagnose the running state of a drainage pipe network, and uses a pipeline closed circuit television (Closed Circuit Television, CCTV), sonar imaging (Sonar), a through-the-earth radar (Ground Penetrating Radar, GPR), an infrared thermal imaging technology (Infrared Thermography), a pipeline scanning and evaluating technology (The Sewer Scanner and Evaluation Technology, SSET), a multiple sensor (Sewer Assessment with Multi-sensors, SAM), a periscope (Quick View) and other visual technologies to analyze and diagnose the leakage, breakage, mixing and other conditions of the drainage pipe network; the method can intuitively provide the internal running state of the drainage pipeline and timely find out running problems such as damage of the drainage pipeline, pipeline blockage, infiltration of underground water and the like. However, the method has high cost, low efficiency and high consumption of manpower and material resources, cannot detect some sewage stealing and discharging behaviors, and lacks continuous monitoring of a pipe network system.
Disclosure of Invention
Aiming at the problems in the prior art, the invention provides a method for diagnosing the running state of a drainage system based on liquid level monitoring, which can accurately diagnose the running problems of the drainage system such as pipeline siltation, groundwater infiltration, rain and sewage mixed connection and the like only by using inspection well liquid level monitoring data, and has the advantages of simple operation and low cost of monitoring equipment; in addition, the method can also improve the diagnosis efficiency of the drainage system and avoid the interference of rainwater on the system.
In order to achieve the technical purpose, the technical scheme of the invention is as follows:
a method for diagnosing an operational state of a drainage system based on liquid level monitoring, comprising the steps of:
(1) Optimally arranging liquid level monitoring points, and dividing the whole drainage system into a plurality of drainage subareas by taking the liquid level monitoring points as boundaries;
(2) After dividing the drainage subareas, collecting and preprocessing collected drainage pipe network liquid level monitoring data;
(3) Taking the preprocessed liquid level monitoring data as an input condition, and calculating a pipeline flow value through a Manning formula or solving a Saint Vigna equation set;
(4) Checking the pipeline flow by comparing the difference between the calculated pipeline flow value and the actually measured pipeline flow value, and determining the fouling degree of the pipeline;
(5) And calculating the sunny day average flow, the underground water infiltration amount and the domestic sewage mixed connection amount of each drainage partition based on the checked pipeline flow, evaluating the underground water infiltration severity and the domestic sewage mixed connection severity of different drainage partitions, and then carrying out detailed investigation on the areas with serious underground water infiltration and serious domestic sewage mixed connection.
Preferably, the specific steps of optimally arranging the liquid level monitoring points in the step (1) are as follows:
(a) Collecting basic information such as topology information of a drainage pipe network, drainage areas, regional sewage sources and the like;
(b) And selecting a monitoring point arrangement method based on the pipe length or the pipe importance of the monitoring area to arrange the liquid level monitoring points, wherein the pipe importance is the product of the pipe length and the pipe cross section area, and the pipe length and/or the pipe importance in each monitoring area and the pipe length/the pipe importance are smaller than the pipe length/the pipe importance when a certain pipe length/a certain pipe importance is selected as a standard for arranging the liquid level monitoring points.
Preferably, the specific steps of collecting the drainage pipe network liquid level monitoring data in the step (2) are as follows: (a) selecting a sunny day of one consecutive week as a typical week; (b) extracting typical intra-week level gauge level monitoring data.
Preferably, the specific steps of preprocessing the drainage pipe network liquid level monitoring data in the step (2) are as follows:
(a) Filling the missing liquid level monitoring value by adopting linear interpolation;
(b) Adopting STL (Seaseal-Trend decomposition procedure based on Loess) time sequence decomposition algorithm to decompose the liquid level monitoring time sequence into three parts of trend term, period term and residual term, setting a residual term threshold, and taking the monitoring value of the residual term exceeding the threshold as an abnormal value;
(c) The liquid level monitoring abnormal value is replaced by a new value, and the replaced new value only comprises a trend term and a period term, and the residual term is 0.
Preferably, in the step (3), when both the inspection wells on the upstream and the downstream of the pipeline have monitoring data, the pipeline flow is calculated by solving the Saint View equation group, when only the inspection well on the upstream of the pipeline has liquid level monitoring data, the water cross-section area A and the hydraulic radius R of the pipeline are calculated by using the liquid level of the inspection well on the upstream, and the pipeline gradient is used for replacing the friction gradient S f Directly through Manning formulaAnd calculating the pipeline flow.
Further preferably, the specific steps of calculating the pipeline flow by solving the san View equation group are as follows:
(a) The san View equation set for solving the open channel unsteady flow problem can be expressed as:
continuity equation:
momentum equation:
wherein A is the water cross-section area of a pipeline, t is the time step, Q is the pipeline flow, x is the pipeline length, Q is the pipeline side flow per unit length, z is the water head, g is the gravity acceleration, S f Is the friction gradient (head loss per unit length);
(b) Solving the san veland equation set by using a finite difference method to obtain a new partial differential equation for calculating the pipeline flow:
friction gradient S f The Manning formula, which simulates constant uniform flow, can be expressed as:
wherein A is the cross-sectional area of water, R is the hydraulic radius, and n is the Manning roughness coefficient; substituting a Manning formula into a partial differential equation, and replacing the differential by utilizing the differential to obtain an implicit differential form of the equation, wherein the implicit differential form is as follows:
will Q t+Δt The equation form for the iterative solution obtained by the term transfer is:
in the method, in the process of the invention,for the average flow rate of the pipeline, +.>L is the length of the pipeline, g is the gravitational acceleration, R is the hydraulic radius, n is the Manning roughness coefficient, z 2 、z 1 At t+Δt respectivelyMarking up and down stream inspection well water head +.>The water cross-sectional areas of the upstream and downstream inspection wells at the time t+delta t are respectively;
and (3) carrying out multiple iterations on the formula I, and solving the pipeline flow at the time t+delta t, wherein the initial iteration value is the pipeline constant flow at the time t=0, and stopping iteration after the two iteration errors are smaller than 0.001 or the iteration times are more than 10, wherein the solved Qtdt is the pipeline flow at the time t+delta t.
Further preferably, when the pipe is pressureless, the water flow section a is calculated from the water depth h upstream and downstream of the pipe, the water flow section a=d 2 And/8 x (θ—sin θ), where D is the pipe diameter, θ is the filling angle, θ=2cos -1 (1-2*h/D); when the pipeline generates pressure flow, a virtual Prussian slit is added to the top end of the pipeline, and the relation between the width of the Prussian slit and the depth of water in the pipeline can be expressed as:
when 1<When the h/D is less than or equal to 1.78,
when h/D>1.78, W slot =0.01D
Wherein W is slot The width of the Prussian slit is shown, D is the pipe diameter, and h is the depth of water in the pipe.
At this time, the water flow cross section a=pi D 2 /4+(h-D)*W slot
Preferably, the specific steps of step (4) are:
(a) Assuming that the difference between the calculated flow and the actually measured flow is only caused by sediment accumulation at the lower end of the pipeline, the pipeline accumulation only affects the elevation of the bottom of the downstream pipeline, selecting a certain moment in a typical period, measuring the instantaneous flow of the pipeline by using the portable flowmeter, and comparing the instantaneous flow with the calculated flow at the moment;
(b) And continuously correcting the elevation of the downstream pipe bottom of the pipeline to ensure that the calculated flow at the monitoring moment is equal to the actually measured flow, wherein the difference value between the elevation of the downstream pipe bottom of the pipeline and the elevation of the pipe bottom in the basic data is the thickness of the pipeline siltation layer.
Preferably, the step (5) of calculating the sunny day average flow, the groundwater infiltration amount and the domestic sewage mixing amount of each drainage partition based on the checked pipeline flow, and evaluating the groundwater infiltration severity and the domestic sewage mixing severity of different drainage partitions comprises the following specific steps:
(a) Taking the pipeline flow of 2:00-4:00 in the early morning as the minimum flow at night, calculating the groundwater infiltration amount of each drainage partition according to a minimum flow method at night, then calculating the infiltration load of unit pipe length in each drainage partition, and then evaluating the groundwater infiltration severity degree of different drainage partitions according to the unit pipe length bearing standard;
(b) Subtracting daily groundwater infiltration amount from typical daily average flow in week to obtain domestic sewage mixing and connecting amounts of different partitions, and evaluating the severity of domestic sewage mixing and connecting of different partitions.
Preferably, the specific steps of the step (5) for performing detailed investigation on the areas with serious groundwater infiltration and serious domestic sewage mixing are as follows: aiming at the subareas with groundwater infiltration amount exceeding the infiltration standard, carrying out CCTV pipeline detection, and checking the integrity degree and connection condition of the pipelines; aiming at the region with the most serious mixed joint of domestic sewage, the source tracing investigation of the domestic sewage source is adopted.
From the above description, it can be seen that the present invention has the following advantages:
(1) The invention can accurately diagnose the running problems of the drainage system such as pipeline siltation, groundwater infiltration, rain and sewage mixed connection and the like only by using inspection well liquid level monitoring data, has high diagnosis efficiency, low cost and wider application range, and is suitable for drainage companies with limited budget and long-term monitoring of the running state of a pipe network.
(2) The invention adopts monitoring data of typical weeks, avoids the interference of rainwater on a system, and reduces the influence of domestic sewage discharge difference on diagnosis results on workdays and rest days.
(3) The invention divides the drainage system into a plurality of subareas by taking the monitoring points as boundaries, identifies the infiltration amount of the underground water and the mixed connection condition of rain and sewage in each area, improves the diagnosis efficiency of the drainage system, and provides references for subsequent mixed connection investigation, pipeline repair and other works.
Drawings
FIG. 1 is a general flow chart for diagnosing an operating condition of a drain system based on fluid level monitoring in accordance with the present invention;
FIG. 2 is a flow chart of iterative computation of pipeline flow;
FIG. 3 is a schematic cross-sectional view of a non-full flow round tube;
FIG. 4 is a schematic diagram of a Prinseman slit;
FIG. 5 is a schematic view of the liquid level monitoring point locations and drain zones in example 2;
FIG. 6 is a graph of pre-treatment fluid level monitoring data for fluid level monitoring Point # 3 of example 2;
FIG. 7 is a graph of liquid level monitoring data after processing at liquid level monitoring point No. 3 in example 2;
Detailed Description
The features of the invention are further illustrated by way of example only and not by way of limitation in the claims.
A method for diagnosing an operational state of a drainage system based on liquid level monitoring, comprising the steps of:
(1) Optimally arranging liquid level monitoring points, and dividing the whole drainage system into a plurality of drainage subareas by taking the liquid level monitoring points as boundaries;
the concrete steps of optimally arranging the liquid level monitoring points are as follows:
(a) Collecting basic information such as topology information of a drainage pipe network, drainage areas, regional sewage sources and the like;
(b) Selecting a monitoring point arrangement method based on the pipe length or pipe importance (pipe importance is the product of the pipe length and the pipe cross-sectional area) of the monitoring area to arrange liquid level monitoring points, wherein when a certain pipe length/certain pipe importance is selected as a standard for arranging liquid level monitoring points, the pipe length and/or pipe importance in each monitoring area is smaller than the pipe length/pipe importance; for example, when the pipeline length is 4km as a standard arrangement liquid level monitoring point, the sum of the pipeline lengths in each monitoring area is smaller than 4km;
(2) After dividing the drainage subareas, collecting and preprocessing collected drainage pipe network liquid level monitoring data;
the specific steps of collecting the liquid level monitoring data of the drainage pipe network are as follows:
(a) Selecting a sunny day of one week as a typical week;
(b) Extracting liquid level monitoring data of a typical pericycle liquid level meter;
the specific steps of preprocessing the drainage pipe network liquid level monitoring data are as follows:
(a) Filling the missing liquid level monitoring value by adopting linear interpolation;
(b) Adopting STL (Seaseal-Trend decomposition procedure based on Loess) time sequence decomposition algorithm to decompose the liquid level monitoring time sequence into three parts of trend term, period term and residual term, setting a residual term threshold, and taking the monitoring value of the residual term exceeding the threshold as an abnormal value;
(c) Replacing the liquid level monitoring abnormal value with a new value, wherein the replaced new value only comprises a trend item and a period item, and the residual item is 0;
(3) Calculating a pipeline flow value by taking the preprocessed liquid level monitoring data as an input condition and using a Manning formula or solving a Santa-Vigna equation set, wherein when two inspection wells on the upstream and downstream of the pipeline have monitoring data, the pipeline flow is calculated by solving the Santa-Vigna equation set, and when only the inspection well on the upstream of the pipeline has liquid level monitoring data, the water cross section area A and the hydraulic radius R of the pipeline are calculated by using the liquid level of the inspection well on the upstream of the pipeline, and the pipeline gradient is used for replacing the friction gradient S f Directly through Manning formulaCalculating the pipeline flow;
the specific steps for calculating the pipeline flow by solving the san View equation group are as follows:
(a) The san View equation set for solving the open channel unsteady flow problem can be expressed as:
continuity equation:
momentum equation:
wherein A is the area of the water section of the pipeline, t is the time step, Q is the pipeline flow, x is the pipeline length, Q is the pipeline side flow per unit length, z is the water head, g is the gravitational acceleration, S f Is the friction gradient (head loss per unit length);
(b) Solving the san veland equation set by using a finite difference method to obtain a new partial differential equation for calculating the pipeline flow:
friction gradient S f The Manning formula, which simulates constant uniform flow, can be expressed as:
wherein A is the cross-sectional area of water, R is the hydraulic radius, and n is the Manning roughness coefficient; substituting a Manning formula into a partial differential equation, and replacing the differential by utilizing the differential to obtain an implicit differential form of the equation, wherein the implicit differential form is as follows:
will Q t+Δt The equation form for the iterative solution obtained by the term transfer is:
in the method, in the process of the invention,for the average flow rate of the pipeline, +.>L is the length of the pipeline, g isGravity acceleration, R is hydraulic radius, n is Manning roughness coefficient, z 2 、z 1 Upstream and downstream inspection well water heads at the time t+delta t respectively, < >>The water cross-sectional areas of the upstream and downstream inspection wells at the time t+delta t are respectively;
wherein, when the pipeline has no pressure flow, the water flow section A is calculated from the water depth h of the upstream and downstream of the pipeline, as shown in fig. 3, the water flow section A=D 2 And/8 x (θ—sin θ), where D is the pipe diameter, θ is the filling angle, θ=2cos -1 (1-2*h/D); when the pressure flow occurs in the pipeline, an imaginary opening slit is added to the top end of the pipeline, which is called a Primann slit, as shown in fig. 4, the Primann slit can keep the pipeline on the free water surface, the slit enables the section of the pipeline to be larger than the actual section, so that the formula of open channel flow is still applicable under the assumption, and the relation between the width of the Primann slit and the depth of water in the pipeline can be expressed as:
when 1<When the h/D is less than or equal to 1.78,
when h/D>1.78, W slot =0.01D
Wherein W is slot The width of the Prussian slit is shown, D is the pipe diameter, and h is the depth of water in the pipe.
At this time, the water flow cross section a=pi D 2 /4+(h-D)*W slot
Carrying out multiple iterations on the formula I, solving the pipeline flow at the time t+delta t, wherein the initial iteration value is the constant flow of the pipeline at the time t=0, stopping iteration after the two iteration errors are smaller than 0.001 or the iteration times are more than 10 times, and obtaining Qtdt which is the pipeline flow at the time t+delta t at the moment, and carrying out iterative calculation on the pipeline flow at the time t+delta t according to an iteration equation, wherein a flow chart is shown in figure 2;
(4) Firstly, selecting a certain moment in a typical week, measuring the instantaneous flow of a pipeline by using a portable flowmeter, and comparing the instantaneous flow with the calculated flow at the moment (assuming that the difference between the calculated flow and the actually measured flow is only caused by sediment accumulation at the lower end of the pipeline, and the pipeline accumulation only affects the elevation of the bottom of a pipeline downstream); then continuously correcting the elevation of the downstream pipe bottom of the pipeline to ensure that the calculated flow at the monitoring moment is equal to the actually measured flow, wherein the difference value between the elevation of the downstream pipe bottom of the pipeline and the elevation of the pipe bottom in the basic data is the thickness of a pipeline siltation layer; for example, in the basic data, the elevation of the downstream pipe bottom is 0.9m, the corrected pipe bottom elevation is 1.1m, and the thickness of the pipeline deposited layer is 0.2m;
(5) Calculating sunny day average flow, underground water infiltration amount and domestic sewage mixed connection amount of each drainage partition based on the checked pipeline flow, evaluating the underground water infiltration severity and the domestic sewage mixed connection severity of different drainage partitions, and then performing detailed investigation on the areas with serious underground water infiltration and serious domestic sewage mixed connection;
the method comprises the specific steps of calculating the sunny day average flow, the underground water infiltration amount and the domestic sewage mixing amount of each drainage partition based on the checked pipeline flow, and evaluating the underground water infiltration severity and the domestic sewage mixing severity of different drainage partitions:
(a) Taking the pipeline flow of 2:00-4:00 in the early morning as the minimum flow at night (the minimum flow at night comprises 10% of daily sewage mixed connection amount and underground water infiltration amount), calculating the underground water infiltration amount of each drainage partition according to a minimum flow method at night, then calculating the infiltration load of unit pipe length in each drainage partition, and taking 50-75m according to the unit pipe length bearing standard (the unit pipe length bearing standard 3 /(km.d)) evaluating the severity of groundwater infiltration for different drainage zones;
(b) Subtracting daily groundwater infiltration amount by using typical daily average flow (calculated from t=0 moment and delta t as calculation time step according to solved t+delta t moment pipeline flow) to obtain domestic sewage mixed connection amounts of different partitions, and evaluating domestic sewage mixed connection severity of different partitions;
the specific steps for carrying out detailed investigation on the areas with serious groundwater infiltration and serious domestic sewage mixed connection are as follows: aiming at the subareas with groundwater infiltration amount exceeding the infiltration standard, carrying out CCTV pipeline detection, and checking the integrity degree and connection condition of the pipelines; aiming at the region with the most serious mixed joint of domestic sewage, the source tracing investigation of the domestic sewage source is adopted.
Example 2
Taking a certain diversion rainwater pipe network in the northeast region T of Shanghai city as an example, the implementation process of the method for diagnosing the running state of the drainage system based on liquid level monitoring is further described.
The pipe network comprises 1102 nodes and 1101 pipes, only one outlet is discharged into the CF municipal pipe network, the total water collecting area is 21ha, the total pipe length is 12.55km, and the pipe diameter range is DN50-DN800. The area is a resident living area, and the clear water and rain water pipeline has obvious water flow.
(1) Optimally arranging liquid level monitoring points, and dividing the whole drainage system into a plurality of drainage subareas by taking the liquid level monitoring points as boundaries;
select pipe importance 150m 3 As a monitoring point arrangement standard, 8 liquid level monitoring points are arranged in the whole area, the whole area is divided into 8 drainage partitions, and a partition schematic diagram is shown in fig. 5. And the No. 6 monitoring point data is abnormal and is an invalid monitoring point, and finally the monitoring data of 7 effective monitoring points are obtained.
(2) Collecting and preprocessing collected drainage pipe network liquid level monitoring data;
the time interval of the monitoring data is 1min, each monitoring point samples 24h/1 min=1440 data every day, about 10-20 data are deleted every day, and the deletion value is filled by a linear interpolation method. The time sequence of the monitoring data is decomposed into a period term, a trend term and a residual term, the larger value of the residual term is marked by using a five-pointed star, as shown in fig. 6, the residual term is removed for the detected abnormal value, only the period term and the trend term are reserved, and the original abnormal value is replaced by a new value. And a complete week from 11 months 11 to 17 months 11 in 2019 is selected as a typical week, and no rainfall occurs in the week, so that the interference of rainwater on a system is avoided.
(3) Calculating a pipeline flow value by taking the preprocessed liquid level monitoring data (shown in fig. 7) as an input condition and using a Manning formula or solving a Saint Vietnam equation set;
(4) Checking the pipeline flow by comparing the difference between the calculated pipeline flow value and the actually measured pipeline flow value, and determining the fouling degree of the pipeline;
through portable flowmeter, measure each monitoring point downstream pipeline flow one by one, comparison calculation flow and actual measurement flow find, and No. 4 monitoring point and No. 8 monitoring point actual measurement flow are great with calculation flow difference, and other point actual measurement flow is comparatively close with calculation flow. And (3) correcting the elevation of the bottom of the downstream pipeline of the pipeline to ensure that the calculated flow at the monitoring moment is equal to the monitored flow, so that the thickness of the deposited layer of the downstream pipeline of the No. 4 monitoring point is about 5.5cm, and the thickness of the deposited layer of the downstream pipeline of the No. 8 monitoring point is about 1.8cm.
(5) And calculating the sunny day average flow, the underground water infiltration amount and the domestic sewage mixed connection amount of each drainage partition based on the checked pipeline flow, evaluating the underground water infiltration severity and the domestic sewage mixed connection severity of different drainage partitions, and then carrying out detailed investigation on the areas with serious underground water infiltration and serious domestic sewage mixed connection.
After checking the pipeline flow, obtaining the data of the sunny flow of each partition, calculating the daily average flow by using the data of the typical weekly flow, wherein the daily average flow of the rainwater system on sunny days is 1938.67m 3 And/d. The minimum flow rate at night was selected from 2:00-4:00 a.m., and the groundwater infiltration amounts of the respective zones were calculated as shown in Table 1. Wherein the groundwater infiltration load of the regions 2, 3, 4 exceeds the infiltration standard, and the infiltration load of the other regions is within the infiltration standard.
TABLE 1 groundwater infiltration capacity for each zone
Based on the total mixed amount, the infiltration amount of the groundwater was subtracted for 24 hours to obtain the mixed amount of domestic sewage in each region, as shown in Table 2. Wherein, region 2 belongs to heavy misce bene, region 4 and region 8 belong to moderate misce bene, and the other regions belong to light misce bene.
Table 2 area domestic sewage mixing and receiving amount
Analysis and diagnosis of the whole system show that in the sunny flow of the rainwater system, groundwater infiltration accounts for 53.61%, and domestic sewage mixed connection accounts for 46.39%. Carrying out CCTV pipeline detection on pipelines with serious groundwater infiltration, and finding out the phenomenon of dislocation of different degrees at the joint of the pipelines, wherein obvious groundwater infiltration marks appear in partial areas; the regional traceability survey of severe domestic sewage mixed connection finds out two domestic sewage mixed connection sources.
It is to be understood that the foregoing detailed description of the invention is merely illustrative of the invention and is not limited to the embodiments of the invention. It will be understood by those of ordinary skill in the art that the present invention may be modified or substituted for elements thereof to achieve the same technical effects; as long as the use requirement is met, the invention is within the protection scope of the invention.

Claims (5)

1. A method for diagnosing an operational condition of a drainage system based on fluid level monitoring, comprising the steps of:
(1) Optimally arranging liquid level monitoring points, and dividing the whole drainage system into a plurality of drainage subareas by taking the liquid level monitoring points as boundaries;
(2) After dividing the drainage subareas, collecting and preprocessing collected drainage pipe network liquid level monitoring data;
(3) Taking the preprocessed liquid level monitoring data as an input condition, and calculating a pipeline flow value through a Manning formula or solving a Saint Vigna equation set;
(4) Checking the pipeline flow by comparing the difference between the calculated pipeline flow value and the actually measured pipeline flow value, and determining the fouling degree of the pipeline;
(5) Calculating sunny day average flow, underground water infiltration amount and domestic sewage mixed connection amount of each drainage partition based on the checked pipeline flow, evaluating the underground water infiltration severity and the domestic sewage mixed connection severity of different drainage partitions, and then performing detailed investigation on the areas with serious underground water infiltration and serious domestic sewage mixed connection;
the specific steps of optimally arranging the liquid level monitoring points in the step (1) are as follows:
(a) Collecting topology information of a drainage pipe network, a drainage area and an area sewage source;
(b) Selecting a monitoring point arrangement method based on the pipe importance of a monitoring area to arrange liquid level monitoring points, wherein the pipe importance is the product of the pipe length and the pipe cross section area, and when a certain pipe importance is selected as a standard to arrange the liquid level monitoring points, the pipe importance in each monitoring area is smaller than the pipe importance;
the specific steps of calculating the sunny day average flow, the underground water infiltration amount and the domestic sewage mixing amount of each drainage partition based on the checked pipeline flow in the step (5) and evaluating the underground water infiltration severity and the domestic sewage mixing severity of different drainage partitions are as follows:
(a) Taking the pipeline flow of 2:00-4:00 in the early morning as the minimum flow at night, calculating the groundwater infiltration amount of each drainage partition according to a minimum flow method at night, then calculating the infiltration load of unit pipe length in each drainage partition, and then evaluating the groundwater infiltration severity degree of different drainage partitions according to the unit pipe length bearing standard;
(b) Subtracting daily groundwater infiltration amount from typical daily average flow in week to obtain domestic sewage mixed connection amounts of different partitions, and evaluating the domestic sewage mixed connection severity of different partitions;
in the step (3), when two inspection wells on the upstream and downstream of the pipeline have monitoring data, calculating the pipeline flow by solving the Saint View equation group, when only the inspection well on the upstream of the pipeline has liquid level monitoring data, calculating the water cross-section area A and the hydraulic radius R of the pipeline by utilizing the liquid level of the inspection well on the upstream, and replacing the friction gradient S with the pipeline gradient f Directly through Manning formulaCalculating the pipeline flow;
the specific steps for calculating the pipeline flow by solving the san View south equation group are as follows:
(a) The san View equation set for solving the open channel unsteady flow problem can be expressed as:
continuity equation:
momentum equation:
wherein A is the water cross-section area of a pipeline, t is the time step, Q is the pipeline flow, x is the pipeline length, Q is the pipeline side flow per unit length, z is the water head, g is the gravity acceleration, S f Is the friction gradient (head loss per unit length);
(b) Solving the san veland equation set by using a finite difference method to obtain a new partial differential equation for calculating the pipeline flow:
friction gradient S f The Manning formula, which simulates constant uniform flow, can be expressed as:
wherein A is the cross-sectional area of water, R is the hydraulic radius, and n is the Manning roughness coefficient; substituting a Manning formula into a partial differential equation, and replacing the differential by utilizing the differential to obtain an implicit differential form of the equation, wherein the implicit differential form is as follows:
will Q t+Δt The equation form for the iterative solution obtained by the term transfer is:
in the method, in the process of the invention,for the average flow rate of the pipeline, +.>L is the length of the pipeline, g is the gravitational acceleration, R is the hydraulic radius, n is the Manning roughness coefficient, z 2 、z 1 Upstream and downstream inspection well water heads at the time t+delta t respectively, < >>The water cross-sectional areas of the upstream and downstream inspection wells at the time t+delta t are respectively;
carrying out multiple iterations on the formula (I), and solving the pipeline flow at the time t+delta t, wherein the initial iteration value is the constant flow of the pipeline at the time t=0, and stopping iteration after the two iteration errors are smaller than 0.001 or the iteration times are more than 10, wherein the solved Qtdt is the pipeline flow at the time t+delta t;
when the pipeline has no pressure flow, the water flow section A is calculated by the water depth h of the upstream and downstream of the pipeline, and the water flow section A=D 2 And/8 x (θ—sin θ), where D is the pipe diameter, θ is the filling angle, θ=2cos -1 (1-2*h/D); when the pipeline generates pressure flow, a virtual Prussian slit is added to the top end of the pipeline, and the relation between the width of the Prussian slit and the depth of water in the pipeline can be expressed as:
when 1<When the h/D is less than or equal to 1.78,
when h/D>1.78, W slot =0.01D
Wherein W is slot The width of the Prussian slit is D, the pipe diameter is D, and h is the depth of water in the pipe;
at this time, the water flow cross section a=pi D 2 /4+(h-D)*W slot
2. The method of claim 1, wherein the specific steps of collecting the drainage network level monitoring data in step (2) are: (a) selecting a sunny day of one consecutive week as a typical week; (b) extracting typical intra-week level gauge level monitoring data.
3. The method of claim 1, wherein the specific steps of preprocessing the drainage network level monitoring data in step (2) are:
(a) Filling the missing liquid level monitoring value by adopting linear interpolation;
(b) Adopting STL (Seaseal-Trend decomposition procedure based on Loess) time sequence decomposition algorithm to decompose the liquid level monitoring time sequence into three parts of trend term, period term and residual term, setting a residual term threshold, and taking the monitoring value of the residual term exceeding the threshold as an abnormal value;
(c) The liquid level monitoring abnormal value is replaced by a new value, and the replaced new value only comprises a trend term and a period term, and the residual term is 0.
4. The method of claim 1, wherein the specific steps of step (4) are:
(a) Assuming that the difference between the calculated flow and the actually measured flow is only caused by sediment accumulation at the lower end of the pipeline, the pipeline accumulation only affects the elevation of the bottom of the downstream pipeline, selecting a certain moment in a typical period, measuring the instantaneous flow of the pipeline by using the portable flowmeter, and comparing the instantaneous flow with the calculated flow at the moment;
(b) And continuously correcting the elevation of the downstream pipe bottom of the pipeline to ensure that the calculated flow at the monitoring moment is equal to the actually measured flow, wherein the difference value between the elevation of the downstream pipe bottom of the pipeline and the elevation of the pipe bottom in the basic data is the thickness of the pipeline siltation layer.
5. The method as set forth in claim 1, wherein the specific steps of the step (5) for examining the areas with serious groundwater infiltration and serious mixing of domestic sewage are as follows: aiming at subareas with groundwater infiltration amount exceeding infiltration standard, carrying out CCTV (closed circuit television) CCTV pipeline detection, and checking the integrity degree and connection condition of the pipelines; aiming at the region with the most serious mixed joint of domestic sewage, the source tracing investigation of the domestic sewage source is adopted.
CN202010497087.3A 2020-06-04 2020-06-04 Method for diagnosing running state of drainage system based on liquid level monitoring Active CN111667168B (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010497087.3A CN111667168B (en) 2020-06-04 2020-06-04 Method for diagnosing running state of drainage system based on liquid level monitoring

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010497087.3A CN111667168B (en) 2020-06-04 2020-06-04 Method for diagnosing running state of drainage system based on liquid level monitoring

Publications (2)

Publication Number Publication Date
CN111667168A CN111667168A (en) 2020-09-15
CN111667168B true CN111667168B (en) 2024-01-23

Family

ID=72385908

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010497087.3A Active CN111667168B (en) 2020-06-04 2020-06-04 Method for diagnosing running state of drainage system based on liquid level monitoring

Country Status (1)

Country Link
CN (1) CN111667168B (en)

Families Citing this family (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112525255B (en) 2020-11-11 2022-08-30 中国电建集团华东勘测设计研究院有限公司 Concealed culvert water outlet inspection process and drainage port tracing method
CN113111478B (en) * 2021-05-13 2024-02-20 北控水务(中国)投资有限公司 Evaluation method and equipment for mixed connection, inflow and infiltration degree of drainage system pipe network
CN113392523B (en) * 2021-06-14 2022-08-16 云南省设计院集团有限公司 Sewage pipe network health condition diagnosis method based on long-duration multi-measuring-point
CN113553686B (en) * 2021-08-03 2024-02-13 中水珠江规划勘测设计有限公司 Method and device for determining fault pipeline and electronic equipment
CN113760880A (en) * 2021-09-07 2021-12-07 天津大学 Pretreatment method of water quality automatic monitoring data
IT202200001325U1 (en) 2022-03-30 2022-06-30 Dyrecta Lab S R L Non-invasive diagnostic system for real-time recognition of water leaks close to the surface in urban and suburban conditions
CN116050037B (en) * 2023-01-13 2024-04-02 三峡智慧水务科技有限公司 Urban drainage system liquid level indirect monitoring and analyzing method based on directed topology network
CN115876409B (en) * 2023-02-16 2023-05-16 广州中工水务信息科技有限公司 Sewage pipeline leakage monitoring and analyzing system and method
CN116821845A (en) * 2023-06-20 2023-09-29 郑州大学 Pipeline siltation condition diagnosis method and device based on multi-sensor data fusion
CN117453800A (en) * 2023-11-09 2024-01-26 哈尔滨航天恒星数据系统科技有限公司 Clogging analysis method based on rainwater pipe network model
CN117890462A (en) * 2024-03-14 2024-04-16 陕西多奇电子科技有限公司 Mining sump liquid mud double-position monitoring system and monitoring method

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039565B1 (en) * 2003-01-03 2006-05-02 Haestad Methods, Inc. Method and system for developing a numerical dynamic sanitary sewer and storm water drainage simulation model
CN102183895A (en) * 2011-03-24 2011-09-14 杭州电子科技大学 Novel modeling and controlling method of drain pipe network hydraulics system
CN104236626A (en) * 2014-09-04 2014-12-24 北京清控人居环境研究院有限公司 Integrated online monitoring system for drainage pipeline liquid level and flow
CN204679098U (en) * 2015-06-17 2015-09-30 上海水顿智能科技有限公司 A kind of experimental provision analyzing drainage pipeline alluvial characteristic
CN105651484A (en) * 2015-12-29 2016-06-08 北京无线电计量测试研究所 Adaptive rainwater pipe network monitoring method and adaptive rainwater pipe network monitoring system
CN106382471A (en) * 2016-11-25 2017-02-08 上海市城市排水有限公司 Municipal drainage pipe network diagnostic assessment method giving consideration to key node
CN108151791A (en) * 2018-01-25 2018-06-12 上海水顿智能科技有限公司 A kind of method of analysis conduit deposting distribution
CN108278491A (en) * 2017-12-18 2018-07-13 苏州航天系统工程有限公司 A kind of method and system finding drainage pipeline networks operation exception
CN108755898A (en) * 2018-05-30 2018-11-06 上海水顿智能科技有限公司 A method of water penetrates under analysis buried pipe entoilage and conduit wastewater exudation is distributed
CN108984921A (en) * 2018-07-24 2018-12-11 中国水利水电科学研究院 A kind of river water level method of calculating flux and device that combination underground water infiltrates

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7039565B1 (en) * 2003-01-03 2006-05-02 Haestad Methods, Inc. Method and system for developing a numerical dynamic sanitary sewer and storm water drainage simulation model
CN102183895A (en) * 2011-03-24 2011-09-14 杭州电子科技大学 Novel modeling and controlling method of drain pipe network hydraulics system
CN104236626A (en) * 2014-09-04 2014-12-24 北京清控人居环境研究院有限公司 Integrated online monitoring system for drainage pipeline liquid level and flow
CN204679098U (en) * 2015-06-17 2015-09-30 上海水顿智能科技有限公司 A kind of experimental provision analyzing drainage pipeline alluvial characteristic
CN105651484A (en) * 2015-12-29 2016-06-08 北京无线电计量测试研究所 Adaptive rainwater pipe network monitoring method and adaptive rainwater pipe network monitoring system
CN106382471A (en) * 2016-11-25 2017-02-08 上海市城市排水有限公司 Municipal drainage pipe network diagnostic assessment method giving consideration to key node
CN108278491A (en) * 2017-12-18 2018-07-13 苏州航天系统工程有限公司 A kind of method and system finding drainage pipeline networks operation exception
CN108151791A (en) * 2018-01-25 2018-06-12 上海水顿智能科技有限公司 A kind of method of analysis conduit deposting distribution
CN108755898A (en) * 2018-05-30 2018-11-06 上海水顿智能科技有限公司 A method of water penetrates under analysis buried pipe entoilage and conduit wastewater exudation is distributed
CN108984921A (en) * 2018-07-24 2018-12-11 中国水利水电科学研究院 A kind of river water level method of calculating flux and device that combination underground water infiltrates

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
文碧岚.排水管道淤积状况模拟分析.给水排水.2015,第41卷(第7期),153-155页. *
沈伟康、滕严婷.一种基于液位监测诊断排水系统运行状态的方法.城市道桥与防洪.2017,84-86页. *

Also Published As

Publication number Publication date
CN111667168A (en) 2020-09-15

Similar Documents

Publication Publication Date Title
CN111667168B (en) Method for diagnosing running state of drainage system based on liquid level monitoring
Zhao et al. Pin-pointing groundwater infiltration into urban sewers using chemical tracer in conjunction with physically based optimization model
Beheshti et al. Infiltration/inflow assessment and detection in urban sewer system
Hoes et al. Locating illicit connections in storm water sewers using fiber-optic distributed temperature sensing
CN111501953A (en) Exception analysis method, system, medium and equipment for drainage pipe network
Xu et al. Quantification of non-stormwater flow entries into storm drains using a water balance approach
Beheshti et al. Detection of extraneous water ingress into the sewer system using tandem methods–a case study in Trondheim city
CN113836622B (en) Drainage pipe network information management and comprehensive application system based on GIS and BIM
CN109242367B (en) Urban sewage treatment rate evaluation and calculation method
Heiderscheidt et al. Stable water isotopes as a tool for assessing groundwater infiltration in sewage networks in cold climate conditions
Mannina et al. Uncertainty in sewer sediment deposit modelling: Detailed vs simplified modelling approaches
Tan et al. Assessment and pathway determination for rainfall-derived inflow and infiltration in sanitary systems: a case study
CN111008357B (en) Method for determining average blocking height of sediment in sewage pipeline
CN115931057B (en) Underground pipeline emission early warning system based on urban water conservancy data
Guo et al. River water intrusion as a source of inflow into the sanitary sewer system
Regueiro-Picallo et al. Towards urban drainage sediment accumulation monitoring using temperature sensors
KR100353391B1 (en) An integrated method for providing monitoring/management data of a sewer system and a system thereof
KR20060083388A (en) A method analyzing infiltration/inflow using chloride ion concentrations and flow and a system monitoring sewer drainage pipes using the same
CN115650327A (en) Method for stripping inflow of dry weather and outdoor water of urban sewage system
Larrarte Suspended solids within sewers: an experimental study
KR102080066B1 (en) System for monitoring illegal waste- water discharge
Panasiuk et al. Identifying sources of infiltration and inflow in sanitary sewers in a northern community: comparative assessment of selected methods
KR20040072277A (en) A sewer pipe maintenance system and construction method thereof
CN106682271B (en) Method for determining SWMM water quality scouring model parameters
Raimondi et al. Uncertainty on flow rate and temperature measurement for the detection of illicit flows in sewers

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