Tidal mixing parameter determination method and device for ocean numerical model and product

Abstract

The invention discloses a tidal mixing parameter determination method and device for an ocean numerical model and a product, and the method comprises the steps: building the ocean numerical model, and obtaining 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 the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value; respectively calculating to obtain 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 ocean numerical model. According to the method, the tide mixing parameters and the seawater temperature value of the ocean numerical model are optimized and estimated through the combined state variable set Kalman filtering, the tide mixing parameters are estimated more accurately through data assimilation, the value error of the tide mixing parameters is reduced, and the simulation prediction precision of the ocean numerical model for the future time period is improved.

Description

Technical Field

[0001] The present invention belongs to the field of ocean data assimilation, and in particular relates to a method, device and product for determining tidal mixing parameters for an ocean numerical model. Background Art

[0002] Tidal mixing is a constant phenomenon in the ocean, controlling the material exchange of seawater and exerting a significant influence on ocean circulation, temperature, and salinity. However, tidal mixing cannot be observed. In ocean numerical models, the strength of tidal mixing is usually quantified by the size of tidal mixing parameters. These parameters are generally given by empirical values ​​or calculated using simulated data from ocean numerical models, and therefore often exhibit large errors compared to actual conditions. Summary of the Invention

[0003] In view of this, the present invention aims to overcome the defects in the prior art and proposes a method, device and product for determining tidal mixing parameters for ocean numerical models.

[0004] To achieve the above object, the technical solution of the present invention is achieved as follows: The present invention discloses a method for determining tidal mixing parameters for an ocean numerical model, comprising: establishing an ocean numerical model, and obtaining a corresponding state vector when the ocean numerical model runs to stability, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value; Determine the observation matrix for Kalman filtering corresponding to the state vector; Based on the state vector, the prediction error covariance matrix for Kalman filtering corresponding to the state vector is obtained; Based on the state vector, observation matrix and 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; Obtaining observation values, which are observation data of seawater temperature values; 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 calculated respectively; The tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model.

[0005] In one implementation of the present invention, establishing an ocean numerical model includes: using a seabed topography dataset as water depth data of the ocean numerical model.

[0006] In one implementation of the present invention, an ocean numerical model is established, including: using a China Ocean Reanalysis dataset as a tidal boundary condition of the ocean numerical model.

[0007] In one implementation of the present invention, the initial fields of the ocean numerical model are obtained from the China Ocean Reanalysis Dataset, including initial temperature, salinity, ocean current, and water level.

[0008] In one implementation of the present invention, the observation value is observation data of seawater temperature, including: obtaining the observation data of seawater temperature in the same period from a data set.

[0009] In one implementation of the present invention, the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model, including: the ocean numerical model performs simulation prediction for future time periods based on the new tidal mixing parameter value.

[0010] In a second aspect, the present invention discloses a device for determining tidal mixing parameters for an ocean numerical model, the device comprising: The ocean numerical model establishment module is used to establish the ocean numerical model. When the ocean numerical model runs to stability, the corresponding state vector is obtained, where the state vector is composed of tidal mixing parameter values ​​and seawater temperature values; A first determination module is used to determine an observation matrix for Kalman filtering corresponding to a state vector; A second determination module is used to obtain a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector; A third determination module is used to calculate 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 determination module is used to obtain an observation value, where the observation value is observation data of the seawater temperature value; a fifth determination module, configured to 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; The sixth determination module, used for analyzing the tidal mixing parameter value, is a new tidal mixing parameter value in the ocean numerical model.

[0011] In a third aspect, the present invention discloses an electronic device comprising: one or more processors; a storage device for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors implement the above method.

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

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

[0014] Compared with the prior art, the present invention has the following advantages: The present invention discloses a method, device, and product for determining tidal mixing parameters for an ocean numerical model, comprising: establishing an ocean numerical model, obtaining a corresponding state vector when the ocean numerical model runs to stability; 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 the Kalman gain of the tidal mixing parameter value and the 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 ocean numerical model. The present invention discloses a method, device, and product for determining tidal mixing parameters for an ocean numerical model, utilizing a joint state variable set Kalman filter to optimize and estimate the tidal mixing parameters and seawater temperature values ​​of the ocean numerical model, and through data assimilation, more accurately estimating the tidal mixing parameters, reducing the error in the tidal mixing parameter value, and improving the simulation and prediction accuracy of the ocean numerical model for future time periods. BRIEF DESCRIPTION OF THE DRAWINGS

[0015] The accompanying drawings, which constitute a part of the present invention, are used to provide a further understanding of the present invention. The exemplary embodiments of the present invention and their descriptions are used to explain the present invention and do not constitute improper limitations on the present invention.

[0016] In the attached figure: Figure 1 This is a schematic diagram of a method for determining tidal mixing parameters for an ocean numerical model according to an embodiment of the present invention; Figure 2 This is a schematic diagram of a device for determining tidal mixing parameters for an ocean numerical model according to an embodiment of the present invention; Figure 3 The figure is a schematic diagram of an electronic device for determining tidal mixing parameters of an ocean numerical model according to an embodiment of the present invention. DETAILED DESCRIPTION

[0017] It should be noted that, in the absence of conflict, the embodiments of the present invention and the features in the embodiments may be combined with each other.

[0018] In the description of the present invention, it should be understood that the terms "center", "longitudinal", "lateral", "up", "down", "front", "back", "left", "right", "vertical", "horizontal", "top", "bottom", "inside", "outside", etc., indicating the orientation or position relationship, are based on the orientation or position relationship shown in the accompanying drawings, and are only for the convenience of describing the present invention and simplifying the description, rather than indicating or implying that the device or element referred to must have a specific orientation, be constructed and operate in a specific orientation, and therefore cannot be understood as limiting the present invention.

[0019] In the description of the present invention, it should be noted that, unless otherwise expressly specified or limited, the terms "mounted," "connected," and "connected" should be understood in a broad sense. For example, they may refer to fixed connections, detachable connections, or integral connections; mechanical connections or electrical connections; direct connections or indirect connections through an intermediate medium; and internal communication between two components. Those skilled in the art will understand the specific meanings of the above terms in the present invention based on specific circumstances.

[0020] In the description of the present invention, it should be further clarified that the terms "first," "second," etc. are used for descriptive purposes only and should not be understood to indicate or imply relative importance or implicitly specify the number of the technical features indicated. Therefore, a feature specified as "first," "second," etc. may explicitly or implicitly include one or more of the features. In the description of the present invention, unless otherwise specified, "plurality" means two or more.

[0021] Under existing technical conditions, the strength of tidal mixing in ocean numerical models is typically quantified using the magnitude of tidal mixing parameters. The values ​​of these parameters are typically given by empirical values ​​or calculated using simulated data from the ocean numerical model, often resulting in significant errors compared to actual conditions. The present invention discloses a method, device, and product for determining tidal mixing parameters for ocean numerical models, including optimizing and estimating the tidal mixing parameters and seawater temperature values ​​of the ocean numerical model using a joint state variable set Kalman filter. Through data assimilation, the tidal mixing parameters are more accurately estimated, reducing errors in the tidal mixing parameter values ​​and improving the simulation and prediction accuracy of the ocean numerical model for future time periods.

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

[0023] In one embodiment of the present invention, the present invention discloses a method for determining tidal mixing parameters for an ocean numerical model, comprising: Step S101: Establish an ocean numerical model. When the ocean numerical model runs to stability, obtain a corresponding state vector, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value. For example, the calculation area of ​​the ocean numerical model is set to the Bohai Sea (37.083°N to 41.033°N, 117.52°E to 122.47°E), the horizontal direction is set to a regular rectangular grid with a resolution of 1 / 20°×1 / 20° and 80×100 grids; six unequally spaced water layers are set vertically, namely 0m, 5m, 15m, 25m, 35m and 65m.

[0024] In this embodiment, an ocean numerical model is established, including: using the GEBCO2014 seafloor topography dataset provided by the British Oceanographic Data Centre (BODC) as the depth data for the ocean numerical model. The China Ocean ReAnalysis (CORA) dataset is used as the tidal boundary conditions for the ocean numerical model. The initial field of the ocean numerical model is obtained from the CORA dataset, including initial temperature, salinity, ocean currents, and water level. Observational seawater temperature data for the same period is obtained from the OISST dataset of the National Oceanic and Atmospheric Administration (NOAA).

[0025] For example, the tidal mixing parameter value is represented by D, the seawater temperature value is represented by T, and the state vector is represented by X, X=[T,D]; Step S102, determining the measurement matrix for Kalman filtering corresponding to the state vector; The observation matrix of X is represented by H p , the observation matrix of T is represented as H, H p =[H,0]; Step S103: Based on the state vector, obtain the prediction error covariance matrix for Kalman filtering corresponding to the state vector; For example, the ocean numerical model is used to make predictions and obtain the predicted value at the next moment: (1); in, represents the i-th predicted value in the set of predicted values ​​at time t+1, represents the i-th predicted value in the analysis value set at time t; M(*) is the prediction model, which is an ocean numerical model; ωi is the prediction error, which conforms to the unbiased Gaussian distribution.

[0026] (2); In formula (2), is the predicted mean of the state vector X at time t+1, n Indicates the number of elements in the predicted value set; is the prediction error covariance matrix of the state vector X at time t+1, expressed as follows: (3);

[0027] (4);

[0028] (5);

[0029] (6);

[0030] (7); in, is the i-th predicted value of the state variable seawater temperature T at time t+1; is the i-th predicted value of the tide 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 sea temperature T and tidal mixing parameter D, is the prediction error covariance matrix of the tide mixing parameter D.

[0031] 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; For example, the prediction error covariance matrix of the state vector X is used , calculate the Kalman gain; (8); R represents the noise of Gaussian distribution, the state vector X=[T,D], the observation matrix H p And by substituting formula (3) into formula (8), we can obtain the Kalman gains of the state variables seawater temperature T and tidal mixing parameter D: (9); (10); in, is the Kalman gain of the seawater temperature T at time t+1, is the Kalman gain of the tide mixing parameter D at time t+1, and R represents the noise of Gaussian distribution; Step S105, obtaining an observation value, which is observation data of the seawater temperature value; Step S106, 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; Exemplarily, the state vector is corrected and updated using the observed data, predicted data, and Kalman gain at time t+1 to obtain the analytical value of the state vector as follows: (11); in, represents the i-th predicted value in the analysis value set at time t; represents the i-th predicted value in the set of predicted values ​​at time t+1, is the observed value at time t+1.

[0032] Substituting formula (9) and formula (10) into formula (11), we get: (12); (13); in, represents the i-th value in the tidal mixing parameter analysis value set at time t+1, represents the i-th value in the set of seawater temperature analysis values ​​at time t+1; Step S107: The tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model.

[0033] Based on the previous embodiment, in another embodiment of the present invention, the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model, including: the ocean numerical model performs simulation predictions for future time periods based on the new tidal mixing parameter value.

[0034] like Figure 2 As shown, the present invention also discloses a device for determining tidal mixing parameters for an ocean numerical model, comprising: The ocean numerical model establishment module 201 is used to establish an ocean numerical model. When the ocean numerical model runs to a stable state, a corresponding state vector is obtained, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value. A first determining module 202 is used to determine an observation matrix for Kalman filtering corresponding to a state vector; A second determining module 203 is configured to obtain a prediction error covariance matrix for Kalman filtering corresponding to the state vector based on the state vector; The third determination module 204 is configured to calculate 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; The fourth determining module 205 is used to obtain an observation value, which is observation data of the seawater temperature value; A fifth determination module 206 is configured to 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. The sixth determination module 207 is used to determine whether the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model.

[0035] The present invention also discloses an electronic device, such as Figure 3 As shown, an embodiment is disclosed, which is a block diagram of an electronic device suitable for determining tidal mixing parameters of the above-mentioned ocean numerical model.

[0036] The electronic device 30 of this embodiment includes a processor 301, which can perform various appropriate actions and processes according to the program stored in the ROM 302 or the program loaded from the storage part 308 into the RAM 303. The processor 301 may include, for example, a general-purpose microprocessor, an instruction set processor and / or a related chipset and / or a dedicated microprocessor, etc. The processor 301 may also include onboard memory for caching purposes. The processor 301 may include a single processing unit or multiple processing units for performing different actions of the method flow according to the embodiment of the present invention.

[0037] RAM 303 stores various programs and data required for the operation of electronic device 30. Processor 301, ROM 302, and RAM 303 are interconnected via bus 304. Processor 301 executes the programs in ROM 302 and / or RAM 303 to perform various operations according to the method flow of the embodiment of the present invention. It should be noted that the programs may also be stored in one or more memories other than ROM 302 and RAM 303, and processor 301 may also execute the programs stored in one or more memories to perform various operations according to the method flow of the embodiment of the present invention.

[0038] According to an embodiment of the present invention, the electronic device 30 may further include an I / O interface 305, which is also connected to the bus 304. The electronic device 30 may further include one or more of the following components connected to the I / O interface 305: an input unit 306 including a keyboard, a mouse, etc.; an output unit 307 including a cathode ray tube, a liquid crystal display, and a speaker; a storage unit 308 including a hard disk; and a communication unit 309 including a network interface card such as a LAN card or a modem. The communication unit 309 performs communication processing via a network such as the Internet. A drive 3010 is also connected to the I / O interface 305 as needed. Removable media 3011, such as a magnetic disk, an optical disk, a magneto-optical disk, a semiconductor memory, etc., is installed in the drive 3010 as needed, so that computer programs read therefrom can be installed into the storage unit 308 as needed.

[0039] The present invention also provides a computer-readable storage medium.

[0040] The computer-readable storage medium may be included in the electronic device / device system described in the above embodiments, or may exist independently and not be incorporated into the electronic device / device. The computer-readable storage medium carries one or more programs that, when executed, implement the method according to the embodiments of the present invention.

[0041] According to an embodiment of the present invention, a computer-readable storage medium may be a non-volatile computer-readable storage medium. Examples include, but are not limited to, a portable computer disk, a hard disk, random access memory (RAM), read-only memory (ROM), erasable programmable read-only memory (EPROM) or flash memory, a portable compact disk read-only memory (CDROM), an optical storage device, a magnetic storage device, or any suitable combination thereof. In the present invention, a computer-readable storage medium may be any tangible medium containing or storing a program that can be used by or in conjunction with an instruction execution system, apparatus, or device.

[0042] Embodiments of the present invention also include a computer program product.

[0043] The computer program product includes a computer program, which contains program code for executing the method provided by the embodiment of the present invention. When the computer program product runs on an electronic device, the program code is used to enable the electronic device to implement the method provided by the embodiment of the present invention.

[0044] In one embodiment, the computer program may be stored on a tangible storage medium such as an optical storage device or a magnetic storage device. In another embodiment, the computer program may be transmitted and distributed in the form of a signal over a network medium. The program code contained in the computer program may be transmitted using any appropriate network medium, including but not limited to wireless, wired, or any suitable combination thereof.

[0045] According to an embodiment of the present invention, the program code for executing the computer program provided by the embodiment of the present invention can be written by any combination of one or more programming languages. Specifically, these computing programs can be implemented using high-level procedural and / or object-oriented programming languages. Programming languages ​​include, but are not limited to, 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 of a remote computing device, the remote computing device can be connected to the user computing device through any type of network, including a local area network or a wide area network, or can be connected to an external computing device.

[0046] The flowcharts and block diagrams in the accompanying drawings illustrate the possible architectures, functions, and operations of the systems, methods, and computer program products according to various embodiments of the present invention. In this regard, each box in the flowchart or block diagram may represent a module, program segment, or portion of code, which contains one or more executable instructions for implementing the specified logical function. It should also be noted that in some alternative implementations, the functions marked in the boxes may occur in an order different from that marked in the accompanying drawings. For example, two boxes shown in succession may actually be executed substantially in parallel, or they may sometimes be executed in the opposite order, depending on the functions involved. It should also be noted that each box in the block diagram or flowchart, as well as the combination of boxes in the block diagram or flowchart, may be implemented using a dedicated hardware-based system that performs the specified function or operation, or may be implemented using a combination of dedicated hardware and computer instructions. It will be understood by those skilled in the art that the features described in the various embodiments of the present invention may be combined and / or coupled in various ways, even if such combinations or couplings are not explicitly described in the present invention. In particular, the features described in the various embodiments of the present invention may be combined and / or coupled in various ways without departing from the spirit and teachings of the present invention. All such combinations and / or combinations fall within the scope of the present invention.

[0047] The embodiments of the present invention have been described above. However, these embodiments are for illustrative purposes only and are not intended to limit the scope of the present invention. Although each embodiment has been described above separately, this does not mean that the measures in each embodiment cannot be advantageously used in combination. The scope of the present invention is defined by the appended claims and their equivalents. Without departing from the scope of the present invention, those skilled in the art may make various substitutions and modifications, which should all fall within the scope of the present invention.

Claims

1. A method for determining tidal mixing parameters for an ocean numerical model, characterized in that: The method includes: establishing an ocean numerical model, and when the ocean numerical model runs to stability, obtaining a corresponding state vector, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value; Determine an observation matrix for Kalman filtering corresponding to the state vector; Based on the state vector, obtaining a prediction error covariance matrix for Kalman filtering corresponding to the state vector; Based on the state vector, the observation matrix and the prediction error covariance matrix, respectively calculating the Kalman gain of the tidal mixing parameter value and the Kalman gain of the seawater temperature value; Obtaining an observation value, wherein the observation value is observation data of the seawater temperature value; 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, a Kalman gain of the tidal mixing parameter value, and a Kalman gain of the seawater temperature value, respectively; The tidal mixing parameter analysis value is the new tidal mixing parameter value in the ocean numerical model.

2. A method for determining tidal mixing parameters for an ocean numerical model according to claim 1, characterized in that: The establishing of the ocean numerical model includes: using a seabed topography data set as the water depth data of the ocean numerical model.

3. The method for determining tidal mixing parameters for an ocean numerical model according to claim 1, wherein: The establishing of the ocean numerical model includes: using the China Ocean Reanalysis Dataset as the tidal boundary condition of the ocean numerical model.

4. The method for determining tidal mixing parameters for an ocean numerical model according to claim 1, wherein: The initial fields of the ocean numerical model are obtained from the China Ocean Reanalysis dataset, including initial temperature, salinity, ocean currents, and water level.

5. The method for determining tidal mixing parameters for an ocean numerical 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 in the same period from the data set.

6. The method for determining tidal mixing parameters for an ocean numerical model according to claim 1, wherein: The tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model, including: the ocean numerical model performs simulation prediction for future time periods based on the new tidal mixing parameter value.

7. A device for determining tidal mixing parameters for an ocean numerical model, characterized by: The device comprises: An ocean numerical model establishment module is used to establish an ocean numerical model. When the ocean numerical model runs to stability, a corresponding state vector is obtained, wherein the state vector is composed of a tidal mixing parameter value and a seawater temperature value; A first determining module, configured to determine an observation matrix for Kalman filtering corresponding to the state vector; A second determining module is configured to obtain, based on the state vector, a prediction error covariance matrix for Kalman filtering corresponding to the state vector; a third determination module, configured to calculate 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, configured to obtain an observation value, wherein the observation value is observation data of the seawater temperature value; a fifth determination module, configured to 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; A sixth determination module is used to determine that the tidal mixing parameter analysis value is a new tidal mixing parameter value in the ocean numerical model.

8. An electronic device, characterized in that: include: one or more processors; A storage device for storing one or more programs, wherein when the one or more programs are executed by the one or more processors, the one or more processors are enabled to implement the method according to any one of claims 1 to 6.

9. A computer-readable storage medium having a computer program stored thereon, characterized in that When the computer program is executed by a processor, it implements the method described in any one of claims 1 to 6.

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

Images