Carbon dioxide column concentration prediction method and system based on generative diffusion and atmospheric diffusion fusion model, and electronic equipment

Abstract

The invention discloses a carbon dioxide column concentration prediction method and system based on a generative diffusion and atmospheric diffusion fusion model, and electronic equipment, and belongs to the technical field of environment prediction, and the method comprises the steps: obtaining the time dynamic characteristics of carbon dioxide column concentration data and meteorological data in a target region, according to the time dynamic characteristics and the time-space dynamic characteristics of the regional gridding meteorological data in the target region, fusing the time dynamic characteristics and the time-space dynamic characteristics; performing forward diffusion operation and reverse denoising on the fusion features to obtain a preliminary prediction result; smooth transition between adjacent moments in the prediction result is constrained through local smoothness loss, and the consistency of the change trend between the preliminary prediction result and the existing data is constrained through trend consistency loss. According to the method, spatio-temporal feature fusion is integrated, forward noise adding and reverse denoising are performed based on the diffusion model to reconstruct the prediction data, the moment sequence smooth continuity and the overall physical rationality of the prediction result are ensured through local and trend consistency loss, and the prediction precision is improved.

Description

Technical Field

[0001] This application belongs to the technical field of environmental prediction, and particularly relates to a method, system, and electronic device for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion. Background Art

[0002] Accurate prediction of carbon dioxide concentration is of great significance for climate research, environmental management, and policy making. Traditional methods mainly rely on atmospheric diffusion models (such as the WRF model), which simulate the diffusion process based on physical equations, but their accuracy is limited by the quality of input parameters, grid resolution, and cumulative errors in long-term prediction. In addition, physical models are difficult to capture complex spatio-temporal patterns and non-linear relationships in the data. At the same time, due to the limited spatial and temporal resolution of satellite XCO2 data, it is impossible to accurately reflect local complex terrains or short-term emission peaks, resulting in low prediction accuracy of carbon concentration in small-scale areas such as cities or industrial parks. Summary of the Invention

[0003] Object of the Invention: This application develops a method, system, and electronic device for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion, aiming to solve the technical problem of low prediction accuracy of carbon concentration in the prior art.

[0004] Technical Solution: In a first aspect, an embodiment of this application provides a method for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion, including:

[0005] Identifying a target area, and obtaining first data and second data at multiple existing moments in the target area; the first data includes carbon dioxide column concentration data and meteorological data, and the second data includes regional grid meteorological data of the target area; Fusing the temporal dynamic features of the first data and the spatio-temporal dynamic features of the second data to obtain fused features; Performing forward diffusion and reverse denoising on the fused features to obtain a preliminary prediction result; Constraining the smooth transition between adjacent moments in the preliminary prediction result based on local smoothness loss, and constraining the consistency of the change trend between the preliminary prediction result and the existing carbon dioxide column concentration data based on trend consistency loss to obtain a final prediction result.

[0006] In some embodiments, the step of obtaining the temporal dynamic features of the first data includes: Converting the first data at each moment into a first vector; the representation formula of the first vector includes: ; Wherein, is the corresponding moment The first vector, , is the number of time instants; is a matrix with dimension ; is the number of parameter types in the first data; Based on the order of the time instants, update the hidden state and memory state of the multi-layer long short-term memory network through multiple of the first vectors; its representation formula includes: ; where is the hidden state of the multi-layer long short-term memory network at the corresponding time instant ; is the memory cell state of the multi-layer long short-term memory network at the corresponding time instant ; is the multi-layer long short-term memory network; is the first vector at the corresponding time instant ; is the hidden state of the multi-layer long short-term memory network at the corresponding time instant ; is the memory cell state of the multi-layer long short-term memory network at the corresponding time instant ; Output the sequence of the updated hidden state as the time dynamic feature of the first data.

[0007] In some embodiments, the steps of obtaining the spatio-temporal dynamic feature of the second data include: Obtain the grid information of the second data, and convert each of the second data into a second vector based on the grid information; the representation formula of the second vector includes: ; where is the second vector at the corresponding time instant , , is the number of time instants; is a matrix with dimension ; is the height of the grid, is the width of the grid, is the number of data channels; Input the second data into a convolutional neural network to obtain the spatial feature of the second data; Input the spatial feature into a multi-layer long short-term memory network, and based on the order of the time instants, extract the spatio-temporal dynamic feature of the second data through multiple of the spatial features; its representation formula includes: ; Among them, is the spatio-temporal dynamic feature of the second data; is the multi-layer long short-term memory network; is the corresponding moment of the spatial feature of the second data; is the number of moments; is of dimension matrix, is the dimension of the moment dynamic feature of the output second data.

[0008] In some embodiments, the step of fusing the moment dynamic feature of the first data and the spatio-temporal dynamic feature of the second data to obtain a fused feature includes: Based on the dimension of the spatio-temporal dynamic feature, splicing the moment dynamic feature and the spatio-temporal dynamic feature to obtain a spliced feature; Performing weighted processing on the spliced feature based on the self-attention mechanism to obtain a weighted feature; Obtaining the similarity between features in the weighted feature based on the self-attention mechanism, and increasing the weight of key features based on the similarity to obtain the fused feature.

[0009] In some embodiments, the step of forward diffusion includes: Obtaining the weight coefficient corresponding to each moment; Adding Gaussian noise to the fused feature based on the weight coefficient to obtain the noise sequence at each moment; its representation formula includes: ; Among them, is the noise sequence corresponding to the fused feature; is the weight coefficient corresponding to the moment t; is the

[0010] In some embodiments, the step of reverse denoising includes: Performing reverse denoising operation on the noise sequence through a noise prediction network based on the weight coefficient to obtain a preliminary prediction result; its representation formula includes: ; Among them, is the result after reverse denoising at the moment, and is the noise intensity parameter in the t-th diffusion process, used to control the degree of Gaussian perturbation, ; is the noise prediction network; is the existing meteorological data, time coding, and carbon dioxide column concentration distribution; is a random noise vector sampled from a standard normal distribution, used to simulate the noise in the diffusion process.

[0011] In some embodiments, the characterization formula of the local smoothness loss includes: ; ; where, is the local smoothness loss; is the time span of future prediction, used to balance the growth amplitude of trends under different time windows; is the serial number of the time; is in the reverse diffusion process the generation result at time; is in the reverse diffusion process the generation result at time; (·) is a dynamic weight function; is a control parameter for adjusting the deviation degree, used to control the decline rate of the weight; is the mean value of; is the standard deviation of.

[0012] In some embodiments, the characterization formula of the trend consistency loss includes: ; where, is the trend consistency loss; is the average change trend of the predicted future data sequence; is the average change trend of the existing data sequence; is an explicit slope adjustment term, is the preset slope of the existing upward trend, is the time span of future prediction, used to balance the growth amplitude of trends under different time windows.

[0013] In a second aspect, an embodiment of the present application further provides a carbon dioxide column concentration prediction system based on a generative diffusion and atmospheric diffusion fusion model, including: A data acquisition module, which is used to confirm the target area and acquire the first data and the second data of multiple existing times in the target area; the first data includes carbon dioxide column concentration data and meteorological data, and the second data includes the regional grid meteorological data of the target area; A spatio-temporal feature fusion module, which is used to fuse the moment dynamic features of the first data and the spatio-temporal dynamic features of the second data to obtain fused features; A concentration prediction module, which is used to perform forward diffusion and reverse denoising on the fused features to obtain a preliminary prediction result; A moment consistency constraint module, which is used to constrain the smooth transition between adjacent moments in the preliminary prediction result based on the local smoothness loss, and to constrain the consistency between the preliminary prediction result and the change trend of the existing carbon dioxide column concentration data based on the trend consistency loss to obtain a final prediction result.

[0014] In a third aspect, an embodiment of the present application further provides an electronic device, including a memory and a processor. The memory stores a computer program executable by the processor. When the processor executes the computer program, the steps of the carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion as described in any one of the first aspects are implemented.

[0015] Beneficial effects: Compared with the prior art, a carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion provided by an embodiment of the present application includes obtaining the moment dynamic features of the carbon dioxide column concentration data and meteorological data in a target area, as well as the spatio-temporal dynamic features of the regional grid meteorological data in the target area, and fusing the moment dynamic features and the spatio-temporal dynamic features to obtain fused features; and performing forward diffusion operation and reverse denoising operation on the fused features to obtain a preliminary prediction result; and constraining the smooth transition between adjacent moments in the prediction result through the local smoothness loss to ensure the smooth continuity of the moment sequence, and constraining the consistency between the preliminary prediction result and the change trend of the existing data through the trend consistency loss to ensure the overall physical rationality. The present application integrates spatio-temporal features to provide sufficient conditional information for the diffusion model, then performs forward noise addition and reverse denoising on the fused features based on the diffusion model to reconstruct future prediction data, and finally ensures the smooth continuity of the moment sequence and the overall physical rationality of the prediction result through local and trend consistency losses, improving the prediction accuracy of carbon concentration. Description of the Drawings

[0016] In order to more clearly illustrate the technical solutions in the embodiments of the present application, the following will briefly introduce the drawings required for the description of the embodiments. Obviously, the following drawings are only some embodiments of the present application. For those skilled in the art, without creative efforts, other drawings can be obtained based on these drawings.

[0017] Figure 1The flowchart of the steps of the carbon dioxide column concentration prediction method based on the generative diffusion and atmospheric diffusion fusion model provided by the embodiments of the present application; Figure 2 The flowchart of the steps of obtaining the temporal dynamic characteristics of the first data in the carbon dioxide column concentration prediction method based on the generative diffusion and atmospheric diffusion fusion model provided by the embodiments of the present application; Figure 3 The flowchart of the steps of obtaining the spatio-temporal dynamic characteristics of the second data in the carbon dioxide column concentration prediction method based on the generative diffusion and atmospheric diffusion fusion model provided by the embodiments of the present application; Figure 4 The flowchart of the steps of fusing the temporal dynamic characteristics of the first data and the spatio-temporal dynamic characteristics of the second data in the carbon dioxide column concentration prediction method based on the generative diffusion and atmospheric diffusion fusion model provided by the embodiments of the present application; Figure 5 The module connection diagram of the carbon dioxide column concentration prediction system based on the generative diffusion and atmospheric diffusion fusion model provided by the embodiments of the present application; Figure 6 The structural diagram of the electronic device provided by the embodiments of the present application; Reference numerals: 10, data acquisition module; 20, spatio-temporal feature fusion module; 30, concentration prediction module; 40, temporal consistency constraint module; 100, memory; 200, processor. Detailed implementation manners

[0018] Next, the technical solutions in the embodiments of the present application will be clearly and completely described in conjunction with the accompanying drawings in the embodiments of the present application. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative efforts belong to the scope of protection of the present application.

[0019] The accurate prediction of carbon dioxide concentration is of great significance for climate research, environmental management and policy making. Traditional methods mainly rely on atmospheric diffusion models (such as the WRF model), which simulate the diffusion process based on physical equations, but their accuracy is limited by the quality of input parameters, grid resolution and cumulative errors in long-term predictions. In addition, physical models are difficult to capture complex spatio-temporal patterns and non-linear relationships in the data. It is worth noting that there is currently a lack of relevant research using deep learning algorithms, especially deep learning algorithms based on the XCO2 dataset, to accurately predict future XCO2. The main advantage of deep learning methods lies in their powerful ability to automatically learn high-level features from a wide range of datasets, which is a key step in bridging the gap between data patterns at different feature levels. Given the excellent feature extraction performance of deep learning neural networks, they have great potential in fusing multi-source data to extract key information.

[0020] There are still the following engineering problems that have not been solved: Problem 1: The spatial and temporal resolution of satellite XCO2 data is limited, and it cannot accurately reflect local complex terrains or short-term emission peaks, affecting the prediction accuracy of carbon concentration in small-scale areas such as cities or industrial parks.

[0021] Problem 2: The lack or imbalance of observational data leads to the problem of data sparsity in the carbon concentration prediction model based on satellite remote sensing, reducing the reliability of the prediction.

[0022] In view of this, embodiments of the present application provide a method for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion, including obtaining the temporal dynamic characteristics of carbon dioxide column concentration data and meteorological data in a target area, as well as the spatio-temporal dynamic characteristics of regional grid meteorological data in the target area, and fusing the temporal dynamic characteristics and spatio-temporal dynamic characteristics to obtain fused characteristics; performing a forward diffusion operation and a reverse denoising operation on the fused characteristics to obtain a preliminary prediction result; and ensuring the smooth transition between adjacent moments in the prediction result through a local smoothness loss constraint, ensuring the smooth continuity of the time series, and ensuring the consistency of the change trend between the preliminary prediction result and the existing data through a trend consistency loss constraint, ensuring the overall physical rationality. The present application integrates spatio-temporal features to provide sufficient conditional information for the diffusion model, then performs forward noise addition and reverse denoising on the fused characteristics based on the diffusion model to reconstruct future prediction data, and finally ensures the smooth continuity of the time series and the overall physical rationality of the prediction result through local and trend consistency losses, improving the prediction accuracy of carbon concentration.

[0023] In some embodiments, please refer to Figure 1 , Figure 1 is a flowchart of the steps of the method for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion provided by embodiments of the present application. The method for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion provided by embodiments of the present application is specifically implemented through steps 100 to 400: Step 100: Confirm the target area and obtain the first data and the second data at multiple moments existing in the target area; the first data includes carbon dioxide column concentration data and meteorological data, and the second data includes regional grid meteorological data of the target area.

[0024] In some embodiments, the historical carbon dioxide column concentration data and meteorological data are normalized to obtain the first data, where the meteorological data includes parameters such as wind speed, temperature, and humidity; the regional grid meteorological data output by the WRF model is subjected to convolutional dimensionality reduction and time series modeling to obtain the second data.

[0025] Step 200: Integrate the moment dynamic features of the first data and the spatio-temporal dynamic features of the second data to obtain integrated features.

[0026] In some embodiments, refer to Figure 2 , Figure 2 , which is the flowchart of the steps for obtaining the moment dynamic features of the first data in the carbon dioxide column concentration prediction method based on the integrated model of generative diffusion and atmospheric diffusion provided by the embodiments of the present application. The steps for obtaining the moment dynamic features of the first data are specifically implemented through steps 210 to 230: Step 210: Convert the first data at each moment into a first vector.

[0027] In some embodiments, the first data includes the existing carbon dioxide column concentration data and meteorological data at [number of moments] moments. The representation formula of the first data at each moment includes: ; where, is the first vector corresponding to the moment , , is the number of moments; is a matrix with a dimension of , is the number of parameter types in the first data. For example, the parameter types include carbon dioxide column concentration, wind speed, temperature, humidity, etc.

[0028] Step 220: Based on the order of moments, update the hidden state and memory state of the multi-layer long short-term memory network through multiple first vectors.

[0029] In some embodiments, to fully capture the time dynamic features, the present application uses a multi-layer long short-term memory network (LSTM) for modeling. The steps are as follows: Obtain a first input sequence based on the first vectors at multiple moments and perform normalization processing. Among them, the representation formula of the first input sequence includes: ; where, is the first input sequence; is the first vector corresponding to the moment , , is the number of moments.

[0030] Based on the order of moments, the representation formula for updating the hidden state and memory state of the multi-layer long short-term memory network through multiple first vectors includes: ; where, is the corresponding moment The hidden state of the multi-layer long short-term memory network; At the corresponding moment The memory cell state of the multi-layer long short-term memory network; Is the multi-layer long short-term memory network; At the corresponding moment The first vector; At the corresponding moment The hidden state of the multi-layer long short-term memory network; At the corresponding moment The memory cell state of the multi-layer long short-term memory network.

[0031] Step 230: Output the updated sequence of hidden states as the temporal dynamic features of the first data.

[0032] Specifically, after iterating through all time steps, the temporal dynamic features of the first data are obtained , Is a matrix with a dimension of , Is The output feature dimension.

[0033] In some embodiments, please refer to Figure 3 , Figure 3 This is the flowchart of the steps for obtaining the spatio-temporal dynamic features of the second data in the carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion provided by the embodiments of the present application. The steps for obtaining the spatio-temporal dynamic features of the second data are specifically implemented through steps 240 to 260: Step 240: Obtain the grid information of the second data, and convert each second data into a second vector based on the grid information.

[0034] In some embodiments, the representation formula of the second vector includes: ; Wherein, Is the second vector at the corresponding moment , , Is the number of moments; Is a matrix with a dimension of , Is the height of the grid, Is the width of the grid, Is the number of data channels.

[0035] Furthermore, the regional gridded meteorological data is the regional gridded data output by WRF , and has obvious spatial correlation.

[0036] Step 250: Input the second data into the convolutional neural network to obtain the spatial features of the second data.

[0037] In some embodiments, for the second data at each moment, a convolutional neural network (CNN) is used to extract spatial features, and the specific process is as follows: For each Perform multi-layer convolutional operations and use the following formula to extract local spatial features: ; Where and are local spatial features; is a convolutional operation; is a matrix with a dimension of ; is the number of feature channels extracted by the convolutional network.

[0038] Step 260: Input the spatial features into a multi-layer long short-term memory network, and based on the order of moments, extract the spatio-temporal dynamic features of the second data through multiple spatial features.

[0039] In some embodiments, the representation formula for extracting the spatio-temporal dynamic features of the second data through multiple spatial features includes: ; Where is the spatio-temporal dynamic feature of the second data; is a multi-layer long short-term memory network; is the spatial feature of the second data corresponding to the moment ; is the number of moments; is a matrix with a dimension of ; is the dimension of the moment dynamic feature of the second data output.

[0040] In some embodiments, please refer to Figure 4 , Figure 4 is the flowchart of the steps for fusing the moment dynamic feature of the first data and the spatio-temporal dynamic feature of the second data in the carbon dioxide column concentration prediction method based on the fusion of the generative diffusion and atmospheric diffusion models provided by the embodiments of the present application. The steps for fusing the moment dynamic feature of the first data and the spatio-temporal dynamic feature of the second data are specifically implemented through steps 270 to 290: Step 270: Based on the dimension of the spatio-temporal dynamic feature, splice the moment dynamic feature and the spatio-temporal dynamic feature to obtain a spliced feature.

[0041] In some embodiments, the representation formula for the spliced feature includes: ; Among them, is a splicing feature; is the moment dynamic feature of the first data; is the spatio-temporal dynamic feature of the second data; is a matrix with a dimension of ).

[0042] Step 280: Perform weighted processing on the splicing feature based on the self-attention mechanism to obtain a weighted feature.

[0043] In some embodiments, in order to model the interaction relationship between features from different sources, the self-attention mechanism is used to weight the splicing feature, and its representation formula includes: ; Among them, is the self-attention mechanism; is the query vector, representing the target position where the attention needs to be calculated currently, = , is the key vector, which interacts with Q and is used to calculate the similarity, = , is the value vector, which stores the information that actually needs to be extracted, = , is a learnable linear transformation matrix, is the splicing feature; is the normalized exponential function; is the transpose of K, is the dimension, is the transpose symbol.

[0044] Step 290: Obtain the similarity between features in the weighted feature based on the self-attention mechanism, and increase the weight of key features based on the similarity to obtain a fused feature.

[0045] It can be understood that the self-attention mechanism assigns higher weights to key features by calculating the similarity between features, thereby generating a fused spatio-temporal feature representation: ; Among them, is the fused feature; is the self-attention mechanism; is the splicing feature; is a matrix with a dimension of .

[0046] Step 300: Perform forward diffusion and reverse denoising on the fused feature based on the diffusion model to obtain a preliminary prediction result.

[0047] In some embodiments, when performing forward diffusion, the weight coefficients corresponding to each moment are obtained, and in the diffusion model, Gaussian noise is added to the fused features based on the weight coefficients to obtain the noise sequences at each moment; the representation formula thereof includes: ; where is the noise sequence corresponding to the -th moment; is the fused feature; is the weight coefficient corresponding to the moment t; is Gaussian noise.

[0048] After steps, basically degenerates into pure noise, laying the foundation for the reverse denoising diffusion process.

[0049] In some embodiments, when performing reverse denoising, based on the weight coefficients, the noise sequences are subjected to reverse denoising operations through the noise prediction network KAN to obtain preliminary prediction results; the representation formula thereof includes: ; where is the result after reverse denoising at the moment; is the noise intensity parameter in the t-th step of the diffusion process, used to control the degree of Gaussian perturbation, ; is the noise prediction network; is the existing meteorological data, time encoding, and carbon dioxide column concentration distribution; is a random noise vector sampled from the standard normal distribution, used to simulate the noise in the diffusion process.

[0050] After steps of iteration, the future XCO2 concentration distribution map predicted by the system is output. This generation process makes full use of the spatio-temporal conditional information extracted by the fusion module and the WRF regional meteorological data, so that the prediction results are significantly improved in terms of spatial resolution and physical rationality.

[0051] In some embodiments, to train the noise prediction network KAN, the present application uses the standard mean squared error (MSE) loss function to fit the predicted noise and the actual sampled noise, and at the same time introduces the reconstruction error as an auxiliary loss. During the training process, the historical data and the obtained by its forward diffusion are correspondingly subjected to supervised learning to gradually optimize the network parameters θ, so that KAN can accurately predict the noise components added at each time step, thereby ensuring the stability and effectiveness of the inverse diffusion process.

[0052] Step 400: Based on the local smoothness loss constraint, make the transition between adjacent moments in the preliminary prediction result smooth, and based on the trend consistency loss constraint, ensure the consistency between the preliminary prediction result and the changing trend of the existing carbon dioxide column concentration data, and obtain the final prediction result.

[0053] In some embodiments, when the preliminary prediction result is within the confidence range, MSE is used to evaluate the consistency. When the preliminary prediction result exceeds the confidence interval, it is considered not to directly penalize sudden events, but to weaken the penalty by reducing the weight, reflecting the inclusiveness of sudden situations. The representation formula of the local smoothness loss includes: ; ; where is the local smoothness loss; is the time span of future prediction, used to balance the growth amplitude of trends under different time windows; is the serial number of the moment; is the generation result at the moment in the diffusion inverse process; is the generation result at the moment in the diffusion inverse process; (·) is the dynamic weight function; is the control parameter for adjusting the deviation degree, used to control the decline rate of the weight; is the mean value of; is the standard deviation of.

[0054] It can be understood that this local smoothness loss formula ensures the smoothness and continuity of the time series generated by the model by minimizing the difference between the prediction results of adjacent time steps. If there are sudden events (such as fires, large-scale emissions, etc.) in the data, the loss function will give greater tolerance to the prediction results of these events. It conforms to the physical characteristics of the carbon concentration prediction task, because the change of XCO2 is usually smooth and gradual, but there are still occasional sudden situations.

[0055] In some embodiments, in addition to local smoothness, the overall future prediction sequence should maintain a long-term change trend similar to the historical data. For this purpose, this embodiment proposes the long-term trend consistency loss , and the representation formula of the trend consistency loss includes: ; where is the trend consistency loss; is the average change trend of the predicted future data sequence; is the average change trend of the existing data sequence; is the explicit slope adjustment term, is the slope of the preset existing upward trend, is the time span of future prediction, which is used to balance the growth amplitude of trends under different time windows.

[0056] Understandably, this formula keeps the trend change amount of future prediction consistent with the historical trend change amount, and at the same time does not ignore a certain growth deviation of XCO2.

[0057] In some embodiments, to balance local smoothness and overall trend consistency, during model training, the two parts of losses are weighted and summed to obtain the final time consistency loss : ; wherein, and are hyperparameters for balancing local and global constraints, and their values can be determined by cross-validation; is the model loss.

[0058] During the training process, and the original loss of the diffusion model are used as the total loss for backpropagation together, so as to prompt the generated XCO2 prediction sequence to be continuously smooth in time and reflect the long-term historical trend change.

[0059] Understandably, the embodiments of the present application provide a method for predicting carbon dioxide column concentration based on the fusion model of generative diffusion and atmospheric diffusion, including obtaining the moment dynamic characteristics of carbon dioxide column concentration data and meteorological data in the target area, as well as the spatio-temporal dynamic characteristics of regional grid meteorological data in the target area, and fusing the moment dynamic characteristics and spatio-temporal dynamic characteristics to obtain the fusion characteristics; performing forward diffusion operation and reverse denoising operation on the fusion characteristics to obtain a preliminary prediction result; and ensuring the smooth transition between adjacent moments in the prediction result through the local smoothness loss constraint, ensuring the smooth continuity of the moment sequence, and ensuring the consistency of the change trend between the preliminary prediction result and the existing data through the trend consistency loss constraint, ensuring the overall physical rationality. The present application integrates spatio-temporal features to provide sufficient conditional information for the diffusion model, then performs forward noise addition and reverse denoising on the fusion characteristics based on the diffusion model to reconstruct future prediction data, and finally ensures the smooth continuity of the moment sequence and the overall physical rationality of the prediction result through local and trend consistency losses, improving the prediction accuracy of carbon concentration.

[0060] From the detailed description of this application, it can be seen that the solution provided by this application realizes high-precision prediction of future XCO2 concentration distribution by integrating spatio-temporal feature fusion, a prediction network based on diffusion models and WRF data, and temporal consistency constraints. The solution of this application uses LSTM and CNN to extract historical and regional meteorological data features respectively, and conducts in-depth fusion through the self-attention mechanism to provide sufficient conditional information for the diffusion model; during the generation process, a forward noise addition and reverse denoising strategy is adopted to reconstruct future data, and the local and trend consistency losses ensure the smooth continuity of the time series and the overall physical rationality. Experimental results show that the method of this application is superior to traditional prediction methods in terms of prediction accuracy, spatio-temporal consistency, and physical constraint satisfaction, and has high engineering application value.

[0061] The implementation details of the above steps can be adjusted and extended according to actual needs to ensure good results under different regions and different data conditions. The solution of this application fully embodies the innovative idea of combining deep generative models with physical meteorological data to achieve complex spatio-temporal dynamic modeling, providing a new technical path for carbon concentration prediction. The following problems can be solved: Problem 1: The spatial and temporal resolutions of satellite XCO2 data are limited, unable to accurately reflect local complex terrains or short-term emission peaks, affecting the prediction accuracy of carbon concentration in small-scale regions such as cities or industrial parks.

[0062] Problem 2: The lack or imbalance of observational data leads to data sparsity problems in carbon concentration prediction models based on satellite remote sensing, reducing the reliability of predictions.

[0063] Therefore, more accurate prediction results can be obtained in the carbon concentration prediction task.

[0064] Correspondingly, the embodiments of this application also provide a carbon dioxide column concentration prediction system based on the fusion model of generative diffusion and atmospheric diffusion. Please refer to Figure 5 , Figure 5 which is the module connection diagram of the carbon dioxide column concentration prediction system based on the fusion model of generative diffusion and atmospheric diffusion provided by the embodiments of this application. The carbon dioxide column concentration prediction system includes: A data acquisition module 10, which is used to confirm the target area and acquire the first data and the second data at multiple existing moments in the target area; the first data includes carbon dioxide column concentration data and meteorological data, and the second data includes the regional grid meteorological data of the target area; A spatio-temporal feature fusion module 20, which is used to fuse the moment dynamic features of the first data and the spatio-temporal dynamic features of the second data to obtain fused features; The concentration prediction module 30 is used to perform forward diffusion and reverse denoising on the fused features to obtain a preliminary prediction result; The temporal consistency constraint module 40 is used to constrain the smooth transition between adjacent times in the preliminary prediction result based on the local smoothness loss, and to constrain the consistency between the preliminary prediction result and the change trend of the existing carbon dioxide column concentration data based on the trend consistency loss, so as to obtain the final prediction result.

[0065] In some embodiments, the spatio-temporal feature fusion module 20 is specifically configured to: Convert the first data at each time into a first vector; the representation formula of the first vector includes: ; Wherein, is the first vector corresponding to the time , , is the number of times; is a matrix with a dimension of , is the number of parameter types in the first data; Based on the order of times, update the hidden state and memory state of the multi-layer long short-term memory network through multiple first vectors; its representation formula includes: ; Wherein, is the hidden state of the multi-layer long short-term memory network corresponding to the time ; is the memory cell state of the multi-layer long short-term memory network corresponding to the time ; is the multi-layer long short-term memory network; is the first vector corresponding to the time ; is the hidden state of the multi-layer long short-term memory network corresponding to the time ; is the memory cell state of the multi-layer long short-term memory network corresponding to the time ; Output the sequence of the updated hidden states as the temporal dynamic features of the first data.

[0066] In some embodiments, the spatio-temporal feature fusion module 20 is specifically configured to: Obtain the grid information of the second data, and convert each second data into a second vector based on the grid information; the representation formula of the second vector includes: ; Wherein, is the corresponding time The second vector, , is the number of time instants; is a matrix with dimension ; is the height of the grid, is the width of the grid, is the number of data channels; Input the second data into the convolutional neural network to obtain the spatial features of the second data; Input the spatial features into the multi-layer long short-term memory network, and based on the order of time instants, extract the spatio-temporal dynamic features of the second data through multiple spatial features; its representation formula includes: ; where, is the spatio-temporal dynamic feature of the second data; is the multi-layer long short-term memory network; is the spatial feature of the second data corresponding to the time instant ; is the number of time instants; is a matrix with dimension ; is the dimension of the time dynamic feature of the output second data.

[0067] In some embodiments, the spatio-temporal feature fusion module 20 is specifically configured to: Based on the dimension of the spatio-temporal dynamic feature, concatenate the time dynamic feature and the spatio-temporal dynamic feature to obtain a concatenated feature; Perform weighted processing on the concatenated feature based on the self-attention mechanism to obtain a weighted feature; Obtain the similarity between features in the weighted feature based on the self-attention mechanism, and based on the similarity, increase the weight of the key features to obtain a fused feature.

[0068] In some embodiments, the concentration prediction module 30 is specifically configured to: Obtain the weight coefficients corresponding to each time instant; In the diffusion model, add Gaussian noise to the fused feature based on the weight coefficients to obtain the noise sequence at each time instant; its representation formula includes: ; where, is the noise sequence corresponding to the th time instant; is the fused feature; is the weight coefficient corresponding to the time instant t; is the Gaussian noise.

[0069] In some embodiments, the concentration prediction module 30 is specifically configured to: Based on the weight coefficient, perform an inverse denoising operation on the noise sequence through the noise prediction network to obtain a preliminary prediction result; its representation formula includes: ; where is the result after inverse denoising at time is the noise intensity parameter in the t-th diffusion process, used to control the degree of Gaussian perturbation, ; is the noise prediction network; is the existing meteorological data, time encoding, and carbon dioxide column concentration distribution; is a random noise vector sampled from the standard normal distribution, used to simulate the noise in the diffusion process.

[0070] Correspondingly, please refer to Figure 6 , Figure 6 which is the structural diagram of the electronic device provided by the embodiment of the present application. The electronic device provided by the embodiment of the present application includes: A memory 100 and a processor 200, which are communicatively connected to each other. The memory 100 stores computer instructions, and the processor 200 executes the computer instructions to perform the carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion provided by the embodiment of the present application. Specifically, the memory 100 and the processor 200 are connected through a communication bus.

[0071] In some embodiments, the processor 200 may be a central processing unit (CPU). The processor 200 may also be other general-purpose processors, digital signal processors (DSPs), application specific integrated circuits (ASICs), field-programmable gate arrays (FPGAs), or other programmable logic devices, discrete gate or transistor logic devices, discrete hardware components, etc. chips, or a combination of the above types of chips.

[0072] In some embodiments, the memory 100, as a non-transitory computer-readable storage medium, can be used to store non-transitory software programs, non-transitory computer-executable programs, and modules, such as the program instructions / modules corresponding to the methods in the method embodiments of the present application. The processor 200 executes various functional applications and data processing of the processor 200 by running the non-transitory software programs, instructions, and modules stored in the memory 100, that is, to implement the methods in the above method embodiments.

[0073] In some embodiments, the memory 100 may include a program storage area and a data storage area. Among them, the program storage area can store an operating device and application programs required for at least one function; the data storage area can store data created by the processor 200, etc. In addition, the memory 100 may include a high-speed random access memory 100, and may also include a non-transitory memory 100, such as at least one disk memory 100 device, a flash memory device, or other non-transitory solid-state memory 100 devices. In some embodiments, the memory 100 may optionally include a memory 100 remotely provided relative to the processor 200, and these remote memories 100 can be connected to the processor 200 through a network. Examples of the above networks include but are not limited to the Internet, enterprise intranets, local area networks, mobile communication networks, and combinations thereof.

[0074] In some embodiments, one or more modules are stored in the memory 100 and, when executed by the processor 200, implement the methods in the above method embodiments.

[0075] In some embodiments, the specific details of the above electronic device can be understood by referring to the corresponding relevant descriptions and effects in the above method embodiments, and will not be elaborated here.

[0076] The above has introduced in detail a carbon dioxide column concentration prediction method, system, and electronic device based on a fusion model of generative diffusion and atmospheric diffusion provided by the embodiments of the present application. Specific examples are used in this article to elaborate on the principle and implementation manner of the present application. The description of the above embodiments is only used to help understand the method and its core idea of the present application; at the same time, for those skilled in the art, according to the idea of the present application, there will be changes in the specific implementation manner and application scope. In summary, the content of this specification should not be construed as a limitation to the present application.

Claims

1. A method for predicting carbon dioxide column concentration based on a fusion model of generative diffusion and atmospheric diffusion, characterized in that Including: Confirm the target area, and obtain the first data and the second data at multiple existing moments in the target area; the first data includes carbon dioxide column concentration data and meteorological data, and the second data includes the regional grid meteorological data of the target area; Fuse the temporal dynamic characteristics of the first data and the spatio-temporal dynamic characteristics of the second data to obtain fused characteristics; Perform forward diffusion and reverse denoising on the fused characteristics based on a diffusion model to obtain a preliminary prediction result; Based on the local smoothness loss, constrain the smooth transition between adjacent moments in the preliminary prediction result, and based on the trend consistency loss, constrain the consistency between the preliminary prediction result and the change trend of the existing carbon dioxide column concentration data to obtain the final prediction result.

2. The carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion according to claim 1, wherein The steps for obtaining the temporal dynamic characteristics of the first data include: Convert the first data at each moment into a first vector; the representation formula of the first vector includes: ； Among them, is the first vector at the corresponding moment , , is the number of moments; is a matrix with dimension , is the number of parameter types in the first data; Based on the order of the moments, update the hidden state and memory state of a multi-layer long short-term memory network through multiple first vectors; its representation formula includes: ； wherein, is the hidden state of the multi-layer long short-term memory network at the corresponding moment; is the memory cell state of the multi-layer long short-term memory network at the corresponding moment; is the corresponding moment; is the memory cell state of the multi-layer long short-term memory network at the corresponding moment; is the multi-layer long short-term memory network; is the corresponding moment; is the first vector at the corresponding moment; is the corresponding moment; is the hidden state of the multi-layer long short-term memory network at the corresponding moment; is the corresponding moment; is the memory cell state of the multi-layer long short-term memory network at the corresponding moment; Output the sequence of the updated hidden state as the temporal dynamic characteristics of the first data.

3. The carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion according to claim 1, wherein The steps for obtaining the spatio-temporal dynamic characteristics of the second data include: Obtain the grid information of the second data, and based on the grid information, convert each second data into a second vector; the representation formula of the second vector includes: ； wherein, is the second vector at the corresponding moment , , is the number of moments; is a matrix with a dimension of , is the height of the grid, is the width of the grid, is the number of data channels; Input the second data into a convolutional neural network to obtain the spatial characteristics of the second data; Input the spatial characteristics into a multi-layer long short-term memory network, and based on the order of the moments, extract the spatio-temporal dynamic characteristics of the second data through multiple spatial characteristics; its representation formula includes: ； Among them, is the spatio-temporal dynamic feature of the second data; is the multi-layer long short-term memory network; is the corresponding moment of the spatial feature of the second data; is the number of moments; is of dimension matrix, is the dimension of the moment dynamic feature of the output second data.

4. The carbon dioxide column concentration prediction method based on the fusion model of generation diffusion and atmospheric diffusion according to claim 1, wherein The steps for fusing the temporal dynamic characteristics of the first data and the spatio-temporal dynamic characteristics of the second data to obtain fused characteristics include: Based on the dimension of the spatio-temporal dynamic characteristics, splice the temporal dynamic characteristics and the spatio-temporal dynamic characteristics to obtain a spliced characteristic; Perform weighted processing on the spliced characteristic based on a self-attention mechanism to obtain a weighted characteristic; Obtain the similarity between features in the weighted characteristic based on the self-attention mechanism, and based on the similarity, increase the weight of key features to obtain the fused characteristic.

5. The carbon dioxide column concentration prediction method based on the fusion model of generative diffusion and atmospheric diffusion according to claim 1, wherein The steps for forward diffusion include: Obtain the weight coefficients corresponding to each moment; In the diffusion model, add Gaussian noise to the fused characteristic based on the weight coefficients to obtain the noise sequence at each moment; its representation formula includes: ； Among them, is the noise sequence corresponding to the th moment; is the fusion feature; is the weight coefficient corresponding to the moment t; is the Gaussian noise.

6. The method for predicting carbon dioxide column concentration based on the fusion model of generative diffusion and atmospheric diffusion according to claim 5, wherein The steps for reverse denoising include: Based on the weight coefficients, perform reverse denoising operations on the noise sequence through a noise prediction network to obtain a preliminary prediction result; its representation formula includes: ； Among them, is the result after reverse denoising at the moment; is the noise intensity parameter in the t-th diffusion process, used to control the degree of Gaussian perturbation, ; is the noise prediction network; is the existing meteorological data, time encoding, and carbon dioxide column concentration distribution; is a random noise vector sampled from the standard normal distribution, used to simulate the noise in the diffusion process.

7. The method for predicting carbon dioxide column concentration based on the fusion model of generative diffusion and atmospheric diffusion according to claim 1, wherein The representation formula of the local smoothness loss includes: ； ； Among them, is the local smoothness loss; is the time span of future prediction, which is used to balance the growth rate of trends under different time windows; is the serial number of the time; in the diffusion inverse process is the generation result at time in the diffusion inverse process is the generation result at time (·) is the dynamic weight function; is the control parameter for adjusting the deviation degree, which is used to control the decline rate of the weight; is the mean value of; is the standard deviation of.

8. The method for predicting carbon dioxide column concentration based on the fusion model of generative diffusion and atmospheric diffusion according to claim 1, wherein The representation formula of the trend consistency loss includes: ； Among them, is the trend consistency loss; is the average change trend of the predicted future data sequence; is the average change trend of the existing data sequence; is an explicit slope adjustment term, is the slope of the preset existing upward trend, is the time span of future prediction, which is used to balance the growth amplitude of trends under different time windows.

9. A carbon dioxide column concentration prediction system based on a fusion model of generative diffusion and atmospheric diffusion, characterized in that, Including: A data acquisition module (10), the data acquisition module (10) is used to confirm the target area and obtain the first data and the second data at multiple existing moments in the target area; the first data includes carbon dioxide column concentration data and meteorological data, and the second data includes the regional grid meteorological data of the target area; A spatio-temporal feature fusion module (20), which is used to fuse the moment dynamic feature of the first data and the spatio-temporal dynamic feature of the second data to obtain a fused feature; A concentration prediction module (30), which is used to perform forward diffusion and reverse denoising on the fused feature to obtain a preliminary prediction result; A moment consistency constraint module (40), which is used to constrain the smooth transition between adjacent moments in the preliminary prediction result based on the local smoothness loss, and to constrain the consistency between the preliminary prediction result and the change trend of the existing carbon dioxide column concentration data based on the trend consistency loss, so as to obtain a final prediction result.

10. An electronic device, characterized in that, It includes a memory (100) and a processor (200). The memory (100) stores a computer program executable by the processor (200). When the processor (200) executes the computer program, the steps of the carbon dioxide column concentration prediction method based on the generative diffusion and atmospheric diffusion fusion model as described in any one of claims 1-8 are implemented.

Images