A sea surface height numerical simulation data correction method based on deep learning

Abstract

The application provides a kind of numerical simulation data correction method of sea surface height based on deep learning, it is related to data processing technical field, specifically includes: preparation data, and data is handled;Based on numerical simulation data, simulation satellite data is constructed;Deep learning assimilation model is constructed;With simulation data, background field data and auxiliary physical quantity as input, with the numerical simulation data of specified time as label, the deep learning assimilation model is trained, and the model capable of outputting the corrected high-resolution SSHA field is obtained;The dataset of real high-resolution wide-swath satellite observation data, low-resolution daily satellite data product, numerical simulation data and auxiliary data is put into the trained model, and the corrected sea surface height anomaly data product is output.The technical scheme of the application overcomes the problem that high-resolution satellite observation and low-resolution satellite product with good continuity cannot be fully integrated in the prior art.

Description

TECHNICAL FIELD

[0001] The present application relates to the field of data processing, in particular to a sea surface height numerical simulation data correction method based on deep learning. BACKGROUND

[0002] Sea surface height anomaly (SSHA) is an important physical quantity for describing the dynamic structure and thermal process of the ocean. The main ways to obtain SSHA include ocean numerical simulation and satellite remote sensing observation. High-resolution ocean models can provide continuous and spatially consistent SSHA simulation data, making up for the shortcomings of satellite observations, such as orbital limitations, data discontinuity, and limited spatial coverage. However, numerical models inevitably have error sources, are highly dependent on initial fields and boundary conditions, and have some simplifications and idealizations in the expression of physical processes. Therefore, systematic biases often occur in dynamic active regions such as boundary currents and tropical circulations, affecting the identification of key ocean phenomena such as mesoscale eddies and thermal content distribution. In addition, the sensitivity of numerical models to forcing terms such as wind stress and heat flux can further amplify errors, limiting their reliability in high-precision ocean prediction and intelligent business applications.

[0003] To reduce model errors, the common method is to introduce multi-source observation data into numerical models for error correction and state reconstruction. The current mainstream technology route is the data assimilation method based on physical mechanisms, such as three-dimensional variation (3D-Var), four-dimensional variation (4D-Var), and ensemble Kalman filter (EnKF). These methods incorporate observation data into the simulation system through the construction of error covariance structure between observation and model state, thereby improving the accuracy and reliability of model results. Although these methods have been successfully applied to atmospheric and oceanic prediction systems and deployed in some business systems, there are still significant limitations. First, physical assimilation algorithms require extremely high computing resources, usually relying on high-performance computing platforms, making it difficult to achieve rapid deployment and real-time operation. Second, these methods are highly dependent on the distribution of observation data, making it difficult to achieve stable assimilation in sparsely observed or discontinuous regions, especially in high-resolution ocean scenarios, where the assimilation effect will decrease significantly with spatial scale refinement. In addition, the system construction of traditional assimilation technology is complex, requiring fine tuning for different observation types, resolutions, and regional characteristics, resulting in poor generality and migration.

[0004] Although assimilation models based on deep learning have made good progress in the field of meteorology, they are still in the initial exploration stage in the field of ocean. There are the following challenges in the SSHA correction problem: firstly, the nonlinear, multiscale coupling of ocean dynamic processes leads to complex error propagation mechanism, and the model is difficult to fully capture; secondly, the limitations of satellite observations in spatiotemporal coverage and quality control restrict the construction of high-precision training samples; thirdly, most of the current assimilation models still rely on observations and historical simulation results for training, and lack of sufficient modeling of physical consistency and causal mechanism, which still has deficiencies in generalization ability, credibility and interpretability.

[0005] Therefore, there is a need for a deep learning-based sea surface height numerical simulation data correction method that fully integrates high-resolution satellite observations and low-resolution satellite products with good continuity to effectively correct high-resolution numerical simulation SSHA, thereby improving the reliability and applicability of ocean models in high-resolution and high-precision application scenarios. SUMMARY

[0006] The main purpose of the present application is to provide a deep learning-based sea surface height numerical simulation data correction method to solve the problem that the prior art cannot fully integrate high-resolution satellite observations and low-resolution satellite products with good continuity.

[0007] To achieve the above-mentioned purpose, the present application provides a deep learning-based sea surface height numerical simulation data correction method, which specifically comprises the following steps:

[0008] S1, preparing data includes: high-resolution wide swath data, low-resolution daily satellite data product, numerical simulation data and auxiliary data, and processing the data.

[0009] S2, constructing simulated satellite data based on numerical simulation data, including: high-resolution wide swath satellite simulation data and low-resolution daily satellite product simulation data.

[0010] S3, constructing a deep learning assimilation model, including an encoder, a processor and a decoder connected to each other.

[0011] S4, taking the simulation data, the background field data and the auxiliary physical quantity as the input, and taking the numerical simulation data at the specified time as the label, training the deep learning assimilation model to obtain a model capable of outputting the corrected high-resolution SSHA field.

[0012] S5, putting the real high-resolution wide swath satellite observation data, low-resolution daily satellite data product, numerical simulation data and auxiliary data into the model trained in step S4, and outputting the corrected sea surface height anomaly data product.

[0013] Further, step S1 specifically comprises the following steps:

[0014] S1.1. Spatio-temporal matching of different source data, and unifying the spatial grid of all data to the high-resolution grid of numerical simulation data through spatial interpolation.

[0015] S1.2. Outlier rejection and consistency processing of data.

[0016] Further, step S2 specifically comprises the following steps:

[0017] S2.1. Construction of high-resolution wide-swath satellite simulation data: selecting numerical simulation SSHA data with the same latitude and longitude range as the high-resolution wide-swath observation data as the high-resolution wide-swath satellite simulation data, and simulating the actual observation error by adding noise.

[0018] S2.2. Construction of low-resolution daily satellite product simulation data: interpolating the numerical simulation SSHA data to the grid of the low-resolution daily satellite data product, and then re-interpolating the interpolated data to the numerical simulation grid, introducing blur processing to approximate the low-resolution observation characteristics, and adding noise to simulate observation uncertainty.

[0019] Further, step S3 specifically comprises the following steps:

[0020] S3.1. Processing of non-gridded observations and data missing through the SetConv layer of the encoder, as shown in formula (1):

[0021] (1) ;

[0022] wherein, represents the estimated value in the grid coordinate, represents a learnable kernel function, represents the target grid point, represents the spatial coordinate of the th observation data point, represents the observation value of the th observation data point.

[0023] Grid representation of multi-source observations is input into the visual Transformer module of the encoder for deep learning feature fusion:

[0024] (2) ;

[0025] wherein, represents the input set of original multi-source observation data, represents the set convolution processing process, The feature extraction and fusion process represented by the encoder visual Transformer, The output feature of the visual Transformer module representing the encoder.

[0026] S3.2, the features are extracted by using the multiple visual Transformers cascaded in the processor, the predicted field is obtained by using formula (3), and the long-distance space-time dependence is captured:

[0027] (3) ;

[0028] Wherein, The global ocean state at the time step is the specific output of the VIT network.

[0029] S3.3, the high-resolution SSHA field is reconstructed by using upsampling in the decoder.

[0030] Further, the background field data in step S4 is the numerical simulation SSHA data near the specified date; the auxiliary physical quantities include: sea surface temperature, sea surface flow field and geopotential height, and the output label is the true numerical simulation SSHA field of the specified date.

[0031] The present application has the following beneficial effects:

[0032] 1. The present application combines high-resolution wide-swath satellite data and low-resolution daily satellite data product organically, uses the spatial fine structure of the former and the time continuity of the latter, and corrects the high-resolution numerical simulation SSHA data, which breaks through the limitation of traditional methods that are difficult to balance spatial resolution and time continuity;

[0033] 2. The present application uses a deep learning assimilation model to correct the numerical simulation results, avoids the dependence on error covariance matrix and large-scale iterative calculation in traditional physical assimilation methods, greatly improves the calculation efficiency, and can realize rapid correction and real-time application;

[0034] 3. The sea surface height numerical simulation data method proposed in the present application introduces background field and auxiliary physical quantity in the model training process, so that the model learns the mode bias law while enhancing the representation ability of complex dynamic process, which can effectively reduce the systematic bias of the model, improve the precision, stability and scalability of the prediction results.

[0035] ​In summary, the present application proposes a numerical simulation data correction method for sea surface height based on deep learning, which combines high-resolution wide-swath satellite observations with low-resolution satellite products with better time continuity, ensuring data accuracy while considering spatial resolution and time continuity. This method effectively reduces the systematic bias of the model in the dynamic active area, significantly improves the reliability of SSHA data, and can generate high-resolution data sets covering a long time series. Compared with traditional assimilation methods, the present application has higher computational efficiency, simpler structure, and good robustness and scalability, providing reliable support for ocean prediction, climate change research, and intelligent ocean systems. BRIEF DESCRIPTION OF DRAWINGS

[0036] In order to more clearly illustrate the technical solutions in the specific embodiments or prior art of the present application, the drawings needed in the specific embodiments or prior art description will be briefly introduced as follows. Obviously, the drawings in the following description are some embodiments of the present application, and those skilled in the art can obtain other drawings according to these drawings without creative labor. In the drawings:

[0037] Figure 1 A flowchart of a numerical simulation data correction method for sea surface height based on deep learning is shown.

[0038] Figure 2 A sea surface height anomaly map of numerical simulation data before correction is shown.

[0039] Figure 3 A real sea surface height anomaly map is shown.

[0040] Figure 4 A sea surface height anomaly map of numerical simulation data after correction is shown.

[0041] Figure 5 A sea surface height anomaly difference distribution map before and after correction is shown.

[0042] Figure 6 A sea surface height anomaly difference distribution map after correction and real sea surface height anomaly data is shown. DETAILED DESCRIPTION

[0043] The technical solutions of the present application will be described in detail below with reference to the drawings. Obviously, the described embodiments are part of the embodiments of the present application, not all embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without creative labor are within the scope of protection of the present application.

[0044] As Figure 1The method for correcting numerical simulation data of sea surface height based on deep learning shown specifically comprises the following steps:

[0045] S1, preparing data comprises: high-resolution wide swath data, low-resolution daily satellite data products, numerical simulation data and auxiliary data, and processing the data.

[0046] S2, constructing simulated satellite data based on numerical simulation data, comprising: high-resolution wide swath satellite simulation data and low-resolution daily satellite product simulation data.

[0047] S3, constructing a deep learning assimilation model, comprising an encoder, a processor and a decoder connected to each other.

[0048] S4, taking the simulation data, the background field data and the auxiliary physical quantity as input, and taking the numerical simulation data at the specified time as label, training the deep learning assimilation model to obtain a model capable of outputting the corrected high-resolution SSHA field; the auxiliary physical quantity is ERA5 auxiliary data.

[0049] S5, putting the real high-resolution wide swath satellite observation data, low-resolution daily satellite data products, numerical simulation data and auxiliary data into the model trained in step S4, and outputting the corrected sea surface height anomaly data product.

[0050] The present application proposes a numerical simulation data correction method based on deep learning in view of the problems of difficult to balance spatial resolution and time continuity in sea surface height anomaly (SSHA) data, systematic deviation of numerical model simulation results and complex calculation of traditional assimilation method. The method provided by the present application fuses high spatial resolution but time sparse wide swath satellite observation and time continuous but low spatial resolution satellite product, significantly improves the spatial resolution and time continuity of SSHA data under the premise of ensuring data accuracy, realizes effective correction of model data, and can generate high-resolution SSHA data set covering long time sequence through data reconstruction, providing more reliable and complete data support for marine dynamic process analysis and business forecast.

[0051] Specifically, step S1 specifically comprises the following steps:

[0052] S1.1, time and space matching is performed on data from different sources, and all data are uniformly interpolated to the high-resolution grid of the numerical simulation data through spatial interpolation to ensure the comparability and fusibility of the data.

[0053] S1.2, outlier rejection and consistency processing are performed on the data to reduce the influence of observation errors and noise on model training.

[0054] Specifically, due to the deficiencies in spatiotemporal coverage and quality control of real high-resolution wide-swath satellite observation data and low-resolution daily satellite products, they are difficult to use directly as large-scale training samples. This invention first constructs two types of simulated satellite data based on numerical simulation data to approximate real observation characteristics and simulate their error distribution, thereby ensuring that the training process can fully learn the differences between observations and models. This simulation strategy can alleviate the problem of sparse and discontinuous observation data and provide a systematic set of high-quality training samples for deep learning models, ensuring that the model can effectively capture the regularity of model bias. The numerical simulation SSHA data is processed to construct two types of satellite simulation data. Step S2 specifically includes the following steps:

[0055] S2.1 Constructing high-resolution wide-swath satellite simulation data: Select numerical simulation SSHA data with the same latitude and longitude range as the high-resolution wide-swath observation data as the high-resolution wide-swath satellite simulation data, and simulate actual observation errors by adding noise.

[0056] S2.2, Constructing Low-Resolution Daily Satellite Product Simulation Data: Numerical simulation SSHA data is interpolated onto the grid of the low-resolution daily satellite data product. Then, the interpolated data is re-interpolated back onto the numerical simulation grid. Fuzzing is introduced to approximate the characteristics of low-resolution observations, and noise is added to simulate observation uncertainties. Specifically, fuzzing involves interpolating 2km resolution model data onto a 0.125-degree grid, and then re-interpolating the data already interpolated to the 0.125-degree grid back onto the model data grid. This results in a lower model data resolution, matching the resolution of the low-resolution daily satellite product, and thus, fuzzing is applied.

[0057] Specifically, step S3 includes the following steps:

[0058] S3.1, the problem of non-mesh observations and missing data is handled by the SetConv layer of the encoder, as shown in formula (1): learnable kernel functions to aggregate sparse observations into regular grids, thereby enhancing the robustness and generalization ability of the model.

[0059] (1);

[0060] in, This represents the estimated value in grid coordinates. This represents a learnable kernel function. Indicates the target grid point. Indicates the first Spatial coordinates of each observation data point Indicates the first The observed values ​​of each observation data point.

[0061] The multi-source observation grid representation obtained by formula (1) is input into the visual Transformer module of the encoder for deep learning feature fusion; the module captures long-range spatiotemporal dependencies through its self-attention mechanism and outputs high-quality grid initial state, i.e., the output feature of the visual Transformer module of the encoder , and the process is as shown in formula (2):

[0062] (2);

[0063] Wherein, represents the input original multi-source observation data set, represents the set convolution processing process, represents the feature extraction and fusion process of the encoder visual Transformer.

[0064] The encoder includes a SetConv layer and a visual Transformer module connected to each other.

[0065] S3.2, a plurality of visual Transformers (i.e., VIT networks) are cascaded in the processor to extract features, i.e., to obtain a predicted field by formula (3), to capture long-range spatiotemporal dependencies: an autoregressive prediction chain is formed by the cascaded visual Transformers, a predicted field is obtained by formula (3), and the accuracy and reliability of the model are improved.

[0066] (3);

[0067] Wherein, is a high-dimensional tensor, representing the global ocean state at time step is the specific output of the VIT network, representing the predicted state change.

[0068] S3.3, in the decoder, the high-resolution SSHA field is reconstructed by upsampling.

[0069] Specifically, step S4 needs to ensure that the deep learning assimilation model can learn the error relationship between the numerical simulation and the observation, and can also identify the deviation characteristics of the model in time evolution when training the deep learning assimilation model. Therefore, the numerical simulation SSHA data of a specified date is set as a label to provide a clear training target, i.e., the idealized true field at that time; at the same time, the numerical simulation results near the date are selected as the background field input, so that the deep learning assimilation model can utilize the model error information contained in the time adjacent data, thereby more effectively realizing the correction. Step S4 includes the following steps:

[0070] ​S4.1, making a training data set

[0071] The numerical simulation SSHA data of a specified date is taken as the label data, i.e., the real sea surface height field. The observation input is constructed by step S2, i.e., two types of simulated satellite data are first generated based on the numerical simulation data: one is a high-resolution wide-swath simulated observation obtained by intercepting the same latitude and longitude region data as the wide-swath satellite coverage range and adding noise; the other is a low-resolution daily satellite product simulated observation obtained by first interpolating the numerical simulation data to the grid of the low-resolution daily satellite product, then interpolating back to the original model grid and performing blurring and noise adding operations. At the same time, in order to provide the bias information of the model itself, the numerical simulation SSHA data near the specified date is also selected as the background field input; and further, auxiliary physical quantities of the same date as these background fields are introduced, including sea surface temperature (SST), sea surface flow field (U, V) and geopotential, as additional input features to enhance the constraint on the model error characteristics. After time matching, spatial alignment and data cleaning processing, the finally formed input data is composed of two types of simulated satellite observations, background field data and auxiliary physical quantities, and the output label is the real numerical simulation SSHA field of the specified date. Based on this process, the complete training set, validation set and test set are constructed.

[0072] S4.2, the deep learning assimilation model constructed by step S3 dynamically fuses low-resolution daily satellite data products, high-resolution wide-swath data and numerical simulation data, efficiently extracts deep semantic features and accurately restores spatial details.

[0073] Specifically, in step S4, the background field data is the numerical simulation SSHA data near the specified date; the auxiliary physical quantities include sea surface temperature, sea surface flow field and geopotential, and the output label is the real numerical simulation SSHA field of the specified date. After the training of the deep learning assimilation model is completed, the data set made of real high-resolution wide-swath satellite observation data, low-resolution daily satellite data products, numerical simulation data and auxiliary data is put into the model trained by step S4. The model corrects and reconstructs the input satellite observation data in spatial and temporal dimensions according to the multi-source data features and error correction rules learned in the training stage, and outputs the corrected high-resolution sea surface height anomaly data product with strong temporal continuity and significantly reduced error.

[0074] The application designs a dynamic standardization preprocessing procedure. Mean and standard deviation of background field, observation field and true value field are independently calculated based on a training set, and fixed standardization parameters are generated and synchronized to a validation set and a test set to eliminate distribution deviation of multi-source data; satellite simulation data and numerical simulation data are used to train a model, sea surface height field of the numerical simulation data is obtained through the satellite simulation data, and after the model is trained, real satellite observation data high-resolution wide swath data and low-resolution daily satellite data products can be used to obtain real sea surface height field. The encoder-processor-decoder cascade architecture is used to realize the whole process prediction, the efficiency is significantly improved compared with the traditional scheme, the end-to-end joint fine-tuning is supported, and the error is reduced. The training loss function only calculates the error of the effective observation area, and the dynamic learning rate scheduling (ReduceLROnPlateau) is combined, and finally the high spatial and temporal resolution numerical simulation data sea surface height anomaly product is output.

[0075] In order to verify the method provided by the application, taking the data of a sea area on January 12, 2004 as an example, the abscissa represents longitude, and the ordinate represents latitude. Figure 2 is a sea surface height anomaly map before the numerical simulation data is corrected; Figure 3 is a real sea surface height anomaly map; Figure 4 is a sea surface height anomaly map after the numerical simulation data is corrected; by comparing the three maps, it can be found that the area near the land has been significantly corrected.

[0076] As shown in Figure 5 , the sea surface height anomaly difference map of the numerical simulation data before and after correction, the deeper the color, the greater the difference, the root mean square difference is 0.0933, indicating that the sea surface height is effectively corrected by deep learning; as shown in Figure 6 , the difference distribution map of the numerical simulation data after correction and the real sea surface height anomaly data, the deeper the color, the greater the difference between the two, although there is a difference, but the difference is small, the root mean square difference is 0.0017, indicating that the numerical simulation data after correction by the deep learning method provided by the application has a small difference with the real sea surface height field.

[0077] Of course, the above description is not a limitation of the application, and the application is not limited to the above examples, and the changes, modifications, additions or replacements made by those skilled in the art within the essential scope of the application should also be within the protection scope of the application.

Claims

1. A method for correcting sea surface height numerical simulation data based on deep learning, characterized in that, Specifically, the steps include the following: S1. Prepare data including: high-resolution wide-swath data, low-resolution daily satellite data products, numerical simulation data and auxiliary data, and process the data; S2, based on numerical simulation data, constructs simulated satellite data, including: high-resolution wide-swath satellite simulation data and low-resolution daily satellite product simulation data; S3 builds deep learning assimilation models, including interconnected encoders, processors, and decoders; S4 takes simulation data, background field data, and auxiliary physical quantities as inputs, and numerical simulation data at a specified time as labels to train a deep learning assimilation model, resulting in a model that can output a corrected high-resolution SSHA field. S5: Input the dataset made from real high-resolution wide-swath satellite observation data, low-resolution daily satellite data products, numerical simulation data and auxiliary data into the model trained in step S4, and output the corrected sea surface height anomaly data product. Step S3 specifically includes the following steps: S3.1, the problem of non-meshable observations and missing data is handled through the SetConv layer of the encoder, as shown in formula (1): （1）； in, This represents the estimated value in grid coordinates. This represents a learnable kernel function. Indicates the target grid point. Indicates the first Spatial coordinates of each observation data point Indicates the first The observed values ​​of each observation data point; The multi-source observations are represented in a gridded format and then input into the visual Transformer module of the encoder for deep learning feature fusion. （2）； in, This represents the original set of multi-source observation data input. This represents the set convolution process. The feature extraction and fusion process of the encoder's visual Transformer. This represents the output characteristics of the encoder's visual Transformer module; S3.2, features are extracted using multiple cascaded visual Transformers in the processor, and the prediction field is obtained using formula (3) to capture long-distance spatiotemporal dependencies: （3）； in, Indicates time step The overall state of the ocean It is the specific output of the VIT network; S3.3 utilizes upsampling in the decoder to reconstruct a high-resolution SSHA field.

2. The method for correcting sea surface height numerical simulation data based on deep learning according to claim 1, characterized in that, Step S1 specifically includes the following steps: S1.1 performs spatiotemporal matching on data from different sources and unifies the spatial grid of all data to the high-resolution grid of the numerical simulation data through spatial interpolation. S1.2, perform outlier removal and data consistency processing.

3. The method for correcting sea surface height numerical simulation data based on deep learning according to claim 1, characterized in that, Step S2 specifically includes the following steps: S2.1, Constructing high-resolution wide-swath satellite simulation data: Select numerical simulation SSHA data with the same latitude and longitude range as the high-resolution wide-swath observation data as the high-resolution wide-swath satellite simulation data, and simulate actual observation errors by adding noise; S2.2, Constructing low-resolution daily satellite product simulation data: Interpolate the numerical simulation SSHA data onto the grid of the low-resolution daily satellite data product, and then re-interpolate the interpolated data onto the numerical simulation grid. Introduce fuzzing to approximate the low-resolution observation characteristics, and add noise to simulate observation uncertainties.

4. The method for correcting sea surface height numerical simulation data based on deep learning according to claim 1, characterized in that, In step S4, the background field data is numerical simulation SSHA data near a specified date; the auxiliary physical quantities include: sea surface temperature, sea surface current field and geopotential height, and the output label is the real numerical simulation SSHA field of the specified date.

Images