Water pollution tracing method and device and electronic equipment

Abstract

The invention relates to the technical field of water pollution management, and discloses a water pollution tracing method and device and electronic equipment, and the method comprises the steps: determining at least one problem point location according to a water quality monitoring value time sequence corresponding to each monitoring point location in a selected range region in a target river channel; according to the position of each problem point location and / or the water quality analysis sequence corresponding to each problem point location, determining at least one to-be-checked river reach in the target river channel; and determining a first pollution source corresponding to the target river channel according to the water quality monitoring value time sequence and the flow monitoring value time sequence corresponding to each discharge port in each to-be-checked river reach. According to the embodiment of the invention, through decomposition and correlation analysis of the water quality index monitoring value time sequence, the change rule of the water quality index at the problem point location can be deeply excavated, reliable data support is provided for water pollution tracing, and the water pollution tracing efficiency is improved in combination with segmented troubleshooting capable of narrowing the water pollution tracing range. And a reliable basis is provided for water pollution treatment.

Description

Technical Field

[0001] The present disclosure relates to the technical field of water pollution management, and in particular to a water pollution source tracing method, device and electronic equipment. Background Art

[0002] In the prior art, common water pollution source tracing research methods include two categories: on-site sampling and measurement methods and mathematical model simulation methods. Compared with on-site sampling and measurement methods, mathematical model simulation methods can predict and simulate pollution sources through probability statistics, optimization theory, and mechanism models, and have higher economic benefits. However, the mathematical model simulation method has problems such as high randomness, large amount of calculation, and susceptibility to various factors such as initial values ​​and boundary conditions, resulting in low accuracy of pollution source tracing results. Summary of the invention

[0003] In view of this, the present disclosure proposes a technical solution for a water pollution source tracing method, device and electronic equipment.

[0004] According to one aspect of the present disclosure, a method for tracing the source of water pollution is provided, comprising: determining at least one problem point according to a time series of water quality monitoring values ​​corresponding to each monitoring point within a selected range in a target river channel, wherein the time series of water quality monitoring values ​​corresponding to any one monitoring point includes a time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point; determining at least one river section to be investigated in the target river channel according to the position of each problem point and / or a water quality analysis sequence corresponding to each problem point, wherein a water quality analysis sequence corresponding to any one problem point is determined by time series decomposition based on a time series of water quality monitoring values ​​and a time series of flow rate monitoring values ​​corresponding to the problem point, and the water quality analysis sequence corresponding to the problem point includes a trend analysis sequence, a periodic analysis sequence and a residual analysis sequence corresponding to the problem point; determining a first pollution source corresponding to the target river channel according to a time series of water quality monitoring values ​​and a time series of flow rate monitoring values ​​corresponding to each outlet in each river section to be investigated, wherein the time series of water quality monitoring values ​​corresponding to any one outlet includes a time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet.

[0005] In a possible implementation, at least one problem point is determined based on the time series of water quality monitoring values ​​corresponding to each monitoring point in the selected range area in the target river channel, including: during a rainy period, for any monitoring point, if the monitoring value at any moment in the time series of monitoring values ​​corresponding to any water quality indicator at the monitoring point exceeds the first preset threshold value corresponding to the water quality indicator, the monitoring point is determined as a problem point; or, during a dry period, for any monitoring point, if the duration for which the monitoring value at a single moment in the time series of monitoring values ​​corresponding to any water quality indicator at the monitoring point exceeds the second preset threshold value corresponding to the water quality indicator meets the first preset duration, and / or the duration for which the time series of monitoring values ​​corresponding to any water quality indicator conforms to the sewage change law meets the second preset duration, the monitoring point is determined as a problem point, wherein the sewage change law is used to reflect the daily cycle law of the change over time of the monitoring value of any water quality indicator that exceeds the pollution standard.

[0006] In a possible implementation, the method further includes: for any problem point, based on the time series of flow velocity monitoring values ​​corresponding to the problem point and the distance between the problem point and the first upstream problem point, determining the time difference between the problem point and the first upstream problem point; for any problem point, based on the time difference between the problem point and the first upstream problem point, aligning the time series of water quality monitoring values ​​corresponding to the problem point with the time series of water quality monitoring values ​​corresponding to the first upstream problem point to determine the aligned time series of water quality monitoring values ​​corresponding to the problem point; for any problem point, performing time series decomposition on the aligned time series of water quality monitoring values ​​corresponding to the problem point to determine the water quality analysis sequence corresponding to the problem point.

[0007] In a possible implementation, at least one river section to be checked in the target river channel is determined based on the position of each problem point and / or the water quality analysis sequence corresponding to each problem point, including: for any problem point, if the monitoring point upstream of the problem point is not a problem point, the river section between the problem point and the upstream adjacent monitoring point is determined as the river section to be checked; or, for any two adjacent problem points, a multi-dimensional correlation analysis is performed based on the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet a first preset condition, and if the correlation analysis results between the two problem points meet the first preset condition, the river section between the two problem points is determined as the river section to be checked, wherein the any two adjacent problem points represent two monitoring points that are adjacent in position and are both problem points.

[0008] In a possible implementation, for any two adjacent problem points, a multi-dimensional correlation analysis is performed according to the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet the first preset condition, including: for any two adjacent problem points, according to the trend analysis sequences corresponding to the two problem points, determining the trend similarity between the two problem points; for any two adjacent problem points, according to the residual analysis sequences corresponding to the two problem points, determining the residual mutation similarity between the two problem points; for any two adjacent problem points, according to the periodic analysis sequences corresponding to the two problem points, determining the similarity of the change rules between the two problem points and the mutual correlation function value between the two problem points. ; For any two adjacent problem points, when the trend similarity corresponding to the two problem points meets the second preset condition, the change law similarity meets the third preset condition, the mutual correlation function value meets the fourth preset condition, and the residual mutation similarity meets the fifth preset condition, the correlation analysis result between the two problem points is determined to meet the first preset condition; For any two adjacent problem points, when the trend similarity corresponding to the two problem points does not meet the second preset condition, or the change law similarity does not meet the third preset condition, or the mutual correlation function value does not meet the fourth preset condition, or the residual mutation similarity does not meet the fifth preset condition, any one of the following is true, the correlation analysis result between the two problem points is determined to not meet the first preset condition.

[0009] In a possible implementation, the first pollution source corresponding to the target river channel is determined according to the time series of water quality monitoring values ​​and the time series of flow monitoring values ​​corresponding to each outlet in each river section to be investigated, including: for any outlet, when the monitoring value at any moment in the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, the pollution flux corresponding to the water quality indicator at the outlet is determined according to the time series of monitoring values ​​and the time series of flow monitoring values ​​corresponding to the water quality indicator at the outlet; for any outlet, when the monitoring value at any moment in the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is greater than or equal to the third preset threshold value corresponding to the water quality indicator, When the first pollution flux threshold is reached, the outlet is determined as the first pollution source corresponding to the target river; for any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is less than or equal to the first pollution flux threshold and greater than or equal to the second pollution flux threshold, determine whether the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition based on the water quality monitoring value time series corresponding to the outlet, and when the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition, determine the outlet as the first pollution source corresponding to the target river.

[0010] In a possible implementation, the method further includes: when the first pollution source corresponding to the target river channel is a rainwater outlet, determining the rainwater pipe network system corresponding to the rainwater outlet; determining the liquid level distribution information, flow velocity distribution information and water quality distribution information of the rainwater pipe network system during the dry period according to the real-time monitoring data corresponding to the rainwater pipe network system, wherein the real-time monitoring data corresponding to the rainwater pipe network system includes a liquid level monitoring value time series, a flow velocity monitoring value time series and a monitoring value time series corresponding to at least one water quality indicator at each monitoring position in the rainwater pipe network system; determining the abnormal liquid level position of the rainwater pipe network system according to the liquid level distribution information; determining the abnormal flow velocity position of the rainwater pipe network system according to the flow velocity distribution information; determining the abnormal water quality position of the rainwater pipe network system according to the water quality distribution information; determining the second pollution source corresponding to the target river channel according to the abnormal liquid level position, abnormal flow velocity position and abnormal water quality position of the rainwater pipe network system.

[0011] In a possible implementation, the method further includes: when the first pollution source corresponding to the target river channel is an overflow outlet, determining the sewage pipe network system corresponding to the overflow outlet; determining the rainwater increment during the rainy day period corresponding to each monitoring location according to the flow monitoring value time series corresponding to each monitoring location in the sewage pipe network system; for any monitoring location, when the rainwater increment during the rainy day period corresponding to the monitoring location is greater than or equal to a preset increment threshold, determining the monitoring location as the third pollution source corresponding to the target river channel.

[0012] According to another aspect of the present disclosure, a water pollution source tracing device is provided, comprising: a problem point determination module, for determining at least one problem point according to a time series of water quality monitoring values ​​corresponding to each monitoring point in a selected range area in a target river channel, wherein the time series of water quality monitoring values ​​corresponding to any one monitoring point includes a time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point; a river section investigation module, for determining at least one river section to be investigated in the target river channel according to the position of each problem point and / or the water quality analysis sequence corresponding to each problem point, wherein the time series of water quality monitoring values ​​corresponding to any one of the problem points includes a time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point; The water quality analysis sequence is determined by decomposing the time series of the water quality monitoring values ​​and the flow rate monitoring values ​​corresponding to the problem point, and the water quality analysis sequence corresponding to the problem point includes the trend analysis sequence, periodic analysis sequence and residual analysis sequence corresponding to the problem point; the first pollution source determination module is used to determine the first pollution source corresponding to the target river channel according to the water quality monitoring value time series and the flow rate monitoring value time series corresponding to each outlet in each river section to be investigated, wherein the water quality monitoring value time series corresponding to any outlet includes the monitoring value time series corresponding to at least one water quality indicator at the outlet.

[0013] According to another aspect of the present disclosure, an electronic device is provided, comprising: a processor; and a memory for storing processor-executable instructions; wherein the processor is configured to implement the above method when executing the instructions stored in the memory.

[0014] In the disclosed embodiments, a variety of monitoring indicators such as water quality indicators, flow velocity and flow rate can be flexibly selected for real-time monitoring in different monitoring scenarios such as monitoring points and outlets in the target river. Based on the monitoring value time series corresponding to at least one water quality indicator at each monitoring point in the selected range area in the target river, real-time monitoring of water quality abnormalities can be achieved and at least one problem point can be determined. According to the location of each problem point and / or the water quality analysis sequence corresponding to each problem point, the target river channel can be checked in sections from two aspects: the location distribution of the problem points and the water flow transmissibility, and at least one river section to be checked in the target river channel can be determined, thereby effectively narrowing the scope of water pollution source tracing and improving the efficiency and accuracy of water pollution source tracing; wherein, the water quality analysis sequence corresponding to any problem point is determined by time series decomposition based on the time series of water quality monitoring values ​​and the time series of flow rate monitoring values ​​corresponding to the problem point. Compared with the method of chemical determination based on on-site sampling or the method of simulation and prediction using mathematical models in the prior art, the present disclosure uses real-time monitoring data presenting a time series to perform time series decomposition, which can deeply explore the long-term changes and periodic fluctuations of at least one water quality indicator at the problem point, thereby providing reliable data support for water pollution source tracing and improving the accuracy and reliability of water pollution source tracing. According to the time series of water quality monitoring values ​​and flow monitoring values ​​corresponding to each outlet in each river section to be investigated, the first pollution source corresponding to the target river can be determined, wherein the time series of water quality monitoring values ​​corresponding to any outlet includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet. Compared with the mathematical model simulation method, the water pollution source tracing method provided by the present invention can make full use of real-time monitoring data, has high accuracy and reliability, and reduces the workload of data processing. In addition, the use of at least one water quality indicator can also dynamically trace the source of water pollution based on the changes in the monitoring values ​​of different water quality indicators, determine the pollution source, and further improve the reliability and accuracy of the tracing results.

[0015] Further features and aspects of the present disclosure will become apparent from the following detailed description of exemplary embodiments with reference to the attached drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0016] The accompanying drawings, which are incorporated in and constitute a part of the specification, illustrate exemplary embodiments, features, and aspects of the disclosure and, together with the description, serve to explain the principles of the disclosure.

[0017] Figure 1 A flow chart of a method for tracing the source of water pollution according to an embodiment of the present disclosure is shown.

[0018] Figure 2 A schematic diagram of a target river channel according to an embodiment of the present disclosure is shown.

[0019] Figure 3A schematic diagram showing time series alignment according to an embodiment of the present disclosure is shown.

[0020] Figure 4 A schematic diagram showing a time series decomposition result according to an embodiment of the present disclosure.

[0021] Figure 5 A schematic diagram of a river section to be investigated according to an embodiment of the present disclosure is shown.

[0022] Figure 6 A block diagram of a water pollution source tracing device according to an embodiment of the present disclosure is shown.

[0023] Figure 7 A block diagram of an electronic device according to an embodiment of the present disclosure is shown. DETAILED DESCRIPTION

[0024] Various exemplary embodiments, features and aspects of the present disclosure will be described in detail below with reference to the accompanying drawings. The same reference numerals in the accompanying drawings represent elements with the same or similar functions. Although various aspects of the embodiments are shown in the accompanying drawings, the drawings are not necessarily drawn to scale unless otherwise specified.

[0025] The word “exemplary” is used exclusively herein to mean “serving as an example, example, or illustration.” Any embodiment described herein as “exemplary” is not necessarily to be construed as preferred or advantageous over other embodiments.

[0026] The term "and / or" in this article is only a description of the association relationship of associated objects, indicating that there may be three relationships. For example, A and / or B can represent the three situations of: A exists alone, A and B exist at the same time, and B exists alone. In addition, the term "at least one" in this article represents any combination of at least two of any one or more of a plurality of. For example, including at least one of A, B, and C can represent any one or more elements selected from the set consisting of A, B, and C. The terms "two adjacent problem points" and "problem point 1 is adjacent to problem point 2" in this article both mean adjacent locations and are two monitoring points of the problem point.

[0027] In addition, in order to better illustrate the present disclosure, numerous specific details are given in the following specific embodiments. It should be understood by those skilled in the art that the present disclosure can also be implemented without certain specific details. In some examples, methods, means, components and circuits well known to those skilled in the art are not described in detail in order to highlight the subject matter of the present disclosure.

[0028] In the prior art, common water pollution source tracing research methods include: on-site sampling and measurement methods and mathematical model simulation methods. Among them, the mathematical model simulation method can include methods based on probability statistics, methods based on optimization theory, and methods based on mechanism models.

[0029] Specifically, methods based on probability statistics are usually based on Bayesian theorem, focusing on evaluating the probability of occurrence of specific events, and then determining the source of pollution through the probability evaluation results; methods based on optimization theory can usually use certain mathematical and physical equations to analyze the movement of pollutants, optimize parameters to make the simulated values ​​and observed values ​​the best match, and use simulated values ​​to determine the source of pollution; methods based on mechanism models are usually based on existing complex water environment models, combined with relevant restrictions or algorithms to trace the source of pollutants.

[0030] Compared with the on-site sampling measurement method, the above-mentioned various mathematical model simulation methods have higher economic benefits, but they also have problems such as high randomness, large amount of calculation, and susceptibility to various factors such as initial values ​​and boundary conditions, resulting in low accuracy of pollution tracing results.

[0031] In view of this, the present disclosure provides a method for tracing the source of water pollution, which can flexibly select multiple monitoring indicators such as water quality indicators, flow rate and flow rate for real-time monitoring in different monitoring scenarios, and deeply explore the long-term and periodic change laws of water quality indicators at problem points through the decomposition and correlation analysis of the time series of water quality indicator monitoring values, which can provide reliable data support for water pollution tracing, thereby improving the accuracy of water pollution tracing; on the other hand, using problem points to conduct segmented inspections of target rivers can effectively narrow the scope of water pollution tracing, thereby improving the efficiency of water pollution tracing and providing a reliable basis for water pollution control. In addition, the present disclosure uses at least one water quality indicator to dynamically determine the pollution source based on the changes in the monitoring values ​​of different water quality indicators, further improving the reliability and accuracy of water pollution tracing results.

[0032] The water pollution source tracing method disclosed in the present invention is described in detail below.

[0033] Figure 1 A flow chart of a method for tracing the source of water pollution according to an embodiment of the present disclosure is shown. The method for tracing the source of water pollution can be executed by an electronic device such as a terminal device or a server. The terminal device can be a user equipment (UE), a mobile device, a user terminal, a terminal, a cellular phone, a cordless phone, a personal digital assistant (PDA), a handheld device, a computing device, a vehicle-mounted device, a wearable device, etc. The method for tracing the source of water pollution can be implemented by a processor calling computer-readable instructions stored in a memory. Alternatively, the method for tracing the source of water pollution can be executed by a server. Figure 1 As shown, the water pollution source tracing method includes:

[0034] In step S11, at least one problem point is determined based on the time series of water quality monitoring values ​​corresponding to each monitoring point in the selected range area in the target river channel, wherein the time series of water quality monitoring values ​​corresponding to any monitoring point includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point.

[0035] The target river channel here can refer to a river channel whose water body is affected by at least one drainage pipe network system, and its specific form can be flexibly set according to actual use requirements, for example, it can be any river channel flowing through an urban built-up area or an industrial park, etc., and the present disclosure does not make specific limitations on this. The specific form of the drainage pipe network system here, for example, can be a rainwater pipe network system, a sewage pipe network system, and a combined pipe network system, etc., which can be flexibly set according to actual use requirements, and the present disclosure does not make specific limitations on this.

[0036] Based on the preset position distribution distance, multiple monitoring points can be set in the target river channel, and water quality monitoring equipment can be deployed at each monitoring point to monitor the water quality of each monitoring point in real time, and determine the time series of water quality monitoring values ​​corresponding to each monitoring point, so as to divide the target river channel into multiple river sections. Among them, the specific value of the preset position distribution distance can be flexibly set according to actual use requirements, for example, it can be set to 5-10 kilometers, etc., and this disclosure does not make specific restrictions on this; the specific form of the water quality monitoring equipment can be flexibly set according to actual use requirements, and this disclosure does not make specific restrictions on this.

[0037] In one example, the water quality monitoring equipment may be a spectral sensor. The spectral sensor may obtain the monitoring value corresponding to at least one water quality indicator at the corresponding monitoring point based on a preset acquisition frequency. Compared with the method commonly used in the prior art of conducting water quality analysis in the laboratory after manual sampling, it may realize in-situ, online, high-frequency, and real-time water quality monitoring and data acquisition. Among them, the specific value of the acquisition frequency of the spectral sensor, for example, may be set to 3-60 minutes / time, 5-30 minutes / time, 8-20 minutes / time, 10-15 minutes / time, 3-5 minutes / time, etc., which is much higher than the traditional detection method, and thus the monitoring value of the water quality indicator may be obtained at a higher frequency, thereby effectively capturing the water quality characteristics of the target river. Preferably, the acquisition frequency of the spectral sensor is 3-5 minutes / time. Preferably, the spectral sensor is a quantum dot spectral sensor.

[0038] The selected range area and the specific value of the selected range here can be flexibly set according to actual usage requirements. For example, it can be set to a section of the target river corresponding to an industrial area, or a section of the target river corresponding to a living area, etc. The present disclosure does not make specific limitations on this.

[0039] The time series of water quality monitoring values ​​corresponding to any monitoring point may include the time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point; the time series of monitoring values ​​corresponding to any water quality indicator may include the monitoring value corresponding to the water quality indicator at the monitoring point at each collection moment.

[0040] The specific types of water quality indicators here can be flexibly set according to actual usage requirements. For example, monitoring items required to be controlled in national or local surface water environmental quality standards can be selected, or they can be set based on common pollutants from urban and rural non-point sources in the target river. The present disclosure does not make specific limitations on this.

[0041] In one possible implementation, the water quality indicators of each monitoring point in the target river may include: pH, turbidity, dissolved oxygen, permanganate index, total organic carbon, five-day biochemical oxygen demand, chemical oxygen demand, water temperature, pH value, total phosphorus, total nitrogen, ammonia nitrogen, organic phosphorus, organic nitrogen, nitrate nitrogen, nitrate nitrogen, phosphate, conductivity, copper, zinc, selenium, arsenic, mercury, cadmium, chromium (hexavalent), cyanide, fluoride, sulfide, volatile phenol, petroleum, anionic surfactant, and at least one of fecal coliform group.

[0042] According to the time series of water quality monitoring values ​​corresponding to each monitoring point in the target river, at least one problem point can be determined. Among them, the problem point can represent the monitoring point where water quality abnormality exists. The specific content of the water quality abnormality here can be flexibly set according to the actual use requirements. For example, it can include the presence of pollutants that are not common in the target river, or the monitoring value of any water quality indicator exceeds the threshold of the water quality indicator in the target river, or the monitoring value time series corresponding to any water quality indicator meets the preset change law, etc. It can be flexibly set according to the actual use requirements, and the present disclosure does not make specific limitations on this.

[0043] The time series of water quality monitoring values ​​corresponding to any monitoring point includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point. Therefore, this solution can simultaneously determine the time series of monitoring values ​​of one or more water quality indicators. If the time series of monitoring values ​​of any water quality indicator exceeds the threshold of the water quality indicator in the target river or the time series of monitoring values ​​corresponding to any water quality indicator meets the preset change law, etc., it means that the monitoring point is a problem point.

[0044] The process of determining at least one problem point according to the time series of water quality monitoring values ​​corresponding to each monitoring point in the target river channel will be described in detail later in combination with possible implementation methods of the present disclosure, which will not be elaborated here.

[0045] In step S12, at least one river section to be inspected in the target river channel is determined based on the location of each problem point and / or the water quality analysis sequence corresponding to each problem point, wherein the water quality analysis sequence corresponding to any problem point is determined by time series decomposition of the water quality monitoring value time series and the flow velocity monitoring value time series corresponding to the problem point, and the water quality analysis sequence corresponding to the problem point includes a trend analysis sequence, a periodic analysis sequence and a residual analysis sequence corresponding to the problem point.

[0046] For any problem point, in order to deeply explore the changing rules of each water quality index at the problem point, so as to facilitate the subsequent analysis of the correlation between different problem points, each problem point can also be provided with a flow rate monitoring device to monitor the flow rate at each problem point in real time, and determine the time series of flow rate monitoring values ​​corresponding to the problem point; then, the time series can be decomposed according to the time series of flow rate monitoring values ​​and water quality monitoring values ​​corresponding to the problem point, and the long-term trend, periodic fluctuations and irregular parts that cannot be explained by long-term trend and periodic fluctuations of at least one water quality index at the problem point can be comprehensively analyzed from the two aspects of long-term trend and periodic fluctuation, so as to determine the water quality analysis sequence corresponding to the problem point.

[0047] The specific form of the flow rate monitoring device can refer to the implementation in the related art, and the present disclosure does not make specific limitations on this. It should be noted that at the same problem point, the water quality monitoring device and the flow rate monitoring device need to be arranged at the same location.

[0048] The specific method of time series decomposition can refer to the implementation methods in the relevant technology, for example, the sliding average method or the time series analysis method based on the internal and external double cycles can be used, and the present disclosure does not make specific limitations on this. The water quality analysis sequence of any problem point can include the trend analysis sequence, periodic analysis sequence and residual analysis sequence corresponding to the problem point.

[0049] In one example, a time series decomposition method can be used to decompose the time series of water quality monitoring values ​​at any problem point into a trend analysis sequence, a periodic analysis sequence, and a residual analysis sequence according to the time series of flow rate monitoring values ​​corresponding to any problem point. The time series decomposition method can be expressed as formula (1):

[0050]

[0051] Among them, PSC i,j,k It represents the monitoring value of the i-th water quality index at the k-th monitoring point in the selected range of the target river at the collection time j; TPSC i,j,kIt represents the trend component of the monitoring value of the i-th water quality index at the k-th monitoring point in the selected range of the target river at the collection time j; SPSC i,j,k Represents the periodic component of the monitoring value of the i-th water quality index at the k-th monitoring point in the selected range of the target river at the collection time j; RPSC i,j,k It represents the residual component of the monitoring value of the i-th water quality index at the k-th monitoring point in the selected range of the target river at the collection time j.

[0052] The following text will combine the possible implementation methods of the present invention to explain in detail the process of performing time series decomposition for any problem point based on the time series of water quality monitoring values ​​and flow rate monitoring values ​​corresponding to the problem point, and determining the water quality analysis sequence corresponding to the problem point, which will not be elaborated here.

[0053] Figure 2 FIG. 1 is a schematic diagram showing a target river channel according to an embodiment of the present disclosure. Figure 2 As shown in the figure, the black solid arrow indicates the target river channel, and the circular marks indicate the six monitoring points set in the target river channel, dividing the target river channel into six river sections. Figure 2 As shown in (a), the three light gray circular marks in the upstream represent three monitoring points that are not problem points, and the three dark gray circular marks in the downstream represent problem points, namely problem point 1, problem point 2, and problem point 3. The river section between the first problem point 1 in the upstream and its adjacent monitoring point in the upstream is problem river section 1, the river section between problem point 1 and problem point 2 is problem river section 2, and the river section between problem point 2 and problem point 3 is problem river section 3. Among them, upstream refers to the direction of the water source corresponding to the target river channel, and downstream refers to the direction where the water flow finally flows after transmission.

[0054] There may be upstream and downstream transfer of sewage between adjacent problematic river sections. Therefore, the factors that lead to water pollution exceeding the standard in the problematic river section located downstream include at least the transfer of sewage from the upstream problematic river section and the existence of independent external sewage input in the problematic river section.

[0055] For any problematic river section, when it is determined that the problematic river section is only affected by the sewage transmission from the upstream problematic river section, it can be judged that there is no pollution source within the scope of the problematic river section, and thus there is no need to conduct further pollution tracing for the problematic river section.

[0056] Based on the above principles, according to the location of each problem point in the target river channel, or the water quality analysis sequence corresponding to each problem point, the problem river sections can be preliminarily screened from the two aspects of the location distribution of the problem points and the sewage transmissibility, so as to determine the river sections to be investigated that may have separate external sewage inputs, thereby reducing the number of river sections that need to be investigated in the process of water pollution source tracing, narrowing the scope of pollution source tracing, and thus improving the efficiency and accuracy of water pollution source tracing.

[0057] The following text will describe in detail the process of determining at least one river section to be investigated in the target river channel according to the location of each problem point or the water quality analysis sequence corresponding to each problem point in combination with possible implementation methods of the present disclosure, which will not be elaborated here.

[0058] In step S13, the first pollution source corresponding to the target river is determined based on the water quality monitoring value time series and flow monitoring value time series corresponding to each outlet in each river section to be inspected, wherein the water quality monitoring value time series corresponding to any outlet includes the monitoring value time series corresponding to at least one water quality indicator at the outlet.

[0059] For any river section to be investigated, a comprehensive analysis can be conducted from the two aspects of water quality change and flow rate change based on the time series of water quality monitoring values ​​and flow rate monitoring values ​​corresponding to each outlet in the river section to be investigated, to determine at least one problematic outlet with abnormal water quality, and to determine the problematic outlet as the first pollution source directly affecting the water quality of the target river. The outlet here can refer to the outlet corresponding to any drainage pipe network system that affects the target river, and can include rainwater outlets and overflow outlets. The time series of water quality monitoring values ​​corresponding to any outlet can include the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet; the specific form of the time series of monitoring values ​​corresponding to each water quality indicator can refer to the time series of water quality monitoring values ​​corresponding to any monitoring point.

[0060] The specific method for determining the time series of water quality monitoring values ​​and flow rate monitoring values ​​corresponding to any outlet can be flexibly set according to actual use requirements, and can refer to the method for determining the time series of water quality monitoring values ​​and flow rate monitoring values ​​corresponding to any monitoring point, and the present disclosure does not make specific limitations on this. The water quality indicators at the outlet may include the following indicators: chemical oxygen demand, conductivity, ammonia nitrogen, five-day biochemical oxygen demand, total phosphorus, total nitrogen, suspended solids, total dissolved solids, petroleum, PH, anionic surfactants, cyanide, sulfide, fluoride, chloride, organic phosphorus, sulfate, mercury, chromium, cadmium, arsenic, lead, nickel, beryllium, silver, selenium, copper, zinc, manganese, iron, volatile phenols, benzene series, aniline, nitrobenzene.

[0061] In one possible implementation, for any outlet, water quality monitoring equipment and flow monitoring equipment can be deployed at the last inspection well upstream of the pipeline of the outlet to obtain the time series of water quality monitoring values ​​and flow monitoring values ​​corresponding to the outlet.

[0062] Optionally, according to the real-time progress of pollution tracing, the monitoring equipment installed at each monitoring point where it is determined that there is no water pollution exceeding the standard can be removed and installed at each outlet of the river section to be investigated. This can be done by refining the equipment step by step and deploying it each time it is needed. This can improve the use efficiency of the monitoring equipment, reduce the investment in monitoring equipment while ensuring the accuracy of the monitoring data, save the equipment investment cost in the process of water pollution tracing, and improve the economy.

[0063] The following text will describe in detail the process of determining the first pollution source corresponding to the target river channel based on the time series of water quality monitoring values ​​and flow monitoring values ​​corresponding to each outlet in each river section to be investigated, in combination with possible implementation methods of the present disclosure, which will not be elaborated here.

[0064] In the disclosed embodiments, a variety of monitoring indicators such as water quality indicators, flow velocity and flow rate can be flexibly selected for real-time monitoring in different monitoring scenarios such as monitoring points and outlets in the target river. Based on the monitoring value time series corresponding to at least one water quality indicator at each monitoring point in the selected range area in the target river, real-time monitoring of water quality abnormalities can be achieved and at least one problem point can be determined. According to the location of each problem point, or the water quality analysis sequence corresponding to each problem point, the target river channel can be checked in sections from two aspects: the location distribution of the problem points and the sewage transmissibility, and at least one river section to be checked in the target river channel can be determined, thereby effectively narrowing the scope of water pollution source tracing and improving the efficiency and accuracy of water pollution source tracing; wherein, the water quality analysis sequence corresponding to any problem point is determined by time series decomposition based on the time series of water quality monitoring values ​​and the time series of flow rate monitoring values ​​corresponding to the problem point. Compared with the method of chemical determination based on on-site sampling or the method of simulation and prediction using mathematical models in the prior art, the present disclosure uses real-time monitoring data presenting a time series to perform time series decomposition, which can deeply explore the long-term changes and periodic fluctuations of at least one water quality indicator at the problem point, thereby providing reliable data support for water pollution source tracing and improving the accuracy and reliability of water pollution source tracing. According to the time series of water quality monitoring values ​​and flow monitoring values ​​corresponding to each outlet in each river section to be investigated, the first pollution source corresponding to the target river can be determined, wherein the time series of water quality monitoring values ​​corresponding to any outlet includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet. Compared with the mathematical model simulation method, the water pollution source tracing method provided by the present invention can make full use of real-time monitoring data, has high accuracy and reliability, and reduces the workload of data processing. In addition, the use of at least one water quality indicator can also dynamically trace the source of water pollution based on the changes in the monitoring values ​​of different water quality indicators, determine the pollution source, and further improve the reliability and accuracy of the tracing results.

[0065] In a possible implementation, at least one problem point is determined based on a time series of water quality monitoring values ​​corresponding to each monitoring point within a selected range in the target river channel, including: during a rainy period, for any monitoring point, if the monitoring value at any moment in the time series of monitoring values ​​corresponding to any water quality indicator at the monitoring point exceeds a first preset threshold value corresponding to the water quality indicator, the monitoring point is determined as a problem point; or, during a dry period, for any monitoring point, if the duration for which the monitoring value at a single moment in the time series of monitoring values ​​corresponding to any water quality indicator at the monitoring point exceeds a second preset threshold value corresponding to the water quality indicator meets the first preset duration, and / or the duration for which the time series of monitoring values ​​corresponding to any water quality indicator complies with a sewage change law meets a second preset duration, the monitoring point is determined as a problem point, wherein the sewage change law is used to reflect the daily cycle law of the change over time of the monitoring value of any water quality indicator that exceeds the pollution standard.

[0066] The dry weather period generally refers to the period when the drainage pipe network system corresponding to the target river channel is not affected by rainfall runoff, and is generally a period of continuous rainfall-free time from the end of rainfall exceeding the first preset time range to the beginning of the next rainfall. The first preset time range here can be flexibly set according to actual usage requirements, for example, 48 hours, 72 hours, etc., and the present disclosure does not make specific limitations on this.

[0067] In contrast, the rainy day period can represent the period when the drainage pipe network system corresponding to the target river channel is affected by rainfall runoff. Usually, when the average hourly rainfall needs to meet the preset rainfall condition, the time period from the beginning of rainfall to the second preset time range after the end of rainfall is defined as the rainy day period. Among them, the specific value of the preset rainfall condition can be flexibly set according to actual use requirements, and this disclosure does not make specific restrictions on this; the specific value of the second preset time range can be flexibly set according to actual use requirements, and this disclosure does not make specific restrictions on this;

[0068] In addition, if the interval between the end of a certain rainfall and the beginning of the next rainfall is less than the third preset time range, the two rainfalls are merged into one rainy day period. The third preset time range here can be flexibly set according to actual usage requirements, such as 3 hours, 6 hours, etc., and the present disclosure does not make specific limitations on this.

[0069] During rainy days, in the absence of external sewage input, the monitoring values ​​corresponding to each water quality indicator in the target river channel will usually decrease due to the dilution effect caused by the increase in water volume in the target river channel.

[0070] Therefore, during rainy days, for any monitoring point, if the monitoring value at any moment in the monitoring value time series corresponding to any water quality indicator at the monitoring point exceeds the first preset threshold value corresponding to the water quality indicator, it indicates that there is external sewage input upstream of the monitoring point, and the monitoring point can be determined as a problem point. Among them, the specific value of the first preset threshold value corresponding to any water quality indicator can be flexibly set according to the actual use requirements of different water quality indicators. For example, it can be set based on the "Surface Water Environmental Quality Standards", etc., and the present disclosure does not make specific limitations on this. In the present disclosure, since there are many kinds of water quality indicators (such as the various water quality indicators described above), the problem point can be dynamically determined based on the situation where different water quality indicators exceed the corresponding first preset threshold value, thereby improving the real-time dynamic tracing effect.

[0071] During dry periods, in the absence of external sewage input, the monitoring values ​​corresponding to each water quality indicator in the target river channel usually remain relatively stable.

[0072] Therefore, during the dry period, for any monitoring point, when the monitoring value at a single moment in the monitoring value time series corresponding to any water quality indicator at the monitoring point exceeds the second preset threshold corresponding to the water quality indicator for a period that satisfies the first preset period, it indicates that external sewage has been input upstream of the monitoring point, and the monitoring point can be identified as a problem point. Among them, the specific value of the second preset threshold corresponding to any water quality indicator can be flexibly set according to the actual use requirements of different water quality indicators, for example, it can be set based on the "Surface Water Environmental Quality Standards", etc., and the present disclosure does not make specific restrictions on this; the specific value of the first preset period can be flexibly set according to the actual use requirements, for example, it can be set to 1 hour, etc., and the present disclosure does not make specific restrictions on this.

[0073] On the other hand, for any monitoring point, if the duration of the monitoring value time series corresponding to any water quality index at the monitoring point exceeds the preset threshold value corresponding to the water quality index and does not meet the first preset duration, if the duration of the monitoring value time series corresponding to the water quality index conforming to the sewage change law meets the second preset duration, it can also be explained that there is an input of external sewage upstream of the monitoring point, but the external sewage may not be continuously and stably input, but intermittently input according to a specific periodic law, so that the monitoring point will not have a long-term water quality index exceeding the standard. In this case, the monitoring point can also be determined as a problem point, so as to improve the accuracy and reliability of water quality abnormality monitoring through a preliminary analysis of the changes in the monitoring values ​​of the water quality index, and avoid the unreasonable setting of time parameters such as monitoring time and acquisition frequency, which affects the accuracy of troubleshooting the problem point. Among them, the sewage change law can reflect the daily cycle law of the monitoring value of any water quality index exceeding the standard in the target river channel over time, and its specific content can be flexibly set according to actual use needs, and the present disclosure does not make specific restrictions on this.

[0074] In one example, the sewage change pattern includes: the monitoring value of any water quality indicator during the peak water usage period during the day is higher than the preset threshold, and the monitoring value during the low water usage period at night is lower than the preset threshold.

[0075] Through the above process, for different weather scenarios, we can conduct a comprehensive and reliable analysis of the water quality monitoring time series corresponding to each monitoring point from two aspects: the monitoring values ​​of water quality indicators at each moment and the changes in water quality indicators, so as to determine the problem points where abnormal water quality occurs.

[0076] In a possible implementation, the method further includes: for any problem point, based on the time series of flow velocity monitoring values ​​corresponding to the problem point and the distance between the problem point and the first problem point upstream, determining the time difference between the problem point and the first problem point upstream; for any problem point, based on the time difference between the problem point and the first problem point upstream, aligning the time series of water quality monitoring values ​​corresponding to the problem point with the time series of water quality monitoring values ​​corresponding to the first problem point upstream to determine the aligned time series of water quality monitoring values ​​corresponding to the problem point; for any problem point, performing time series decomposition on the aligned time series of water quality monitoring values ​​corresponding to the problem point to determine the water quality analysis sequence corresponding to the problem point.

[0077] In order to improve the accuracy and reliability of subsequent correlation analysis based on the water quality analysis sequences corresponding to any two adjacent problem points (i.e., two monitoring points that are adjacent in location and are both problem points), for any problem point, the time series of the flow velocity monitoring value corresponding to the problem point can be used to align the time series of the water quality monitoring value corresponding to the problem point, thereby eliminating the influence of the transmission time required for the same water flow to flow from the upstream problem point to the problem point.

[0078] Specifically, the first problem point upstream (i.e., the first problem point upstream) can be used as the benchmark for time series alignment. For any problem point other than the first problem point upstream, the flow velocities of the two problem points at the same collection time can be determined based on the flow velocity monitoring value time series corresponding to the problem point and the flow velocity monitoring value time series corresponding to the first problem point upstream, and the average value can be solved as the flow velocity of the same stream of water flowing from the most upstream stream to the problem point at the collection time. Furthermore, the required transfer time for the same stream of water to flow from the most upstream stream to the problem point can be determined in combination with the distance between the problem point and the first problem point upstream, as the time difference between the problem point and the first problem point upstream.

[0079] According to the time difference between the problem point and the first problem point upstream, the time series of water quality monitoring values ​​corresponding to the problem point can be aligned with the time series of water quality monitoring values ​​corresponding to the first problem point upstream.

[0080] Figure 3 FIG. 2 is a schematic diagram showing a time series alignment according to an embodiment of the present disclosure. Figure 3 As shown in the figure, problem point 1 is the first problem point upstream, and problem point 2 is the problem point downstream. For the monitoring value time series corresponding to the chemical oxygen demand at problem point 1, the concentration peak and peak time of the chemical oxygen demand at problem point 1 on November 11 can be determined based on the first black dot mark. The time required for the same stream of water to flow to problem point 2, that is, the time difference, is determined through the above process. The time difference can be used to move the monitoring value time series corresponding to the chemical oxygen demand at problem point 2 forward, so that the monitoring value time series corresponding to the chemical oxygen demand at problem point 2 is aligned with the monitoring value time series corresponding to the chemical oxygen demand at problem point 1.

[0081] After completing the alignment of the time series of different problem points, for any problem point, the aligned time series of water quality monitoring values ​​corresponding to the problem point can be decomposed to determine the water quality analysis sequence corresponding to the problem point; then, through the water quality analysis sequence corresponding to each problem point, combined with the positional relationship of the problem point in the target river channel, the problematic river sections that may have external sewage input can be screened to determine the river sections to be investigated.

[0082] In one possible implementation, at least one river section to be investigated in the target river channel is determined based on the location of each problem point and / or the water quality analysis sequence corresponding to each problem point, including: for any problem point, if the monitoring point upstream of the problem point is not a problem point, the river section between the problem point and the upstream adjacent monitoring point is determined as the river section to be investigated; or, for any two adjacent problem points, a multi-dimensional correlation analysis is performed based on the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet a first preset condition, and if the correlation analysis results between the two problem points meet the first preset condition, the river section between the two problem points is determined as the river section to be investigated, wherein any two adjacent problem points represent two monitoring points that are adjacent in location and are both problem points.

[0083] Taking into account the upstream and downstream transmissibility of sewage in the target river channel, as well as the dilution and degradation effects of sewage during flow, for any problem point, if the adjacent monitoring point upstream of the problem point is not a problem point, it can be directly determined that there is exogenous sewage input between the problem point and the adjacent monitoring point upstream of the problem point, thereby determining that the river section between the problem point and the adjacent monitoring point upstream of the problem point is to be investigated.

[0084] With the above Figure 2 For example, Figure 2 As shown in the figure, the black solid line represents the target river channel, and the circular marks represent the six monitoring points set up in the target river channel, dividing the target river channel into six river sections. Figure 2 As shown in (b), the dark circular mark at the far downstream represents the only problem point in the target river channel, and the five light gray circular marks upstream represent monitoring points that are not problem points. At this time, it can be directly determined that there is external sewage input between the problem point and its upstream adjacent monitoring points, and the river section between the problem point and its upstream adjacent monitoring points is determined as the river section to be investigated.

[0085] like Figure 2As shown in (c) in the figure, the dark gray circular marks indicate the problem points in the target river, namely, problem point 1, problem point 2, problem point 3 and problem point 4; the light gray circular marks indicate the monitoring points that are not problem points, among which monitoring point A is located between problem point 1 and problem point 2. At this time, it can be directly determined that there is external sewage input between problem point 2 and its upstream adjacent monitoring point A, and the river section between problem point 2 and monitoring point A is determined as the river section to be investigated. In contrast, since the water flow from problem point 1 does not cause abnormal water quality at monitoring point A when it flows to monitoring point A, it means that the natural dilution or degradation during the flow process is sufficient to eliminate the water quality abnormality. Therefore, there is no possibility of external sewage input between problem point 1 and monitoring point A, and the river section between problem point 1 and monitoring point A is not the river section to be investigated.

[0086] For any two adjacent problem points (i.e., two monitoring points that are adjacent in location and both are problem points), it is necessary to perform a multi-dimensional correlation analysis based on the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet the first preset condition, so as to judge whether the problem point located downstream is only affected by the problem point located upstream, and then determine that the river section between the two problem points is the river section to be investigated. Among them, the specific content of the first preset condition can be flexibly set according to actual use requirements, for example, including the correlation analysis results between any two adjacent problem points indicating that there is upstream and downstream transferability of sewage between the two adjacent problem points, etc., and the present disclosure does not make specific limitations on this. The specific method of correlation analysis can be flexibly set according to actual use requirements, depending on the specific form of the water quality analysis sequence corresponding to the problem point, and the present disclosure does not make specific limitations on this.

[0087] In a possible implementation, for any two adjacent problem points, a multi-dimensional correlation analysis is performed based on the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet the first preset condition, including: for any two adjacent problem points, based on the trend analysis sequences corresponding to the two problem points, determining the trend similarity between the two problem points; for any two adjacent problem points, based on the residual analysis sequences corresponding to the two problem points, determining the residual mutation similarity between the two problem points; for any two adjacent problem points, based on the periodic analysis sequences corresponding to the two problem points, determining the similarity of the change rules between the two problem points, as well as the mutual similarity between the two problem points. Correlation function value; for any two adjacent problem points, when the trend similarity corresponding to the two problem points meets the second preset condition, the change law similarity meets the third preset condition, the mutual correlation function value meets the fourth preset condition, and the residual mutation similarity meets the fifth preset condition, the correlation analysis result between the two problem points is determined to meet the first preset condition; for any two adjacent problem points, when the trend similarity corresponding to the two problem points does not meet the second preset condition, or the change law similarity does not meet the third preset condition, or the mutual correlation function value does not meet the fourth preset condition, or the residual mutation similarity does not meet the fifth preset condition, any one of the following is true, the correlation analysis result between the two problem points is determined to not meet the first preset condition.

[0088] Figure 4 A schematic diagram showing a time series decomposition result according to an embodiment of the present disclosure is shown. Figure 4 As shown in the figure, the time series of the monitoring values ​​of the chemical oxygen demand at any problem point is decomposed into a trend analysis series, a periodic analysis series and a residual analysis series. The trend analysis series can reflect the long-term changes in the monitoring values ​​of the chemical oxygen demand at the problem point; the periodic analysis series can reflect the periodic fluctuations in the monitoring values ​​of the chemical oxygen demand at the problem point; and the residual analysis series can generally represent the part of the monitoring values ​​of the chemical oxygen demand corresponding to the problem point that cannot be explained by the long-term trend and periodic fluctuations.

[0089] For any two adjacent problem points, a similarity analysis of long-term changes can be performed based on the trend analysis sequences corresponding to the two problem points to obtain the trend similarity between the two problem points.

[0090] Preferably, when the trend analysis sequences corresponding to the two problem points keep changing synchronously (including rising and falling simultaneously), the trend similarity between the two problem points can be determined as trend similarity. When the trend analysis sequences corresponding to the two problem points cannot keep changing synchronously, the trend similarity between the two problem points can be determined as trend dissimilarity.

[0091] In one example, problem point 1 is adjacent to problem point 2, and the trend analysis sequence corresponding to problem point 1 shows an upward trend from time point A to time point B, and remains relatively stable at other times. If the trend analysis sequence corresponding to problem point 2 also only shows an upward trend from time point A to time point B, and remains relatively stable at other times, it can be determined that the trend similarity between problem point 1 and problem point 2 is trend similarity. If the trend analysis sequence corresponding to problem point 2 does not maintain an upward trend from time point A to time point B, but shows an upward trend from time point A to time point C, and a downward trend from time point C to time B, it can be determined that the trend similarity between problem point 1 and problem point 2 is trend dissimilar.

[0092] In addition, for any two adjacent problem points, the Euclidean distance, Manhattan distance, cosine similarity and other data indicators for quantitatively measuring similarity can be calculated based on the trend analysis sequences corresponding to the two problem points, or the dynamic time warping method and other methods can be used to determine the trend similarity between the two problem points. The specific method for determining trend similarity can be flexibly set according to actual usage requirements, and the present disclosure does not make specific limitations on this.

[0093] It should be noted that when the trend similarity between any two adjacent problem points is trend similarity, the change rates of the trend analysis sequences corresponding to the two problem points may be different and do not need to be completely consistent.

[0094] For any two adjacent problem points, the similarity analysis of residual mutation time can be performed according to the residual analysis sequences corresponding to the two problem points to obtain the residual mutation similarity between the two problem points.

[0095] Preferably, when the difference between the time points corresponding to the mutation positions in the residual analysis sequence corresponding to the two problem points is less than or equal to a preset threshold, the residual mutation similarity between the two problem points can be determined to be mutation similarity. The mutation position can represent the position where the residual value changes sharply and significantly, and the definition of the mutation point of the statistical image in the relevant technology can be referred to, and the present disclosure does not make specific limitations on this.

[0096] With the above Figure 4 For example, Figure 4As shown, the residual value at time point A changes sharply from the lowest trough in the entire residual analysis sequence to the highest peak in the entire residual analysis sequence. Therefore, the position corresponding to time point A can be determined as the mutation position of the residual analysis sequence.

[0097] When the difference between the time points corresponding to the mutation positions in the residual analysis sequence corresponding to the two problem points is greater than the preset threshold, the residual mutation similarity between the two problem points can be determined as mutation dissimilarity. The specific value of the preset threshold here can be flexibly set according to actual usage requirements, and this disclosure does not make specific limitations on this.

[0098] In one example, problem point 1 is adjacent to problem point 2, the time point corresponding to the mutation position in the residual analysis sequence corresponding to problem point 1 is time point A, and the time point corresponding to the mutation position in the residual analysis sequence corresponding to problem point 2 is time point B. If the difference between time point A and time point B is less than or equal to a preset threshold, it can be determined that the residual mutation similarity between problem point 1 and problem point 2 is mutation similarity. If the difference between time point A and time point B is greater than a preset threshold, it can be determined that the residual mutation similarity between problem point 1 and problem point 2 is mutation dissimilar.

[0099] In addition, for any two adjacent problem points, the edit distance, Jaccard similarity coefficient and other data indicators for quantitatively measuring similarity can be calculated based on the residual analysis sequences corresponding to the two problem points, or the dynamic time warping method and other methods can be used to determine the residual mutation similarity between the two problem points. The specific method for determining the residual mutation similarity can be flexibly set according to actual usage requirements, and the present disclosure does not make specific limitations on this.

[0100] For any two adjacent problem points, we can perform similarity analysis on the periodic variation rules based on the periodic analysis sequences corresponding to the two problem points to obtain the similarity of the variation rules between the two problem points.

[0101] Preferably, when the time points corresponding to each zero position (i.e., the position where the data value is 0), each peak position, and each trough position in the periodic analysis sequence corresponding to the two problem points are the same, the similarity of the change rules between the two problem points can be determined as similar change rules. When there is a difference in the time points corresponding to at least one zero position, at least one peak position, or at least one trough position in the periodic analysis sequence corresponding to the two problem points, the similarity of the change rules between the two problem points can be determined as dissimilar change rules.

[0102] In addition, for any two adjacent problem points, the Pearson correlation coefficient, Spearman rank correlation coefficient and other data indicators for quantitatively measuring similarity can be calculated based on the periodic analysis sequences corresponding to the two problem points, or the mutual information analysis method, the longest common subsequence comparison method and other methods can be used to judge the similarity of the change patterns between the two problem points. The specific method for judging the similarity of the change patterns can be flexibly set according to actual usage requirements, and the present disclosure does not make specific limitations on this.

[0103] In one example, problem point 1 is adjacent to problem point 2. In the periodic analysis sequence corresponding to problem point 1, for any two adjacent problem points, a cross-correlation function can be constructed based on the periodic analysis sequences corresponding to the two problem points, and a cross-correlation analysis can be performed to obtain the cross-correlation function value between the two problem points, thereby quantifying and accurately expressing the cross-correlation between the two problem points. The specific form of the cross-correlation function here can refer to the implementation in the relevant technology, and the present disclosure does not specifically limit this.

[0104] In one example, the cross-correlation function can be expressed as formula (2):

[0105]

[0106] Among them, x represents any problem point; x(t) represents the periodic analysis sequence corresponding to the problem point x; y represents the problem point adjacent to the problem point x, among which the problem point y and the problem point x are two monitoring points that are adjacent in position and are both problem points; y(t) represents the periodic analysis sequence corresponding to the problem point y; t represents time.

[0107] By using the cross-correlation function shown in formula (2), the time dependence and similarity of the periodic analysis sequences corresponding to two adjacent problem points can be accurately reflected. Moreover, compared with other methods for determining the cross-correlation function value, this method can reduce the influence of the time difference between two adjacent problem points and more accurately reflect the correlation of the periodic analysis sequences corresponding to the two adjacent problem points at the time difference corresponding to the two adjacent problem points.

[0108] For any two problem points, when the trend similarity corresponding to the two problem points meets the second preset condition, the change law similarity meets the third preset condition, the mutual correlation function value meets the fourth preset condition, and the residual mutation similarity meets the fifth preset condition, it can be determined that the correlation analysis result between the two problem points meets the first preset condition.

[0109] Among them, the second preset condition may include that the trend similarity between any two problem points is trend similarity; the third preset condition may include that the change law similarity between any two problem points is change law similarity; the fourth preset condition may include that the mutual correlation function value between any two problem points is greater than or equal to a preset mutual correlation threshold. The specific value of the mutual correlation threshold here can be flexibly set according to actual usage requirements, and the present disclosure does not make specific restrictions on this; the fifth preset condition may include that the residual mutation similarity between any two problem points is mutation similarity.

[0110] Specifically, for any two problem points, when the trend similarity corresponding to the two problem points is trend similarity, the change law similarity is change law similarity, the mutual correlation function value is greater than the preset mutual correlation threshold, and the residual mutation similarity is mutation similarity, it can be determined that the correlation analysis result between the two problem points has upstream and downstream transmission, that is, it meets the first preset condition.

[0111] For any two problem points, if the trend similarity corresponding to the two problem points does not meet the second preset condition, or the change law similarity does not meet the third preset condition, or the mutual correlation function value does not meet the fourth preset condition, or the residual mutation similarity does not meet any one of the fifth preset conditions, it can be determined that the correlation analysis result between the two problem points does not meet the first preset condition.

[0112] Specifically, for any two problem points, if the trend similarity corresponding to the two problem points is trend dissimilarity, or the change law similarity is change law dissimilarity, or the correlation function value is less than or equal to the preset correlation threshold, or the residual mutation similarity is mutation dissimilarity, any one of the following is true, it can be determined that the correlation analysis result between the two problem points has no upstream and downstream transmission, that is, the first preset condition is not met.

[0113] When the correlation analysis results between any two problem points meet the first preset condition, the river section between the two problem points is determined as the river section to be checked.

[0114] Through the above process, based on the flow direction of the target river channel and utilizing the upstream and downstream transmissibility of sewage, we can exclude river sections that are only affected by upstream sewage and lock in river sections to be investigated with external sewage input, thereby narrowing the scope of water pollution tracing and improving the accuracy of water pollution tracing.

[0115] After determining the river section to be investigated, the outlet can be traced for the river section to be investigated. Specifically, for any river section to be investigated, if there is only one outlet in the river section to be investigated, the outlet can be directly determined as the first pollution source causing water pollution in the target river; if there are multiple outlets in the river section to be investigated, the first pollution source corresponding to the target river can be determined based on the water quality monitoring value time series and flow monitoring value time series corresponding to each outlet in the river section to be investigated.

[0116] In a possible implementation, the first pollution source corresponding to the target river is determined based on the time series of water quality monitoring values ​​and the time series of flow monitoring values ​​corresponding to each outlet in each river section to be inspected, including: for any outlet, when the monitoring value at any moment in the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, the pollution flux corresponding to the water quality indicator at the outlet is determined based on the time series of monitoring values ​​and the time series of flow monitoring values ​​corresponding to the water quality indicator at the outlet; for any outlet, when the monitoring value at any moment in the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator When the flux is greater than the first pollution flux threshold, the outlet is determined as the first pollution source corresponding to the target river; for any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is less than or equal to the first pollution flux threshold and greater than or equal to the second pollution flux threshold, determine whether the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition based on the water quality monitoring value time series corresponding to the outlet, and when the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition, the outlet is determined as the first pollution source corresponding to the target river.

[0117] For any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, the pollution flux corresponding to the water quality indicator at the outlet can be determined based on the water quality monitoring value time series and flow monitoring value time series corresponding to the outlet. Among them, the specific method for determining the pollution flux can refer to the implementation method in the relevant technology, and the present disclosure does not make specific limitations on this. The specific value of the third preset threshold value corresponding to any water quality indicator here can be flexibly set according to actual usage requirements, and the present disclosure does not make specific limitations on this.

[0118] For any outlet, if the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is greater than the first pollution flux threshold, the outlet can be directly determined as the first pollution source corresponding to the target river. Among them, the specific value of the first pollution flux threshold can be flexibly set according to actual use requirements. For example, the water quality indicator limit of Class V water * the average flow of the outlet can be selected as the first pollution flux threshold, etc., based on the functional classification of Class I, II, III, IV, V, and Class V surface water grades in the "Surface Water Environmental Quality Standards". This disclosure does not make specific restrictions on this.

[0119] For any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is less than or equal to the first pollution flux threshold value and greater than or equal to the second pollution flux threshold value, it is necessary to conduct a linkage analysis on the outlet and the river section to be inspected where it is located, and determine the correlation analysis result between the outlet and the problem point adjacent to its downstream; and when the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition, the outlet can be determined as a problem outlet.

[0120] Among them, the specific method for determining the correlation analysis result between any one outlet and its downstream adjacent problem point can refer to the process of the correlation analysis method for any two adjacent problem points recorded in the previous text; the specific content of the sixth preset condition can be flexibly set according to the actual use requirements, and can include the correlation analysis result between any one outlet and its downstream adjacent problem point, indicating that there is upstream and downstream transferability of sewage between the outlet and its downstream adjacent problem point, etc., and the present disclosure does not make specific restrictions on this. The specific value of the second pollution flux threshold can be flexibly set according to the actual use requirements. For example, you can refer to the functional classification of Class I, II, III, IV, V, and Class V surface water grades in the "Surface Water Environmental Quality Standards", and select the water quality index limit of Class III water * the average outlet flow as the second pollution flux threshold, etc. The present disclosure does not make specific restrictions on this.

[0121] For any outlet, if the monitoring value at each moment in the monitoring value time series corresponding to each water quality indicator at the outlet does not exceed the third preset threshold value corresponding to each water quality indicator, or the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, but the pollution flux corresponding to the water quality indicator is less than the second pollution flux threshold value, it can be directly determined that the outlet is not the main factor causing abnormal water quality at the adjacent problem point downstream, and there is no need to analyze the drainage network system upstream of the outlet.

[0122] Figure 5 FIG. 1 is a schematic diagram of a river section to be investigated according to an embodiment of the present disclosure. Figure 5 As shown, each of the river sections to be investigated 1, 2 and 3 includes two outlets, namely outlet 1 to outlet 6. Using the above analysis method, outlet 1, outlet 3 and outlet 5 can be determined as the first pollution source corresponding to the target river.

[0123] Through the above process, we can use the comprehensive linkage analysis between the outlets and the problem points to determine the outlets in each river section to be investigated that cause abnormal water quality at adjacent problem points downstream, determine the first pollution source that directly affects the water quality of the target river, improve the efficiency and accuracy of water pollution tracing, and improve the timeliness of water pollution control and management.

[0124] In a possible implementation, the method also includes: when the first pollution source corresponding to the target river channel is a rainwater outlet, determining the rainwater pipe network system corresponding to the rainwater outlet; determining the liquid level distribution information, flow velocity distribution information and water quality distribution information of the rainwater pipe network system during the dry period according to the real-time monitoring data corresponding to the rainwater pipe network system, wherein the real-time monitoring data corresponding to the rainwater pipe network system includes a liquid level monitoring value time series, a flow velocity monitoring value time series and a monitoring value time series corresponding to at least one water quality indicator at each monitoring position in the rainwater pipe network system; determining the abnormal liquid level position of the rainwater pipe network system according to the liquid level distribution information; determining the abnormal flow velocity position of the rainwater pipe network system according to the flow velocity distribution information; determining the abnormal water quality position of the rainwater pipe network system according to the water quality distribution information; determining the second pollution source corresponding to the target river channel according to the abnormal liquid level position, abnormal flow velocity position and abnormal water quality position of the rainwater pipe network system.

[0125] When the first pollution source corresponding to the target river is a rainwater outlet, the rainwater pipe network system upstream of the rainwater outlet can be determined, and the real-time monitoring data corresponding to the rainwater pipe network system can be divided into real-time monitoring data during the dry period and real-time monitoring data during the rainy period according to the weather conditions. Then, according to the real-time monitoring data during the dry period, the liquid level distribution information, flow velocity distribution information and water quality distribution information of the rainwater pipe network system during the dry period can be drawn with high and low values. Among them, the real-time monitoring data corresponding to the rainwater pipe network system include a time series of liquid level monitoring values, a time series of flow rate monitoring values ​​and a time series of monitoring values ​​corresponding to at least one water quality indicator corresponding to each monitoring position in the rainwater pipe network system. The specific form of each time series can refer to the time series corresponding to any monitoring point; the monitoring position of the rainwater pipe network system can be flexibly set according to actual usage requirements, and the present disclosure does not make specific limitations on this; among them, the specific forms of liquid level distribution information, flow rate distribution information and water quality distribution information can refer to the implementation methods in the relevant technology. For example, the liquid level distribution information can be set as a liquid level thermodynamic map, the flow rate distribution information can be set as a flow rate thermodynamic map, the water quality distribution information can be set as a water quality thermodynamic map, etc. The present disclosure does not make specific limitations on this.

[0126] According to the liquid level distribution information of the rainwater pipe network system during the dry period, the abnormal liquid level position of the rainwater pipe network system can be determined. Specifically, during the dry period, the rainwater pipe network system is not affected by rainfall runoff, and the amount of water in the rainwater pipe network system is usually small, or there is no water at all, and the corresponding liquid level is usually low; when the liquid level corresponding to any monitoring position in the rainwater pipe network system is high, it means that this section of the rainwater pipe may have problems such as mixed or misconnected sewage into the rainwater pipe (rainwater and sewage mixing), or groundwater infiltration (external water intrusion). Therefore, when there is a location where the liquid level is significantly higher than other areas in the liquid level distribution information of the rainwater pipe network system during the dry period, this location can be determined as an abnormal liquid level location in the rainwater pipe network system, which can be used as a basis for in-depth investigation of the fundamental factors affecting the water quality of the target river.

[0127] According to the flow velocity distribution information of the rainwater pipe network system during the dry period, the abnormal flow velocity position of the rainwater pipe network system can be determined. Specifically, during the dry period, the rainwater pipe network system is not affected by rainfall runoff, the amount of water in the rainwater pipe network system is usually small, or there is no water at all, and the corresponding flow velocity is usually low; when the flow velocity corresponding to any section of the rainwater pipe in the rainwater pipe network system is high, it means that this section of the rainwater pipe may have problems such as rainwater and sewage mixing, external water intrusion, etc. Therefore, when there is a location with a flow velocity significantly higher than other areas in the liquid level distribution information of the rainwater pipe network system during the dry period, this location can be determined as the abnormal flow velocity location in the rainwater pipe network system, which can be used as the basis for in-depth investigation of the fundamental factors affecting the water quality of the target river.

[0128] According to the water quality distribution information of the rainwater pipe network system during the dry period, the abnormal water quality location of the rainwater pipe network system can be determined. When the monitoring value of any water quality index in the water quality distribution information of the rainwater pipe network system during the dry period is significantly higher than that of other areas, it indicates that there may be external sewage input at this location, and this location can be determined as the abnormal water quality location in the rainwater pipe network system as a basis for in-depth investigation of the fundamental factors affecting the water quality of the target river. Among them, the water quality indicators in the water quality distribution information corresponding to the rainwater pipe network system may include the following indicators: chemical oxygen demand, conductivity, ammonia nitrogen, five-day biochemical oxygen demand, total phosphorus, total nitrogen, suspended solids, total dissolved solids, petroleum, PH, anionic surfactants, cyanide, sulfide, fluoride, chloride, organic phosphorus, sulfate, mercury, chromium, cadmium, arsenic, lead, nickel, beryllium, silver, selenium, copper, zinc, manganese, iron, volatile phenols, benzene series, aniline, nitrobenzene.

[0129] According to the abnormal liquid level positions, abnormal flow velocity positions and abnormal water quality positions of the rainwater pipe network system, it is possible to determine the overlapping positions where the abnormal liquid level, abnormal flow velocity and abnormal water quality exist at the same time. These overlapping positions obviously have abnormal sewage input, which is transmitted through the sewage in the rainwater pipe network system and finally passes through the rainwater outlet, that is, the first pollution source corresponding to the target river, affecting the water quality of the target river. Therefore, these overlapping positions can be determined as the second pollution source corresponding to the target river.

[0130] Through the above process, the abnormal liquid level position, abnormal flow velocity position and abnormal water quality position in the upstream rainwater pipe network system of each problematic outlet can be comprehensively analyzed from the three aspects of liquid level, flow velocity and water quality, and the overlapping position of these three abnormal positions can be determined as the specific mixing or inflow and infiltration position, which has high accuracy and reliability, and can provide a reliable basis for determining the target pollution source.

[0131] In a possible implementation, the method also includes: when the first pollution source corresponding to the target river channel is an overflow outlet, determining the sewage pipe network system corresponding to the overflow outlet; determining the rainwater increment during the rainy day period corresponding to each monitoring location according to the flow monitoring value time series corresponding to each monitoring location in the sewage pipe network system; for any monitoring location, when the rainwater increment during the rainy day period corresponding to the monitoring location is greater than or equal to a preset increment threshold, determining the monitoring location as the third pollution source corresponding to the target river channel.

[0132] When the first pollution source corresponding to the target river channel is an overflow port, the sewage pipe network system upstream of the overflow port can be determined, and the flow monitoring value time series corresponding to the sewage pipe network system can be classified according to the weather conditions, and the data in the dry day period and the data in the rainy day period in the flow monitoring value time series can be distinguished, and then the rainwater increment in the rainy day period corresponding to each monitoring position can be determined. Among them, the monitoring position of the sewage pipe network system can be flexibly set according to actual use requirements, and this disclosure does not make specific restrictions on this; the specific method for determining the rainwater increment in the rainy day period corresponding to each monitoring position can be flexibly set according to actual use requirements, and this disclosure does not make specific restrictions on this.

[0133] In one example, for any monitoring location, the average flow rate of the monitoring location during the dry day period can be calculated based on the flow rate monitoring value time series of the monitoring location in at least one dry day before the rainy day period; the average flow rate of the monitoring location during the rainy day period can be calculated based on the flow rate monitoring value time series of the monitoring location during the rainy day period; and then the rainwater increment of the rainy day period corresponding to the monitoring location can be determined based on the average flow rate of the monitoring location during the rainy day period and the average flow rate during the dry day period. Specifically, the rainwater increment of the rainy day period corresponding to any monitoring location can be expressed as formula (3):

[0134]

[0135] Among them, Q R Indicates the rainfall increment during the rainy period; q rain It represents the average flow rate of the monitoring location during the rainy period; q dry It indicates the average flow value of the monitoring location during the dry period.

[0136] During rainy days, the sewage pipe network system is usually less affected by rainfall runoff. That is to say, in the sewage pipe network system, the flow rate of any section of sewage pipe during rainy days is usually smaller than the flow rate change during dry days. Therefore, for any monitoring location, when the rainwater increment during the rainy day period corresponding to the monitoring location is greater than or equal to the preset increment threshold, it can be judged that there are problems such as rainwater and sewage mixing and external water intrusion at the monitoring location, which will cause abnormal sewage discharge in the sewage pipe network system, and further affect the water quality of the target river. At this time, the monitoring location can be determined as the third pollution source corresponding to the target river.

[0137] Through the above process, different types of monitoring indicators such as liquid level, flow rate, water quality, and flow rate can be comprehensively utilized according to the type of drainage pipe network system to analyze the locations of problems such as rainwater and sewage mixing, external water intrusion, etc. in the drainage pipe network system with high accuracy and reliability, thereby providing a reliable basis for determining the second or third pollution source corresponding to the target river; and constructing a complete pollution tracing chain of monitoring point-river section-outlet-abnormal position of the pipe network, progressively exploring the root causes affecting the water quality of the target river, rather than just stopping at tracing the direct factors on the surface such as the outlet, which can further improve the effect of water pollution tracing.

[0138] In the embodiments of the present disclosure, a variety of monitoring indicators such as water quality indicators, flow rate and flow rate can be flexibly selected for real-time monitoring in different monitoring scenarios such as monitoring points and outlets. According to the monitoring value time series corresponding to at least one water quality indicator at each monitoring point in the selected range area in the target river, real-time monitoring of abnormal water quality phenomena can be achieved. Determine at least one problem point. According to the location of each problem point, or the water quality analysis sequence corresponding to each problem point, the target river can be segmented and checked from two aspects: the location distribution of the problem point and the sewage transmissibility, and at least one river section to be checked in the target river can be determined, thereby effectively narrowing the scope of water pollution tracing and improving the efficiency and accuracy of water pollution tracing; wherein, the water quality analysis sequence corresponding to any problem point is determined by time series decomposition based on the water quality monitoring value time series and the flow rate monitoring value time series corresponding to the problem point, which is compared with the prior art method of chemical determination based on on-site sampling or the use of mathematical methods. The method of model simulation and prediction disclosed in the present invention uses real-time monitoring data presenting a time series to perform time series decomposition, which can deeply explore the long-term changes and periodic fluctuations of at least one water quality indicator at the problem point, thereby providing reliable data support for water pollution source tracing, and improving the accuracy and reliability of water pollution source tracing. According to the water quality monitoring value time series and flow monitoring value time series corresponding to each outlet in each river section to be investigated, the first pollution source corresponding to the target river channel can be determined, wherein the water quality monitoring value time series corresponding to any outlet includes the monitoring value time series corresponding to at least one water quality indicator at the outlet. Furthermore, it is also possible to combine the type of drainage pipe network system corresponding to the outlet, and comprehensively utilize different monitoring data such as liquid level, flow rate, water quality, and flow rate to analyze the location of problems such as rainwater and sewage mixing and external water intrusion in the drainage pipe network system, rather than just tracing the source of direct factors such as the outlet. This can further improve the effect of water pollution control, have high accuracy and reliability, and thus provide a reliable basis for determining the second or third pollution source corresponding to the target river; and construct a complete pollution tracing chain of monitoring point-river section-outlet-abnormal position of the pipe network, and achieve water pollution tracing efficiently and quickly by progressively narrowing the scope of pollution tracing layer by layer. Compared with the mathematical model simulation method, the water pollution tracing method provided by the present disclosure can make full use of real-time monitoring data, has high accuracy and reliability, and reduces the workload of data processing.

[0139] In addition, the present disclosure also provides a water pollution source tracing device, an electronic device and a computer-readable storage medium, all of which can be used to implement any water pollution source tracing method provided by the present disclosure. The corresponding technical solutions and descriptions can be found in the corresponding records in the method part and will not be repeated here.

[0140] Figure 6 FIG. 1 is a block diagram of a water pollution source tracing device according to an embodiment of the present disclosure. Figure 6As shown, the device 600 includes:

[0141] The problem point determination module 601 is used to determine at least one problem point according to the time series of water quality monitoring values ​​corresponding to each monitoring point in the selected range area in the target river channel, wherein the time series of water quality monitoring values ​​corresponding to any monitoring point includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point;

[0142] The river section inspection module 602 is used to determine at least one river section to be inspected in the target river channel according to the location of each problem point and / or the water quality analysis sequence corresponding to each problem point, wherein the water quality analysis sequence corresponding to any problem point is determined by time series decomposition based on the water quality monitoring value time series and the flow rate monitoring value time series corresponding to the problem point, and the water quality analysis sequence corresponding to the problem point includes a trend analysis sequence, a periodicity analysis sequence and a residual analysis sequence corresponding to the problem point;

[0143] The first pollution source determination module 603 is used to determine the first pollution source corresponding to the target river channel according to the water quality monitoring value time series and flow monitoring value time series corresponding to each outlet in each river section to be inspected, wherein the water quality monitoring value time series corresponding to any outlet includes the monitoring value time series corresponding to at least one water quality indicator at the outlet.

[0144] In one possible implementation, the problem point determination module 601 is also used for: during a rainy period, for any monitoring point, if the monitoring value at any moment in the monitoring value time series corresponding to any water quality indicator at the monitoring point exceeds the first preset threshold value corresponding to the water quality indicator, the monitoring point is determined as a problem point; or, during a dry period, for any monitoring point, if the duration for which the monitoring value at a single moment in the monitoring value time series corresponding to any water quality indicator at the monitoring point exceeds the second preset threshold value corresponding to the water quality indicator meets the first preset duration, and / or the duration for which the monitoring value time series corresponding to any water quality indicator conforms to the sewage change law meets the second preset duration, the monitoring point is determined as a problem point, wherein the sewage change law is used to reflect the daily cycle law of the change over time of the monitoring value of any water quality indicator that exceeds the pollution standard.

[0145] In one possible implementation, the device 600 also includes a time series decomposition module, which is used to: for any problem point, determine the time difference between the problem point and the first problem point upstream based on the flow velocity monitoring value time series corresponding to the problem point and the distance between the problem point and the first problem point upstream; for any problem point, align the water quality monitoring value time series corresponding to the problem point with the water quality monitoring value time series corresponding to the first problem point upstream based on the time difference between the problem point and the first problem point upstream, and determine the aligned water quality monitoring value time series corresponding to the problem point; for any problem point, perform time series decomposition on the aligned water quality monitoring value time series corresponding to the problem point to determine the water quality analysis sequence corresponding to the problem point.

[0146] In one possible implementation, the river section investigation module 602 is also used for: for any problem point, if the monitoring point upstream of the problem point is not a problem point, determining the river section between the problem point and the upstream adjacent monitoring point as the river section to be investigated; for any two adjacent problem points, performing a multi-dimensional correlation analysis based on the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet the first preset condition, and if the correlation analysis results between the two problem points meet the first preset condition, determining the river section between the two problem points as the river section to be investigated.

[0147] In a possible implementation, the river section inspection module 602 is also used to: for any two adjacent problem points, determine the trend similarity between the two problem points according to the trend analysis sequences corresponding to the two problem points; for any two adjacent problem points, determine the similarity of the change rules between the two problem points and the cross-correlation function value between the two problem points according to the periodic analysis sequences corresponding to the two problem points; for any two adjacent problem points, determine the residual mutation similarity between the two problem points according to the residual analysis sequences corresponding to the two problem points; for any two adjacent problem points, determine the similarity of the residual mutation between the two problem points according to the trend analysis sequences corresponding to the two problem points. When the trend similarity meets the second preset condition, the change law similarity meets the third preset condition, the mutual correlation function value meets the fourth preset condition, and the residual mutation similarity meets the fifth preset condition, the correlation analysis result between the two problem points is determined to meet the first preset condition; for any two adjacent problem points, when the trend similarity corresponding to the two problem points does not meet the second preset condition, or the change law similarity does not meet the third preset condition, or the mutual correlation function value does not meet the fourth preset condition, or the residual mutation similarity does not meet the fifth preset condition, any one of the following is true, the correlation analysis result between the two problem points is determined to not meet the first preset condition.

[0148] In a possible implementation, the first pollution source determination module 603 is also used to: for any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, determine the pollution flux corresponding to the water quality indicator at the outlet according to the water quality monitoring value time series and the flow monitoring value time series corresponding to the outlet; for any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is greater than the first pollution flux threshold value, determine the outlet. It is the first pollution source corresponding to the target river; for any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is less than or equal to the first pollution flux threshold value and greater than or equal to the second pollution flux threshold value, determine whether the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition according to the water quality monitoring value time series corresponding to the outlet, and when the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition, determine the outlet as the first pollution source corresponding to the target river.

[0149] In one possible implementation, the device 600 also includes a second pollution source determination module, which is used to: when the first pollution source corresponding to the target river channel is a rainwater outlet, determine the rainwater pipe network system corresponding to the rainwater outlet; determine the liquid level distribution information, flow velocity distribution information and water quality distribution information of the rainwater pipe network system during the dry period according to the real-time monitoring data corresponding to the rainwater pipe network system, wherein the real-time monitoring data corresponding to the rainwater pipe network system includes a liquid level monitoring value time series, a flow velocity monitoring value time series and a monitoring value time series corresponding to at least one water quality indicator at each monitoring position in the rainwater pipe network system; determine the abnormal liquid level position of the rainwater pipe network system according to the liquid level distribution information; determine the abnormal flow velocity position of the rainwater pipe network system according to the flow velocity distribution information; determine the abnormal water quality position of the rainwater pipe network system according to the water quality distribution information; determine the second pollution source corresponding to the target river channel according to the abnormal liquid level position, abnormal flow velocity position and abnormal water quality position of the rainwater pipe network system.

[0150] In one possible implementation, the second pollution source determination module is also used for: when the first pollution source corresponding to the target river channel is the overflow outlet, the sewage pipe network system corresponding to the overflow outlet; according to the flow monitoring value time series corresponding to each monitoring location in the sewage pipe network system, respectively determine the rainwater increment during the rainy day period corresponding to each monitoring location; for any monitoring location, when the rainwater increment during the rainy day period corresponding to the monitoring location is greater than or equal to a preset increment threshold, determine the monitoring location as the third pollution source corresponding to the target river channel.

[0151] In some embodiments, the functions or modules included in the device provided by the embodiments of the present disclosure can be used to execute the method described in the above method embodiments. The specific implementation can refer to the description of the above method embodiments, and for the sake of brevity, it will not be repeated here.

[0152] The embodiment of the present disclosure also provides a computer-readable storage medium on which computer program instructions are stored, and the computer program instructions implement the above method when executed by a processor. The computer-readable storage medium can be a volatile or non-volatile computer-readable storage medium.

[0153] An embodiment of the present disclosure further proposes an electronic device, comprising: a processor; and a memory for storing instructions executable by the processor; wherein the processor is configured to implement the above method when executing the instructions stored in the memory.

[0154] Figure 7 1900 is a block diagram of an electronic device according to an embodiment of the present disclosure. For example, the apparatus 1900 may be provided as a server or a terminal device. Figure 7, the apparatus 1900 includes a processing component 1922, which further includes one or more processors, and a memory resource represented by a memory 1932 for storing instructions, such as an application, that can be executed by the processing component 1922. The application stored in the memory 1932 may include one or more modules, each corresponding to a set of instructions. In addition, the processing component 1922 is configured to execute instructions to perform the above method.

[0155] The device 1900 may also include a power supply component 1926 configured to perform power management of the device 1900, a wired or wireless network interface 1950 configured to connect the device 1900 to a network, and an input / output interface 1958 (I / O interface). The device 1900 may operate based on an operating system stored in the memory 1932, such as Windows Server 2000. TM , MacOS X TM , Unix TM , Linux TM , FreeBSD TM or similar.

[0156] In an exemplary embodiment, a non-volatile computer-readable storage medium is also provided, such as a memory 1932 including computer program instructions, which can be executed by the processing component 1922 of the device 1900 to perform the above method.

[0157] The present disclosure may be a system, a method and / or a computer program product. The computer program product may include a computer-readable storage medium carrying computer-readable program instructions for causing a processor to implement various aspects of the present disclosure.

[0158] A computer-readable storage medium may be a tangible device that can hold and store instructions used by an instruction execution device. A computer-readable storage medium may be, for example, but not limited to, an electrical storage device, a magnetic storage device, an optical storage device, an electromagnetic storage device, a semiconductor storage device, or any suitable combination of the above. More specific examples (a non-exhaustive list) of computer-readable storage media include: a portable computer disk, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or flash memory), a static random access memory (SRAM), a portable compact disk read-only memory (CD-ROM), a digital versatile disk (DVD), a memory stick, a floppy disk, a mechanical encoding device, such as a punch card or a raised structure in a groove on which instructions are stored, and any suitable combination of the above. The computer-readable storage medium used herein is not to be interpreted as a transient signal itself, such as a radio wave or other freely propagating electromagnetic wave, an electromagnetic wave propagating through a waveguide or other transmission medium (e.g., a light pulse through a fiber optic cable), or an electrical signal transmitted through a wire.

[0159] The computer-readable program instructions described herein can be downloaded from a computer-readable storage medium to each computing / processing device, or downloaded to an external computer or external storage device via a network, such as the Internet, a local area network, a wide area network, and / or a wireless network. The network can include copper transmission cables, optical fiber transmissions, wireless transmissions, routers, firewalls, switches, gateway computers, and / or edge servers. The network adapter card or network interface in each computing / processing device receives the computer-readable program instructions from the network and forwards the computer-readable program instructions for storage in the computer-readable storage medium in each computing / processing device.

[0160] The computer program instructions for performing the operations of the present disclosure may be assembly instructions, instruction set architecture (ISA) instructions, machine instructions, machine-dependent instructions, microcode, firmware instructions, state setting data, or source code or object code written in any combination of one or more programming languages, including object-oriented programming languages, such as Smalltalk, C++, etc., and conventional procedural programming languages, such as "C" language or similar programming languages. The computer-readable program instructions may be executed entirely on the user's computer, partially on the user's computer, as a separate software package, partially on the user's computer, partially on a remote computer, or entirely on a remote computer or server. In the case of a remote computer, the remote computer may be connected to the user's computer through any type of network, including a local area network (LAN) or a wide area network (WAN), or may be connected to an external computer (e.g., using an Internet service provider to connect through the Internet). In some embodiments, by using the state information of the computer-readable program instructions to personalize an electronic circuit, such as a programmable logic circuit, a field programmable gate array (FPGA), or a programmable logic array (PLA), the electronic circuit may execute the computer-readable program instructions, thereby implementing various aspects of the present disclosure.

[0161] Various aspects of the present disclosure are described herein with reference to flowcharts and / or block diagrams of methods, devices (systems) and computer program products according to embodiments of the present disclosure. It should be understood that each box in the flowchart and / or block diagram and the combination of boxes in the flowchart and / or block diagram can be implemented by computer-readable program instructions.

[0162] These computer-readable program instructions can be provided to a processor of a general-purpose computer, a special-purpose computer, or other programmable data processing device, thereby producing a machine, so that when these instructions are executed by the processor of the computer or other programmable data processing device, a device that implements the functions / actions specified in one or more boxes in the flowchart and / or block diagram is generated. These computer-readable program instructions can also be stored in a computer-readable storage medium, and these instructions cause the computer, programmable data processing device, and / or other equipment to work in a specific manner, so that the computer-readable medium storing the instructions includes a manufactured product, which includes instructions for implementing various aspects of the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0163] Computer-readable program instructions may also be loaded onto a computer, other programmable data processing apparatus, or other device so that a series of operating steps are performed on the computer, other programmable data processing apparatus, or other device to produce a computer-implemented process, thereby causing the instructions executed on the computer, other programmable data processing apparatus, or other device to implement the functions / actions specified in one or more boxes in the flowchart and / or block diagram.

[0164] The flow chart and block diagram in the accompanying drawings show the possible architecture, function and operation of the system, method and computer program product according to multiple embodiments of the present disclosure. In this regard, each square box in the flow chart or block diagram can represent a part of a module, program segment or instruction, and the part of the module, program segment or instruction contains one or more executable instructions for realizing the specified logical function. In some alternative implementations, the function marked in the square box can also occur in a sequence different from that marked in the accompanying drawings. For example, two continuous square boxes can actually be executed substantially in parallel, and they can sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each square box in the block diagram and / or flow chart, and the combination of the square boxes in the block diagram and / or flow chart can be implemented with a dedicated hardware-based system that performs the specified function or action, or can be implemented with a combination of special hardware and computer instructions.

[0165] The embodiments of the present disclosure have been described above, and the above description is exemplary, not exhaustive, and is not limited to the disclosed embodiments. Many modifications and changes will be apparent to those of ordinary skill in the art without departing from the scope and spirit of the described embodiments. The selection of terms used herein is intended to best explain the principles of the embodiments, practical applications, or technical improvements in the market, or to enable other persons of ordinary skill in the art to understand the embodiments disclosed herein.

Claims

1. A method for tracing the source of water pollution, characterized in that: include: Determine at least one problem point according to the time series of water quality monitoring values ​​corresponding to each monitoring point in the selected range area in the target river channel, wherein the time series of water quality monitoring values ​​corresponding to any one monitoring point includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point; According to the location of each problem point and / or the water quality analysis sequence corresponding to each problem point, at least one river section to be checked in the target river channel is determined, wherein the water quality analysis sequence corresponding to any problem point is determined by time series decomposition based on the time series of water quality monitoring values ​​and flow rate monitoring values ​​corresponding to the problem point, and the water quality analysis sequence corresponding to the problem point includes a trend analysis sequence, a periodicity analysis sequence and a residual analysis sequence corresponding to the problem point; According to the time series of water quality monitoring values ​​and flow monitoring values ​​corresponding to each outlet in each river section to be inspected, the first pollution source corresponding to the target river is determined, wherein the time series of water quality monitoring values ​​corresponding to any outlet includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the outlet.

2. The method according to claim 1, characterized in that The step of determining at least one problem point according to the time series of water quality monitoring values ​​corresponding to each monitoring point within the selected range area in the target river channel comprises: During the rainy period, for any monitoring point, if the monitoring value at any time in the monitoring value time series corresponding to any water quality indicator at the monitoring point exceeds the first preset threshold value corresponding to the water quality indicator, the monitoring point is determined as a problem point; or, During the drought period, for any monitoring point, if the duration that the monitoring value at a single moment in the monitoring value time series corresponding to any water quality indicator at the monitoring point exceeds the second preset threshold value corresponding to the water quality indicator meets the first preset duration, and / or the duration that the monitoring value time series corresponding to any water quality indicator conforms to the sewage change law meets the second preset duration, the monitoring point is determined as a problem point, wherein the sewage change law is used to reflect the daily cycle law of the monitoring value of any water quality indicator that exceeds the pollution standard changing over time.

3. The method according to claim 1 or 2, characterized in that: The method further comprises: For any problem point, determine the time difference between the problem point and the first problem point upstream according to the flow velocity monitoring value time series corresponding to the problem point and the distance between the problem point and the first problem point upstream; For any problem point, according to the time difference between the problem point and the first problem point upstream, the time series of water quality monitoring values ​​corresponding to the problem point is aligned with the time series of water quality monitoring values ​​corresponding to the first problem point upstream, and the aligned time series of water quality monitoring values ​​corresponding to the problem point is determined; For any problem point, the aligned time series of water quality monitoring values ​​corresponding to the problem point is decomposed to determine the water quality analysis sequence corresponding to the problem point.

4. The method according to claim 1 or 2, characterized in that: The step of determining at least one river section to be investigated in the target river channel according to the location of each problem point and / or the water quality analysis sequence corresponding to each problem point includes: For any problem point, if the monitoring point upstream of the problem point is not a problem point, the river section between the problem point and the monitoring point upstream is determined as the river section to be checked; or, For any two adjacent problem points, a multi-dimensional correlation analysis is performed according to the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet the first preset condition, and when the correlation analysis results between the two problem points meet the first preset condition, the river section between the two problem points is determined as the river section to be checked, wherein the any two adjacent problem points represent two monitoring points that are adjacent in location and are both problem points.

5. The method according to claim 4, characterized in that The method of performing a multi-dimensional correlation analysis on any two adjacent problem points according to the water quality analysis sequences corresponding to the two problem points to determine whether the correlation analysis results between the two problem points meet the first preset condition includes: For any two adjacent problem points, determine the trend similarity between the two problem points according to the trend analysis sequences corresponding to the two problem points; For any two adjacent problem points, the residual mutation similarity between the two problem points is determined based on the residual analysis sequences corresponding to the two problem points; For any two adjacent problem points, determine the similarity of the change rules between the two problem points and the mutual correlation function value between the two problem points according to the periodic analysis sequences corresponding to the two problem points; For any two adjacent problem points, when the trend similarity corresponding to the two problem points meets the second preset condition, the change law similarity meets the third preset condition, the mutual correlation function value meets the fourth preset condition, and the residual mutation similarity meets the fifth preset condition, it is determined that the correlation analysis result between the two problem points meets the first preset condition; For any two adjacent problem points, if the trend similarity corresponding to the two problem points does not satisfy the second preset condition, or the change law similarity does not satisfy the third preset condition, or the mutual correlation function value does not satisfy the fourth preset condition, or the residual mutation similarity does not satisfy any one of the fifth preset conditions, it is determined that the correlation analysis result between the two problem points does not satisfy the first preset condition.

6. The method according to claim 1 or 2, characterized in that: The determining the first pollution source corresponding to the target river channel according to the water quality monitoring value time series and the flow monitoring value time series corresponding to each outlet in each river section to be investigated includes: For any outlet, when a monitoring value at any moment in a monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds a third preset threshold value corresponding to the water quality indicator, a pollution flux corresponding to the water quality indicator at the outlet is determined based on the monitoring value time series corresponding to the water quality indicator at the outlet and the flow monitoring value time series; For any outlet, if the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is greater than the first pollution flux threshold value, the outlet is determined as the first pollution source corresponding to the target river; For any outlet, when the monitoring value at any moment in the monitoring value time series corresponding to at least one water quality indicator at the outlet exceeds the third preset threshold value corresponding to the water quality indicator, and the pollution flux corresponding to the water quality indicator is less than or equal to the first pollution flux threshold value and greater than or equal to the second pollution flux threshold value, determine whether the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition based on the water quality monitoring value time series corresponding to the outlet, and when the correlation analysis result between the outlet and the problem point adjacent to its downstream meets the sixth preset condition, determine the outlet as the first pollution source corresponding to the target river.

7. The method according to claim 1 or 2, characterized in that: The method further comprises: In a case where the first pollution source corresponding to the target river channel is a rainwater outlet, determining a rainwater pipe network system corresponding to the rainwater outlet; According to the real-time monitoring data corresponding to the rainwater pipe network system, determine the liquid level distribution information, flow rate distribution information and water quality distribution information of the rainwater pipe network system during the dry period, wherein the real-time monitoring data corresponding to the rainwater pipe network system includes a time series of liquid level monitoring values, a time series of flow rate monitoring values ​​and a time series of monitoring values ​​corresponding to at least one water quality indicator corresponding to each monitoring position in the rainwater pipe network system; Determining the abnormal liquid level position of the rainwater pipe network system according to the liquid level distribution information; Determining the abnormal flow velocity position of the rainwater pipe network system according to the flow velocity distribution information; Determining the location of water quality anomaly in the rainwater pipe network system according to the water quality distribution information; According to the abnormal liquid level position, abnormal flow velocity position and abnormal water quality position of the rainwater pipe network system, the second pollution source corresponding to the target river channel is determined.

8. The method according to claim 1 or 2, characterized in that: The method further comprises: In a case where the first pollution source corresponding to the target river channel is an overflow outlet, determining a sewage pipe network system corresponding to the overflow outlet; Determine the rainwater increment in the rainy day period corresponding to each monitoring position according to the flow monitoring value time series corresponding to each monitoring position in the sewage pipe network system; For any monitoring location, when the rainwater increment during the rainy day period corresponding to the monitoring location is greater than or equal to a preset increment threshold, the monitoring location is determined as the third pollution source corresponding to the target river.

9. A water pollution tracing device, characterized in that: include: A problem point determination module is used to determine at least one problem point according to the time series of water quality monitoring values ​​corresponding to each monitoring point in the selected range area of ​​the target river channel, wherein the time series of water quality monitoring values ​​corresponding to any one monitoring point includes the time series of monitoring values ​​corresponding to at least one water quality indicator at the monitoring point; A river section inspection module is used to determine at least one river section to be inspected in the target river channel according to the location of each problem point and / or the water quality analysis sequence corresponding to each problem point, wherein the water quality analysis sequence corresponding to any problem point is determined by time series decomposition based on the water quality monitoring value time series and the flow rate monitoring value time series corresponding to the problem point, and the water quality analysis sequence corresponding to the problem point includes a trend analysis sequence, a periodic analysis sequence and a residual analysis sequence corresponding to the problem point; The first pollution source determination module is used to determine the first pollution source corresponding to the target river channel according to the water quality monitoring value time series and flow monitoring value time series corresponding to each outlet in each river section to be inspected, wherein the water quality monitoring value time series corresponding to any outlet includes the monitoring value time series corresponding to at least one water quality indicator at the outlet.

10. An electronic device, characterized in that: include: processor; a memory for storing processor-executable instructions; The processor is configured to implement the method described in any one of claims 1 to 8 when executing the instructions stored in the memory.

Images