A method and device for tracing the pollution source of a rainwater outlet

Abstract

A method and device for tracing the pollution source of a rainwater outlet, relating to the technical field of computer-aided pollution tracing, a method for tracing the pollution source of a rainwater outlet, comprising the following steps: obtaining pipe network data and first outlet data; determining monitoring factors according to the first outlet data; determining encrypted monitoring areas according to the pipe network data and monitoring factors, and obtaining pollution source coordinates; obtaining pollution source data and second outlet data corresponding to the pollution source coordinates; calculating the pollution load of the pollution source and the pollution load of the outlet according to the pollution source data and the second outlet data; calculating the determination coefficient of the pollution source and the outlet according to the pollution load of the pollution source and the pollution load of the outlet, and judging whether the tracing is completed according to the determination coefficient. The present application solves the problem that the current sewage tracing method cannot verify the accuracy of the results by calculating the determination coefficient of the pollution source and the outlet to determine the correlation between the pollution load of the pollution source and the pollution load of the rainwater outlet.

Description

Technical Field

[0001] The invention belongs to the technical field of computer-aided pollution source tracing, and in particular to a pollution source tracing method and device for a rainwater outlet. Background Art

[0002] The urban rainwater pipe network is an important part of the urban drainage system. It is mainly used to collect and discharge rainwater runoff in the city. It includes rainwater inlets for collecting rainwater on the ground, branch pipes and main pipes for transporting rainwater, and rainwater outlets for discharging the collected rainwater into rivers and lakes. Its operating status directly affects the city’s flood control and drainage capacity, residents’ quality of life and the ecological environment.

[0003] However, the current urban rainwater pipe network faces many challenges. For example, due to serious aging, insufficient design standards, inadequate maintenance and management, and the combination of rainwater and sewage, there are still a large number of functional defects such as siltation, mud accumulation, and full pipe flow in the drainage network, as well as structural defects such as mixing of rainwater and sewage, which cause a large amount of pollutants to be directly discharged into the water body, leading to frequent black and smelly problems in urban water bodies.

[0004] Existing pollution source tracing methods for rainwater outlets mainly include: manual investigation and tracing, that is, tracing the source by technical personnel opening the manhole covers one by one on site for observation or sampling and analysis; tracing the source with pipeline detection equipment, that is, mainly using closed-circuit television monitoring systems, pipeline periscopes, sonar, ground-penetrating radar and other equipment to conduct comprehensive inspections of pipelines; and computer-assisted tracing.

[0005] For example, a Chinese patent with publication number CN113947033B discloses a system and method for tracing pollutants in a drainage network based on artificial intelligence, and the method steps include: distributing monitoring points in the drainage network and obtaining pollutant concentration value data monitored by the monitoring points; a pollutant simulation operation module, which simulates the state of pollutants in the drainage pipe through the SWMM simulator according to the pollutant concentration value data to obtain simulation results; the drainage network pollution analysis module processes the data in the drainage network to obtain the flow direction with the highest pollution concentration in the drainage network; the sewage flow direction processing module obtains the concentration value data of pollutants closest to the sewage treatment plant, and predicts the concentration value data when the sewage flows into the sewage treatment plant; it uses the machine learning decision tree method to classify the flow direction of the largest pollutants in the drainage network, realizing fully automated prediction and classification analysis of the drainage network without human intervention.

[0006] However, current sewage source tracing methods do not verify the accuracy of the results. For example, when there are multiple pollutants in the pipe network, the pollution source is easily missed, resulting in poor accuracy of the tracing results and poor reliability of the tracing method.

[0007] Therefore, it is urgent to develop a pollution source tracing method for rainwater outlets to solve the problems in the existing technology. Summary of the invention

[0008] The purpose of the present invention is to provide a pollution source tracing method for rainwater outlets, by calculating the determination coefficient between the pollution source and the outlet to compare the correlation between the pollution source and the rainwater outlet, so as to solve the problem of poor tracing accuracy and reliability caused by the current sewage source tracing method proposed in the above background technology not verifying the accuracy of the results.

[0009] In order to solve the above technical problems, the specific technical solutions of the present invention are as follows:

[0010] A method for tracing the pollution source of a rainwater outlet comprises the following steps:

[0011] Obtain pipe network data and first outlet data;

[0012] Determine the monitoring factor based on the first row of data;

[0013] According to the pipe network data and monitoring factors, determine the encrypted monitoring area and obtain the coordinates of the pollution source;

[0014] Obtain pollution source data and second outlet data corresponding to the pollution source coordinates;

[0015] Calculate the pollution load of the pollution source and the pollution load of the outlet according to the pollution source data and the second outlet data;

[0016] According to the pollution load of the pollution source and the pollution load of the outlet, the determination coefficient of the pollution source and the outlet is calculated, and the determination coefficient is used to determine whether the source tracing is completed.

[0017] Furthermore, the pollution load is calculated by the following formula:

[0018] ;

[0019] in, is the pollution load, in t / d; is the discharge flow rate, in m³ / s; is the monitoring factor concentration, in mg / L;

[0020] The determination coefficient is calculated by the following formula:

[0021] ;

[0022] in, is the coefficient of determination, is the pollution load of the pollution source in the i-th time period; is the outlet pollution load in the time period corresponding to the i-th time period of the pollution source; is the average value of the pollution load of the pollution source in n time periods.

[0023] Furthermore, it also includes:

[0024] Before calculating the coefficient of determination of pollution sources and outlets, it also includes calculating relative errors;

[0025] The calculation formula of the relative error is as follows:

[0026] ;

[0027] in, is the relative error, is the pollution load of the pollution source, It is the pollution load at the outlet.

[0028] Further, the step of determining the monitoring factor according to the first row of data comprises the following steps:

[0029] Determine the pollution factor based on the first outlet data;

[0030] When there is only one pollution factor, the monitoring factor is the pollution factor;

[0031] When there are multiple pollution factors, the monitoring factor is determined by the optimization selection index of the pollution factor;

[0032] The optimization selection index is calculated by the following formula:

[0033] ;

[0034] in, is the optimal selection index of the ith pollution factor, is the concentration of the ith pollution factor, is the emission flow of the ith pollution factor, is the river discharge or related normalization factor, is the relevant adjustment coefficient.

[0035] Furthermore, the method of determining the encrypted monitoring area and obtaining the coordinates of the pollution source according to the pipe network data and monitoring factors includes the following steps:

[0036] Obtain monitoring factor data of key nodes;

[0037] Calculate the concentration change rate of the monitoring factor;

[0038] Determine the intensified monitoring area based on the concentration change rate of the monitoring factor;

[0039] Get regional monitoring data for the encrypted monitoring area.

[0040] Furthermore, the method of determining the encrypted monitoring area and obtaining the coordinates of the pollution source according to the pipe network data and the monitoring factors also includes the following steps:

[0041] Based on the regional monitoring data, Kriging interpolation analysis is performed on several coordinate positions in the encrypted monitoring area to obtain the predicted concentration of monitoring factors at several coordinate positions;

[0042] Predict concentrations based on monitoring factors and determine the predicted coordinates of pollution sources.

[0043] A pollution source tracing device for a rainwater outlet, comprising:

[0044] A data acquisition module is used to acquire pipe network data and first outlet data;

[0045] A monitoring factor determination module, used to determine the monitoring factor according to the first row of port data;

[0046] The pollution source coordinate determination module is used to determine the encrypted monitoring area and obtain the pollution source coordinates based on the pipe network data and monitoring factors;

[0047] The source tracing completion judgment module is used to obtain the pollution source data and the second outlet data corresponding to the pollution source coordinates, and calculate the pollution load of the pollution source and the pollution load of the outlet based on the pollution source data and the second outlet data, and calculate the determination coefficient of the pollution source and the outlet based on the pollution load of the pollution source and the pollution load of the outlet, and judge whether the source tracing is completed based on the determination coefficient.

[0048] Furthermore, the source tracing completion judgment module is also used to calculate the relative error before calculating the determination coefficient of the pollution source and the outlet.

[0049] A computer program product comprises a computer program, which implements the steps of the method when executed by a processor.

[0050] A computer device comprises a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method.

[0051] The present invention has the following advantages:

[0052] (1) This application calculates the pollution load of pollution sources and rainwater outlets, and measures the degree of fit between the pollution load of pollution sources and the pollution load of rainwater outlets through the coefficient of determination, so as to evaluate the contribution of pollution sources to rainwater outlet pollution, thereby determining the correlation between the pollution load of pollution sources and the pollution load of rainwater outlets, more efficiently screening out possible pollution sources, avoiding misjudgment and omission, and improving the accuracy and reliability of pollution source tracing, thereby providing reliable data support for subsequent illegal discharge evidence and pollution control. At the same time, the calculation of pollution load can also provide data support for subsequent pollution control.

[0053] (2) This application applies the optimization selection index formula to evaluate the relative importance of each pollution factor in the source tracing process. The optimization selection index formula can screen out the most representative characteristic pollution factors by comprehensively considering the concentration of pollution factors, discharge flow, and river flow or normalization factors, and combining relevant adjustment coefficients, thereby improving the pertinence and efficiency of monitoring. Under limited monitoring resources, priority monitoring of pollution factors with high optimization selection indexes can lock in pollution sources more quickly and reduce unnecessary monitoring work.

[0054] (3) This application uses Kriging interpolation to analyze the predicted concentration of monitoring factors at coordinate locations to predict water quality in unmonitored areas, thereby assisting in finding the actual coordinate locations of pollution sources. At the same time, the application of Kriging interpolation can also optimize the layout of encrypted monitoring to more scientifically divide monitoring blocks and improve monitoring coverage and accuracy.

[0055] Other features and advantages of the present invention will be disclosed in detail in the following specific embodiments and drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0056] Figure 1 It is a schematic diagram of the process of the present invention;

[0057] Figure 2 This is a schematic diagram of the concentration prediction of the monitoring factors of the present invention. DETAILED DESCRIPTION

[0058] In order to better understand the purpose, structure and function of the present invention, the present invention is further described in detail below in conjunction with the accompanying drawings.

[0059] A method for tracing the pollution source of a rainwater outlet, such as Figure 1 As shown, the following steps are included:

[0060] S1: Obtain the pipe network data and the first outlet data;

[0061] S2: Determine the monitoring factor based on the first row of data;

[0062] S3: Determine the encrypted monitoring area and obtain the coordinates of the pollution source based on the pipe network data and monitoring factors;

[0063] S4: Obtain pollution source data and second outlet data corresponding to the pollution source coordinates;

[0064] S5: Calculate the pollution load of the pollution source and the pollution load of the outlet according to the pollution source data and the second outlet data;

[0065] S6: Calculate the determination coefficient of the pollution source and the outlet based on the pollution load of the pollution source and the pollution load of the outlet, and determine whether the tracing is completed based on the determination coefficient.

[0066] In this embodiment, the pipe network data in S1 includes a drainage pipe network topology diagram, etc., which is used to understand the overall layout and connection relationship of the pipe network. Specifically, the pipe network data can be obtained by consulting relevant pipe network archives to obtain information such as the distribution of pipes of different diameters, the location of inspection wells, and the connection direction of each pipe section.

[0067] The first outlet data in S1 includes the detection data of sewage at the outlet site and the river flow and outlet discharge at that time.

[0068] In this embodiment, at the rainwater outlet site, a portable water quality multi-parameter tester was used to quickly test the water sample, and the pH value was measured to be 6.5, the conductivity was 1200μS / cm, and the dissolved oxygen was 5mg / L. The water sample was sent to the laboratory, and the chemical analysis methods such as GC-MS or ICP-MS were used to accurately determine the composition and concentration of pollutants. It was found that the water sample contained an ammonia nitrogen concentration of 15mg / L and a total phosphorus concentration of 2mg / L. The outlet discharge when the water sample was taken was 0.5 m³ / s, and the river flow was 100m³ / s.

[0069] Combining the results of on-site rapid testing and laboratory analysis, it was found that the concentrations of ammonia nitrogen and total phosphorus were high. Combined with the characteristics of known pollution sources, it can be preliminarily judged that the pollution source may be domestic sewage or agricultural wastewater.

[0070] Determining the monitoring factor according to the first row of data in S2 includes the following steps:

[0071] Determine the pollution factor based on the first outlet data;

[0072] When there is a pollution factor, the monitoring factor is the pollution factor;

[0073] When there are multiple pollution factors, the monitoring factor is determined by optimizing the selection index of the pollution factors.

[0074] In this embodiment, the pollution factors in the water sample include ammonia nitrogen (NH3-N) and total phosphorus (TP).

[0075] The optimization selection index is obtained by the following formula:

[0076] ;

[0077] in, is the optimal selection index of the ith pollution factor, which is used to indicate the relative importance of the ith pollution factor in the process of tracing the source of the river outlet;

[0078] is the concentration of the ith pollution factor, which is used to reflect the actual content of the pollution factor in the outlet;

[0079] is the discharge flow of the ith pollution factor, which is used to indicate the rate at which the pollution factor enters the river;

[0080] is the river flow or related normalization factor, which is used to compare and normalize the selected index of pollution factors at different outlets or in different time periods;

[0081] is the relevant adjustment factor, which is an optional adjustment factor used to consider toxicity or other factors other than toxicity that may affect the selection of pollution factors, such as the persistence, detectability, and regulatory restrictions of the pollution factors. In practical applications, if these factors have little effect on the comparison results or can be ignored, then Can be set to 1 or omitted.

[0082] Optionally, the relevant adjustment coefficient may also be directly set to 1 or adjusted based on experience.

[0083] In this example, the concentration of ammonia nitrogen (NH3-N) in the water sample is 15 mg / L, and the concentration of total phosphorus (TP) is 2 mg / L. The discharge volume at the outlet is 0.5 m³ / s, and the river flow is 100 m³ / s. Since the difference between ammonia nitrogen pollution and total phosphorus pollution is not large, the relevant adjustment coefficients of the two are the same, both of which are 1.

[0084] The calculation process of the optimal selection index of ammonia nitrogen (NH3-N) is as follows:

[0085] ;

[0086] The optimization selection index calculation process of total phosphorus (TP) is as follows:

[0087] ;

[0088] It can be seen that the optimization selection index of ammonia nitrogen is significantly higher than that of total phosphorus, which means that under the current monitoring conditions, the relative importance of ammonia nitrogen as a pollution factor is higher. Therefore, ammonia nitrogen is used as a monitoring factor to more effectively track the source of pollution.

[0089] In S3, based on the pipe network data and monitoring factors, the encrypted monitoring area is determined and the coordinates of the pollution source are obtained, including the following steps:

[0090] S31: Obtain monitoring factor data of key nodes;

[0091] The key node may be an intersection of a rainwater pipe network, such as an intersection between a main outlet pipe and a branch pipe, and the monitoring factor data may include a concentration of the monitoring factor and a flow rate of the monitoring factor.

[0092] In this embodiment, integrated flow and water quality monitoring equipment is installed in the rainwater pipe network, and ammonia nitrogen is selected as the monitoring factor. The baseline value of ammonia nitrogen under a stable and pollution-free state is measured to be 0.5 mg / L. Through real-time monitoring, the ammonia nitrogen concentration at a key node is obtained to be 15 mg / L, and the flow rate is 0.5 m³ / s.

[0093] S32: Calculate the concentration change rate of the monitoring factor;

[0094] In this embodiment, the concentration change rate of the monitoring factor is calculated by the following formula:

[0095] ;

[0096] in, is the concentration change rate of the monitoring factor at the kth key node, is the water quality parameter value of the kth key node, It is the reference value of this parameter under stable and pollution-free conditions.

[0097] In this embodiment, the concentration change rate of ammonia nitrogen at this location is calculated as follows:

[0098] ;

[0099] It can be seen that the concentration change rate of ammonia nitrogen here reached 29 times.

[0100] The concentration change rate calculation formula is used to calculate the change rate of water quality parameters at key nodes relative to the stable and pollution-free state, so as to quickly determine the source of pollution. This calculation formula can quickly narrow the scope of investigation and locate possible pollution source areas. This has an important guiding role in subsequent encrypted monitoring and pipeline inspections.

[0101] S33: Determine the encrypted monitoring area based on the concentration change rate of the monitoring factor.

[0102] In this example, the concentration change rate of ammonia nitrogen at the front reached 29 times, which is a very significant increase, strongly indicating that the pollution may come from the upstream direction of the key node. This result provides a clear indication for subsequent intensified monitoring, that is, the rainwater pipe network area upstream of the node should be focused on and the rainwater pipe network area upstream of the node should be determined as the intensified monitoring area.

[0103] Optionally, if the determined encrypted monitoring area is small, the pollution source can be identified by manual investigation directly in the encrypted monitoring area.

[0104] S35: Acquire regional monitoring data of the encrypted monitoring area, wherein the regional monitoring data includes the coordinate position, flow data and water quality parameters of each monitoring point in the encrypted monitoring area.

[0105] In this embodiment, the key nodes in the encrypted monitoring area are used as monitoring points, and the integrated flow and water quality monitoring equipment is further arranged to obtain the coordinate position, flow data and water quality parameters of each monitoring point for subsequent spatial interpolation and prediction analysis. Among them, the water quality parameter is the concentration of the monitoring factor.

[0106] S36: Based on the regional monitoring data, Kriging interpolation analysis is performed to obtain the coordinates of the pollution source.

[0107] Specifically, the steps include:

[0108] Based on the regional monitoring data, Kriging interpolation analysis is performed on several coordinate positions in the encrypted monitoring area to obtain the predicted concentration of monitoring factors at several coordinate positions;

[0109] Predict concentrations based on monitoring factors and determine the predicted coordinates of pollution sources.

[0110] In this embodiment, the pollution source coordinates refer to the pollution source predicted coordinates.

[0111] In this embodiment, the regional monitoring data includes the following Table 1, and the monitoring factor prediction concentrations of several coordinate positions after Kriging interpolation analysis are as follows: Figure 2 shown.

[0112] Table 1:

[0113]

[0114] According to Table 1, if you want to predict the ammonia nitrogen concentration at a point with longitude 116.41 and latitude 39.93, the following steps are included:

[0115] According to the Kriging equations, the weight coefficients of each monitoring point are calculated as follows:

[0116] ;

[0117] According to the weight coefficient of each monitoring point, the predicted concentration of the monitoring factor at the coordinate position is calculated as follows:

[0118] .

[0119] The predicted ammonia nitrogen concentration at the coordinate location of longitude 116.41 and latitude 39.93 is 15.45 mg / L.

[0120] The pollution source coordinates are the predicted coordinates of possible pollution sources. When there are one or more pollution sources, the predicted concentration of the monitoring factor at the coordinate position of the pollution source is generally higher than the predicted concentration of the monitoring factor at the coordinate position around the pollution source. Based on this, one or more pollution source coordinates can be obtained.

[0121] Optionally, if the encrypted monitoring area in S33 is large, making it difficult to determine the coordinates of the pollution source in S36, the concentration can be predicted based on the monitoring factors, the encrypted monitoring area can be re-determined, and S35 and S36 can be repeated.

[0122] If the encrypted monitoring area is small, the problematic pipe section can be directly inspected, the pollution source can be identified by manual investigation, and S36 can be skipped.

[0123] After obtaining the coordinates of the pollution source, the problematic pipe section can be inspected and the pollution source can be locked by manual investigation, such as using CCTV, QV, sonar, ground penetrating radar and other equipment to conduct a comprehensive inspection of the pipeline. CCTV equipment can place the camera inside the pipeline and observe the internal conditions of the pipeline through video transmission, such as whether the pipeline is damaged, leaking, or whether there are foreign objects. Sonar equipment can detect whether there are cavities inside the pipeline, and ground penetrating radar can detect the soil conditions around the pipeline to determine whether there are external influences caused by pipeline problems.

[0124] By observing the internal conditions of the pipeline, the relationship between pipeline problems and pollution can be analyzed. For example, if a pipeline is found to have a broken part and there are signs of sewage seepage around the broken part, it can be determined that the pollution source may have entered the pipeline through the broken part, thereby locating the specific location of the pollution source.

[0125] After locking in the pollution source, it is necessary to calibrate the pollution load of the drainage household and the pollution load of the rainwater outlet.

[0126] In this embodiment, the pollution source data corresponding to the pollution source coordinates in S4 is the pollution source data of the actual coordinates of the pollution source obtained according to the pollution source prediction coordinates, and the pollution source data includes the concentration and emission flow of the monitoring factor at the pollution source. The second outlet data includes the concentration and emission flow of the monitoring factor at the outlet.

[0127] The pollution load in S5 is calculated by the following formula:

[0128] ;

[0129] in, is the pollution load at the coordinate position, in t / d; is the discharge flow rate at the coordinate position, in m³ / s; is the concentration of the monitoring factor at that coordinate position, in mg / L.

[0130] The coefficient of determination in S6 is calculated by the following formula:

[0131] ;

[0132] in, is the coefficient of determination, is the pollution load of the pollution source in the i-th time period; is the outlet pollution load in the time period corresponding to the i-th time period of the pollution source; is the average value of the pollution load of the pollution source in n time periods. Wherein, n is determined as needed. Optionally, the time period corresponding to the i-th time period of the pollution source can be obtained by means such as water flow velocity prediction in the pipe network, which is a prior art and will not be described in detail in this application.

[0133] In this embodiment, the pollution load of the rainwater outlet is: 8t / h, 10t / h and 12t / h; the pollution load of the pollution source is: 7.5t / h, 10.2t / h and 11.8t / h.

[0134] The calculation process of the coefficient of determination is as follows:

[0135] ;

[0136] ;

[0137] ;

[0138] .

[0139] In this embodiment, the determination coefficient threshold range is 0.9 to 1. If it is within the threshold range, the traceability is completed.

[0140] In this embodiment, when there are multiple pollution sources, the pollution load of the pollution source is the sum of the pollution loads of all pollution sources when calculating the determination coefficient.

[0141] This application calculates the pollution load of pollution sources and rainwater outlets, and measures the degree of fit between the pollution load of pollution sources and the pollution load of rainwater outlets through the coefficient of determination, so as to evaluate the contribution of pollution sources to rainwater outlet pollution, thereby determining the correlation between the pollution load of pollution sources and the pollution load of rainwater outlets, more efficiently screening out possible pollution sources, avoiding misjudgment and missed judgment, and improving the accuracy and reliability of pollution source tracing, thereby providing reliable data support for subsequent illegal discharge evidence and pollution control. At the same time, the calculation of pollution load can also provide data support for subsequent pollution control.

[0142] A pollution source tracing device for a rainwater outlet, comprising:

[0143] A data acquisition module is used to acquire pipe network data and first outlet data;

[0144] A monitoring factor determination module, used to determine the monitoring factor according to the first row of port data;

[0145] The pollution source coordinate determination module is used to determine the encrypted monitoring area and obtain the pollution source coordinates based on the pipe network data and monitoring factors;

[0146] The source tracing completion judgment module is used to obtain the pollution source data and the second outlet data corresponding to the pollution source coordinates, and calculate the pollution load of the pollution source and the pollution load of the outlet based on the pollution source data and the second outlet data, and calculate the determination coefficient of the pollution source and the outlet based on the pollution load of the pollution source and the pollution load of the outlet, and judge whether the source tracing is completed based on the determination coefficient.

[0147] A computer program product comprises a computer program, which implements the steps of the method when executed by a processor.

[0148] A computer device comprises a memory, a processor and a computer program stored in the memory, wherein the processor executes the computer program to implement the steps of the method.

[0149] Embodiment 2

[0150] The difference between this embodiment and the first embodiment is that before calculating the determination coefficient of the pollution source and the outlet, the relative error is also calculated. If the relative error is less than the threshold error, the determination coefficient of the pollution source and the outlet is calculated.

[0151] The calculation formula of the relative error is as follows:

[0152] ;

[0153] in, is the relative error, Pollution load for pollution sources; It is the pollution load at the outlet.

[0154] The data shows that the drainage flow rate of the pollution source is 2m³ / s. In this embodiment, the relative error calculation is performed based on the pollution load of 10 hours. The ammonia nitrogen concentration of the drainage in the first 5 hours is measured by water quality monitoring equipment to be 10 mg / L, and increases linearly to 20 mg / L in the next 5 hours.

[0155] ;

[0156] .

[0157] At the same time, the rainwater outlet corresponding to the pollution source that had been locked was monitored for 10 consecutive hours using flow and water quality monitoring equipment, and the outlet pollution load for the time period corresponding to the pollution load of the above-mentioned pollution source was calculated. The 10-hour rainwater outlet pollution load was found to be 255 tons, where the calculation formula is the same as above.

[0158] Calculate the relative error based on the pollution load of the pollution source and the pollution load of the outlet ,as follows:

[0159] ;

[0160] The relative error is only 1.96%, indicating that the pollution load at the pollution source is very close to the pollution load at the rainwater outlet, and the traceability result is accurate and reliable. Optionally, the error threshold can be less than 5%.

[0161] In the process of tracing the source, the difference between the pollution load of the pollution source and the pollution load of the rainwater outlet can be quantified by calculating the RE between the two. The smaller the RE value, the smaller the difference between the pollution source and the monitoring value of the rainwater outlet, and the higher the accuracy of the tracing result, which helps to identify and correct possible prediction deviations.

[0162] In this embodiment, there is only one pollution source and one outlet. If there are multiple pollution sources, It is the sum of the pollution loads of all pollution sources corresponding to the outlet, and the same applies to the outlet.

[0163] In the complex and ever-changing urban drainage system, a single indicator may not accurately reflect the true situation of the pollution source. This application also uses relative error The data were verified by using the coefficient of determination R² and the relative error The traceability results are verified from the error aspect and from the correlation aspect by using the coefficient of determination R² to reduce the risk of traceability errors caused by misjudgment of a single indicator and ensure that the traceability results perform well in different dimensions.

[0164] At the same time, the calculation process of relative error is relatively simple, and it can quickly compare pollution loads over a long period of time or multiple time periods. Setting the calculation of relative error before calculating the determination coefficient can help to quickly verify the connection between the pollution source and the outlet, quickly obtain the calculation results, and save computing resources.

[0165] Optionally, the calculation of the relative error may also be performed after the calculation of the determination coefficient or simultaneously with the calculation of the determination coefficient.

[0166] A pollution source tracing device for a rainwater outlet, comprising:

[0167] A data acquisition module is used to acquire pipe network data and first outlet data;

[0168] A monitoring factor determination module, used to determine the monitoring factor according to the first row of port data;

[0169] The pollution source coordinate determination module is used to determine the encrypted monitoring area and obtain the pollution source coordinates based on the pipe network data and monitoring factors;

[0170] The source tracing completion judgment module is used to obtain the pollution source data and the second outlet data corresponding to the pollution source coordinates, and calculate the pollution load of the pollution source and the pollution load of the outlet based on the pollution source data and the second outlet data, and calculate the determination coefficient of the pollution source and the outlet based on the pollution load of the pollution source and the pollution load of the outlet, and judge whether the source tracing is completed based on the determination coefficient.

[0171] Wherein, the source tracing completion judgment module is also used to calculate the relative error before calculating the determination coefficient of the pollution source and the outlet.

[0172] It is to be understood that the present invention is described by some embodiments, and it is known to those skilled in the art that various changes or equivalent substitutions may be made to these features and embodiments without departing from the spirit and scope of the present invention. In addition, under the teachings of the present invention, these features and embodiments may be modified to adapt to specific circumstances and materials without departing from the spirit and scope of the present invention. Therefore, the present invention is not limited by the specific embodiments disclosed herein, and all embodiments falling within the scope of the claims of this application are within the scope of protection of the present invention.

[0173] In addition, it should be understood that although the present specification is described according to implementation modes, not every implementation mode contains only one independent technical solution. This description of the specification is only for the sake of clarity. Those skilled in the art should regard the specification as a whole. The technical solutions in each embodiment may also be appropriately combined to form other implementation modes that can be understood by those skilled in the art.

Claims

1. A method for tracing the pollution source of a rainwater outlet, characterized in that: The steps include: Obtain pipe network data and first outlet data; Determine the monitoring factor based on the first row of data; According to the pipe network data and monitoring factors, determine the encrypted monitoring area and obtain the coordinates of the pollution source; Obtain pollution source data and second outlet data corresponding to the pollution source coordinates; Calculate the pollution load of the pollution source and the pollution load of the outlet according to the pollution source data and the second outlet data; According to the pollution load of the pollution source and the pollution load of the outlet, the determination coefficient of the pollution source and the outlet is calculated, and the determination coefficient is used to determine whether the source tracing is completed.

2. The method for tracing the pollution source of a rainwater outlet according to claim 1, characterized in that: The pollution load of the pollution source and the pollution load of the discharge outlet are calculated by the following formula: ； in, is the pollution load, in t / d; is the discharge flow rate, in m³ / s; is the monitoring factor concentration, in mg / L; The determination coefficient is calculated by the following formula: ； in, is the coefficient of determination, is the pollution load of the pollution source in the i-th time period; is the outlet pollution load in the time period corresponding to the i-th time period of the pollution source; is the average value of the pollution load of the pollution source in n time periods.

3. The method for tracing the pollution source of a rainwater outlet according to claim 2, characterized in that: Also includes: Before calculating the coefficient of determination of pollution sources and outlets, it also includes calculating relative errors; The calculation formula of the relative error is as follows: ； in, is the relative error, is the pollution load of the pollution source, It is the pollution load at the outlet.

4. The method for tracing the pollution source of a rainwater outlet according to claim 1, characterized in that: Determining the monitoring factor according to the first row of data comprises the following steps: Determine the pollution factor based on the first outlet data; When there is only one pollution factor, the monitoring factor is the pollution factor; When there are multiple pollution factors, the monitoring factor is determined by the optimization selection index of the pollution factor; The optimization selection index is calculated by the following formula: ； in, is the optimal selection index of the ith pollution factor, is the concentration of the ith pollution factor, is the emission flow of the ith pollution factor, is the river discharge or related normalization factor, is the relevant adjustment coefficient.

5. The method for tracing the pollution source of a rainwater outlet according to claim 4, characterized in that: The method of determining the encrypted monitoring area and obtaining the coordinates of the pollution source according to the pipe network data and monitoring factors includes the following steps: Obtain monitoring factor data of key nodes; Calculate the concentration change rate of the monitoring factor; Determine the intensified monitoring area based on the concentration change rate of the monitoring factor; Get regional monitoring data for the encrypted monitoring area.

6. The method for tracing the pollution source of a rainwater outlet according to claim 5, characterized in that: The method of determining the encrypted monitoring area and obtaining the coordinates of the pollution source according to the pipe network data and monitoring factors also includes the following steps: Based on the regional monitoring data, Kriging interpolation analysis is performed on several coordinate positions in the encrypted monitoring area to obtain the predicted concentration of monitoring factors at several coordinate positions; Predict concentrations based on monitoring factors and determine the predicted coordinates of pollution sources.

7. A pollution source tracing device for a rainwater outlet, characterized in that: include: A data acquisition module is used to acquire pipe network data and first outlet data; A monitoring factor determination module, used to determine the monitoring factor according to the first row of port data; The pollution source coordinate determination module is used to determine the encrypted monitoring area and obtain the pollution source coordinates based on the pipe network data and monitoring factors; The source tracing completion judgment module is used to obtain the pollution source data and the second outlet data corresponding to the pollution source coordinates, and calculate the pollution load of the pollution source and the pollution load of the outlet based on the pollution source data and the second outlet data, and calculate the determination coefficient of the pollution source and the outlet based on the pollution load of the pollution source and the pollution load of the outlet, and judge whether the source tracing is completed based on the determination coefficient.

8. The pollution source tracing device for a rainwater outlet according to claim 7, characterized in that: The source tracing completion judgment module is also used to calculate the relative error before calculating the determination coefficient of the pollution source and the outlet.

9. A computer program product, comprising a computer program, characterized in that When the computer program is executed by a processor, the steps of the method according to any one of claims 1 to 6 are implemented.

10. A computer device comprising a memory, a processor and a computer program stored in the memory, characterized in that: The processor executes the computer program to implement the steps of the method according to any one of claims 1 to 6.

Images