Universe pollutant diffusion field simulation method and system based on sparse data

Abstract

The invention provides a global pollutant diffusion field simulation method and system based on sparse data, and relates to the technical field of pollutant diffusion simulation, and the method comprises the steps: S1, collecting meteorological data, pollution source data and site monitoring concentration data at each natural moment; s2, generating an initial particle set of the natural moment T according to the natural moment T and the meteorological data and pollution source data before the natural moment T; s3, combining the particle set at the natural moment T and the particle moment t with the meteorological data and pollution source data at the natural moment T, estimating to obtain the state of the particle set at the natural moment T and the particle moment t + 1, performing error evaluation, and when an error meets an acceptable condition, executing the step S3; according to the state estimation of the particle set at the natural moment T and the particle moment t + 1, obtaining an RF estimation field at the natural moment T, and turning to step S4; s5, a predicted pollutant concentration field at the natural moment T + 1 is obtained through prediction of the generator; and S6, evaluating the predicted pollutant concentration field at the natural moment T + 1 through a discriminator.

Description

Technical Field

[0001] The present invention relates to the technical field of pollutant diffusion simulation, and in particular to a method and system for simulating a global pollutant diffusion field based on sparse data. Background Art

[0002] At present, the methods that can be used to simulate the diffusion field of atmospheric pollutants can be divided into the following three categories: 1. Classical methods: including Gaussian model, Lagrangian method, CMAQ, CALPUFF, etc. These methods are usually based on theoretical models of physics and chemistry, assuming that emission sources, meteorological conditions, etc. are relatively known, and are solved through specific mathematical formulas or equations; 2. Machine learning / deep learning algorithms: such as attention networks, gated recurrent units, and attention mechanisms. These methods rely on data-driven models and can automatically adjust model parameters for prediction and simulation by learning the complex relationship between input and output in historical data; 3. Statistical algorithms: such as particle filtering, genetic algorithms, simulated annealing, etc. These methods are optimization algorithms, usually used to solve parameter tuning problems in complex systems, and can find suitable solutions through search and optimization. In diffusion simulation, they can be used to optimize model parameters or sample and infer different hypotheses.

[0003] However, each of these three methods has certain limitations, especially for the global pollutant diffusion field simulation scenario. Because the monitoring stations in a large area are unevenly distributed, there is no emission information or the emission information is incomplete in some areas. In this case, pollutant diffusion simulation and estimation will face certain challenges. The lack or incomplete emission information will affect the accuracy and reliability of the model. Therefore, a single method often has limitations in accuracy, adaptability, and real-time performance. Summary of the invention

[0004] The purpose of the present invention is to provide a method and system for simulating the global pollutant diffusion field based on sparse data, which mainly combines particle filtering and generative adversarial network algorithms. It can also solve the problem of simulating the global pollutant diffusion field well in sparse data and complex environments.

[0005] In order to solve the above technical problems, the present invention adopts the following solutions: A global pollutant diffusion field simulation method based on sparse data, the global pollutant diffusion field simulation method comprising the following steps: S1. Collect meteorological data, pollution source data and site monitoring concentration data at each natural moment, and package the meteorological data, pollution source data and site monitoring concentration data at the same natural moment into one moment data; S2. Generate an initial particle set at natural time T according to meteorological data and pollution source data before natural time T, wherein the initial particle set refers to a particle set whose particle time is t=0; S3, combining the particle set at natural time T and particle time t with the meteorological data and pollution source data at natural time T, estimating the state of the particle set at natural time T and particle time t+1 through the diffusion model, and performing error evaluation. When the error meets the acceptable condition, the RF estimation field at natural time T is obtained according to the state estimation of the particle set at natural time T and particle time t+1, and going to step S4; S4, according to the concentration distribution field at natural time T and before natural time T, combined with the meteorological data and pollution source data at natural time T and before natural time T, the predicted pollutant concentration field at natural time T+1 is obtained through the generator prediction; S5. Obtain the real pollutant concentration field at natural time T+1, combine it with the predicted pollutant concentration field at natural time T+1, and evaluate the predicted pollutant concentration field at natural time T+1 through the discriminator.

[0006] Furthermore, the meteorological data includes pressure, wind speed, wind direction, temperature, precipitation, relative humidity, solar radiation, cloud cover, and grid ID; the pollution source data includes the location of the pollution source, emission intensity, emission height, and type of pollutants emitted; the site monitoring concentration data is the pollutant concentration data monitored by the monitoring site at each moment, including PM2.5, PM10, NO2, SO2, and CO.

[0007] Furthermore, in S2, the process of generating the initial particle set at natural time T according to the meteorological data and pollution source data at natural time T and before natural time T is: Generate a particle set at natural time T, initialize the particle set at natural time T according to meteorological data and pollution source data at natural time T and before natural time T, randomly initialize each particle through noise or initialize with prior information to obtain the initial state of the particle, and then obtain the initial particle set at natural time T.

[0008] Further, in S3, the error evaluation includes the following steps: SA1, error evaluation of the particle set at natural time T and particle time t+1; SA2. After obtaining the error evaluation result, the particle set at natural time T and particle time t+1 is updated according to the site monitoring concentration data at natural time T to obtain an updated particle set at natural time T and particle time t+1, and then go to step SA3; SA3. Judge the error evaluation results: When the error meets the acceptable conditions, the RF estimation field of the natural time T is obtained according to the updated particle set at the natural time T and the particle time t+1, and go to step S4; when the error does not meet the acceptable conditions, update t to t', t'=t+1, and go to step SA1.

[0009] Furthermore, in SA2, the process of obtaining the updated particle set at natural time T and particle time t+1 is specifically as follows: Update the weight of each particle of the particle set at natural time T and particle time t+1 according to the site monitoring concentration data at natural time T; Particle resampling is performed according to the updated weights of each particle to update the weight of the particle set at the natural time T and the particle time t+1, thereby obtaining the updated particle set at the natural time T and the particle time t+1.

[0010] Furthermore, the process of updating the weight of each particle of the particle set at natural time T and particle time t+1 according to the site monitoring concentration data at natural time T is as follows: Obtain the site monitoring concentration data at natural time T and the state of the particle set at natural time T and particle time t, compare the site monitoring concentration data at natural time T with the state of the particle set at natural time T and particle time t, and update the weight of each particle according to the comparison result.

[0011] Furthermore, the particle resampling is used to remove particles with lower weights in the particle set, and the weight of the particle set at the natural time T and the particle time t+1 is updated.

[0012] Furthermore, the process of estimating the RF estimated field at natural time T according to the updated particle set at natural time T and particle time t+1 is: The weighted average calculation is performed on the particles in the particle set at the updated natural time T and the particle time t+1, and the obtained weighted average value is used as the final concentration field, that is, the RF estimation field at the natural time T.

[0013] A global pollutant diffusion field simulation system based on sparse data, using a global pollutant diffusion field simulation method based on sparse data, comprising: Data collection module: collects meteorological data, pollution source data and site monitoring concentration data at each natural moment, and packages the meteorological data, pollution source data and site monitoring concentration data at the same natural moment into one moment data; Particle set initialization module: Generates an initial particle set at natural time T according to meteorological data and pollution source data before natural time T, wherein the initial particle set refers to a particle set whose particle time is t=0; Concentration distribution field estimation module: The particle set at natural time T and particle time t is combined with the meteorological data and pollution source data at natural time T, and the state of the particle set at natural time T and particle time t+1 is estimated through the diffusion model, and the error is evaluated. When the error meets the acceptable conditions, the RF estimation field at natural time T is obtained according to the state of the particle set at natural time T and particle time t+1; Generator prediction module: Based on the concentration distribution field at natural time T and before natural time T, combined with the meteorological data and pollution source data at natural time T and before natural time T, the generator predicts the pollutant concentration field at natural time T+1; Discriminator evaluation module: obtain the real pollutant concentration field at natural time T+1, combine it with the predicted pollutant concentration field at natural time T+1, and evaluate the predicted pollutant concentration field at natural time T+1 through the discriminator.

[0014] Beneficial effects of the present invention: The present invention provides a method and system for simulating a global pollutant diffusion field based on sparse data. In the case of sparse data and uneven distribution of monitoring sites, PF-GAN can effectively make up for the data gaps between monitoring sites and fill in blank areas through its flexible generation model, nonlinear modeling ability, adversarial training process and combination with particle filtering, and generate a reasonable pollutant concentration distribution field under irregular data input, and provide more accurate diffusion estimation results. In addition, by alternately updating particle filtering and GAN, the estimation of the pollutant concentration field is gradually optimized, so that the model can make full use of existing data to accurately simulate the pollutant diffusion field in the sparse data scenario. After each particle filtering step, the pollutant concentration field is predicted by the generator, and then adversarial training is performed by the discriminator. The GAN training process can improve the accuracy of the generated pollutant concentration field, while the particle filter maintains the state estimation of the pollutant concentration field. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] Figure 1 It is a flowchart of a method for simulating a global pollutant diffusion field based on sparse data in Example 1 of the present invention.

[0016] Figure 2 It is a schematic diagram of the process of particle filtering in Example 1 of the present invention.

[0017] Figure 3 Schematic diagram of the process of adversarial network training in Example 1 of the present invention. DETAILED DESCRIPTION

[0018] 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.

[0019] 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.

[0020] 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.

[0021] 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.

[0022] 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.

[0023] 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.

[0024] The present invention will be described in detail below with reference to the accompanying drawings and in conjunction with embodiments: Example 1 At present, monitoring stations over a large area are usually unevenly distributed, resulting in some blank areas with no emission information or incomplete emission information. In this case, pollutant diffusion simulation and estimation will face certain challenges. The lack or incomplete emission information will affect the accuracy and reliability of the model. Therefore, a single method often has limitations in accuracy, adaptability, and real-time performance.

[0025] In order to solve the above problems, a global pollutant diffusion field simulation method based on sparse data is proposed in this embodiment. Based on the nested algorithm of statistical algorithm and machine / deep learning, the particle filter-generative adversarial algorithm (PF-GAN) is adopted, which can better solve the diffusion simulation problem of sparse data, has higher adaptability and better real-time performance, and can take into account both accuracy and flexibility to adapt to rapidly changing environments. The PF-GAN combines particle filtering and generative adversarial networks, and can handle the challenges brought by sparse data well without being restricted by gridding. Therefore, in this specific scenario, it has obvious advantages, as follows: (1) No gridding requirement, flexible response to uneven data distribution: PF-GAN does not rely on a regular grid structure and can directly process data with sparse data and uneven distribution of monitoring sites. Traditional grid-based methods, such as PF-CNN, require that data must be regularly gridded, which will cause many difficulties in a sparse data environment. PF-GAN, by generating a model, can generate a reasonable pollutant concentration field under irregular data input and fill in the data gaps between sites. (2) Using adversarial training to fill data gaps: The core advantage of the generative adversarial network (GAN) lies in its generative ability. The generator network can generate concentration fields based on existing sparse data, meteorological and pollution source information. The adversarial training process between the generator and the discriminator can make the generated data closer and closer to the actual pollutant concentration distribution. Even if the original data points are sparse, GAN can learn the global characteristics of the data and generate a smooth and reasonable concentration field to fill the gaps in the sparse data area. (3) Efficient modeling of nonlinear relationships: The pollutant diffusion process is usually highly nonlinear and is affected by a variety of complex factors (such as meteorology, geographical environment, emission sources, etc.). Especially in the generator part, this nonlinear relationship can be captured through a complex neural network structure. It is more flexible than traditional physical model-based methods (such as Gaussian diffusion model) and can automatically learn the complex patterns in the pollutant diffusion process. (4) Enhanced adaptability to observed data through particle filtering: Particle filtering (PF) can dynamically estimate the state of the pollutant concentration field and update it using observed data. By updating the weights of particles, PF can enhance the adaptability to the uncertainty and dynamic changes of observed data, especially in environments with sparse data and high noise. When combined with GAN, particle filtering can help the generator generate pollutant concentration fields more accurately, ensuring that the generated results conform to the distribution characteristics of the observed data. (5) Improving simulation accuracy and generating diverse simulation results: PF-GAN can generate multiple different concentration field simulation results, which provides diverse perspectives and uncertainty assessments for the diffusion prediction of atmospheric pollutants. This is very important in practical applications because the diffusion of atmospheric pollutants is highly uncertain. Generating multiple simulation results helps provide more comprehensive decision support. (6) Enhance the robustness and generalization ability of the model: PF-GAN can effectively avoid overfitting and enhance the generalization ability of the model in the face of sparse data or different types of pollutant diffusion fields through generative adversarial training. Especially when the site distribution is uneven and the data points are scarce, GAN can learn the potential spatial structure extracted from the uneven data, so as to better apply it to the diffusion simulation of other regions and time periods.

[0026] In summary, in the case of sparse data and uneven distribution of monitoring sites, PF-GAN can effectively make up for the data gaps between monitoring sites, fill in the blank areas, generate reasonable pollutant concentration fields under irregular data input, and provide more accurate diffusion prediction results through its flexible generative model, nonlinear modeling ability, adversarial training process, and combination with particle filtering.

[0027] Specifically, a global pollutant diffusion field simulation method based on sparse data, such as Figure 1 As shown, the method comprises the following steps: S1. Collect meteorological data, pollution source data and site monitoring concentration data at each natural moment, and package the meteorological data, pollution source data and site monitoring concentration data at the same natural moment into one moment data; S2. Generate an initial particle set at natural time T according to meteorological data and pollution source data before natural time T, wherein the initial particle set refers to a particle set whose particle time is t=0; S3, combining the particle set at natural time T and particle time t with the meteorological data and pollution source data at natural time T, estimating the state of the particle set at natural time T and particle time t+1 through the diffusion model, and performing error evaluation. When the error meets the acceptable condition, the RF estimation field at natural time T is obtained according to the state estimation of the particle set at natural time T and particle time t+1, and going to step S4; S4, according to the concentration distribution field at natural time T and before natural time T, combined with the meteorological data and pollution source data at natural time T and before natural time T, the predicted pollutant concentration field at natural time T+1 is obtained through the generator prediction; S5. Obtain the real pollutant concentration field at natural time T+1, combine it with the predicted pollutant concentration field at natural time T+1, and evaluate the predicted pollutant concentration field at natural time T+1 through the discriminator.

[0028] The meteorological data are: , where the dimension of meteorological data is: meteorological data at time T, including pressure, wind speed, wind direction, temperature, precipitation, relative humidity, solar radiation, cloud cover, grid ID and other D-dimensional features.

[0029] The pollution source data are: , where the dimension of pollution source data is: pollution source information at time T, including location, emission intensity, etc., with a total of S-dimensional features.

[0030] The monitoring concentration data of the site are: , where the dimension of the site monitoring concentration data is: pollutant concentration data of M sites at time T, including PM2.5, PM10, NO2, SO2, CO and other parameters.

[0031] The time series data is: T time data { , , }, the meteorological data, pollution source data and site monitoring concentration data at the same time are packaged into one time data.

[0032] After the data was collected, it was found that the distribution of monitoring sites was often uneven under the global pollutant diffusion field simulation, and the data in some areas was scarce or completely missing. Traditional grid-based methods, such as PF-CNN, require that the data must be regularly gridded, which will cause many difficulties in this sparse data environment. Therefore, the present invention first uses particle filtering (PF) to dynamically estimate the state of the pollutant concentration field and uses the observed data to update it. The process of dynamically estimating the state of the pollutant concentration field using particle filtering (PF), such as Figure 2 As shown in Figure 1, the first step is to initialize the particle set. First, generate a particle set at natural time T. Initialize the particle set at natural time T according to the meteorological data and pollution source data at natural time T and before natural time T. Randomly initialize each particle through noise or initialize it with prior information to obtain the initial state of the particle. Then, obtain the initial particle set at natural time T, where each particle represents a possible pollutant concentration field. After initializing the particle set, multiple possible states of the pollutant concentration field can be obtained. The specific process is as follows: Input data: meteorological data at the initial time t=0 and pollution source data , that is, the particle set is initialized with the meteorological data and pollution source data at the initial time t=0; Get the initialized particle set: ; Each particle , represents a possible pollutant concentration field, is the initial state of the ith particle, and N is the total number of particles.

[0033] The second step is particle prediction, which mainly predicts the pollutant concentration field of particles at the next moment based on meteorological data and pollution source data. Here, it is necessary to first obtain the particles at the current moment T from the particle set. And its corresponding meteorological data and pollution source data , and then the state of the particle set at particle time t=1 is predicted by the diffusion model. The specific process is as follows: Input data: Particles at the current moment , meteorological data and pollution source data , that is, obtaining the meteorological data and pollution source data of particles from the particle set; Prediction step: For each particle , update the concentration field of the particles according to the diffusion model f(·), that is, input the meteorological data and pollution source data of the particles in the initial particle set at natural time T and particle time t=0 into the diffusion model in sequence, and predict the state of the particles at the next time through the diffusion model to obtain the state of the particle set at particle time t=1, that is, the concentration field of the particle set at particle time t=1. Specifically, the prediction is: ,in, Represents the updated particle concentration field.

[0034] The third step is particle weighting. By comparing the concentration field predicted by the particle with the actual observation data, the particle weight is updated. The update of the particle weight reflects the degree of fit of the particle to the observation data. The particle with a higher weight indicates that the prediction of the concentration field is more accurate. Specifically, the weight of each particle in the particle set at natural time T and particle time t+1 is updated according to the site monitoring concentration data at natural time T. The process is as follows: Input data: The monitoring concentration data of the site at the current time T , according to the concentration field predicted at the current time T ; Weighting step: The weight of the particle is calculated by the likelihood function: ; Among them, the weight of each particle , indicating the degree of match between each particle state and the monitoring concentration data at the site.

[0035] The fourth step is the resampling step, which mainly resamples according to the weight of the particles to generate a new set of particles. This can not only remove particles with lower weights, but also ensure that particles with higher weights are more likely to be selected, which has a greater impact on the prediction of the next moment. Specifically, the particles are resampled according to the updated weights of each particle to update the weights of the particle set at natural time T and particle time t+1. The process is as follows: Input data: current particle set , particle weight ; Resampling step: Resample according to the particle weights to generate a new particle set: ; Among them, the dimension of the new particle set is the same as the dimension of the particle set, both are R M .

[0036] The fifth step is the result output, which is mainly to calculate the concentration field prediction result at the current time based on the weighted average of the particles. Specifically, the RF estimation field at the natural time T is obtained based on the updated particle set at the natural time T and the particle time t+1. The concentration distribution field at the current time t can be estimated through particle filtering to provide input for the subsequent GAN generation process. The process is as follows: Output data: The weight of each particle in the updated particle set is calculated by weighted average, and the estimated particle weighted average value is used as the final pollutant concentration field, that is, the concentration distribution field of the particle set at the current time t. The formula is: ; Among them, the concentration distribution field at the current time T, Represents the prediction result of the pollutant concentration field at the current time T.

[0037] After obtaining the predicted results of the pollutant concentration field at the current time T through particle filtering, the generative adversarial network (GAN) training is selected to generate the future pollutant concentration field. The adversarial network is trained adversarially through the generator and the discriminator, which can make the generated pollutant concentration field closer to the true value.

[0038] like Figure 3 As shown, first, a generator is constructed, through which the pollutant concentration field at the next moment can be generated according to the pollutant concentration field, meteorological data and pollution source data at the current moment. Specifically, the concentration distribution field of the particle set at the next moment T+1 is generated by the generator according to the concentration distribution field, meteorological data and pollution source data of the updated particle set at the current moment T. The generator generates a future pollutant concentration field close to the real one through adversarial training, and the process is as follows: Input data: prediction results of pollutant concentration field output by particle filter , meteorological data , pollution source data ; Generator model: Generator G generates the predicted pollutant concentration field at the next time T+1 based on the input data: ;in, It represents the generated pollutant concentration field, that is, the predicted pollutant concentration field at the next time T+1.

[0039] Then, a discriminator is constructed, which can evaluate whether the generated pollutant concentration field is realistic and whether it matches the real observed data. Specifically, the predicted pollutant concentration field of the particle set at the next time T+1 is obtained, and it is judged whether the predicted pollutant concentration field at the next time T+1 matches the real pollutant concentration field at the next time T+1. The discriminator makes the pollutant concentration field generated by the generator closer to the real value through adversarial training with the generator. The process is as follows: Input data: The predicted pollutant concentration field generated at the next moment T+1 And the actual pollutant concentration field at the next moment T+1 ; Discriminator model: The discriminator D outputs a value indicating whether the input concentration field is real ,in, Represents the discriminator's evaluation of the authenticity of the generated results.

[0040] Next, the generator and discriminator are optimized according to the matching results. Specifically, by optimizing the loss functions of the generator and the discriminator, the pollutant concentration field generated by the generator is made closer to the true value. The process is as follows: Loss function (generator): ; Among them, L G ∈R, represents the loss of the generator; Loss function (discriminator): ; Among them, L D ∈R, represents the loss of the discriminator.

[0041] Finally, the accuracy of the model can be evaluated by comparing the actual monitoring data and simulation results based on the error between the calculated pollutant concentration field and the actual observed data, such as the mean square error (MSE). For example: , the pollutant concentration field that can be simulated The actual pollutant concentration field Error assessment.

[0042] Example 2 A global pollutant diffusion field simulation system based on sparse data, using a global pollutant diffusion field simulation method based on sparse data, comprising: Data collection module: collects meteorological data, pollution source data and site monitoring concentration data at each natural moment, and packages the meteorological data, pollution source data and site monitoring concentration data at the same natural moment into one moment data; Particle set initialization module: Generates an initial particle set at natural time T according to meteorological data and pollution source data before natural time T, wherein the initial particle set refers to a particle set whose particle time is t=0; Concentration distribution field estimation module: The particle set at natural time T and particle time t is combined with the meteorological data and pollution source data at natural time T, and the state of the particle set at natural time T and particle time t+1 is estimated through the diffusion model, and the error is evaluated. When the error meets the acceptable conditions, the RF estimation field at natural time T is obtained according to the state of the particle set at natural time T and particle time t+1; Generator prediction module: Based on the concentration distribution field at natural time T and before natural time T, combined with the meteorological data and pollution source data at natural time T and before natural time T, the generator predicts the pollutant concentration field at natural time T+1; Discriminator evaluation module: obtain the real pollutant concentration field at natural time T+1, combine it with the predicted pollutant concentration field at natural time T+1, and evaluate the predicted pollutant concentration field at natural time T+1 through the discriminator.

[0043] 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 global pollutant diffusion field simulation method based on sparse data, characterized in that: The global pollutant diffusion field simulation method comprises the following steps: S1. Collect meteorological data, pollution source data and site monitoring concentration data at each natural moment, and package the meteorological data, pollution source data and site monitoring concentration data at the same natural moment into one moment data; S2. Generate an initial particle set at natural time T according to meteorological data and pollution source data before natural time T, wherein the initial particle set refers to a particle set whose particle time is t=0; S3, combining the particle set at natural time T and particle time t with the meteorological data and pollution source data at natural time T, estimating the state of the particle set at natural time T and particle time t+1, and performing error evaluation. When the error meets the acceptable condition, the RF estimation field at natural time T is obtained according to the state estimation of the particle set at natural time T and particle time t+1, and going to step S4; S4, according to the concentration distribution field at natural time T and before natural time T, combined with the meteorological data and pollution source data at natural time T and before natural time T, the predicted pollutant concentration field at natural time T+1 is obtained through the generator prediction; S5. Obtain the real pollutant concentration field at natural time T+1, combine it with the predicted pollutant concentration field at natural time T+1, and evaluate the predicted pollutant concentration field at natural time T+1 through the discriminator.

2. The method for simulating global pollutant diffusion field based on sparse data according to claim 1 is characterized in that: The meteorological data include pressure, wind speed, wind direction, temperature, precipitation, relative humidity, solar radiation, cloud cover, and grid ID; the pollution source data include the location of the pollution source, emission intensity, emission height, and type of pollutants emitted; the site monitoring concentration data is the pollutant concentration data monitored by the monitoring site at each moment, including PM2.5, PM10, NO2, SO2, and CO.

3. The method for simulating global pollutant diffusion field based on sparse data according to claim 1, characterized in that: In S2, the process of generating the initial particle set at natural time T according to the meteorological data and pollution source data at natural time T and before natural time T is: Generate a particle set at natural time T, initialize the particle set at natural time T according to meteorological data and pollution source data at natural time T and before natural time T, randomly initialize each particle through noise or initialize with prior information to obtain the initial state of the particle, and then obtain the initial particle set at natural time T.

4. The method for simulating global pollutant diffusion field based on sparse data according to claim 1 is characterized in that: In S3, the error evaluation includes the following steps: SA1, error evaluation of the particle set at natural time T and particle time t+1; SA2. After obtaining the error evaluation result, the particle set at natural time T and particle time t+1 is updated according to the site monitoring concentration data at natural time T to obtain an updated particle set at natural time T and particle time t+1, and then go to step SA3; SA3. Judge the error evaluation results: When the error meets the acceptable conditions, the RF estimation field of the natural time T is obtained according to the updated particle set at the natural time T and the particle time t+1, and go to step S4; when the error does not meet the acceptable conditions, update t to t', t'=t+1, and go to step SA1.

5. The method for simulating global pollutant diffusion field based on sparse data according to claim 4 is characterized in that: In SA2, the process of obtaining the updated particle set at natural time T and particle time t+1 is as follows: Update the weight of each particle of the particle set at natural time T and particle time t+1 according to the site monitoring concentration data at natural time T; Particle resampling is performed according to the updated weights of each particle to update the weight of the particle set at the natural time T and the particle time t+1, thereby obtaining the updated particle set at the natural time T and the particle time t+1.

6. The method for simulating global pollutant diffusion field based on sparse data according to claim 5 is characterized in that: The process of updating the weight of each particle in the particle set at natural time T and particle time t+1 according to the site monitoring concentration data at natural time T is: Obtain the site monitoring concentration data at natural time T and the state of the particle set at natural time T and particle time t, compare the site monitoring concentration data at natural time T with the state of the particle set at natural time T and particle time t, and update the weight of each particle according to the comparison result.

7. The method for simulating global pollutant diffusion field based on sparse data according to claim 5, characterized in that: The particle resampling is used to remove particles with lower weights in the particle set, and the weight of the particle set at the natural time T and the particle time t+1 is updated.

8. The method for simulating global pollutant diffusion field based on sparse data according to claim 4 is characterized in that: The process of estimating the RF estimated field at natural time T based on the updated particle set at natural time T and particle time t+1 is: The weighted average calculation is performed on the particles in the particle set at the updated natural time T and the particle time t+1, and the obtained weighted average value is used as the final concentration field, that is, the RF estimation field at the natural time T.

9. A global pollutant diffusion field simulation system based on sparse data, characterized in that: A method for simulating a global pollutant diffusion field based on sparse data as described in any one of claims 1 to 8 is applied, comprising: Data collection module: collects meteorological data, pollution source data and site monitoring concentration data at each natural moment, and packages the meteorological data, pollution source data and site monitoring concentration data at the same natural moment into one moment data; Particle set initialization module: Generates an initial particle set at natural time T according to meteorological data and pollution source data before natural time T, wherein the initial particle set refers to a particle set whose particle time is t=0; Concentration distribution field estimation module: The particle set at natural time T and particle time t is combined with the meteorological data and pollution source data at natural time T, and the state of the particle set at natural time T and particle time t+1 is estimated through the diffusion model, and the error is evaluated. When the error meets the acceptable conditions, the RF estimation field at natural time T is obtained according to the state of the particle set at natural time T and particle time t+1; Generator prediction module: Based on the concentration distribution field at natural time T and before natural time T, combined with the meteorological data and pollution source data at natural time T and before natural time T, the generator predicts the pollutant concentration field at natural time T+1; Discriminator evaluation module: obtain the real pollutant concentration field at natural time T+1, combine it with the predicted pollutant concentration field at natural time T+1, and evaluate the predicted pollutant concentration field at natural time T+1 through the discriminator.

Images