A pollution center tracing method and system based on air quality monitoring network

Abstract

The present invention provides a pollution center source tracing method and system based on an air quality monitoring network, and relates to the technical field of atmospheric monitoring. The pollution center source tracing method mainly lags nearby monitoring stations based on the monitoring sequence of the monitoring station, forms a direction vector at the monitoring station according to the causal test result between the lagged monitoring sequence of the nearby monitoring station and the monitoring sequence of the monitoring station, and extends in the direction of the direction vector, thereby tracing the pollutant center. Not only can the data of the existing atmospheric monitoring network be reasonably and effectively utilized, but the data of the atmospheric monitoring network can also be fully combined with the information of the monitoring station to effectively infer the location of the pollution source.

Description

Technical Field

[0001] The present invention relates to the technical field of atmosphere monitoring, and in particular to a pollution center source tracing method and system based on an air quality monitoring network. Background Art

[0002] The existing air quality monitoring system is mainly based on the monitoring values ​​of fixed national and provincial control stations. However, the construction cost of fixed stations such as national and provincial control stations is high and their number is small. They are unable to achieve full-area monitoring and local pollution monitoring. They can only perceive the pollution situation within the monitoring range of the equipment, but cannot perceive the source of the pollution.

[0003] In view of the above problems, the present invention utilizes an air quality monitoring network built with existing basic equipment and data mining technology to achieve tracing of pollution sources. Summary of the invention

[0004] The purpose of the present invention is to provide a pollution center tracing method and system based on an air quality monitoring network. The pollution center tracing method mainly lags the nearby monitoring stations based on the monitoring sequence of the monitoring stations, forms a direction vector at the monitoring station according to the causal test result between the lagged monitoring sequence of the nearby monitoring stations and the monitoring sequence of the monitoring station, and extends in the direction of the direction vector, thereby tracing the pollutant center. Not only can the data of the existing atmospheric monitoring network be reasonably and effectively utilized, but the data of the atmospheric monitoring network can also be fully combined with the monitoring station information to effectively infer the location of the pollution source.

[0005] In order to solve the above technical problems, the present invention adopts the following solutions:

[0006] A pollution center tracing method based on an air quality monitoring network, the method comprising the following steps:

[0007] S1: Real-time detection of whether the pollutant monitoring value of each monitoring station exceeds the preset peak value. If so, the monitoring station is used as the monitoring station to be traced, and the process goes to step S2;

[0008] S2: Obtain the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced according to the current time, and find the monitoring site to be inspected near the monitoring site to be traced on the map according to the latitude and longitude information of the monitoring site to be traced, and obtain the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be inspected;

[0009] Among them, the monitoring sequence {t1, t2, …, tn} and the lag monitoring sequence {t1-j, t2-j, …, tn-j} are composed of pollutant monitoring values ​​sorted in chronological order, n represents the current time, and j represents the preset lag order;

[0010] S3: Obtain the causal test results based on the monitoring sequence {t1-j, t2-j, …, tn-j} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j}, and connect them on the map according to the causal test results to form a direction vector between the monitoring site to be traced and the monitoring site to be tested; obtain the direction vectors between other monitoring sites based on the monitoring site to be traced and the monitoring site to be tested; extend different direction vectors until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source.

[0011] Furthermore, in S2, in the monitoring sequence {t1, t2, ..., tn}, tn represents the pollutant monitoring value at the current moment, wherein ti represents the pollutant monitoring value of the monitoring station every hour, i∈[1,n].

[0012] Furthermore, in S2, the preset lag order is related to the site distance between the current monitoring site to be traced and the current monitoring site to be inspected, and the wind speed at the current moment.

[0013] Furthermore, S3 specifically includes the following steps:

[0014] S31: Obtain a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and connect them on a map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causal test result;

[0015] S32: Traverse other monitoring sites near the monitoring site to be traced according to the pollutant monitoring value of the monitoring site to be traced, obtain other monitoring sites that match the pollutant monitoring value as the first auxiliary monitoring site, and obtain the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site according to the current time;

[0016] S33: Find a second auxiliary monitoring site near the first auxiliary monitoring site on the map according to the latitude and longitude information of the first auxiliary monitoring site, and obtain a delayed monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site;

[0017] S34: Obtain a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connect the first auxiliary monitoring site and the second auxiliary monitoring site on a map according to the causal test result to form a direction vector between them;

[0018] S35: Extend the direction vector between the monitoring site to be traced and the monitoring site to be inspected and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source.

[0019] Further, in S32, the process of traversing other monitoring stations near the monitoring station to be traced according to the pollutant monitoring value of the monitoring station to be traced is specifically as follows:

[0020] First, the corresponding pollutant type is obtained according to the pollutant monitoring value of the monitoring site to be traced, and other monitoring sites near the monitoring site to be traced are traversed with the current pollutant type, and it is determined in turn whether the pollutant type corresponding to the pollutant monitoring value in other monitoring sites near the monitoring site to be traced is the same as the current pollutant type. If so, it is determined that the other monitoring sites match the pollutant monitoring value.

[0021] Further, in S32, the process of traversing other monitoring stations near the monitoring station to be traced according to the pollutant monitoring value of the monitoring station to be traced is specifically as follows:

[0022] First, the corresponding pollutant type is obtained according to the pollutant monitoring value of the monitoring site to be traced, and other monitoring sites near the monitoring site to be traced are traversed with the current pollutant type, and the pollutant monitoring values ​​of other monitoring sites near the monitoring site to be traced are judged in turn whether they exceed the preset peak value. If so, it is judged whether the pollutant type corresponding to the pollutant monitoring value is the same as the current pollutant type. If they are the same, it is determined that the other monitoring sites match the pollutant monitoring value.

[0023] Further, in S31, the causal test result is obtained according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and the process of connecting the direction vector between the monitoring site to be traced and the monitoring site to be tested on the map according to the causal test result is specifically as follows:

[0024] A Granger causality test is performed on the monitoring sequence {t1, t2, …, tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the monitoring site to be tested. If the causality test result is less than the preset value, they are connected on the map to form a direction vector from the monitoring site to be tested to the monitoring site to be traced.

[0025] Further, in S34, a causal test result is obtained according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and a process of connecting the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site on a map according to the causal test result is specifically as follows:

[0026] A Granger causality test is performed on the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site. If the causality test result is less than the preset value, they are connected on the map to form a direction vector from the second auxiliary monitoring site to the first auxiliary monitoring site.

[0027] A pollution center tracing system based on an air quality monitoring network, comprising:

[0028] Pollutant monitoring value detection module: detect in real time whether the pollutant monitoring value of each monitoring site exceeds the preset peak value. If so, the monitoring site is used as the monitoring site to be traced and go to step S2;

[0029] Monitoring sequence acquisition module: obtain the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced according to the current time, and find the monitoring site to be inspected near the monitoring site to be traced on the map according to the latitude and longitude information of the monitoring site to be traced, and obtain the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be inspected;

[0030] Among them, the monitoring sequence {t1, t2, …, tn} and the lag monitoring sequence {t1-j, t2-j, …, tn-j} are composed of pollutant monitoring values ​​sorted in chronological order, n represents the current time, and j represents the preset lag order;

[0031] The pollution center source tracing module is based on causal testing: the causal testing results are obtained according to the monitoring sequence {t1-j, t2-j, …, tn-j} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j}, and the direction vector between the monitoring site to be traced and the monitoring site to be tested is connected on the map according to the causal testing results; the direction vectors between other monitoring sites are obtained according to the monitoring site to be traced and the monitoring site to be tested; different direction vectors are extended until they intersect, and the longitude and latitude information of the intersection point is used as the longitude and latitude information of the pollution source.

[0032] Furthermore, the pollution center tracing module based on causal verification includes the following submodules:

[0033] The first causality test module: obtains the causality test result according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and connects them on the map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causality test result;

[0034] The first auxiliary monitoring site data acquisition module: traverses other monitoring sites near the monitoring site to be traced according to the pollutant monitoring value of the monitoring site to be traced, obtains other monitoring sites that match the pollutant monitoring value as the first auxiliary monitoring site, and obtains the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site according to the current time;

[0035] The second auxiliary monitoring site data acquisition module: finds the second auxiliary monitoring site near the first auxiliary monitoring site on the map according to the latitude and longitude information of the first auxiliary monitoring site, and obtains the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site;

[0036] The second causality test module: obtains a causality test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connects the first auxiliary monitoring site and the second auxiliary monitoring site on the map according to the causality test result to form a direction vector;

[0037] Tracing module: Extend the direction vector between the monitoring site to be traced and the monitoring site to be inspected, and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source.

[0038] Beneficial effects of the present invention:

[0039] The present invention provides a pollution center tracing method and system based on an air quality monitoring network. The pollution center tracing method mainly lags the monitoring sequence of a nearby monitoring station based on the monitoring sequence of the monitoring station, and obtains the corresponding change relationship between the two monitoring stations through the causal test results between the lagged monitoring sequence of the nearby monitoring station and the monitoring sequence of the monitoring station, so as to reasonably and effectively use the data of the existing atmospheric monitoring network.

[0040] Moreover, a direction vector is formed according to the corresponding change relationship between the two monitoring stations and the positional relationship between the two monitoring stations, and the direction of the direction vector is extended to trace the pollutant center. This can not only reasonably and effectively utilize the data of the existing atmospheric monitoring network, but also fully combine the data of the atmospheric monitoring network with the monitoring station information to effectively infer the location of the pollution source. BRIEF DESCRIPTION OF THE DRAWINGS

[0041] Figure 1 This is a schematic diagram of the process of the pollution center source tracing method in Example 1 of the present invention;

[0042] Figure 2 Schematic diagram of monitoring stations in the atmosphere monitoring network in Example 1 of the present invention;

[0043] Figure 3 This is a schematic diagram of a structure in which a radius r is expanded outward from the monitoring site to be inspected in Example 1 of the present invention;

[0044] Figure 4 is the direction vector V between the traceability monitoring site A and the monitoring site B to be inspected in Example 1 of the present invention BA Schematic diagram of

[0045] Figure 5 is the direction vector V between the first auxiliary monitoring site C and the second auxiliary monitoring site D in Example 1 of the present invention DC Schematic diagram of

[0046] Figure 6 According to the direction vector V in Embodiment 1 of the present invention BA and direction vector V DC Schematic diagram of pollution sources for tracing pollution centers. DETAILED DESCRIPTION

[0047] The following will be combined with the drawings in the embodiments of the present invention to clearly and completely describe the technical solutions in the embodiments of the present invention. Obviously, the described embodiments are only part of the embodiments of the present invention, rather than all the embodiments. The following description of at least one exemplary embodiment is actually only illustrative and is by no means intended to limit the present invention and its application or use. Based on the embodiments of the present invention, all other embodiments obtained by ordinary technicians in this field without creative work are within the scope of protection of the present invention.

[0048] The relative arrangement of components and steps, the numerical expressions and numerical values ​​set forth in these embodiments do not limit the scope of the present invention unless specifically stated otherwise.

[0049] At the same time, it should be understood that for the convenience of description, the sizes of the various parts shown in the drawings are not drawn according to the actual proportional relationship.

[0050] Additionally, descriptions of well-known structures, functions, and configurations may be omitted for clarity and conciseness.One of ordinary skill in the art will recognize that various changes and modifications may be made to the examples described herein without departing from the spirit and scope of the present disclosure.

[0051] Technologies, methods, and apparatus known to ordinary technicians in the relevant field may not be discussed in detail, but where appropriate, such technologies, methods, and apparatus should be considered part of the authorization specification.

[0052] In all examples shown and discussed herein, any specific values ​​should be interpreted as merely exemplary and not as limiting. Therefore, other examples of the exemplary embodiments may have different values.

[0053] The present invention will be described in detail below with reference to the accompanying drawings and in conjunction with embodiments:

[0054] Example 1

[0055] In this embodiment, if Figure 1 As shown, a pollution center tracing method based on an air quality monitoring network is provided, and the pollution center tracing method comprises the following steps:

[0056] S1: Real-time detection of whether the pollutant monitoring value of each monitoring station exceeds the preset peak value. If so, the monitoring station is used as the monitoring station to be traced, and the process goes to step S2;

[0057] S2: Obtain the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced according to the current time, and find the monitoring site to be inspected near the monitoring site to be traced on the map according to the latitude and longitude information of the monitoring site to be traced, and obtain the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be inspected;

[0058] Among them, the monitoring sequence {t1, t2, …, tn} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} are composed of pollutant monitoring values ​​sorted in chronological order, ti represents the pollutant monitoring value of the monitoring station every hour, i∈[1,n], n represents the current time, and tn represents the pollutant monitoring value at the current time;

[0059] j represents the preset hysteresis order, which is related to the site distance between the current monitoring site to be traced and the current monitoring site to be inspected, and the wind speed at the current moment, and is usually 2 or 3, which can be adjusted according to actual conditions;

[0060] S3: Obtain the causal test results based on the monitoring sequence {t1-j, t2-j, …, tn-j} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j}, and connect them on the map according to the causal test results to form a direction vector between the monitoring site to be traced and the monitoring site to be tested; obtain the direction vectors between other monitoring sites based on the monitoring site to be traced and the monitoring site to be tested; extend different direction vectors until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source.

[0061] Preferably, S3 specifically includes the following steps:

[0062] S31: Obtain a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and connect them on a map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causal test result;

[0063] S32: Traverse other monitoring sites near the monitoring site to be traced according to the pollutant monitoring value of the monitoring site to be traced, obtain other monitoring sites that match the pollutant monitoring value as the first auxiliary monitoring site, and obtain the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site according to the current time;

[0064] S33: Find a second auxiliary monitoring site near the first auxiliary monitoring site on the map according to the latitude and longitude information of the first auxiliary monitoring site, and obtain a delayed monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site;

[0065] S34: Obtain a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connect the first auxiliary monitoring site and the second auxiliary monitoring site on a map according to the causal test result to form a direction vector between them;

[0066] S35: Extend the direction vector between the monitoring site to be traced and the monitoring site to be inspected and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source.

[0067] Based on the above principles, the present invention is further described:

[0068] In this embodiment, the pollutant monitoring values ​​of all national control stations and provincial control stations in the region within a period of time are first obtained. The peak or high value of the pollutant monitoring value can be marked, and the corresponding time and longitude and latitude of the station can be recorded; and the peak or high value of the pollutant is determined by the actual local business and control needs. For example, if PM2.5>50ug / m3, the peak value can be preset in advance based on the actual business and control needs, and the pollutant monitoring value of each monitoring station can be detected in real time to detect whether it exceeds the preset peak value. If so, the current monitoring station is recorded as the monitoring station to be traced, so that the pollutant center can be traced to the monitoring station to be traced.

[0069] Specifically, the atmospheric monitoring network can obtain the site information of the current monitoring site, such as longitude and latitude information, and pollutant monitoring values. Based on the longitude and latitude information, the following can be obtained: Figure 2 The monitoring sites A, B, C, D, E and F are shown, so as to obtain the pollutant monitoring values ​​of the six monitoring sites and perform real-time detection on the pollutant monitoring values ​​of the six monitoring sites.

[0070] When it is detected that the pollutant monitoring value of monitoring station A exceeds the preset peak value, monitoring station A is used as the monitoring station to be traced, and the longitude and latitude information of the monitoring station to be traced can be obtained (X A ,Y A ) and monitoring sequence {t1, t2, …, tn}, where ti represents the pollutant monitoring value of the monitoring station every hour, i∈[1,n]; n represents the current moment, tn represents the pollutant monitoring value of the current moment; n-1 represents the previous hour, tn-1 represents the pollutant monitoring value of the previous hour, and so on. n can be 24, that is, the hourly pollutant monitoring value of the monitoring station A at 24 moments can be taken, and it can also be adjusted according to the data collection situation.

[0071] Next, according to the latitude and longitude information of the monitoring site to be traced (X A ,Y A ) Find the monitoring site to be inspected near the monitoring site to be traced, and a distance value r can be preset for the longitude and latitude information. Take the monitoring site to be traced as the center point and expand the radius r outward to find the monitoring site within the circle. It can be used as the monitoring site to be inspected, such as Figure 3 As shown, with monitoring site A as the center point, the radius r is expanded outward to find monitoring site B located within the circle, and monitoring site B is used as the monitoring site to be tested.

[0072] There may be more than one monitoring site to be tested near the monitoring site to be traced, and there may be multiple monitoring sites to be tested. The lag monitoring sequence {t1-j, t2-j, ..., tn-j} of each monitoring site to be tested is obtained respectively, j∈N+, j represents the preset lag order, wherein the preset lag order is related to the site distance between the current monitoring site to be traced and the current monitoring site to be tested, and the wind speed at the current moment, and is usually 2 or 3, which can be adjusted according to actual conditions. Then, the monitoring site to be traced is sequentially subjected to causal tests with each monitoring site to be tested. In this embodiment, the Granger causality test is adopted, and the specific process is as follows:

[0073] When Granger causality test is performed on the monitoring sequence {t1, t2, ..., tn} of the traceability monitoring site A and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site B to be tested, the monitoring sequence {t1, t2, ..., tn} is used as the sequence T A , the lagged monitoring sequence {t1-j, t2-j, …, tn-j} is used as the sequence T B , at this time, if P value If it is less than the preset value, then the sequence T B The change of sequence T A The direction vector between the monitoring site A to be traced and the monitoring site B to be inspected is V BA , the direction is from the monitoring site B to be inspected to the monitoring site A to be traced, V BA The direction is used to indicate the direction of the pollution source propagation corresponding to the current pollutant monitoring value of the monitoring station to be traced, such as Figure 4 shown.

[0074] When there is more than one monitoring site to be inspected near the monitoring site to be traced, there is also more than one direction vector, then all direction vectors are retained, that is, all direction vectors related to the monitoring site A to be traced are retained.

[0075] After the direction vector related to the monitoring site to be traced is determined, V BA After that, at this time, the first auxiliary monitoring site is found near the monitoring site to be traced, and other monitoring sites near the monitoring site to be traced are traversed according to the pollutant monitoring value of the monitoring site to be traced, and the first auxiliary monitoring site matching the pollutant monitoring value is obtained. The specific process is as follows:

[0076] Look for other monitoring stations (excluding monitoring stations A and B) with the same pollutant peak around the monitoring station A to be traced. It can be determined in turn whether the pollutant monitoring values ​​of other monitoring stations near the monitoring station to be traced exceed the preset peak value. If so, it is determined whether the pollutant type corresponding to the pollutant monitoring value is the same as the current pollutant type. If they are the same, it is determined that the other monitoring station matches the pollutant monitoring value and has the same pollutant peak value, and it is used as the first auxiliary monitoring station.

[0077] When the first auxiliary monitoring site found around the monitoring site A to be traced is monitoring site C, a causal test is performed on the monitoring sites near monitoring site C. Specifically, the second auxiliary monitoring site is found on the map as monitoring site D according to the latitude and longitude information of the first auxiliary monitoring site C, and the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site C and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site D are obtained, and a Granger causality test is performed. At this time, if Pvalue is less than the preset value, the same as above, then the direction vector V between the first auxiliary monitoring site C and the second auxiliary monitoring site D is obtained DC , V DC The direction of is also used to indicate the direction of the pollution source propagation corresponding to the current pollutant monitoring value of the first auxiliary monitoring station C, such as Figure 5 shown.

[0078] When there is a direction vector V that satisfies the conditions BA and direction vector V DC When the direction vector between the monitoring site to be traced and the monitoring site to be inspected and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site are both available, the direction vector between the monitoring site to be traced and the monitoring site to be inspected and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site are extended until they intersect, and the longitude and latitude information of the intersection point is used as the longitude and latitude information of the pollution source. Specifically, Figure 6 As shown, the direction vector V BA and direction vector V DC Extend in the opposite direction of propagation until they intersect. The intersection is the longitude and latitude information of the pollution source that caused the pollution peak.

[0079] In summary, a pollution center source tracing method based on an air quality monitoring network is provided. A direction vector is formed based on the causal test result and the positional relationship between the monitoring sites. The direction vector points from one monitoring site to another. Specifically, when a Granger causality test is performed on the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, if the causal test result is less than a preset value, the current direction vector points from the monitoring site to be tested to the monitoring site to be traced, indicating that the air quality change of the monitoring site to be tested can cause the air quality change of the monitoring site to be traced, that is, the change of the pollutant monitoring value. At this time, the direction vector is used to represent a direction of pollution source propagation corresponding to the current pollutant monitoring value of the monitoring site to be traced, and is extended in the direction of the direction vector, thereby realizing the tracing of the pollutant center. Not only can the data of the existing atmospheric monitoring network be reasonably and effectively utilized, but the data of the atmospheric monitoring network can also be fully combined with the monitoring site information to effectively infer the location of the pollution source.

[0080] Example 2

[0081] A pollution center tracing system based on an air quality monitoring network, comprising:

[0082] Pollutant monitoring value detection module: detect in real time whether the pollutant monitoring value of each monitoring site exceeds the preset peak value. If so, the monitoring site is used as the monitoring site to be traced and go to step S2;

[0083] Monitoring sequence acquisition module: obtain the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced according to the current time, and find the monitoring site to be inspected near the monitoring site to be traced on the map according to the latitude and longitude information of the monitoring site to be traced, and obtain the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be inspected;

[0084] Among them, the monitoring sequence {t1, t2, …, tn} and the lag monitoring sequence {t1-j, t2-j, …, tn-j} are composed of pollutant monitoring values ​​sorted in chronological order, n represents the current time, and j represents the preset lag order;

[0085] The pollution center source tracing module is based on causal testing: the causal testing results are obtained according to the monitoring sequence {t1-j, t2-j, …, tn-j} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j}, and the direction vector between the monitoring site to be traced and the monitoring site to be tested is connected on the map according to the causal testing results; the direction vectors between other monitoring sites are obtained according to the monitoring site to be traced and the monitoring site to be tested; different direction vectors are extended until they intersect, and the longitude and latitude information of the intersection point is used as the longitude and latitude information of the pollution source.

[0086] Preferably, the pollution center source tracing module based on causal verification includes the following submodules:

[0087] The first causality test module: obtains the causality test result according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and connects them on the map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causality test result;

[0088] The first auxiliary monitoring site data acquisition module: traverses other monitoring sites near the monitoring site to be traced according to the pollutant monitoring value of the monitoring site to be traced, obtains other monitoring sites that match the pollutant monitoring value as the first auxiliary monitoring site, and obtains the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site according to the current time;

[0089] The second auxiliary monitoring site data acquisition module: finds the second auxiliary monitoring site near the first auxiliary monitoring site on the map according to the latitude and longitude information of the first auxiliary monitoring site, and obtains the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site;

[0090] The second causality test module: obtains a causality test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connects the first auxiliary monitoring site and the second auxiliary monitoring site on the map according to the causality test result to form a direction vector;

[0091] Tracing module: Extend the direction vector between the monitoring site to be traced and the monitoring site to be inspected, and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source.

[0092] The above description is only a preferred embodiment of the present invention and does not limit the present invention in any form. According to the technical essence of the present invention, within the spirit and principles of the present invention, any simple modification, equivalent replacement and improvement made to the above embodiment still falls within the protection scope of the technical solution of the present invention.

Claims

1. A pollution center tracing method based on an air quality monitoring network, characterized in that: The method comprises the following steps: S1: Real-time detection of whether the pollutant monitoring value of each monitoring station exceeds the preset peak value. If so, the monitoring station is used as the monitoring station to be traced, and the process goes to step S2; S2: Obtain the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced according to the current time, and find the monitoring site to be inspected near the monitoring site to be traced on the map according to the latitude and longitude information of the monitoring site to be traced, and obtain the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be inspected; Among them, the monitoring sequence {t1, t2, …, tn} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} are composed of pollutant monitoring values ​​sorted in chronological order, n represents the current moment, tn represents the pollutant monitoring value at the current moment, and ti represents the pollutant monitoring value of the monitoring station every hour, i∈[1,n]; j represents the preset lag order; S3: Obtain causal test results based on the monitoring sequence {t1-j, t2-j, …, tn-j} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j}, and connect on the map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causal test results; obtain direction vectors between other monitoring sites based on the monitoring site to be traced and the monitoring site to be tested; extend different direction vectors until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source; The S3 specifically includes the following steps: S31: Obtain a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and connect them on a map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causal test result; S32: Traverse other monitoring sites near the monitoring site to be traced according to the pollutant monitoring value of the monitoring site to be traced, obtain other monitoring sites that match the pollutant monitoring value as the first auxiliary monitoring site, and obtain the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site according to the current time; S33: Find a second auxiliary monitoring site near the first auxiliary monitoring site on the map according to the latitude and longitude information of the first auxiliary monitoring site, and obtain a delayed monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site; S34: Obtain a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connect the first auxiliary monitoring site and the second auxiliary monitoring site on a map according to the causal test result to form a direction vector between them; S35: Extend the direction vector between the monitoring site to be traced and the monitoring site to be inspected and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source; In S32, the process of traversing other monitoring stations near the monitoring station to be traced according to the pollutant monitoring value of the monitoring station to be traced is specifically as follows: First, the corresponding pollutant type is obtained according to the pollutant monitoring value of the monitoring station to be traced, and the other monitoring stations near the monitoring station to be traced are traversed with the current pollutant type, and it is determined in turn whether the pollutant type corresponding to the pollutant monitoring value of other monitoring stations near the monitoring station to be traced is the same as the current pollutant type. If so, it is determined that the other monitoring stations match the pollutant monitoring value; Alternatively, first obtain the corresponding pollutant type based on the pollutant monitoring value of the monitoring site to be traced, traverse other monitoring sites near the monitoring site to be traced with the current pollutant type, and determine in turn whether the pollutant monitoring values ​​of other monitoring sites near the monitoring site to be traced exceed the preset peak value. If exceeded, determine whether the pollutant type corresponding to the pollutant monitoring value is the same as the current pollutant type. If the same, determine that the other monitoring sites match the pollutant monitoring value.

2. A pollution center source tracing method based on an air quality monitoring network according to claim 1, characterized in that: In S2, the preset hysteresis order is related to the site distance between the current monitoring site to be traced and the current monitoring site to be inspected, and the wind speed at the current moment.

3. The pollution center tracing method based on the air quality monitoring network according to claim 1 is characterized in that: In S31, the causal test result is obtained according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and the process of connecting the direction vector between the monitoring site to be traced and the monitoring site to be tested on the map according to the causal test result is specifically as follows: A Granger causality test is performed on the monitoring sequence {t1, t2, …, tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the monitoring site to be tested. If the causality test result is less than the preset value, they are connected on the map to form a direction vector from the monitoring site to be tested to the monitoring site to be traced.

4. The pollution center tracing method based on the air quality monitoring network according to claim 1 is characterized in that: In S34, the process of obtaining a causal test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connecting the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site on the map according to the causal test result is specifically as follows: A Granger causality test is performed on the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site. If the causality test result is less than the preset value, they are connected on the map to form a direction vector from the second auxiliary monitoring site to the first auxiliary monitoring site.

5. A pollution center tracing system based on an air quality monitoring network, characterized in that: include: Pollutant monitoring value detection module: Real-time detection of whether the pollutant monitoring value of each monitoring site exceeds the preset peak value. If so, the monitoring site is used as a monitoring site to be traced; Monitoring sequence acquisition module: obtain the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced according to the current time, and find the monitoring site to be inspected near the monitoring site to be traced on the map according to the latitude and longitude information of the monitoring site to be traced, and obtain the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be inspected; Among them, the monitoring sequence {t1, t2, …, tn} and the lagged monitoring sequence {t1-j, t2-j, …, tn-j} are composed of pollutant monitoring values ​​sorted in chronological order, n represents the current moment, tn represents the pollutant monitoring value at the current moment, and ti represents the pollutant monitoring value of the monitoring station every hour, i∈[1,n]; j represents the preset lag order; Based on the causal test, the pollution center source tracing module is carried out: the causal test results are obtained according to the monitoring sequence {t1-j, t2-j, ..., tn-j} and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j}, and the direction vector between the monitoring site to be traced and the monitoring site to be tested is connected on the map according to the causal test results; the direction vector between other monitoring sites is obtained according to the monitoring site to be traced and the monitoring site to be tested; different direction vectors are extended until they intersect, and the longitude and latitude information of the intersection point is used as the longitude and latitude information of the pollution source; The pollution center tracing module based on causal testing includes the following submodules: The first causality test module: obtains the causality test result according to the monitoring sequence {t1, t2, ..., tn} of the monitoring site to be traced and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the monitoring site to be tested, and connects them on the map to form a direction vector between the monitoring site to be traced and the monitoring site to be tested according to the causality test result; The first auxiliary monitoring site data acquisition module: traverses other monitoring sites near the monitoring site to be traced according to the pollutant monitoring value of the monitoring site to be traced, obtains other monitoring sites that match the pollutant monitoring value as the first auxiliary monitoring site, and obtains the monitoring sequence {t1, t2, …, tn} of the first auxiliary monitoring site according to the current time; The second auxiliary monitoring site data acquisition module: finds the second auxiliary monitoring site near the first auxiliary monitoring site on the map according to the latitude and longitude information of the first auxiliary monitoring site, and obtains the lagged monitoring sequence {t1-j, t2-j, …, tn-j} of the second auxiliary monitoring site; The second causality test module: obtains a causality test result according to the monitoring sequence {t1, t2, ..., tn} of the first auxiliary monitoring site and the lagged monitoring sequence {t1-j, t2-j, ..., tn-j} of the second auxiliary monitoring site, and connects the first auxiliary monitoring site and the second auxiliary monitoring site on the map according to the causality test result to form a direction vector; Source tracing module: Extend the direction vector between the monitoring site to be traced and the monitoring site to be inspected and the direction vector between the first auxiliary monitoring site and the second auxiliary monitoring site until they intersect, and use the longitude and latitude information of the intersection point as the longitude and latitude information of the pollution source; The specific process of traversing other monitoring stations near the monitoring station to be traced is as follows: First, the corresponding pollutant type is obtained according to the pollutant monitoring value of the monitoring station to be traced, and the other monitoring stations near the monitoring station to be traced are traversed with the current pollutant type, and it is determined in turn whether the pollutant type corresponding to the pollutant monitoring value of other monitoring stations near the monitoring station to be traced is the same as the current pollutant type. If so, it is determined that the other monitoring stations match the pollutant monitoring value; Alternatively, first obtain the corresponding pollutant type based on the pollutant monitoring value of the monitoring site to be traced, traverse other monitoring sites near the monitoring site to be traced with the current pollutant type, and determine in turn whether the pollutant monitoring values ​​of other monitoring sites near the monitoring site to be traced exceed the preset peak value. If exceeded, determine whether the pollutant type corresponding to the pollutant monitoring value is the same as the current pollutant type. If the same, determine that the other monitoring sites match the pollutant monitoring value.

Images