Method, device and product for determining tidal mixing parameters for ocean numerical models

Abstract

The application discloses a method and device for determining a tidal mixing parameter of a marine numerical model and a product, and comprises the following steps: establishing a marine numerical model to obtain a corresponding state vector; determining an observation matrix corresponding to the state vector; obtaining a prediction error covariance matrix corresponding to the state vector based on the state vector; respectively calculating a Kalman gain of a tidal mixing parameter value and a Kalman gain of a seawater temperature value; respectively calculating a tidal mixing parameter analysis value and a seawater temperature analysis value; and the tidal mixing parameter analysis value is a new tidal mixing parameter value in the marine numerical model. The application optimizes and estimates the tidal mixing parameter and the seawater temperature value of the marine numerical model by using joint state variable set Kalman filtering, more accurately estimates the tidal mixing parameter through data assimilation, reduces the value error of the tidal mixing parameter, and improves the simulation and prediction accuracy of the marine numerical model on a future period.

Description

TECHNICAL FIELD

[0001] The application belongs to the field of ocean data assimilation, and particularly relates to a method and device for determining a tidal mixing parameter of an ocean numerical model and a product. BACKGROUND

[0002] Tidal mixing phenomena exist in the ocean at all times, control the exchange of matter in seawater, and have a great influence on ocean circulation and seawater temperature and salinity. Tidal mixing phenomena cannot be observed, and the strength of tidal mixing phenomena is usually quantified by the size of a tidal mixing parameter in an ocean numerical model. The value of the tidal mixing parameter is generally given by an empirical value or calculated from simulation data of the ocean numerical model, and thus there is often a large error from the actual situation. SUMMARY

[0003] Therefore, the application aims to overcome the defects in the prior art and provides a method and device for determining a tidal mixing parameter of an ocean numerical model and a product.

[0004] To achieve the above object, the technical scheme of the application is as follows:

[0005] The application discloses a method for determining a tidal mixing parameter of an ocean numerical model, comprising: establishing an ocean numerical model, obtaining a corresponding state vector when the ocean numerical model is run to be stable, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value;

[0006] determining an observation matrix for Kalman filtering corresponding to the state vector;

[0007] obtaining a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector;

[0008] respectively calculating a Kalman gain of the tidal mixing parameter value and a Kalman gain of the seawater temperature value based on the state vector, the observation matrix and the prediction error covariance matrix;

[0009] obtaining an observation value, which is observation data of the seawater temperature value;

[0010] respectively calculating a tidal mixing parameter analysis value and a seawater temperature analysis value based on the observation value, the state vector, an observation matrix of the seawater temperature value in the state vector, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value;

[0011] The tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model.

[0012] In an implementation manner of the application, the establishment of the ocean numerical model comprises: using seabed topographic data sets for water depth data of the ocean numerical model.

[0013] In an implementation form of the present application, the ocean numerical model is established, comprising: the tidal boundary condition of the ocean numerical model adopts a Chinese ocean reanalysis dataset.

[0014] In an implementation form of the present application, the initial field of the ocean numerical model is obtained from the Chinese ocean reanalysis dataset, comprising initial temperature, salinity, sea current and water level.

[0015] In an implementation form of the present application, the observation value is observation data of seawater temperature value, comprising: obtaining observation data of seawater temperature of the same period from the dataset.

[0016] In an implementation form of the present application, the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model, comprising: the ocean numerical model performs simulation prediction for a future period based on the new tidal mixing parameter value.

[0017] In a second aspect, the present application discloses a tidal mixing parameter determination device for an ocean numerical model, the device comprising:

[0018] An ocean numerical model establishing module is configured to establish an ocean numerical model, and obtain a corresponding state vector when the ocean numerical model runs to be stable, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value;

[0019] A first determining module is configured to determine an observation matrix for Kalman filtering corresponding to the state vector;

[0020] A second determining module is configured to obtain a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector;

[0021] A third determining module is configured to respectively calculate a Kalman gain of the tidal mixing parameter value and a Kalman gain of the seawater temperature value based on the state vector, the observation matrix and the prediction error covariance matrix;

[0022] A fourth determining module is configured to obtain an observation value, which is observation data of the seawater temperature value;

[0023] A fifth determining module is configured to respectively calculate a tidal mixing parameter analysis value and a seawater temperature analysis value based on the observation value, the state vector, the observation matrix of the seawater temperature value in the state vector, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value;

[0024] A sixth determining module is configured to make the tidal mixing parameter analysis value a new tidal mixing parameter value in the ocean numerical model.

[0025] In a third aspect, the present application discloses an electronic device, comprising: one or more processors; a storage device configured to store one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the method described above.

[0026] In a fourth aspect, the present application discloses a computer-readable storage medium having stored thereon a computer program, which, when executed by a processor, implements the method described above.

[0027] In a fifth aspect, the present application discloses a computer program product comprising a computer program, which, when executed by a processor, implements the method described above.

[0028] Compared with the prior art, the present application has the following advantages:

[0029] The present application discloses a method, device and product for determining a tidal mixing parameter of a marine numerical model, comprising: establishing a marine numerical model, obtaining a corresponding state vector when the marine numerical model is running to be stable; determining an observation matrix for Kalman filtering corresponding to the state vector; obtaining a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector; respectively calculating a Kalman gain of the tidal mixing parameter value and a Kalman gain of the seawater temperature value; obtaining an observation value, which is observation data of the seawater temperature value; respectively calculating a tidal mixing parameter analysis value and a seawater temperature analysis value; the tidal mixing parameter analysis value is a new tidal mixing parameter value in the marine numerical model. The present application discloses a method, device and product for determining a tidal mixing parameter of a marine numerical model, which optimizes and estimates the tidal mixing parameter and the seawater temperature value of the marine numerical model by using joint state variable set Kalman filtering, estimates the tidal mixing parameter more accurately through data assimilation, reduces the value error of the tidal mixing parameter, and improves the simulation and prediction accuracy of the marine numerical model for future periods. BRIEF DESCRIPTION OF DRAWINGS

[0030] The accompanying drawings, which form a part of the present application, are intended to provide further understanding of the present application and are incorporated herein for illustrative purposes. The schematic embodiments of the present application and their descriptions are used to explain the present application and do not constitute improper limitations on the present application.

[0031] In the drawings:

[0032] Figure 1 A schematic diagram of a method for determining a tidal mixing parameter of a marine numerical model according to an embodiment of the present application;

[0033] Figure 2 A schematic diagram of a device for determining a tidal mixing parameter of a marine numerical model according to an embodiment of the present application;

[0034] Figure 3A schematic diagram of an electronic device for determining a tidal mixing parameter of a marine numerical model is provided. DETAILED DESCRIPTION

[0035] It should be noted that the embodiments in the present application and the features in the embodiments can be combined with each other without conflict.

[0036] In the description of the present application, it should be understood that the terms "center", "longitudinal", "transverse", "upper", "lower", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inner", "outer" and the like indicate the orientation or positional relationship based on the orientation or positional relationship shown in the drawings, and are only for the convenience of describing the present application and simplifying the description, and therefore cannot be understood as indicating or implying that the device or element referred to must have a particular orientation, be constructed and operated in a particular orientation, and therefore cannot be understood as a limitation on the present application.

[0037] In the description of the present application, it should be noted that unless otherwise explicitly specified and limited, the terms "mounting", "connection", "connection" should be understood broadly, for example, it can be fixed connection, or detachable connection, or integrally connected; it can be mechanical connection, or electrical connection; it can be directly connected, or indirectly connected through an intermediate medium, or it can be the communication inside two elements. For those skilled in the art, the specific meaning of the above terms in the present application can be understood through specific circumstances.

[0038] In the description of the present application, it should be further noted that the terms "first", "second" and the like are only for descriptive purposes and cannot be understood as indicating or implying relative importance or implicitly indicating the number of the technical features indicated. Therefore, the features limited by "first", "second" and the like can explicitly or implicitly include one or more features. In the description of the present application, unless otherwise specified, the meaning of "multiple" is two or more.

[0039] Under the existing technical conditions, in the marine numerical model, the strength of the tidal mixing phenomenon is usually quantified by the size of the tidal mixing parameter, and the value of the tidal mixing parameter is generally given by the experience value or calculated through the simulation data of the marine numerical model, so there is often a problem of large error from the real situation. The present application discloses a tidal mixing parameter determination method, device and product for a marine numerical model, comprising: using joint state variable set Kalman filtering to optimize and estimate the tidal mixing parameter and seawater temperature value of the marine numerical model, more accurately estimating the tidal mixing parameter through data assimilation, reducing the value error of the tidal mixing parameter, and improving the simulation prediction accuracy of the marine numerical model for future period.

[0040] The present application will be described in detail below with reference to the accompanying drawings and in conjunction with the embodiments.

[0041] In one embodiment of the present application, a method for determining tidal mixing parameters for a marine numerical model is disclosed, comprising:

[0042] Step S101, establishing a marine numerical model, and obtaining a corresponding state vector when the marine numerical model runs to stable, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value;

[0043] For example, the calculation area of the marine numerical model is set as the Bohai Sea (37.083°N to 41.033°N, 117.52°E to 122.47°E), the horizontal direction is set as a regular rectangular grid with a resolution of 1 / 20°x1 / 20° and a grid number of 80x100; and the vertical direction is set as 6 unequal interval water layers, i.e. 0m, 5m, 15m, 25m, 35m and 65m.

[0044] In this embodiment, the marine numerical model is established, including: the water depth data of the marine numerical model adopts the GEBCO2014 seafloor topography dataset provided by the British Oceanographic Data Centre (BODC); the tidal boundary condition of the marine numerical model adopts the China Ocean ReAnalysis (CORA); the initial field of the marine numerical model is obtained from the China Ocean ReAnalysis (CORA), including the initial temperature, salinity, current and water level; and the observation data of seawater temperature at the same period is obtained from the OISST dataset of the National Oceanic and Atmospheric Administration (NOAA).

[0045] For example, the tidal mixing parameter value is represented as D, the seawater temperature value is represented as T, and the state vector is represented as X, X=[T,D];

[0046] Step S102, determining an observation matrix for Kalman filtering corresponding to the state vector;

[0047] The observation matrix of X is represented as H p , the observation matrix of T is represented as H, and H p =[H,0];

[0048] Step S103, obtaining a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector;

[0049] For example, the marine numerical model is used for prediction to obtain a predicted value at the next time:

[0050] (1);

[0051] wherein, denotes the i-th prediction value in the prediction value set at time t+1, denotes the i-th prediction value in the analysis value set at time t; M(·) is a prediction model, the prediction model is a marine numerical model; ωi is a prediction error, which conforms to an unbiased Gaussian distribution.

[0052] (2);

[0053] In formula (2), is the predicted mean value of the state vector X at time t+1, n denotes the number of elements in the prediction value set;

[0054] is the prediction error covariance matrix of the state vector X at time t+1, and is expressed as follows:

[0055] (3);

[0056] (4);

[0057] (5);

[0058] (6);

[0059] (7);

[0060] wherein, is the i-th prediction value of the state variable seawater temperature T at time t+1; is the i-th prediction value of the tidal mixing parameter D at time t+1, is the predicted mean value of the state variable seawater temperature T at time t+1, is the predicted mean value of the tidal mixing parameter D at time t+1; is the prediction error covariance matrix of the state variable seawater temperature T, and is the cross error covariance matrix of the state variables seawater temperature T and tidal mixing parameter D, is the prediction error covariance matrix of the tidal mixing parameter D.

[0061] Step S104, based on the state vector, the observation matrix and the prediction error covariance matrix, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value are calculated respectively;

[0062] For example, the prediction error covariance matrix of the state vector X is , calculate the Kalman gain;

[0063] (8);

[0064] R represents the noise of Gaussian distribution, the state vector X = [T, D], the observation matrix H p and formula (3) into formula (8), the Kalman gain of the state variable seawater temperature value T and the tidal mixing parameter D is respectively:

[0065] (9);

[0066] (10);

[0067] wherein, is the Kalman gain of the seawater temperature value T at t+1, is the Kalman gain of the tidal mixing parameter D at t+1, R represents the noise of Gaussian distribution;

[0068] Step S105, obtaining the observation value, the observation value is the observation data of the seawater temperature value;

[0069] Step S106, based on the observation value, the state vector, the observation matrix of the seawater temperature value in the state vector, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value, the tidal mixing parameter analysis value and the seawater temperature analysis value are respectively calculated;

[0070] Exemplarily, the observation data, the predicted data and the Kalman gain at t+1 are used to correct and update the state vector, and the analysis value of the state vector is obtained, as follows:

[0071] (11);

[0072] wherein, represents the i th predicted value in the analysis value set at t; represents the i th predicted value in the predicted value set at t+1, is the observation value at t+1.

[0073] Formula (9) and formula (10) are brought into formula (11), and the following formula is obtained:

[0074] (12);

[0075] (13);

[0076] wherein, represents the i th value in the analysis value set of the tidal mixing parameter at t+1, denotes the i-th value in the set of sea water temperature analysis values at time t+1;

[0077] In step S107, the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model.

[0078] Based on the above embodiment, in another embodiment of the present application, the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model, comprising: the ocean numerical model, based on the new tidal mixing parameter value, performing simulation prediction for a future period.

[0079] As shown in Figure 2 The present application also discloses a tidal mixing parameter determination device for an ocean numerical model, comprising:

[0080] The ocean numerical model establishing module 201 is configured to establish an ocean numerical model, and obtain a corresponding state vector when the ocean numerical model runs to be stable, wherein the state vector is composed of a tidal mixing parameter value and a sea water temperature value;

[0081] The first determination module 202 is configured to determine an observation matrix for Kalman filtering corresponding to the state vector;

[0082] The second determination module 203 is configured to obtain a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector;

[0083] The third determination module 204 is configured to respectively calculate a Kalman gain of the tidal mixing parameter value and a Kalman gain of the sea water temperature value based on the state vector, the observation matrix and the prediction error covariance matrix;

[0084] The fourth determination module 205 is configured to obtain an observation value, wherein the observation value is observation data of the sea water temperature value;

[0085] The fifth determination module 206 is configured to respectively calculate a tidal mixing parameter analysis value and a sea water temperature analysis value based on the observation value, the state vector, the observation matrix of the sea water temperature value in the state vector, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the sea water temperature value;

[0086] The sixth determination module 207 is configured to make the tidal mixing parameter analysis value a new tidal mixing parameter value in the ocean numerical model.

[0087] The present application also discloses an electronic device, as shown in Figure 3 An embodiment is disclosed, which is applicable to a block diagram of the electronic device for determining the tidal mixing parameter of the ocean numerical model.

[0088] The electronic device 30 of this embodiment includes a processor 301 that can perform various appropriate actions and processes in accordance with a program stored in a ROM 302 or a program loaded into a RAM 303 from a storage section 308. The processor 301 can include, for example, a general-purpose microprocessor, an instruction set processor, and / or a related chipset, and / or a special-purpose microprocessor, and so on. The processor 301 can also include an on-board memory for cache use. The processor 301 can include a single processing unit or multiple processing units for performing different actions of the method flows according to embodiments of the present application.

[0089] In the RAM 303, various programs and data required for the operation of the electronic device 30 are stored. The processor 301, the ROM 302, and the RAM 303 are connected to each other through a bus 304, and the processor 301 performs various operations of the method flows according to embodiments of the present application by executing the programs in the ROM 302 and / or the RAM 303. Note that the programs can also be stored in one or more memories other than the ROM 302 and the RAM 303, and the processor 301 can perform various operations of the method flows according to embodiments of the present application by executing the programs stored in the one or more memories.

[0090] According to embodiments of the present application, the electronic device 30 can also include an I / O interface 305 that is also connected to the bus 304. The electronic device 30 can also include one or more of the following components connected to the I / O interface 305: an input section 306 including a keyboard, a mouse, and so on; an output section 307 including a display such as a cathode ray tube, a liquid crystal display, and a speaker, and so on; a storage section 308 including a hard disk, and so on; and a communication section 309 including a network interface card such as a LAN card, a modem, and so on. The communication section 309 performs communication processing via a network such as the Internet. A drive 3010 is also connected to the I / O interface 305 as necessary. A removable recording medium 3011 such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, and so on is attached to the drive 3010 as necessary, so that a computer program read therefrom is installed in the storage section 308 as necessary.

[0091] The present application also provides a computer-readable storage medium.

[0092] The computer-readable storage medium can be included in the electronic device / apparatus system described in the above embodiments; or can exist separately from the electronic device / apparatus and not be assembled in the electronic device / apparatus. The above computer-readable storage medium carries one or more programs that, when executed, implement the method according to embodiments of the present application.

[0093] According to an embodiment of the present application, the computer readable storage medium can be a non-transitory computer readable storage medium. Examples can include, but are not limited to, a portable computer diskette, a hard disk, a random access memory (RAM), a read-only memory (ROM), an erasable programmable read-only memory (EPROM or Flash memory), a portable compact disc read-only memory (CD-ROM), an optical storage device, a magnetic storage device, or any suitable combination of the foregoing. In the present application, a computer readable storage medium can be any tangible medium that can contain, or store a program for use by or in connection with an instruction execution system, apparatus, or device.

[0094] An embodiment of the present application also includes a computer program product.

[0095] The computer program product includes a computer program containing program code for executing the method provided by the embodiments of the present application, which is used to enable the electronic device to implement the method provided by the embodiments of the present application when the computer program product is run on the electronic device.

[0096] In one embodiment, the computer program can rely on tangible storage media such as optical storage devices, magnetic storage devices, etc. In another embodiment, the computer program can also be transmitted, distributed, etc. in the form of signals on network media. The program code contained in the computer program can be transmitted in any suitable network media, including, but not limited to, wireless, wired, etc., or any suitable combination of the foregoing.

[0097] According to an embodiment of the present application, the program code for executing the computer program provided by the embodiments of the present application can be written in any combination of one or more programming languages, and specifically, these computer programs can be implemented using high-level procedural and / or object-oriented programming languages. The programming language includes, but is not limited to, such as Java, C++, python, C language or similar programming languages. The program code can be executed entirely on the user computing device, partially on the user device, partially on a remote computing device, or entirely on a remote computing device or server. In the case involving a remote computing device, the remote computing device can be connected to the user computing device through any kind of network, including a local area network or a wide area network, or can be connected to an external computing device.

[0098] The computer program product of the present application can be a computer program product comprising a computer readable storage medium and a computer program mechanism embedded in the computer readable storage medium. Such computer program product can further include a computer readable storage medium and program means for causing a processor or other programmable processing apparatus to function in a particular manner, such that the computer program mechanism that can be executed by such a processor or processing apparatus cause the processor or processing apparatus to implement the functions of the embodiments of the present application. The computer program product can be intended for use in one or more systems or a processor-controlled device.

[0099] The above-described embodiments of the application have been described in order to allow easy understanding of the application. Those of ordinary skill in the art will understand that various modifications and changes can be devised without departing from the scope of the application. It is intended that the scope of the application sought to be protected not be limited by the above description, but instead be determined by the following claims, their equivalents, and the permissible scope of the following claims.

Claims

1. A method for determining a tidal mixing parameter for a numerical ocean model, characterized in that, The method comprises: establishing a marine numerical model, obtaining a corresponding state vector when the marine numerical model runs to stability, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value; determining an observation matrix corresponding to the state vector for Kalman filtering; obtaining a prediction error covariance matrix corresponding to the state vector for Kalman filtering based on the state vector; based on the state vector, the observation matrix and the prediction error covariance matrix, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value are respectively calculated; obtaining an observation value, which is the observation data of the seawater temperature value; based on the observation value, the state vector, the observation matrix of the seawater temperature value in the state vector, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value, the tidal mixing parameter analysis value and the seawater temperature analysis value are respectively calculated; the tidal mixing parameter analysis value is the new tidal mixing parameter value in the marine numerical model.

2. The method for determining a tidal mixing parameter for a numerical ocean model according to claim 1, wherein, The water depth data of the marine numerical model adopts a seabed topography data set.

3. The method for determining the tidal mixing parameter for a numerical ocean model according to claim 1, wherein, The tidal boundary condition of the marine numerical model adopts a Chinese marine reanalysis data set.

4. The method for determining a tidal mixing parameter for a numerical ocean model according to claim 1, wherein, The initial field of the marine numerical model is obtained from the Chinese marine reanalysis data set, including initial temperature, salinity, current and water level.

5. The method for determining the tidal mixing parameter for a numerical ocean model according to claim 1, wherein, The observation value is the observation data of the seawater temperature value, including: obtaining the observation data of the seawater temperature of the same period from the data set.

6. The method for determining a tidal mixing parameter for a numerical ocean model according to claim 1, wherein, The tidal mixing parameter analysis value is the new tidal mixing parameter value in the marine numerical model, including: the marine numerical model, based on the new tidal mixing parameter value, simulates and predicts the future period.

7. A device for determining a tidal mixing parameter for a numerical ocean model, characterized in that: The device comprises: a marine numerical model establishing module for establishing a marine numerical model, obtaining a corresponding state vector when the marine numerical model runs to stability, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value; a first determining module for determining an observation matrix corresponding to the state vector for Kalman filtering; a second determining module for obtaining a prediction error covariance matrix corresponding to the state vector for Kalman filtering based on the state vector; a third determining module for calculating the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value based on the state vector, the observation matrix and the prediction error covariance matrix; a fourth determining module for obtaining an observation value, which is the observation data of the seawater temperature value; a fifth determining module for calculating the tidal mixing parameter analysis value and the seawater temperature analysis value based on the observation value, the state vector, the observation matrix of the seawater temperature value in the state vector, the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value; a sixth determining module for the tidal mixing parameter analysis value being the new tidal mixing parameter value in the marine numerical model.

8. An electronic device, comprising: comprise: one or more processors; a memory device for storing one or more programs, wherein the one or more programs, when executed by the one or more processors, cause the one or more processors to implement the method according to any one of claims 1 to 6.

9. A computer readable storage medium having stored thereon a computer program, characterized in that , the computer program, which when executed by a processor, implements the method according to any one of claims 1 to 6.

10. A computer program product comprising a computer program, characterized in that, , the computer program, which when executed by a processor, implements the method according to any one of claims 1 to 6.

Images