JP2005121520A - Marine environment analysis method/program/program recording medium/device - Google Patents

Abstract

<P>PROBLEM TO BE SOLVED: To provide a marine environment analysis device for clarifying an influence on the field of a flow exerted from a submarine topography, and resultantly showing clearly a physical process of vortex generation. <P>SOLUTION: In this marine environment analysis device, optionally set submarine topographical data are calculated by a submarine topographical data input part 1, and space-time data of the marine environment such as a water temperature, a salt content, or a flow velocity are calculated by a numerical value marine model part 2. Potential energy and kinetic energy are calculated by an energy calculation part 3, and a conversion rate between the energies is calculated by an energy conversion rate calculation part 4. Differential values between the energies and the energy conversion rates calculated from the two different submarine topographical data are calculated by an energy parameter difference calculation part 5. The energy parameter is displayed by an energy parameter display part 6. Thus, an energy distribution and an energy conversion process by the submarine topography are grasped relatively, to thereby clarify the physical process of the vortex generation. <P>COPYRIGHT: (C)2005,JPO&NCIPI

Description

  The present invention relates to a marine environment analysis technique for clarifying the interaction between a medium-scale eddy generated around Japan and seafloor topography.

  In ocean physics research, it is said that the interaction between medium-scale vortices and submarine irregularities that occur in the sea area around Japan contributes to the fluctuations of the large meandering observed in the western boundary current such as the Kuroshio Current (non-patent literature) 1).

Conventionally, when analyzing the interaction between a medium-scale eddy and seafloor topography, various seafloor topography of a numerical ocean model (see Non-Patent Documents 2 and 3) can be changed in various ways, and the flow field can be calculated using the numerical ocean model. Trial-and-error methods such as discussing the similarity between the results and observations were taken.
T.A. Endoh and T.W. Hibiya, “Numerical simulation of the transient response of the Kuroshio leading to the large format of south of Japan”, J. et al. Geophys. Res. 106, 26833-26850 (2001). G. L. Mellor, “Users guide for a three-dimensional, primitive equation, numerical ocean model (Jully 1998 version)”, Rep. , Program in Atomospheric and Ocean Science, Princeton University, pp. 41 (1998). G. L. Mellor, T .; Ezer, and L.L. -Y. Oey, “The pres sure gradient condrum of sigma coordinated ocean models”, J. Am. Atoms. Oceanic Technol. 11, 1126-1134 (1994).

  However, the conventional method clarifies how the flow field is affected by the seabed topography, and as a result, there is a problem that the physical process of vortex generation cannot be clearly shown.

  Accordingly, an object of the present invention is to clarify how the flow field is affected by the seafloor topography, and as a result, to provide a marine environment analysis method capable of clearly showing the physical process of vortex generation. is there.

  In order to achieve the above object, the marine environment analysis method of the present invention arbitrarily sets the bottom shape of a numerical ocean model, and calculates spatiotemporal data of the marine environment such as water temperature, salinity, and flow velocity based on the numerical ocean model. To do. Based on the spatio-temporal data of the marine environment, potential energy, kinetic energy, and energy conversion rate are calculated. For these parameters, the difference between the different seafloor topography is taken and the difference result is displayed.

  In the present invention, the potential energy, the kinetic energy, the energy conversion rate, and the difference value of these energy parameters regarding the two seafloor topologies are calculated for the calculation results of the numerical ocean model executed based on two arbitrarily set seafloor topologies. By displaying as above, it is possible to relatively grasp the energy distribution and the physical process of energy conversion due to the seabed topography.

  As described above, according to the present invention, it is possible to relatively grasp the energy distribution and conversion process due to the seabed topography and clarify the physical process of vortex generation. Specifically, it is as follows.

  Since the time change of the spatial change tendency of the energy parameter can be displayed, a region where the energy and the energy conversion rate are remarkable can be confirmed.

  Since the temporal change of the spatial average value of the energy parameter can be displayed, the cause of the energy change due to the season can be estimated.

  Since the time change of the energy balance can be confirmed, it can be confirmed which conversion from which energy to which energy is remarkable.

  Regarding the energy parameters, the difference between the two different seafloor topography by the numerical ocean model can be confirmed, so that the characteristics of the seafloor topography contributing to energy conversion can be estimated.

  Next, embodiments of the present invention will be described in detail with reference to the drawings.

  Referring to FIG. 1, a marine environment analysis apparatus according to an embodiment of the present invention includes a seafloor topography data input unit 1, a numerical ocean model unit 2, an energy calculation unit 3, an energy conversion rate calculation unit 4, an energy parameter. A difference calculation unit 5 and an energy parameter display unit 6 are included.

  The seafloor terrain data input unit 1 inputs / smooths an operator input or existing numerical seafloor terrain data. Then, the seafloor topography data 11 is output to the numerical ocean model unit 2.

  The numerical ocean model unit 2 receives the seafloor topography data 11 and calculates spatiotemporal data of the marine environment such as water temperature, salinity, and flow velocity. Then, the marine environment data 12 is output to the energy calculation unit 3 and the energy conversion rate calculation unit 4.

  The energy calculation unit 3 calculates potential energy and kinetic energy in the mean field and vortex field based on the marine environment data 12 calculated by the numerical ocean model unit 2. Then, the calculated energy 13 is output to the energy conversion rate calculation unit 4, the energy parameter difference calculation unit 5, and the energy parameter display unit 6.

  The energy conversion rate calculation unit 4 calculates an energy conversion rate related parameter based on the marine environment data 12 calculated by the numerical ocean model unit 2 and the energy 13 by the energy calculation unit 3. Then, the calculated energy conversion rate 14 is output to the energy parameter difference calculation unit 5 and the energy parameter display unit 6.

  The energy parameter difference calculation unit 5 calculates a difference value between the parameters calculated by the energy calculation unit 3 and the energy conversion rate calculation unit 4 based on two different seabed topography data. Then, the calculated energy difference value 15 and energy conversion rate difference value 16 are output to the energy parameter display unit 6.

  The energy parameter display unit 6 displays the time series data of the spatial distribution and the spatial average value of each parameter calculated by the energy calculation unit 3, the energy conversion rate calculation unit 4, and the energy parameter difference calculation unit 5.

  Next, the operation of the marine environment analyzing apparatus of this embodiment will be described in detail with reference to the drawings.

  The submarine terrain data input unit 1 inputs / smooths an operator input or existing submarine terrain data. FIG. 2 shows a specific method of inputting seafloor topographic data by the operator. There are a method of inputting depth data by designating an arbitrary point on a calculation grid on a two-dimensional screen, and a method of inputting depth data under a line segment connecting two arbitrarily designated points. Also, existing numerical seafloor topographic data is input electronically and smoothed. In smoothing, Non-Patent Document 3 (GL Leller, T. Ezer, and L.-Y. -1134 (1994).) Is set so as to satisfy the criteria of the expression (1), and the calculation by the numerical ocean model unit 2 does not hinder the calculation.

ここで、Ｈ_i+1およびＨ_iは、近接する２点間における深度データ、αは傾斜ファクターである。

Here, H _{i + 1} and H _i are depth data between two adjacent points, and α is a slope factor.

  The numerical ocean model unit 2 calculates marine environment data such as water temperature, salinity, and flow velocity for three-dimensional space and arbitrary time based on numerical ocean models such as Princeton Ocean Model and Modular Ocean Model. Set above. The principle of calculation is to numerically calculate the equation of motion of seawater on the rotating earth, the mass conservation equation, the calorie conservation equation, and the salinity conservation equation. For details, see Non-Patent Document 2 (GL Mellor, “Users Guide for a three-dimensional, primitive equation, numeric ocean model (July 1998 version) ”, Rep., Program in Atomic and Oceanic 41. ing.

Princeton Ocean Model
http: // www. aos. Princeton. edu / WWWPUBLIC / htdocs. pom
Modular Ocean Model
http: // www. gfdl. gov / ˜smg / MOM / MON. html.

  The energy calculation unit 3 calculates the next energy based on the spatiotemporal data calculated by the numerical ocean model unit 2.

Mean field potential energy (Pm)
Potential energy of vortex field (Pe)
Mean field kinetic energy (km)
Vortex kinetic energy (Ke).

  Moreover, the energy calculation part 3 calculates the space average value of each said energy related parameter for every time by Formula (2).

ここで、ｉは経度方向の格子点番号、ｊは緯度方向の格子点番号、Ｘ（ｉ，ｊ）はエネルギーパラメータ、Ｍは経度方向の格子点数、Ｎは緯度方向の格子点数を示す。

Here, i is a lattice point number in the longitude direction, j is a lattice point number in the latitude direction, X (i, j) is an energy parameter, M is the number of lattice points in the longitude direction, and N is the number of lattice points in the latitude direction.

  The energy conversion rate calculation unit 4 calculates the following energy conversion rate related parameters for each latitude, longitude, and time based on the spatio-temporal data by the numerical ocean model unit 2 and each energy item by the energy calculation unit 3.

Pm → Km conversion rate (PmKm)
Pm → Pe conversion rate (PmPe)
Km → Ke conversion rate (KmKe)
Pe → Ke conversion rate (PeKe)
Average field energy dissipation (Dm)
Energy dissipation of vortex field (De)
Mean field energy source by external force (Fm)
Energy source of vortex field by external force (Fe)
Other conversion rate related parameters.

  Moreover, the energy conversion rate calculation part 4 calculates the space average value of said each energy conversion rate related parameter for every time. Calculation is similar to equation (2).

  The energy parameter difference calculation unit 5 calculates a difference value between the parameters calculated by the energy calculation unit 3 and the energy conversion rate calculation unit 4 based on two different seabed topography data. In calculating the difference value, the energy parameter is calculated for each horizontal grid according to the equation (3).

ここで、Ｘ₁（ｉ，ｊ）およびＸ₂（ｉ，ｊ）は、異なる２つの海底地形データにもとづいて算出されたエネルギーパラメータを示す。

Here, X ₁ (i, j) and X ₂ (i, j) represent energy parameters calculated based on two different seafloor topographic data.

  The energy parameter display unit 6 displays the time series data of the spatial distribution and the spatial average value of each parameter calculated by the energy calculation unit 3, the energy conversion rate calculation unit 4, and the energy parameter difference calculation unit 5.

  FIG. 3 shows a display example of a spatial distribution of energy parameters or energy parameter difference values. By setting the energy parameter name, longitude, latitude, and time point by operator operation, the spatial distribution of the energy parameter is displayed as longitude on the horizontal axis and latitude on the vertical axis at the set time point. When there are a plurality of set time points, it is also possible to display the spatial distribution of energy parameters at the plurality of time points by animation. Thus, since the temporal change of the spatial change tendency of the energy parameter can be displayed, a region where the energy and the energy conversion rate are remarkable can be confirmed.

  FIG. 4 shows a time-series display example of the spatial average value of energy parameters or energy parameter difference values. By setting the energy parameter name and the time point by the operator operation, the time series data of the spatial average value of the energy parameter at the set time point is displayed. Thus, since the temporal change of the spatial average value of the energy parameter can be displayed, the cause of the energy change due to the season can be estimated.

  FIG. 5 shows a block diagram of the energy balance. By setting the time point by the operator's operation, a block diagram of the energy balance at the set time point is displayed. Thus, since the time change of the energy balance can be confirmed, it can be confirmed which conversion from which energy to which energy is remarkable.

  The ocean environment analysis apparatus records a program for realizing its function on a computer-readable recording medium in addition to being realized by dedicated hardware, and the program recorded on this recording medium is recorded on the ocean. It may be realized by being read and executed by a computer system to be an environment analysis apparatus. The computer-readable recording medium refers to a recording medium such as a floppy disk, a magneto-optical disk, a CD-ROM, or a storage device such as a hard disk device built in a computer system. Furthermore, a computer-readable recording medium is a medium that dynamically holds a program (transmission medium or transmission wave) in a short time, as in the case of transmitting a program via the Internet. Such as a volatile memory, which holds a program for a certain period of time.

It is a figure which shows the structure of the marine environment analysis apparatus of one Embodiment of this invention. It is an example of the input of the seabed topography data by an operator. It is an example of a display of the spatial distribution of an energy parameter or an energy parameter difference value. It is an example of a time series display of the spatial average value of an energy parameter or an energy parameter difference value. It is a block diagram of energy balance.

Explanation of symbols

DESCRIPTION OF SYMBOLS 1 Submarine topography data input part 2 Numerical ocean model part 3 Energy calculation part 4 Energy conversion rate calculation part 5 Energy parameter difference calculation part 6 Energy parameter display part 11 Submarine landform data 12 Marine environment data 13 Energy 14 Energy conversion rate 15 Energy difference value 16 Energy conversion rate difference value 21 Input arbitrary point of seafloor topographic data 22 Line segment input of seafloor topographic data 31 Spatial distribution of energy parameter 41 Time series of spatial average value of energy parameter 51 Kinetic energy of mean field Km
52 Potential field potential energy Pm
53 Kinetic energy Ke of vortex field
54 Potential energy Pe of vortex field
501 Pm → Km conversion rate 502 Km → Ke conversion rate 503 Pm → Pe conversion rate 504 Pe → Ke conversion rate 505 Average field energy dissipation Dm
506 Energy dissipation De of vortex field
507 Energy source of mean field by external force 508 Energy source of vortex field by external force

Claims (8)

A computer clarification of how the flow field is affected by seafloor topography, and as a result, a marine environment analysis method that clearly shows the physical process of eddy generation,
A first step wherein the computer accepts input of two different seafloor terrain data;
A second step in which the computer applies marine topography data received in the first step to a numerical ocean model to generate marine environment data for an arbitrary spatial position and an arbitrary time;
A third step in which the computer calculates various energies for an arbitrary spatial position and an arbitrary time based on the marine environment data generated in the second step;
A fourth step in which the computer calculates difference data for the two different seafloor topologies for the various energies calculated in the third step;
A method comprising a fifth step in which the computer causes the display means to display the various energies calculated in the third step and the difference data calculated in the fourth step. The computer calculates a conversion rate between energies for an arbitrary spatial position and an arbitrary time based on the marine environment data generated in the second step and the various energies calculated in the third step. A sixth step,
A seventh step in which the computer calculates difference data for the two different submarine topography of the conversion rate between the energies calculated in the sixth step;
The computer further includes an eighth step of causing the display means to display the conversion rate between the energy calculated in the sixth step and the difference data calculated in the seventh step. The method described in 1. A ninth step in which the computer calculates spatial average data at any time of the various energies calculated in the third step and / or the conversion rate between the energies calculated in the sixth step;
A tenth step in which the computer calculates difference data of the spatial average data calculated in the ninth step with respect to the two different seabed topography;
The computer further comprises an eleventh step of causing the display means to display the spatial average data calculated in the ninth step and the difference data calculated in the tenth step. The method described in 1. The program which makes a computer perform each step of any one of Claim 1 to 3. A computer-readable recording medium in which a program for causing a computer to execute each step according to any one of claims 1 to 3 is recorded. It is an ocean environment analysis device that clarifies how the flow field is affected by the seafloor topography, and as a result, clearly shows the physical process of vortex generation,
A first means for accepting input of two different seafloor terrain data;
Applying the seafloor topography data received by the first means to a numerical ocean model to generate marine environment data for an arbitrary spatial position and an arbitrary time;
Third means for calculating various energies for an arbitrary spatial position and an arbitrary time based on the marine environment data generated by the second means;
A fourth means for calculating differential data relating to the two different seafloor topography for the various energies calculated by the third means;
An apparatus having a fifth means for causing the display means to display various energies calculated by the third means and difference data calculated by the fourth means. A sixth means for calculating a conversion rate between the energy for an arbitrary spatial position and an arbitrary time based on the marine environment data generated by the second means and the various energies calculated by the third means; In addition,
The fourth means further calculates difference data regarding the two different seafloor topography, of the conversion rate between the energy calculated by the sixth means,
The apparatus according to claim 6, wherein the fifth means further causes the display means to display the conversion rate between the energy calculated by the sixth means and the difference data calculated by the fourth means. The third means further calculates spatial average data at various times of the calculated various energies,
The sixth means further calculates spatial average data at an arbitrary time of the conversion rate between the calculated energies,
The fourth means further calculates difference data for the two different seafloor topologies of the spatial average data calculated by the third means and / or the sixth means;
The fifth means further causes the display means to display the spatial average data calculated by the third means and / or the sixth means and the difference data calculated by the fourth means. The apparatus according to 6 or 7.

Images