Method and system for monitoring multi-source pollution at estuary based on dynamic traceability

Abstract

The invention provides an estuary multi-source pollution monitoring method and system based on dynamic traceability. The method comprises the following steps: acquiring pollution information of an estuary of a target water area; acquiring pollution information of a pollution point of the target water area; constructing a multi-source pollution convection diffusion model; the pollution information of the estuary of the target water area and the pollution information of the pollution points of the target water area are input into the multi-source pollution convection diffusion model, and the pollution concentration value of each pollution point diffused to the estuary is obtained; according to the pollution concentration value corresponding to each pollution point, obtaining the pollution contribution degree of each pollution point to the estuary; and according to the pollution contribution degree corresponding to each pollution point, in combination with a preset water area monitoring criterion, carrying out pollution point monitoring on the target water area. Therefore, by applying the method provided by the invention, the traceability, monitoring and treatment efficiency of the estuary pollution can be improved.

Description

Technical Field

[0001] The present application relates to the field of water pollution monitoring, and in particular to a method and system for monitoring multi-source pollution at estuaries based on dynamic tracing. Background Art

[0002] With the intensification of coastal economic activities, estuaries have become key nodes for water pollution control due to the interaction between concentrated pollutant emissions from pollution sources and ocean dynamic conditions.

[0003] Traditional monitoring methods are mostly based on static diffusion models, which infer pollutant migration by fixing the coordinates of pollution sources. However, they ignore the spatiotemporal dynamics of multi-source emissions (such as moving ships and seasonal runoff), causing the tracing results to deviate from the actual contribution ratio. At the same time, the global high-resolution grid calculations they rely on are inefficient under complex hydrodynamic conditions and cannot meet real-time monitoring needs. Moreover, they fail to effectively solve problems such as dynamic source coordinate drift and multi-source data fusion lag, resulting in insufficient adaptability of the model under the influence of tides and ocean currents, which in turn leads to inefficient tracing, monitoring and control of estuary pollution. Summary of the Invention

[0004] Based on this, the purpose of this application is to provide a method and system for monitoring multi-source pollution at estuaries based on dynamic tracing, which can effectively improve the efficiency of tracing, monitoring and control of estuary pollution.

[0005] The purpose of this application can be achieved through the following technical solutions:

[0006] A method for monitoring multi-source pollution at an estuary based on dynamic tracing, comprising the following steps: obtaining pollution information at the estuary of a target water area, wherein the pollution information at the estuary at least includes the geographical coordinate information of the estuary and the pollutant concentration at the estuary; obtaining pollution information at pollution points in the target water area, wherein the pollution information at the pollution points at least includes the geographical coordinate information of several pollution points and corresponding pollution emission intensity information; constructing a multi-source pollution convection diffusion model; inputting the pollution information at the estuary of the target water area and the pollution information of the pollution points in the target water area into the multi-source pollution convection diffusion model, and obtaining pollution concentration values of each of the pollution points diffused to the estuary; obtaining the pollution contribution of each of the pollution points to the pollution at the estuary according to the pollution concentration values corresponding to each of the pollution points; and conducting pollution point monitoring on the target water area according to the pollution contribution corresponding to each of the pollution points in combination with preset water area monitoring criteria.

[0007] A multi-source pollution monitoring system for an estuary based on dynamic tracing, wherein the multi-source pollution monitoring system for an estuary based on dynamic tracing comprises: an estuary pollution information acquisition unit, used to acquire pollution information of the estuary of a target water area, wherein the pollution information of the estuary at least includes the geographical coordinate information of the estuary and the concentration of pollutants at the estuary; a pollution point pollution information acquisition unit, used to acquire pollution information of pollution points in the target water area, wherein the pollution information of the pollution points at least includes the geographical coordinate information of several pollution points and the corresponding pollution emission intensity information; a multi-source pollution convection diffusion model construction unit, used to A multi-source pollution convection diffusion model is constructed; a pollution concentration value calculation unit is used to input the pollution information of the estuary of the target water area and the pollution information of the pollution points in the target water area into the multi-source pollution convection diffusion model, and obtain the pollution concentration value of each pollution point diffused to the estuary; a pollution contribution calculation unit is used to obtain the pollution contribution of each pollution point to the estuary according to the pollution concentration value corresponding to each pollution point; a pollution point monitoring unit is used to monitor the pollution points of the target water area according to the pollution contribution corresponding to each pollution point in combination with the preset water area monitoring criteria.

[0008] Compared with the prior art, the method described in the present application first obtains the pollution information of the estuary of the target water area and the pollution information of the pollution points of the target water area, then constructs a multi-source pollution convection diffusion model, and inputs the pollution information of the estuary of the target water area and the pollution information of the pollution points of the target water area into the multi-source pollution convection diffusion model, obtains the pollution concentration value of each of the pollution points diffused to the estuary, and finally obtains the pollution contribution of each of the pollution points to the pollution at the estuary based on the pollution concentration value corresponding to each of the pollution points, and based on the pollution contribution corresponding to each of the pollution points, combined with the preset water area monitoring criteria, performs pollution point monitoring on the target water area. Therefore, the method described in this application, by obtaining the real-time geographic coordinates and pollution data of the estuary and pollution points, constructs an adaptive multi-source diffusion model, realizes the direct correlation between the emission intensity of the pollution points and the concentration changes at the estuary, and reversely analyzes the dynamic contribution of each pollution point through the concentration value output by the model, and realizes responsibility classification in combination with the preset monitoring criteria, thereby improving the implementation effect of pollution control decisions, and improving the matching efficiency of the monitoring results of estuary pollution and control actions, thereby effectively improving the traceability, monitoring and control efficiency of estuary pollution.

[0009] For better understanding and implementation, the present application is described in detail below with reference to the accompanying drawings. BRIEF DESCRIPTION OF THE DRAWINGS

[0010] Figure 1 A flowchart of the steps of a method for monitoring multi-source pollution at an estuary based on dynamic tracing provided in this application;

[0011] Figure 2 A flowchart of the steps for obtaining pollution information of the estuary of the target waters in a multi-source pollution monitoring method of the estuary based on dynamic tracing provided in this application;

[0012] Figure 3 A flowchart of the steps for obtaining a first observation sensor combination in a multi-source pollution monitoring method for an estuary based on dynamic tracing provided by this application;

[0013] Figure 4 A flowchart of the steps for obtaining pollution concentration values of each pollution point that spreads to the estuary in a multi-source pollution monitoring method at the estuary based on dynamic source tracing provided by this application;

[0014] Figure 5 A flowchart of the steps for obtaining the distance between the estuary of the target water area and the pollution point in the target water area in a multi-source pollution monitoring method of the estuary based on dynamic tracing provided by this application;

[0015] Figure 6 This application provides a structural principle diagram of a multi-source pollution monitoring system for estuaries based on dynamic tracing. DETAILED DESCRIPTION

[0016] This application provides a method and system for monitoring multi-source pollution at an estuary based on dynamic source tracing. To make the purpose, technical solution, and effects of this application more clear and explicit, the application is further described in detail below with reference to the accompanying drawings and examples. It should be understood that the specific examples described herein are intended only to explain this application and are not intended to limit this application.

[0017] It will be understood by those skilled in the art that, unless expressly stated otherwise, the singular forms "a", "an", "said" and "the" used herein may also include the plural forms. It should be further understood that the term "comprising" used in the specification of the present application refers to the presence of the features, integers, steps, operations, elements and / or components, but does not exclude the presence or addition of one or more other features, integers, steps, operations, elements, components and / or groups thereof. It should be understood that when we refer to an element as being "connected" or "coupled" to another element, it may be directly connected or coupled to the other element, or there may be intermediate elements. In addition, "connected" or "coupled" as used herein may include wireless connections or wireless couplings. The term "and / or" used herein includes all or any units and all combinations of one or more associated listed items.

[0018] It will be understood by those skilled in the art that, unless otherwise defined, all terms (including technical and scientific terms) used herein have the same meaning as commonly understood by those skilled in the art to which this application belongs. It should also be understood that terms such as those defined in common dictionaries should be understood to have meanings consistent with their meanings in the context of the prior art and will not be interpreted in an idealized or overly formal sense unless specifically defined as herein.

[0019] The invention will be further explained below through description of embodiments in conjunction with the accompanying drawings.

[0020] Example 1

[0021] Please refer to Figure 1 , Figure 1 This application provides a flowchart of a method for monitoring multi-source pollution at an estuary based on dynamic tracing. The method comprises the following steps:

[0022] S10, obtaining pollution information of the estuary of the target water area;

[0023] S20, obtaining pollution information of pollution points in the target water area;

[0024] S30. Construct a multi-source pollution convection and diffusion model;

[0025] S40, inputting the pollution information of the estuary of the target waters and the pollution information of the pollution points of the target waters into the multi-source pollution convection diffusion model, and obtaining the pollution concentration value of each pollution point that diffuses to the estuary;

[0026] S50, obtaining the pollution contribution of each pollution point to the estuary according to the pollution concentration value corresponding to each pollution point;

[0027] S60: Monitoring pollution points in the target water area according to the pollution contribution corresponding to each pollution point and in combination with preset water area monitoring criteria.

[0028] Compared with the existing technology, this application has significantly improved the accuracy and operability of estuary pollution monitoring by systematically integrating dynamic pollution source tracking and multi-source contribution grading mechanism. In the technical solution of this application, by dynamically obtaining the real-time geographic coordinates and pollution data of the estuary and pollution points, an adaptive multi-source diffusion model is constructed to directly link the emission intensity of the pollution point with the concentration change at the estuary, thereby avoiding the traceability distortion problem caused by the data update lag in the traditional model. In addition, in order to solve the problem that the contribution of multi-source pollution is difficult to quantify, the existing technology mostly adopts a single concentration threshold or empirical weight allocation, which cannot distinguish the actual impact ratio of different pollution sources (such as industrial emissions, agricultural runoff, etc.) in complex scenarios; in this application, the dynamic contribution of each pollution point is reversely analyzed through the concentration value output by the model, and the responsibility classification is realized in combination with the preset monitoring criteria, so that the pollution control decision-making is transformed from "fuzzy attribution" to "precise positioning". Therefore, the technical solution of this application improves the efficiency of tracing, monitoring and governance of estuary pollution through dynamic process monitoring and data-driven mechanisms, and improves the matching efficiency of estuary pollution monitoring results and governance actions, providing a feasible technical path for the monitoring, prevention and control and governance of estuary pollution.

[0029] In step S10, pollution information of the estuary of the target water area is obtained.

[0030] Among them, the target waters refer to waters with an estuary used to implement the method described in this application; the pollution information of the estuary includes at least the geographical coordinate information of the estuary and the pollutant concentration at the estuary, such as the latitude and longitude coordinates of the estuary and the concentration of pollutants such as lead, cadmium or crude oil at the estuary.

[0031] Please refer to Figure 2 In one embodiment, step S10 includes the following steps:

[0032] S101. Acquire several different sensors from the sensor pool to construct a target water area observation sensor combination.

[0033] The sensor pool includes a plurality of sensors for observing ground objects, and the target water area observation sensor combination includes at least two different sensors for observing pollution information of the target water area.

[0034] S102, calculating the subject consistency T of the target water area observation sensor combination according to the subject consistency calculation formula of the target water area observation sensor combination r :

[0035]

[0036] Among them, T iis the subject matching degree of the i-th sensor in the target water area observation sensor combination, n is the number of sensors in the target water area observation sensor combination, and,

[0037]

[0038] In one embodiment, if the water body of the estuary is polluted by the proliferation of green algae, the chlorophyll content in the water body will increase significantly, and the chlorophyll a absorption peak of 440nm needs to match the sensor blue light band. For example, the characteristic band of chlorophyll a is 440nm. If the sensor covers 430-450nm, then T i =100%; if sensor A covers 80% of the characteristic band of chlorophyll a and sensor B covers 60%, then the subject consistency of the target water area observation sensor combination is T r =(80%+60%) / 2=70%.

[0039] S103, calculating the spatiotemporal coverage C of the target water area observation sensor combination according to the spatiotemporal coverage calculation formula of the target water area observation sensor combination v :

[0040]

[0041] Among them, S i is the coverage area of the target water area by the i-th sensor in the target water area observation sensor combination within the preset time window, R is the spatial range of the target water area; ∪ is the union of time and space, that is, coverage at any time at the same location is counted as coverage; ∩ is the spatial intersection, that is, only the coverage within the target water area is calculated.

[0042] In this embodiment, S i Can be visualized as a space-time cube.

[0043] In one embodiment, technicians can use Google Earth Engine or ArcGIS Pro to perform spatiotemporal raster operations to generate coverage heat maps to calculate the spatiotemporal coverage. For example, the study area is 100 km 2 The target water area observation sensor combination covers 80km 2 , then C v =80%.

[0044] S104: Calculate the overall accuracy suitability R of the target water area observation sensor combination according to the overall accuracy suitability calculation formula of the target water area observation sensor combination. s :

[0045]

[0046] Among them, p(a i ) is the accuracy suitability of the i-th sensor in the target water area observation sensor combination, that is, the inverse of the resolution of the sensor; m i is the grid number of the resolution corresponding to the i-th sensor in the target water area observation sensor combination; M is the total number of grids in the target water area; and n1 is the number of resolution types.

[0047] In one embodiment, p(a i ) is a score related to the resolution of the sensor. In this embodiment, it can be directly equivalent to the inverse of the resolution of the sensor.

[0048] In one embodiment, a resolution of 10 m corresponds to 100 grids; when the sensor has three resolutions of 10 m, 30 m, and 50 m, n1=3.

[0049] In one embodiment, a weighted average of the accuracy suitability is calculated by taking the proportion of the number of grids of sensors with different resolutions as a weight.

[0050] In one embodiment, if the total number of grids in the target water area is M=1000, there are sensor A: resolution 10m, coverage grid number m1=500; sensor B: resolution 30m, coverage grid number m2=300; sensor C: resolution 50m, coverage grid number m3=200; R S =0.064.

[0051] S105. According to preset evaluation criteria, combined with the subject consistency of the target water area observation sensor combination, the spatiotemporal coverage of the target water area observation sensor combination, and the overall accuracy suitability of the target water area observation sensor combination, the overall observation capability satisfaction evaluation of the target water area observation sensor combination is performed to obtain a first observation sensor combination.

[0052] Among them, the first observation sensor combination is the target water area observation sensor combination with the highest overall observation capability satisfaction.

[0053] Please refer to Figure 3 In one embodiment, step S105 includes the following steps:

[0054] S1051. Calculate the overall observation capability satisfaction OCSI of the target water area observation sensor combination according to the calculation formula for the overall observation capability satisfaction of the target water area observation sensor combination:

[0055]

[0056] Among them, C iis the coverage area of the target water area by the i-th sensor in the target water area observation sensor combination within the preset time window, C a is the total area of the target water area.

[0057] S1052: Set the target water area observation sensor combination with the highest overall observation capability satisfaction as the first observation sensor combination.

[0058] In one embodiment, when the method described in this application is applied to the multi-source pollution monitoring at the estuary of a river tributary, C a =100km 2 , there is sensor A (satellite) with a resolution of 10m, C i =80km 2 , T i =90%, p(a i )=0.1; sensor B (drone), resolution is 0.1m, C i =20km 2 , T i =100%, p(a i )=10; sensor C (mobile vehicle), resolution is 1m, C i =5km 2 , T i =80%, p(a i )=1, then the overall observation capability satisfaction OCSI of the target water area observation sensor combination is 211.2%.

[0059] In some embodiments, the OCSI needs to be constrained between 0% and 100%. A value exceeding 100% indicates that the combined capability of the sensors far exceeds the requirement, and technicians can normalize it to 100%.

[0060] S106: Use the first observation sensor combination to observe the estuary of the target water area, obtain observation data of the estuary of the target water area, and acquire pollution information of the estuary of the target water area.

[0061] The observation data of the estuary of the target water area includes first observation data and second observation data.

[0062] In one embodiment, step S106 includes the following steps:

[0063] S1061: Use the first observation sensor combination to observe the estuary of the target water area at several consecutive time nodes to obtain first observation data.

[0064] In this embodiment, the first observation sensor combination is used to observe the estuary of the target water area at several consecutive time nodes to obtain the time stamps (t1, t2, ..., t m ) sensor time series data, that is, the first observation data.

[0065] S1062: Set a sliding window for data fusion and assign a Gaussian weight w(t j ):

[0066]

[0067] The sliding window is a time period centered on the transit time of the remote sensing satellite over the target waters, t j is the observation timestamp of the first observation data, t sat is the transit time of the remote sensing satellite over the target water area, σ is the control weight attenuation speed, and σ=Δt / 3, Δt is the width of the preset sliding window, and e is the natural logarithm.

[0068] In this embodiment, the transit time t of the remote sensing satellite to the target water area is used. sat As the center, set the window width of the dynamic window △t, for example, ±6 hours; at the same time, use the Gaussian decay weight function, from t sat The more recent the data, the higher the weight.

[0069] S1063: Fusing the first observation data in the sliding window according to the Gaussian weight to generate second observation data Dn(t j ):

[0070]

[0071] Among them, D(t j ) is at time t j The first observation data observed at the time is n, and n is the number of the first observation data in the sliding window.

[0072] In one embodiment, t sat The time is 10:00, the window width Δt = 6 hours, and the first observation data measured by the sensor combination at 8:00 (Gaussian weight 0.8), 12:00 (Gaussian weight 0.6), and 14:00 (Gaussian weight 0.3) are 20, 25, and 30. The second observation data obtained is Dn(t j )=23.2.

[0073] In addition, this application also provides some steps for monitoring multi-source pollution at estuaries. When the steps are applied to the multi-source pollution monitoring method at estuaries, abnormal observation data can be eliminated. The steps include:

[0074] S107, obtaining historical data of observation data of the estuary of the target water area;

[0075] S108, constructing an isolation forest using the historical data, and setting the tree upper limit of the isolation forest to a preset upper limit value;

[0076] In one embodiment, the preset upper limit value is 100.

[0077] S109, calculating an anomaly score of the observation data of the estuary of the target water area based on the isolated forest and the preset upper limit;

[0078] S110: If the abnormality score of the observation data of the estuary of the target water area is greater than a preset value, the observation data of the estuary of the target water area is eliminated.

[0079] In one embodiment, if the anomaly score of the observation data of the estuary of the target water area is greater than 0.65, the observation data of the estuary of the target water area is determined to be an outlier and is removed.

[0080] In step S20, pollution information of pollution points in the target water area is obtained.

[0081] The pollution information of the pollution points includes at least the geographical coordinate information of several pollution points and the corresponding pollution emission intensity information, such as the latitude and longitude coordinates of the pollution points and the concentration of pollutants such as lead, cadmium or crude oil at the pollution points.

[0082] In one embodiment, a technician may make adaptive modifications to the technical solutions described in steps S101-S110 and apply them to step S20 to achieve the same or similar technical effects.

[0083] In one embodiment, the target water area includes at least one pollution point. It can be understood that when the number of pollution points in the target water area is any number, the method described in this application can be applied.

[0084] For step S30, a multi-source pollution convection diffusion model is constructed.

[0085] In one embodiment, step S30 includes the following steps:

[0086] S301: Acquire water area information of the target water area.

[0087] The water area information of the target water area at least includes the diffusion coefficient and water flow velocity field of the target water area;

[0088] S302: Construct a multi-source pollution convection diffusion model based on the water area information of the target water area:

[0089]

[0090] Among them, C is the pollutant concentration field of the target water area, t is time, D is the diffusion coefficient of the target water area, u is the water flow velocity field of the target water area, and S is the amount of pollutants released per second by the pollution point of the target water area.

[0091] In this embodiment, C is a spatially continuous function that describes the concentration distribution of pollutants within the calculation domain (including the vicinity of the pollution source, the propagation path, the estuary, etc.). For example, C may be higher at the pollution source; at the estuary, C is the comprehensive result of all upstream pollution sources after water flow and diffusion.

[0092] In this embodiment, the value of C at the estuary is the value of the multi-source pollution convection diffusion model at a specific location (x 入海口 ,y 入海口 ) solution.

[0093] In step S40, the pollution information of the estuary of the target water area and the pollution information of the pollution points of the target water area are input into the multi-source pollution convection diffusion model to obtain the pollution concentration value of each pollution point that diffuses to the estuary.

[0094] Please refer to Figure 4 In one embodiment, step S40 includes the following steps:

[0095] S401, obtaining the distance between the estuary of the target water area and the pollution point of the target water area;

[0096] In one embodiment, the pollution points of the target waters include point source pollution, line source pollution and non-point source pollution; and, please refer to Figure 5 , the step S401 includes the following steps:

[0097] S4011. When the pollution point of the target water area is point source pollution, obtain the Euclidean distance between the estuary of the target water area and the point source pollution.

[0098] In one embodiment, when the pollution point of the target water area is point source pollution, S=Qδ(x-x0,y-y0), where Q is the emission rate of the point source pollution, and δ is the Dirac function, indicating that there is emission only at (x0,y0).

[0099] S4012: When the pollution point of the target water area is line source pollution, obtain the Euclidean distance between the estuary of the target water area and the end point of the line source pollution.

[0100] In other embodiments, when the pollution point of the target water area is a line source pollution, it can also be regarded as the pollutant discharged along a continuous path (such as a river and an oil pipeline), that is, x and y are regarded as the path coordinate set of the line source. For example, the coordinates of the center line of a tributary river are {(x1, y1), (x2, y2), ..., (x n ,y n )},

[0101]

[0102] Where L is the total length of the line source, Qtotal is the total emission rate, and it is discretized into multiple point source emissions.

[0103] S4013: When the pollution point of the target water area is non-point source pollution, the pollution point of the target water area is rasterized and converted into distributed point source pollution to obtain distance.

[0104] In other embodiments, when the pollution point of the target water area is non-point source pollution, it can be regarded as the pollutant being discharged widely in an area (such as farmland or urban surface runoff), that is, x and y are regarded as the boundary range of the non-point source or the coordinates of the polygon vertices. For example, the boundary coordinates of a farmland area are {(x1, y1), (x2, y2), ..., (x n ,y n )}, forming a polygon,

[0105] (When (x,y)∈surface source area),

[0106] Where A is the area of the source, and Qtotal is the total emission rate, evenly distributed over each unit area.

[0107] S402: Input the pollution information of the estuary of the target waters and the pollution information of the pollution points of the target waters into the multi-source pollution convection diffusion model, and calculate the pollution concentration value C of each pollution point at the estuary based on the distance. i ':

[0108]

[0109] Wherein, Q is the river flow of the target water area, and x is the distance.

[0110] In these embodiments, x and y represent the position of any point within the target water area, which can generally correspond to a geographic coordinate system (such as longitude and latitude) or a plane projection coordinate system (such as metric coordinates in the east-west and north-south directions); x0 and y0 are the precise location coordinates of the pollution source (such as a sewage outlet or a factory chimney).

[0111] In step S50, the pollution contribution of each pollution point to the estuary is obtained according to the pollution concentration value corresponding to each pollution point.

[0112] In one embodiment, step S50 includes the following steps:

[0113] S501. According to the pollution concentration value corresponding to each pollution point and the pollution contribution calculation formula, the pollution contribution ηi of each pollution point to the estuary is obtained:

[0114]

[0115] Among them, C i is the pollution concentration value diffused from the i-th pollution point to the estuary, and n2 is the number of pollution points.

[0116] In step S60, pollution point monitoring is performed on the target water area according to the pollution contribution corresponding to each pollution point and in combination with the preset water area monitoring criteria.

[0117] In one embodiment, technicians can adopt the following practices to monitor pollution points in target waters: first, hierarchical monitoring, listing pollution points with contributions higher than a preset threshold (such as 20%) as key monitoring targets, deploying online sensors to monitor emission intensity and water quality changes in real time, and regular inspections of low-contribution points; second, path tracking, adding mobile monitoring buoys on key diffusion paths between estuaries and high-contribution pollution points (such as river confluences and mainstream ocean currents) to capture the migration dynamics of pollution clusters; third, dynamic response, when the contribution of a pollution point suddenly increases in a short period of time (such as exceeding a threshold of 50%), automatically triggering drone patrols and satellite remote sensing collaborative verification to quickly locate abnormal emission sources; fourth, governance linkage, linking contribution rankings with pollutant discharge permit review and environmental law enforcement priorities, for example, mandatory installation of smart sewage gates for pollution points that have contributed to the top 10% for three consecutive months to achieve closed-loop management of pollution control; in addition, the layout of monitoring points can be optimized based on historical contribution data, redundant monitoring equipment with long-term low contributions can be eliminated, and operation and maintenance costs can be reduced.

[0118] Example 2

[0119] Please refer to Figure 6The present application also provides a multi-source pollution monitoring system for an estuary based on dynamic source tracing, to implement the steps of the multi-source pollution monitoring method for an estuary based on dynamic source tracing in the above embodiment. The multi-source pollution monitoring system for an estuary based on dynamic source tracing includes: an estuary pollution information acquisition unit 1001, a pollution point pollution information acquisition unit 1002, a multi-source pollution convection and diffusion model construction unit 1003, a pollution concentration value calculation unit 1004, a pollution contribution calculation unit 1005, and a pollution point monitoring unit 1006.

[0120] The estuary pollution information acquisition unit 1001 is used to acquire pollution information of the estuary of the target water area, wherein the estuary pollution information at least includes geographical coordinate information of the estuary and pollutant concentrations of the estuary;

[0121] The pollution point pollution information acquisition unit 1002 is used to acquire pollution information of pollution points in the target water area, wherein the pollution information of the pollution points at least includes geographical coordinate information of several pollution points and corresponding pollution emission intensity information;

[0122] The multi-source pollution convection diffusion model construction unit 1003 is used to construct a multi-source pollution convection diffusion model;

[0123] The pollution concentration value calculation unit 1004 is used to input the pollution information of the estuary of the target water area and the pollution information of the pollution points of the target water area into the multi-source pollution convection diffusion model to obtain the pollution concentration value of each pollution point that diffuses to the estuary;

[0124] The pollution contribution calculation unit 1005 is used to obtain the pollution contribution of each pollution point to the estuary according to the pollution concentration value corresponding to each pollution point;

[0125] The pollution point monitoring unit 1006 is used to monitor the pollution points of the target water area according to the pollution contribution corresponding to each pollution point in combination with the preset water area monitoring criteria.

[0126] It should be noted that the above embodiment provides a multi-source pollution monitoring system for an estuary based on dynamic tracing. When implementing a multi-source pollution monitoring method for an estuary based on dynamic tracing, the division of the above-mentioned functional modules is only used as an example. In actual applications, the above-mentioned functions can be assigned to different functional modules as needed, that is, the internal structure of the equipment can be divided into different functional modules to complete all or part of the functions described above.

[0127] In addition, the multi-source pollution monitoring system for an estuary based on dynamic tracing provided in the above embodiment and the multi-source pollution monitoring method for an estuary based on dynamic tracing provided in Example 1 belong to the same concept. The implementation process thereof is detailed in the method embodiment, i.e., Example 1, and will not be repeated here.

[0128] The above-described embodiments merely represent several implementation methods of the present application. While the descriptions are relatively specific and detailed, they should not be construed as limiting the scope of the present invention. It should be noted that a person skilled in the art may make various modifications and improvements without departing from the spirit of the present application, and the present application is intended to encompass such modifications and variations.

Claims

1. A method for monitoring multi-source pollution at an estuary based on dynamic source tracing, comprising the following steps: Obtaining pollution information at the estuary of the target waters, wherein the pollution information at the estuary includes at least geographical coordinate information of the estuary and pollutant concentrations at the estuary; Obtaining pollution information of pollution points in the target waters, wherein the pollution information of the pollution points includes at least geographic coordinate information of several pollution points and corresponding pollution emission intensity information; Construct a multi-source pollution convection diffusion model; Inputting the pollution information of the estuary of the target waters and the pollution information of the pollution points of the target waters into the multi-source pollution convection diffusion model to obtain the pollution concentration value of each pollution point that diffuses to the estuary; Obtaining the pollution contribution of each pollution point to the estuary according to the pollution concentration value corresponding to each pollution point; According to the pollution contribution corresponding to each of the pollution points, combined with the preset water area monitoring criteria, pollution point monitoring is carried out on the target water area.

2. According to the method for monitoring multi-source pollution at an estuary based on dynamic source tracing according to claim 1, the step of obtaining pollution information at the estuary of the target water area comprises: Acquire a plurality of different sensors from a sensor pool to construct a target water area observation sensor combination, wherein the sensor pool includes a plurality of sensors for observing ground objects, and the target water area observation sensor combination includes at least two different sensors for observing pollution information of the target water area; According to the subject consistency calculation formula of the target water area observation sensor combination, the subject consistency T of the target water area observation sensor combination is calculated. r : Among them, T i is the subject matching degree of the i-th sensor in the target water area observation sensor combination, n is the number of sensors in the target water area observation sensor combination, and, According to the calculation formula of the spatiotemporal coverage of the target water area observation sensor combination, the spatiotemporal coverage C of the target water area observation sensor combination is calculated. v : Among them, S i is the coverage area of the target water area by the i-th sensor in the target water area observation sensor combination within the preset time window, R is the spatial range of the target water area; ∪ is the spatiotemporal union, that is, coverage at any time at the same location is counted as coverage; ∩ is the spatial intersection, that is, only coverage within the target water area is counted; According to the calculation formula of the overall accuracy suitability of the target water area observation sensor combination, the overall accuracy suitability R of the target water area observation sensor combination is calculated. s : Among them, p(a i ) is the accuracy suitability of the i-th sensor in the target water area observation sensor combination, that is, the inverse of the resolution of the sensor; m i is the grid number of the resolution corresponding to the i-th sensor in the target water area observation sensor combination; M is the total number of grids in the target water area; n1 is the number of resolution types; According to the preset evaluation criteria, combined with the subject consistency of the target water area observation sensor combination, the spatiotemporal coverage of the target water area observation sensor combination, and the overall accuracy suitability of the target water area observation sensor combination, the overall observation capability satisfaction evaluation of the target water area observation sensor combination is performed to obtain a first observation sensor combination; The first observation sensor combination is used to observe the estuary of the target water area, obtain observation data of the estuary of the target water area, and acquire pollution information of the estuary of the target water area.

3. According to the method for monitoring multi-source pollution in estuaries based on dynamic source tracing according to claim 2, the step of performing an overall observation capability satisfaction evaluation on the target water area observation sensor combination and obtaining the first observation sensor combination comprises: According to the calculation formula of the overall observation capability satisfaction of the target water area observation sensor combination, the overall observation capability satisfaction OCSI of the target water area observation sensor combination is calculated: Among them, C i is the coverage area of the target water area by the i-th sensor in the target water area observation sensor combination within the preset time window, C a is the total area of the target waters; The target water area observation sensor combination with the highest overall observation capability satisfaction is set as the first observation sensor combination.

4. The method for monitoring multi-source pollution in estuaries based on dynamic source tracing according to claim 2 is characterized in that: The observation data of the estuary of the target waters includes first observation data and second observation data; The step of using the first observation sensor combination to observe the estuary of the target water area to obtain observation data of the estuary of the target water area includes: Using the first observation sensor combination to observe the estuary of the target water area at a plurality of consecutive time nodes to obtain first observation data; Set a sliding window for data fusion and assign Gaussian weights w(t j ): The sliding window is a time period centered on the transit time of the remote sensing satellite over the target waters, t j is the observation timestamp of the first observation data, t sat is the transit time of the remote sensing satellite over the target water area, σ is the control weight decay rate, and σ=Δt / 3, Δt is the width of the preset sliding window, and e is the natural logarithm; The first observation data in the sliding window are fused according to the Gaussian weights to generate the second observation data Dn(t j ): Among them, D(t j ) is at time t j The first observation data observed at the time is n, and n is the number of the first observation data in the sliding window.

5. The method for monitoring multi-source pollution at an estuary based on dynamic source tracing according to any one of claims 2 to 4, further comprising: Obtaining historical data of observation data of the estuary of the target water area; constructing an isolation forest using the historical data, and setting the tree upper limit of the isolation forest to a preset upper limit value; Calculating an anomaly score of the observation data of the estuary of the target water area according to the isolation forest and the preset upper limit value; If the abnormal score of the observation data of the estuary of the target water area is greater than a preset value, the observation data of the estuary of the target water area is eliminated.

6. According to the method for monitoring multi-source pollution at an estuary based on dynamic source tracing according to claim 1, the step of constructing a multi-source pollution convection and diffusion model comprises: Acquiring water area information of the target water area, wherein the water area information of the target water area at least includes a diffusion coefficient and a water flow velocity field of the target water area; Based on the water area information of the target water area, a multi-source pollution convection diffusion model is constructed: Among them, C is the pollutant concentration field of the target water area, t is time, D is the diffusion coefficient of the target water area, u is the water flow velocity field of the target water area, and S is the amount of pollutants released per second by the pollution point of the target water area.

7. According to the method for monitoring multi-source pollution at an estuary based on dynamic source tracing according to claim 6, the steps of inputting pollution information at the estuary of the target waters and pollution information of pollution points in the target waters into the multi-source pollution convection diffusion model, and obtaining pollution concentration values of each pollution point that diffuses to the estuary include: Obtaining the distance between the estuary of the target water area and the pollution point in the target water area; The pollution information of the estuary of the target waters and the pollution information of the pollution points of the target waters are input into the multi-source pollution convection diffusion model, and the pollution concentration value C of each pollution point diffused to the estuary is calculated based on the distance. i ': Wherein, Q is the river flow of the target water area, and x is the distance.

8. The method for monitoring multi-source pollution in estuaries based on dynamic source tracing according to claim 7 is characterized by: The pollution points of the target waters include point source pollution, line source pollution and non-point source pollution; The step of obtaining the distance between the estuary of the target water area and the pollution point of the target water area comprises: When the pollution point of the target water area is point source pollution, obtaining the Euclidean distance between the estuary of the target water area and the point source pollution; When the pollution point of the target water area is line source pollution, obtaining the Euclidean distance between the estuary of the target water area and the end point of the line source pollution; When the pollution point of the target water area is non-point source pollution, the pollution point of the target water area is rasterized and converted into distributed point source pollution for distance acquisition.

9. The method for monitoring multi-source pollution at an estuary based on dynamic source tracing according to any one of claims 7 or 8, wherein the step of obtaining the pollution contribution of each pollution point to the pollution at the estuary based on the pollution concentration value corresponding to each pollution point comprises: According to the pollution concentration value corresponding to each pollution point, combined with the pollution contribution calculation formula, the pollution contribution ηi of each pollution point to the estuary is obtained: Among them, C i is the pollution concentration value diffused from the i-th pollution point to the estuary, and n2 is the number of pollution points.

10. A multi-source pollution monitoring system for estuaries based on dynamic tracing, characterized in that: The multi-source pollution monitoring system for estuaries based on dynamic source tracing includes: An estuary pollution information acquisition unit, configured to acquire pollution information of an estuary of a target water area, wherein the estuary pollution information includes at least geographical coordinate information of the estuary and pollutant concentrations of the estuary; A pollution point pollution information acquisition unit, configured to acquire pollution information of pollution points in a target water area, wherein the pollution information of the pollution points at least includes geographical coordinate information of a plurality of pollution points and corresponding pollution emission intensity information; Multi-source pollution convection diffusion model construction unit, used to construct a multi-source pollution convection diffusion model; a pollution concentration value calculation unit, configured to input the pollution information of the estuary of the target waters and the pollution information of the pollution points of the target waters into the multi-source pollution convection diffusion model, and obtain the pollution concentration value of each of the pollution points that diffuses to the estuary; a pollution contribution calculation unit, configured to obtain the pollution contribution of each pollution point to the estuary according to the pollution concentration value corresponding to each pollution point; The pollution point monitoring unit is used to monitor the pollution points of the target water area according to the pollution contribution corresponding to each pollution point and in combination with the preset water area monitoring criteria.

Images