JP7116852B1 - Parallel coupled ocean acoustic prediction system and method of operation - Google Patents
Parallel coupled ocean acoustic prediction system and method of operation Download PDFInfo
- Publication number
- JP7116852B1 JP7116852B1 JP2022013517A JP2022013517A JP7116852B1 JP 7116852 B1 JP7116852 B1 JP 7116852B1 JP 2022013517 A JP2022013517 A JP 2022013517A JP 2022013517 A JP2022013517 A JP 2022013517A JP 7116852 B1 JP7116852 B1 JP 7116852B1
- Authority
- JP
- Japan
- Prior art keywords
- ocean
- sound field
- acoustic
- module
- environment
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 46
- 238000004364 calculation method Methods 0.000 claims abstract description 54
- 230000008878 coupling Effects 0.000 claims abstract description 16
- 238000010168 coupling process Methods 0.000 claims abstract description 16
- 238000005859 coupling reaction Methods 0.000 claims abstract description 16
- 230000003746 surface roughness Effects 0.000 claims abstract description 11
- 238000000605 extraction Methods 0.000 claims abstract description 7
- 230000001131 transforming effect Effects 0.000 claims abstract description 4
- 230000009471 action Effects 0.000 claims abstract description 3
- 230000008569 process Effects 0.000 claims description 23
- 230000007613 environmental effect Effects 0.000 claims description 6
- 238000012935 Averaging Methods 0.000 claims description 3
- 238000012545 processing Methods 0.000 claims description 2
- 239000000284 extract Substances 0.000 claims 1
- 230000006870 function Effects 0.000 abstract description 4
- 230000000694 effects Effects 0.000 abstract description 3
- 230000004927 fusion Effects 0.000 abstract description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 11
- 239000011159 matrix material Substances 0.000 description 5
- 238000011160 research Methods 0.000 description 4
- 238000010586 diagram Methods 0.000 description 2
- 230000010354 integration Effects 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 230000008859 change Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000004088 simulation Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B17/00—Monitoring; Testing
- H04B17/30—Monitoring; Testing of propagation channels
- H04B17/391—Modelling the propagation channel
- H04B17/3912—Simulation models, e.g. distribution of spectral power density or received signal strength indicator [RSSI] for a given geographic region
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
- G06F17/11—Complex mathematical operations for solving equations, e.g. nonlinear equations, general mathematical optimization problems
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B11/00—Transmission systems employing sonic, ultrasonic or infrasonic waves
-
- H—ELECTRICITY
- H04—ELECTRIC COMMUNICATION TECHNIQUE
- H04B—TRANSMISSION
- H04B13/00—Transmission systems characterised by the medium used for transmission, not provided for in groups H04B3/00 - H04B11/00
- H04B13/02—Transmission systems in which the medium consists of the earth or a large mass of water thereon, e.g. earth telegraphy
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- Signal Processing (AREA)
- General Physics & Mathematics (AREA)
- Mathematical Physics (AREA)
- Pure & Applied Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Data Mining & Analysis (AREA)
- Electromagnetism (AREA)
- Algebra (AREA)
- Operations Research (AREA)
- Databases & Information Systems (AREA)
- Software Systems (AREA)
- General Engineering & Computer Science (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
【解決手段】物理的な海洋と水中の音場モデリングの融合の分野に属し、並列結合海洋音響予測システムは、海洋環境予測モジュールと、海洋-音響結合モジュールと、海洋音響伝播ソリューションモジュールと、を含む。海洋環境予測モジュールの前記システム全体における役割は、システムへ初期環境の入力を提供することであり、海洋音響結合モジュールの機能には、パラメータの初期化、地形の抽出、及び海洋環境のパラメータ断面の抽出が含まれ、温度、塩分及び海面のラフネスの出力を含む海洋モジュールのパラメータを音響モデルの入力に変換することにより、海洋モデルと音響モデルが組み合わされ、海洋音響伝播ソリューションモジュールの役割は、海洋音場の計算、海洋音場及びソナーの作用範囲を含む結果の出力である。【効果】水中音場のリアルタイム予測を実現し、三次元音場の正確な計算と世界の海域のソナー作用範囲の推定を迅速に実現できる。【選択図】図３A parallel coupling ocean acoustic prediction system, belonging to the field of fusion of physical ocean and underwater sound field modeling, comprises an ocean environment prediction module, an ocean-acoustic coupling module, and an ocean acoustic propagation solution module. include. The role of the Marine Environment Prediction Module in the overall system is to provide the initial environment input to the system, and the functions of the Ocean Acoustic Coupling Module include parameter initialization, terrain extraction, and parametric cross-section of the ocean environment. The ocean and acoustic models are combined by extracting and transforming the ocean module parameters, including temperature, salinity and sea surface roughness outputs, into acoustic model inputs, and the role of the Ocean Acoustic Propagation Solution Module is to Results output including sound field calculations, ocean sound field and sonar coverage. [Effect] Real-time prediction of the underwater sound field can be realized, and accurate calculation of the three-dimensional sound field and estimation of the sonar action range of the world's oceans can be realized quickly. [Selection drawing] Fig. 3
Description
本発明は、物理的な海洋と水中の音場モデリングの融合の分野に属し、並列結合海洋音響予測システム及び操作方法に関する。 The present invention is in the field of fusion of physical ocean and underwater sound field modeling, and relates to a parallel joint ocean acoustic prediction system and method of operation.
現在、国内外でさまざまな音の伝播モデルが確立されており、波動方程式の解の方法の違いにより、ノーマルモードモデル、放射線モデル、放物型方程式モデル、いくつかの結合モデル。 音場モデルは、独立して計算の面でより成熟してきた。 しかしながら、音場の既存の計算のほとんどはシリアルであり、並列アルゴリズムに関する研究は少ないである。小規模な二次元音の伝播問題を扱う場合、既存の音場計算モデルは実際のニーズを満たすことができるものの、近年、三次元音の伝播問題、大時間、大面積の音場計算の必要性があり、既存の音場計算効率は実際の科学研究のニーズを満たすことが困難になっている。 At present, various sound propagation models have been established in Japan and overseas, and due to differences in the method of solving the wave equation, there are normal mode models, radiation models, parabolic equation models, and several coupled models. Sound field models have become more mature in terms of independent computation. However, most of the existing computations of sound fields are serial, and there are few studies on parallel algorithms. When dealing with small-scale two-dimensional sound propagation problems, existing sound field calculation models can meet the actual needs, but in recent years, three-dimensional sound propagation problems, the need for long-time, large-area sound field calculations However, the existing sound field computational efficiency is difficult to meet the needs of practical scientific research.
海洋学では、天然資源省の第一海洋学研究所の研究者Qiao Fangli研究員が波動誘導混合理論を使用して、その研究チームは、世界的な高解像度の海波-潮流-還流の結合数値モデルFIO-COMを確立し、夏の海洋循環モデルによってシミュレートされたSSTが高くなる傾向、表面下層の温度が低くなる傾向、及び混合層が浅くなる傾向などの一般的な問題を克服し、表面下層の海洋の温度のシミュレーション偏差は約90％減少し、海洋数値モデルに通用する。 海洋モデルは、温度、塩分、海面の粗さなどのより正確な水文環境パラメータを提供できる。 In oceanography, Qiao Fangli, a researcher at the Department of Natural Resources' 1st Institute of Oceanography, used wave-induced mixing theory and his research team produced global high-resolution ocean wave-tidal current-backflow coupling figures. establish a model FIO-COM to overcome common problems such as higher SST, lower subsurface temperatures, and shallower mixed layers simulated by summer ocean circulation models; The simulated temperature anomalies in the subsurface ocean are reduced by about 90% and are compatible with numerical ocean models. Ocean models can provide more accurate hydrological environmental parameters such as temperature, salinity and surface roughness.
現在、海洋学と水中音響学の両方において、各研究分野は多くの理論と実践を蓄積してきた。 しかしながら、海洋学と音響学に焦点を当てたクロススタディはまだ少ない。既存の音場モデルはシミュレーション計算機能のみであり、水中音場予測機能はない。また、既存の研究の結果は計算規模が小さく、国内外の全球海洋音場に適した大規模なタイムスパン予測システムはない。 At present, both in oceanography and hydroacoustics, each research field has accumulated a great deal of theory and practice. However, there are still few cross-studies focusing on oceanography and acoustics. The existing sound field model has only a simulation calculation function, and does not have an underwater sound field prediction function. In addition, the results of existing research have a small computational scale, and there is no large-scale time-span prediction system suitable for the global ocean sound field in Japan and overseas.
本願発明で解決する技術的課題は、並列結合海洋音響予測システム及びその操作方法を提供することである。本発明のシステムは、動的海洋モデルと水中音響伝播モデルの協調の操作により、海洋モデルは、音響モデルの予測に必要な初期環境変数を生成し、さまざまな変化を伴うさまざまな海洋現象の周囲の音場、ソナー作用範囲などのタイムリーな予測を提供する。 The technical problem solved by the present invention is to provide a parallel joint ocean acoustic prediction system and its operation method. The system of the present invention is based on the coordinated manipulation of a dynamic ocean model and an underwater acoustic propagation model, whereby the ocean model generates the initial environmental variables necessary for prediction of the acoustic model, and the surroundings of various ocean phenomena with various changes. Provides timely predictions of sound field, sonar coverage, etc.
本願発明は、以下の技術案により、実現される。
海洋環境予測モジュールと、海洋-音響結合モジュールと、海洋音響伝播ソリューションモジュールと、を含む並列結合海洋音響予測システムであって、海洋環境予測モジュールの前記システム全体における役割は、システムへ初期環境の入力を提供することであり、海洋音響結合モジュールの機能には、パラメータの初期化、地形の抽出、及び海洋環境のパラメータ断面の抽出が含まれ、温度、塩分及び海面のラフネス（roughness）の出力を含む海洋モジュールのパラメータを音響モデルの入力に変換することにより、海洋モデルと音響モデルが組み合わされ、前記海洋音響伝播ソリューションモジュールの役割は、海洋音場の計算、海洋音場及びソナーの作用範囲を含む結果の出力であり、
The present invention is realized by the following technical proposals.
A parallel coupled ocean acoustic prediction system comprising an ocean environment prediction module, an ocean-acoustic coupling module, and an ocean acoustic propagation solution module, wherein the role of the ocean environment prediction module in the overall system is to input an initial environment to the system The functions of the Ocean Acoustic Coupling Module include parameter initialization, terrain extraction, and parametric cross-section extraction of the ocean environment, and outputs temperature, salinity and sea surface roughness. An ocean model and an acoustic model are combined by transforming the parameters of the ocean module, including the ocean module, into the input of the acoustic model, and the role of the ocean sound propagation solution module is to calculate the ocean sound field, to convert the ocean sound field and the sonar coverage into is the resulting output containing
さらに、前記海洋モデルは、海波-潮流-還流の結合数値モデルを用いる。
Furthermore, the ocean model uses a coupled numerical model of ocean waves-tidal currents-backflow.
また、本発明は、並列結合海洋音響予測システムの操作方法を提供する。前記操作方法は、具体的には、海洋環境予測モジュールは、初期環境の入力を提供するステップ１と、 海洋音響結合モジュールを通じて、初期化情報を読み取り、海洋環境予測モジュールにより提供された海洋水文環境予測結果を取得するステップ２であって、初期化情報により提供された計算タイプパラメータは、地形及び環境のパラメータを読み取る各プロセスの領域範囲を決定し、下記3つの計算タイプがあり、第１のタイプは、音源点と受信点の間の垂直断面における音場を計算し、このタイプは、音源点と受信点の間の垂直断面の地形と海洋環境パラメータを抽出することのみが必要となり、第２のタイプは、音源を円の中心として、 360°方位角内の水平断面における音場を計算し、この場合、円の中心として音源を抽出し、360°方位角の範囲内の地形および海洋環境のパラメータを抽出することが必要となり、第３のタイプは、音場計算に基づいて、指定されたソナー品質係数に従って、ソナーの最大作用範囲を計算し、この場合、必要な計算領域の地形および海洋環境のパラメータを抽出する必要となる、ステップ２と、海洋音響伝播ソリューションモジュールでは、システムの並列計算を行い、上記３つのタイプの音場計算には、マルチプロセスの並列計算が必要となり、第１のタイプ及び第2のタイプの計算では、異なる周波数ポイントと異なる方位角の音場計算は、周波数ポイントと方位角を各プロセスに均等に分布する必要であり、三次元音場波数積分の計算では、波数を各プロセスに均等に分布する必要があり、第３のタイプの計算では、規則的な長方形グリッドを使用して地理的空間を分割し、海域を経度と緯度の2次元で分割し、メインプロセスは各サブプロセスにタスクを割り当て、各サブプロセスは割り当てられたタスクに従って音場計算を実行し、各周波数ポイントでの各音源の音場計算結果を取得し、その後、メインプロセスは計算タイプの要求に従って、サブプロセスの音場計算結果を、合計及び平均化を含む処理をさらに行い、音場とソナー作用範囲を含む計算結果を出力するステップ３と、を含む。
The present invention also provides a method of operating a parallel joint ocean acoustic prediction system. Specifically, the operation method includes
従来技術と比較して、本願発明は以下の有利な効果を有する。
本願発明のシステムは、物理的な海洋モデルと音場モデルを組み合わせて、水中音場のリアルタイム予測を実現する。このシステムは、大量の水中音場計算の問題を解決するために、任意の高性能コンピュータシステムにも適用できる。大量の水中音場計算と計算速度の遅さの問題を解決し、三次元音場の正確な計算を迅速に実現し、地球の海でのソナーの距離を推定する。
Compared with the prior art, the present invention has the following advantageous effects.
The system of the present invention combines a physical ocean model and a sound field model to achieve real-time prediction of the underwater sound field. The system can be applied to any high performance computer system to solve a large number of underwater sound field calculation problems. Solve the problem of large amount of underwater sound field calculation and slow calculation speed, quickly realize accurate calculation of three-dimensional sound field, and estimate the sonar distance in the earth's ocean.
本願発明の音場モデルは、グローバル行列結合法線波モデルを使用し、グローバル行列（ＤＧＭ）法を使用して結合係数を解く。従来の方法と比較して、グローバルマトリックス結合のノーマルモードモデルは、行列変換や行列乗算などの多くの操作が回避され、モデルの安定性と計算効率が向上する。 The sound field model of the present invention uses a global matrix-coupled normal wave model and uses the global matrix (DGM) method to solve for the coupling coefficients. Compared with traditional methods, the normal mode model of global matrix combination avoids many operations such as matrix transformation and matrix multiplication, which improves the model stability and computational efficiency.
本願発明のシステムで使用される音場モデルは、波数積分技術を使用して二次元モデルを三次元の音の伝播の課題に拡張し、各波数の音場計算の並列化を実現するため、このモデルは、三次元の音の伝播の課題にも迅速に対処できる。音場計算に基づいて、本願発明のシステムはソナー品質係数を設定し、音場結果をさらに処理するため、本願発明のシステムはソナーの距離を予測する機能を備えている。 The sound field model used in the system of the present invention uses wavenumber integration technology to extend the two-dimensional model to the problem of three-dimensional sound propagation, and realizes parallelization of sound field calculations for each wavenumber. The model can also quickly address the problem of three-dimensional sound propagation. Based on the sound field calculations, the system of the present invention sets a sonar quality factor and further processes the sound field results, so that the system of the present invention has the ability to predict sonar distances.
本願発明の技術案について、図面を参照しながら、下記実施例によりさらに説明するものの、本願発明の保護範囲は、実施例のいかなる形態に限定されない。 The technical solution of the present invention will be further described by the following examples with reference to the drawings, but the protection scope of the present invention is not limited to any form of the examples.
実施例１
並列結合海洋音響予測システムは、海洋環境予測モジュールと、海洋-音響結合モジュールと、海洋音響伝播ソリューションモジュールと、を含む。
Example 1
A parallel coupled ocean acoustic prediction system includes an ocean environment prediction module, an ocean-acoustic coupling module, and an ocean acoustic propagation solution module.
海洋環境予測モジュールの前記システム全体における役割は、システムへ初期環境の入力を提供することである。 The role of the Marine Environmental Prediction Module in the overall system is to provide initial environmental input to the system.
海洋音響結合モジュールの機能には、パラメータの初期化、地形の抽出、及び海洋環境のパラメータ断面の抽出が含まれ、温度、塩分及び海面のラフネス（roughness）の出力を含む海洋モジュールのパラメータを音響モデルの入力に変換することにより、海洋モデルと音響モデルが組み合わされる。 The functionality of the Ocean Acoustic Coupling Module includes parameter initialization, terrain extraction, and parametric cross-section extraction of the ocean environment, and the ocean module parameters, including output of temperature, salinity and sea surface roughness. The ocean and acoustic models are combined by transforming them into model inputs.
前記海洋音響伝播ソリューションモジュールの役割は、海洋音場の計算、海洋音場及びソナーの作用範囲を含む結果の出力である。 The role of the Ocean Acoustic Propagation Solution Module is to calculate the ocean sound field and output results including the ocean sound field and sonar coverage.
本発明の結合海洋音響予測システムは、ＦＯＲＴＲＡＮ－９０プログラミング言語のモジュール化を採用しており、全体的な設計フレームワークが図１に示されている。 システム全体は並列化によって実現され、各音源、周波数点（ポイント）、地理空間の計算タスクは各プロセスに均等に分布される。地理空間は、図2に示すように、規則的な長方形グリッドの並列分割方式を採用している。 本願発明のシステムは最終的に3つのタイプの計算結果を提供できる。1つのタイプは距離-深さ垂直断面の音場計算結果、もう1つのタイプは距離-方位角水平断面の音場計算結果、最後の1つのタイプは、世界または地域の任意の深度でのソナー作用範囲の推定結果である。 The joint ocean acoustic prediction system of the present invention employs the modularization of the FORTRAN-90 programming language, and the overall design framework is shown in FIG. The whole system is realized by parallelization, and each sound source, frequency point (point) and geospatial computation task are evenly distributed to each process. Geospatial employs a regular rectangular grid parallel division scheme, as shown in Figure 2. The system of the present invention can finally provide three types of calculation results. One type is range-depth vertical section sound field calculation results, another type is range-azimuth horizontal section sound field calculation results, and the last type is sonar at any depth in the world or region. This is the result of estimation of the range of action.
システム実行中具体的な計算プロセスを図3に示す。
まず、 海洋音響結合モジュールを通じて、音源の位置、周波数、計算タイプなどを含む初期化情報を読み取り、海洋環境予測モジュールにより提供された、温度、塩分及び海面のラフネスなどの海洋水文環境予測結果を取得した。
Figure 3 shows the specific calculation process during system execution.
First, through the ocean acoustic coupling module, read the initialization information including the position of the sound source, frequency, calculation type, etc., and obtain the marine hydrological environment prediction results such as temperature, salinity and sea surface roughness provided by the marine environment prediction module. did.
次に、海洋音響結合モジュールは、初期化情報により提供された計算タイプパラメータは、地形及び環境のパラメータを読み取る各プロセスの領域範囲を決定し、下記3つの計算タイプがあり、第１のタイプは、音源点と受信点の間の垂直断面の地形と海洋環境パラメータを抽出する必要となり、第２のタイプは、音源を円の中心として、 360°方位角の範囲内の地形および海洋環境のパラメータを抽出する必要となり、第３のタイプは、必要な計算領域の地形および海洋環境のパラメータを抽出する必要となる。 Next, the ocean acoustic coupling module uses the calculation type parameter provided by the initialization information to determine the area range of each process that reads the terrain and environment parameters. , it is necessary to extract topographical and oceanic environment parameters in a vertical section between the source point and the receiving point, and the second type is to extract the topographical and oceanic environment parameters within a 360° azimuth angle with the source as the center of a circle. and the third type requires extracting the topographical and marine environmental parameters of the required computational domain.
次に、海洋音響伝播ソリューションモジュールでは、 メインプロセスは各サブプロセスにタスクを割り当て、各サブプロセスは割り当てられたタスクに従って音場計算を実行し、各サブプロセスの計算プロセスは次のとおりである。１）音波伝播経路の周りの経度と緯度のグリッドポイントの水深、温度、塩分、およびその他のデータを読み取った。２）グリッドポイントのデータを補間して、水深、温度、および塩分を取得し、音波の伝播経路経験式を使用して音速断面を計算した。３）音源の位置、深さ、周波数、音速、海面のラフネス、海底の地形およびパラメータに応じて、単一の周波数ポイントで単一の音源の音場を計算した。 Next, in the ocean sound propagation solution module, the main process assigns a task to each sub-process, each sub-process performs sound field calculation according to the assigned task, and the calculation process of each sub-process is as follows. 1) Read depth, temperature, salinity, and other data at longitude and latitude grid points around the sound propagation path. 2) The grid point data was interpolated to obtain water depth, temperature and salinity, and the sound velocity cross-section was calculated using the empirical formula for sound wave propagation paths. 3) Calculated the sound field of a single source at a single frequency point according to the source position, depth, frequency, sound velocity, sea surface roughness, seafloor topography and parameters.
最後に、計算タイプのメインプロセスに従って、音場計算結果をさらに処理した。 Finally, the sound field calculation results were further processed according to the calculation type main process.
単一グリッド点のソナー作用距離の計算フローチャートは、図4に示されるように、音源の進行距離での単一点音伝播損失を順次計算し、音場伝播損失をソナーの品質係数と比較し、ソナーの最大作用距離を得た。 The calculation flowchart of the sonar working distance of a single grid point is shown in Figure 4, sequentially calculating the single-point sound propagation loss in the traveling distance of the sound source, comparing the sound field propagation loss with the sonar quality factor, Obtained the maximum working distance of sonar.
本願発明について、下記具体的の実施形態と併せて詳細に説明する。
実施例2
2つのポイントA（124.1°E、26.38°N）とB（124.8°E、24.63°N）を選択し、地形図が図5に示されるように、AからBへの垂直断面音場伝播を計算する。 図6に沖縄海溝域の海面の高さ異常を示す。ポイントAとポイントBの音の伝播も黒潮前線の渦の影響を受けていることがわかる。ポイントAは海面の低温域にあり、ポイントBは海の高温域にあり、音の速度断面は、水平方向に変化する。ポイントAとポイントBの垂直断面の音場予測を実行する。最初に、1km間隔でポイントAとポイントBの間の補間（Interpolation）を実行し、各補間ポイントの緯度と経度を取得した。次に、海洋モデルの水文予測結果から、AからBまでの線の近くの長方形の領域の温度と塩分値を読み取る。Global Land and Sea Database（GEBCO）から地形データを読み取り、各補間ポイント周辺の4つの緯度と経度のグリッドの水深、温度、塩分データに基づいて2次元線形補間を実行し、A、Bと補間ポイントの水深、温度、塩分などを取得した。実験式に従って音速を求め、1km間隔の各ポイントでのレベル変化の音速を求めた。最後に、音場モデルを音場伝播経路の水深、音速、音源位置、深さ、周波数、海面粗さ、海底パラメータなどの音場モデルに入力し、音場計算を実行した。単一の垂直断面の音場計算結果については、複数の周波数点の音場に対して並列計算を実行し、最終的にすべてのプロセスをまとめて、最終的に周波数ポイントの平均距離-深さ垂直断面の音場予測結果を取得した。たとえば、選択された2つのポイントAとBについて、この経路に沿った距離-深さ垂直断面の音場予測を提供した。音源はポイントAの水深50mのところにあった。音源の中心周波数は200Hzで、帯域幅は1/3オクターブである。図7は、結合海洋音響予測システムによる沖縄の溝を通る垂直断面音場予測の結果を示す。沖縄の溝の側面にある音源から放出される音波が伝播していることがわかる。溝を通って、音波は沈み、海底に沿って移動し、溝の反対側を通過する音響エネルギーはわずかであった。
The present invention will be described in detail in conjunction with the following specific embodiments.
Example 2
We choose two points A (124.1°E, 26.38°N) and B (124.8°E, 24.63°N) and set the vertical cross-section sound field propagation from A to B as the topographic map is shown in Fig. 5. calculate. Figure 6 shows the sea level anomaly in the Okinawa trench. It can be seen that the sound propagation at points A and B is also affected by the vortex of the Kuroshio front. Point A is in the low-temperature region of the sea surface, point B is in the high-temperature region of the sea, and the cross section of sound velocity varies horizontally. Perform a sound field prediction of the vertical section of points A and B. First, we performed interpolation between points A and B at 1km intervals and obtained the latitude and longitude of each interpolation point. Next, read the temperature and salinity values of the rectangular area near the line from A to B from the hydrological prediction results of the ocean model. It reads terrain data from the Global Land and Sea Database (GEBCO) and performs two-dimensional linear interpolation based on the depth, temperature, and salinity data in four latitude and longitude grids around each interpolation point, and calculates A, B and the interpolation point. We acquired the water depth, temperature, salinity, etc. The speed of sound was obtained according to the empirical formula, and the sound speed of the level change at each point with an interval of 1 km was obtained. Finally, the sound field model was input into the sound field model such as water depth, sound velocity, sound source position, depth, frequency, sea surface roughness, and seabed parameters of the sound field propagation path, and the sound field calculation was executed. For the sound field calculation result of a single vertical section, parallel calculation is performed for the sound field of multiple frequency points, and finally all the processes are put together to finally obtain the average distance-depth of the frequency points The sound field prediction results for the vertical cross section were obtained. For example, for two selected points A and B, we provide sound field predictions for the distance-depth perpendicular cross-section along this path. The sound source was at point A at a depth of 50m. The center frequency of the sound source is 200Hz and the bandwidth is 1/3 octave. Figure 7 shows the results of vertical section sound field prediction through the Okinawa trench by the coupled ocean acoustic prediction system. It can be seen that the sound waves emitted from the sound source on the side of the trench in Okinawa are propagating. Through the trench, the sound waves subsided and traveled along the seafloor, with little acoustic energy passing through the other side of the trench.
実施例3
音源ポイントA（124.1°E、26.38°N）について、Aを中心とした1回転の距離-方位角の水平断面の音場を計算した。ポイントA周辺の水平断面の音場予測については、まず、海洋モデルの予測結果から、A周辺の円形領域の水深、温度、塩分、海面粗さなどの予測結果を読み取った。平均値は方位角360°間隔で取得され、各方位角伝播経路の水深、温度、塩分などは1km間隔で取得された。各補間点のデータは周囲の4つの緯度と経度のポイントの2次線形補間データによって取得された。
Example 3
For the sound source point A (124.1°E, 26.38°N), the sound field of the distance-azimuth horizontal cross-section of one rotation about A was calculated. Regarding the sound field prediction of the horizontal cross section around Point A, first, the prediction results of the water depth, temperature, salinity, sea surface roughness, etc. of the circular area around A were read from the prediction results of the ocean model. Average values were obtained at azimuth angle intervals of 360°, and water depth, temperature, salinity, etc. along each azimuth propagation path were obtained at 1 km intervals. The data of each interpolated point were obtained by quadratic linearly interpolated data of the surrounding four latitude and longitude points.
実験式によれば、音速は水深、温度、塩分などから計算された。水深や音速、音源情報、海面粗さ、海底パラメータなどのデータを音場モデルに入力し、音場計算を行った。各方位角と各周波数点の音場は、並列かつ同期して計算された。沖縄の溝の端にあるA点（東経124.1度、北緯26.38度）を選択し、地形図5を参照し、この点から200km以内の深さ50mの水平断面の音場予測を行った。音源はA点（東経124.1度、北緯26.38度）にあり、水深は50mである。音源の中心周波数は200Hzで、帯域幅は1/3オクターブである。図8は、結合海洋音響予測システムによる沖縄溝海域のA点を中心とした受信深度50mの水平断面の音場予測結果である。水深50mでは、溝の上の音響エネルギーが非常に弱く、溝が音響エネルギーを遮蔽する効果があり、溝を通過する音響エネルギーが非常に小さいことがわかる。
According to the empirical formula, the speed of sound was calculated from water depth, temperature, salinity, etc. Data such as water depth, sound velocity, sound source information, sea surface roughness, seafloor parameters, etc. were input into the sound field model, and the sound field was calculated. The sound field for each azimuth angle and each frequency point was calculated in parallel and synchronously. Selecting point A (124.1 degrees east longitude, 26.38 degrees north latitude) at the edge of the trench in Okinawa, referring to
実施例4
世界の海域または特定の地域でのソナー作用範囲の予測は、垂直および水平断面での音場予測よりも計算量が多くなる。 計算プロセスでは、図2に示すように、計算タスクを均等に分割する必要がある。 本実施例では、世界の海域でのソナー作用範囲予測が実行された。 全地球の海域で音源点が1°間隔で配置されていると仮定して、周辺域での音源の伝播損失を計算し、相反性の原理と所定のソナーに基づいてソナー作用範囲を求めることができ。海域の計算タスクは、プロセスの数に応じて規則的な長方形のグリッドに分割され、各プロセスは、地形、温度、塩分、海面の粗さ、およびその他のパラメータの独自の計算範囲を読み取った。 計算では、まず、音源が海域にあるかどうかを判断し、音源が陸地にある場合、直接スキップする。 単一のターゲット音源に対応するソナー作用範囲の計算プロセスは次のとおりである。１）360°の範囲でターゲット音源の周りに一定の角度間隔でn個の方位角を決定した。２）単一の方位角の水平距離伝播損失を計算した。３）計算結果をソナー品質係数と比較して、単一の方位角のソナー範囲を決定した。４）同じ方法を使用して残りの方位角の動作距離を計算し、すべての方位角の動作距離を合計して平均し、システムに対応するターゲット音源の動作距離を取得した。図9は、全球海洋におけるソナー作用範囲の予測図である。（a）ソナーは水深50mに設置され、（b）ソナーは海底に設置され、計算されたメッシュ間隔は1°×1°であり、各グリッドポイントには、周波数ポイントの数が7であり、方位角の数が16であり、ソナー品質係数が100dBであると想定されている。単一のネットワークポイント、単一の周波数ポイント、および単一の方位角により音の伝播計算を行い、必要なコアの最大数は180（緯度グリッドポイントの数）×360（経度グリッドポイントの数）× 7（周波数ポイントの数）×16（方位角の数）×70％（地球の海の比率）= 5080320（個）。 計算タスクは、プロセスの数に応じて均等に分割されます。 ソナーが水深50mにある場合、北極と南極はちょうど地表音響チャネルの深さにあることがわかる。この場合、ソナーの作用距離は長くなる。ソナーは海底に位置する場合、深海ソナーの作用距離は浅瀬海域のソナーの作用範囲よりも優れている。
Example 4
Predicting the sonar operating range over the world's oceans or specific regions is more computationally intensive than sound field prediction in vertical and horizontal planes. In the computational process, the computational tasks should be divided evenly, as shown in Figure 2. In this example, a sonar operating range prediction was performed for the world's oceans. Assuming that the sound source points are spaced at 1° intervals in the global sea area, calculate the propagation loss of the sound source in the surrounding area, and obtain the sonar coverage based on the principle of reciprocity and a given sonar. can. The sea area computation task was divided into a regular rectangular grid according to the number of processes, with each process reading its own computational range of topography, temperature, salinity, sea surface roughness, and other parameters. The calculation first determines whether the sound source is in the sea area, and skips directly if the sound source is in land. The sonar coverage calculation process for a single target source is as follows. 1) We determined n azimuth angles at regular angular intervals around the target sound source over a range of 360°. 2) Calculated the horizontal distance propagation loss for a single azimuth. 3) The calculated results were compared to the sonar quality factor to determine the sonar range for a single azimuth. 4) Using the same method to calculate the working distances of the remaining azimuth angles, summing and averaging the working distances of all the azimuth angles to obtain the working distance of the target sound source corresponding to the system. FIG. 9 is a projection diagram of the sonar operating range in the global ocean. (a) the sonar was installed at a depth of 50 m, (b) the sonar was installed at the seafloor, the calculated mesh spacing was 1° × 1°, each grid point had 7 frequency points, It is assumed that the number of azimuth angles is 16 and the sonar quality factor is 100 dB. A single network point, a single frequency point, and a single azimuth for sound propagation calculations, the maximum number of cores required is 180 (number of latitude grid points) x 360 (number of longitude grid points) × 7 (number of frequency points) × 16 (number of azimuth angles) × 70% (percentage of the Earth's oceans) = 5080320 (pieces). Computational tasks are divided evenly according to the number of processes. If the sonar is at a depth of 50m, it can be seen that the north and south poles are just at the depth of the surface acoustic channel. In this case, the working distance of the sonar is increased. When the sonar is located on the ocean floor, the working range of deep water sonar is superior to that of shallow water sonar.
Claims (3)
海洋環境予測モジュールの前記システム全体における役割は、海洋音響結合モジュールへ初期環境の入力を提供することであり、
海洋音響結合モジュールの機能には、パラメータの初期化、地形の抽出、及び海洋環境のパラメータ断面の抽出が含まれ、温度、塩分及び海面のラフネス（roughness）の出力を含む海洋モジュールのパラメータを音響モデルの入力に変換することにより、海洋モデルと音響モデルが組み合わされ、
前記海洋音響伝播ソリューションモジュールの役割は、海洋モデルと音響モデルの組み合わせに基づいて、海洋音場の計算、海洋音場及びソナーの作用範囲を含む結果の出力であり、
the role of the ocean environment prediction module in the overall system is to provide initial environmental input to the ocean acoustic coupling module ;
The functionality of the Ocean Acoustic Coupling Module includes parameter initialization, terrain extraction, and parametric cross-section extraction of the ocean environment, and the ocean module parameters, including output of temperature, salinity and sea surface roughness. The ocean and acoustic models are combined by transforming them into model inputs,
The role of the Ocean Acoustic Propagation Solution Module is to calculate the ocean sound field, output the results including the ocean sound field and sonar coverage , based on the combination of the ocean model and the acoustic model ;
前記操作方法は、具体的には、
海洋環境予測モジュールは、海洋音響結合モジュールへ初期環境の入力を提供するステップ１と、
海洋音響結合モジュールを通じて、初期化情報を読み取り、海洋環境予測モジュールにより提供された海洋水文環境予測結果を取得するステップ２であって、
初期化情報により提供された計算タイプパラメータは、地形及び環境のパラメータを読み取る各プロセスの領域範囲を決定し、下記３つの計算タイプがあり、第１のタイプは、音源点と受信点の間の垂直断面における音場を計算し、このタイプは、音源点と受信点の間の垂直断面の地形と海洋環境パラメータを抽出し、第２のタイプは、音源を円の中心として、 ３６０°方位角内の水平断面における音場を計算し、この場合、円の中心として音源を抽出し、３６０°方位角の範囲内の地形および海洋環境のパラメータを抽出し、第３のタイプは、音場計算に基づいて、指定されたソナー品質係数に従って、ソナーの最大作用範囲を計算し、この場合、必要な計算領域の地形および海洋環境のパラメータを抽出する、ステップ２と、
海洋音響伝播ソリューションモジュールでは、上記３つのタイプの音場計算でマルチプロセスの並列計算を行い、第１のタイプ及び第２のタイプの計算では、異なる周波数ポイントと異なる方位角の音場計算は、周波数ポイントと方位角を各プロセスに均等に分布させ、三次元音場波数積分の計算では、波数を各プロセスに均等に分布させ、第３のタイプの計算では、規則的な長方形グリッドを使用して地理的空間を分割し、海域を経度と緯度の２次元で分割し、
メインプロセスは各サブプロセスにタスクを割り当て、各サブプロセスは割り当てられたタスクに従って音場計算を実行し、各周波数ポイントでの各音源の音場計算結果を取得し、その後、メインプロセスは計算タイプの要求に従って、サブプロセスの音場計算結果を、合計及び平均化を含む処理をさらに行い、音場とソナー作用範囲を含む計算結果を出力するステップ３と、
を含む、請求項１または２に記載の並列結合海洋音響予測システムの操作方法。 A method of operating the parallel coupled ocean acoustic prediction system, comprising:
Specifically, the operation method is
step 1, wherein the ocean environment prediction module provides input of the initial environment to the ocean acoustic coupling module ;
Step 2 of reading the initialization information through the ocean acoustic coupling module and obtaining the marine hydrological environment prediction result provided by the marine environment prediction module,
The computation type parameter provided by the initialization information determines the area extent of each process that reads the terrain and environment parameters.There are three computation types: Calculate the sound field in the vertical section, this type extracts the terrain and ocean environment parameters in the vertical section between the source point and the receiving point, the second type uses the sound source as the center of the circle, and the 360° azimuth angle Calculate the sound field in the horizontal section within, in this case extracting the sound source as the center of the circle, and extracting the parameters of the terrain and marine environment within a 360° azimuth angle, the third type is the sound field calculation calculating the maximum working range of the sonar according to a specified sonar quality factor based on, where extracting the topographical and marine environment parameters of the required computational domain, step 2;
In the ocean sound propagation solution module, multi-process parallel calculation is performed for the above three types of sound field calculations . , the frequency points and azimuth angles are distributed evenly over each process, the calculation of the three-dimensional sound field wavenumber integral distributes the wavenumbers evenly over each process, and the third type of calculation uses a regular rectangular grid to divide the geographic space, divide the sea area into two dimensions of longitude and latitude,
The main process assigns a task to each sub-process, each sub-process performs sound field calculation according to the assigned task, obtains the sound field calculation result of each sound source at each frequency point, and then the main process enters the calculation type step 3 of further processing the sound field calculation results of the sub-process including summing and averaging according to the request of , and outputting the calculation results including the sound field and sonar action range;
3. A method of operating a parallel joint ocean acoustic prediction system according to claim 1 or 2, comprising:
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202110769864.X | 2021-07-08 | ||
| CN202110769864.XA CN113259034B (en) | 2021-07-08 | 2021-07-08 | Parallel coupled marine acoustic forecasting system and operation method |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP7116852B1 true JP7116852B1 (en) | 2022-08-10 |
| JP2023010540A JP2023010540A (en) | 2023-01-20 |
Family
ID=77190939
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2022013517A Active JP7116852B1 (en) | 2021-07-08 | 2022-01-31 | Parallel coupled ocean acoustic prediction system and method of operation |
Country Status (2)
| Country | Link |
|---|---|
| JP (1) | JP7116852B1 (en) |
| CN (1) | CN113259034B (en) |
Cited By (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116973899A (en) * | 2023-07-17 | 2023-10-31 | 哈尔滨工程大学 | Under-ice multi-target positioning method and device based on convex clustering |
| CN117574272A (en) * | 2023-10-20 | 2024-02-20 | 青岛海洋地质研究所 | Ocean data processing and classifying method |
| CN117910279A (en) * | 2024-03-18 | 2024-04-19 | 西北工业大学 | Three-dimensional sound field rapid processing method based on space coordinate rotation transformation |
Families Citing this family (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN114707327B (en) * | 2022-04-01 | 2022-11-18 | 自然资源部第一海洋研究所 | Parallel marine environment underwater sound characteristic diagnosis method, storage medium and device |
| CN114925496B (en) * | 2022-04-19 | 2023-05-12 | 中国科学院声学研究所 | Ocean environment noise forecasting method |
| CN114510848B (en) * | 2022-04-20 | 2022-08-02 | 自然资源部第一海洋研究所 | Offshore wind farm underwater noise calculation method, software and measurement device |
| CN115391069B (en) * | 2022-10-27 | 2023-02-03 | 山东省计算中心（国家超级计算济南中心） | Parallel communication method and system based on ocean mode ROMS |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003004848A (en) | 2001-06-25 | 2003-01-08 | Nec Corp | System for conjecturing and visualizing sonar searching range, method using the same and program for the system |
| JP2007040734A (en) | 2005-08-01 | 2007-02-15 | Nec Corp | Transmission mode recommend system, sonar system, transmission mode recommend method, and its program and storage medium |
| WO2014083790A1 (en) | 2012-11-27 | 2014-06-05 | 日本電気株式会社 | Environment measuring system and environment measuring method |
| CN107576388A (en) | 2017-08-22 | 2018-01-12 | 哈尔滨工程大学 | Three-dimensional structure sound source radiation sound field forecasting procedure under a kind of shallow sea channel |
| CN113010836A (en) | 2020-12-04 | 2021-06-22 | 中国船舶重工集团公司第七一五研究所 | Parabolic equation model method for forecasting underwater sound field generated by sound source in air |
Family Cites Families (7)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US6876598B1 (en) * | 2003-06-12 | 2005-04-05 | Ronald Francis Pannatoni | Method for predicting propagation of sound in an ocean |
| CN102156817B (en) * | 2011-04-11 | 2013-04-03 | 国家海洋信息中心 | Rapid prediction method for drifting path and diffusion of marine spilled oil |
| US10168253B2 (en) * | 2014-05-30 | 2019-01-01 | General Electric Company | Marine riser management system including subsea acoustic monitoring platform and an associated method |
| CN109657257A (en) * | 2017-10-12 | 2019-04-19 | 中国船舶重工集团公司第七六研究所 | Labyrinth radiated sound field calculation method in a kind of ocean acoustic channel |
| CN107885934B (en) * | 2017-11-07 | 2021-04-02 | 哈尔滨工程大学 | Elastic structure acoustic radiation forecasting method under ocean channel based on coupled FEM-PE |
| CN108216492B (en) * | 2017-12-22 | 2022-09-20 | 自然资源部第一海洋研究所 | High-precision submerged buoy array system for realizing marine data area monitoring |
| CN109000779B (en) * | 2018-07-12 | 2019-12-20 | 青岛中海潮科技有限公司 | Rapid matching system for ocean sound propagation reliability model |
-
2021
- 2021-07-08 CN CN202110769864.XA patent/CN113259034B/en active Active
-
2022
- 2022-01-31 JP JP2022013517A patent/JP7116852B1/en active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2003004848A (en) | 2001-06-25 | 2003-01-08 | Nec Corp | System for conjecturing and visualizing sonar searching range, method using the same and program for the system |
| JP2007040734A (en) | 2005-08-01 | 2007-02-15 | Nec Corp | Transmission mode recommend system, sonar system, transmission mode recommend method, and its program and storage medium |
| WO2014083790A1 (en) | 2012-11-27 | 2014-06-05 | 日本電気株式会社 | Environment measuring system and environment measuring method |
| CN107576388A (en) | 2017-08-22 | 2018-01-12 | 哈尔滨工程大学 | Three-dimensional structure sound source radiation sound field forecasting procedure under a kind of shallow sea channel |
| CN113010836A (en) | 2020-12-04 | 2021-06-22 | 中国船舶重工集团公司第七一五研究所 | Parabolic equation model method for forecasting underwater sound field generated by sound source in air |
Cited By (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116973899A (en) * | 2023-07-17 | 2023-10-31 | 哈尔滨工程大学 | Under-ice multi-target positioning method and device based on convex clustering |
| CN116973899B (en) * | 2023-07-17 | 2024-04-02 | 哈尔滨工程大学 | Under-ice multi-target positioning method and device based on convex clustering |
| CN117574272A (en) * | 2023-10-20 | 2024-02-20 | 青岛海洋地质研究所 | Ocean data processing and classifying method |
| CN117910279A (en) * | 2024-03-18 | 2024-04-19 | 西北工业大学 | Three-dimensional sound field rapid processing method based on space coordinate rotation transformation |
| CN117910279B (en) * | 2024-03-18 | 2024-05-28 | 西北工业大学 | Three-dimensional sound field rapid processing method based on space coordinate rotation transformation |
Also Published As
| Publication number | Publication date |
|---|---|
| JP2023010540A (en) | 2023-01-20 |
| CN113259034B (en) | 2021-09-28 |
| CN113259034A (en) | 2021-08-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP7116852B1 (en) | Parallel coupled ocean acoustic prediction system and method of operation | |
| Bilskie et al. | Data and numerical analysis of astronomic tides, wind‐waves, and hurricane storm surge along the northern Gulf of Mexico | |
| Simmons et al. | Modeling and prediction of internal waves in the South China Sea | |
| Baptista et al. | A cross-scale model for 3D baroclinic circulation in estuary–plume–shelf systems: II. Application to the Columbia River | |
| KR101613186B1 (en) | Method and system for predicting ocean circulation and wave in coastal areas | |
| Lermusiaux et al. | Coupled ocean–acoustic prediction of transmission loss in a continental shelfbreak region: Predictive skill, uncertainty quantification, and dynamical sensitivities | |
| Tsujino et al. | Simulating present climate of the global ocean–ice system using the Meteorological Research Institute Community Ocean Model (MRI. COM): Simulation characteristics and variability in the Pacific sector | |
| KR101563488B1 (en) | Method and system for predicting wave in coastal areas | |
| Zhang et al. | Method of establishing an underwater digital elevation terrain based on kriging interpolation | |
| Bohluly et al. | Evaluation of wind induced currents modeling along the Southern Caspian Sea | |
| Zhao et al. | A modification of HASM for interpolating precipitation in China | |
| CN115221721A (en) | ROMS-based regional sea area tidal current numerical simulation method | |
| Zieger et al. | Hindcasting of tropical cyclone winds and waves | |
| KR102492075B1 (en) | Apparatus for predicting of ocean and method thereof | |
| Wornom et al. | On coupling the SWAN and WAM wave models for accurate nearshore wave predictions | |
| Deng et al. | A method for assessing the coastline recession due to the sea level rise by assuming stationary wind-wave climate | |
| CN110287620B (en) | Spherical coordinate system density interface forward modeling method and system suitable for earth surface observation surface | |
| Liu et al. | Response of sound propagation characteristics to Luzon cold eddy coupled with tide in the Northern South China Sea | |
| Lin et al. | Characteristics of gravity and magnetic fields and deep structural responses in the southern part of the Kyushu-Palau Ridge | |
| Alcaras et al. | Kriging interpolation of bathymetric data for 3D model of the Bay of Pozzuoli (Italy) | |
| Jakacki et al. | eBalticGrid–an interactive platform for the visualisation of results from a high-resolution operational Baltic Sea model | |
| Renkl et al. | The Alongshore Tilt of Mean Dynamic Topography and its Implications for Model Validation and Ocean Monitoring | |
| Li et al. | Application of a 3D parabolic equation model with higher-order operator approximation | |
| Huang et al. | Two-Dimensional Tidal Modeling Based on Synchronously Optimized Boundary Conditions: Case Study of Haizhou Bay of the Yellow Sea in China | |
| Kelly et al. | Validation of the FAST forecast model for the storm surges due to hurricanes Wilma and Ike |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20220131 |
|
| A871 | Explanation of circumstances concerning accelerated examination |
Free format text: JAPANESE INTERMEDIATE CODE: A871 Effective date: 20220131 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220217 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20220419 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20220615 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20220705 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20220729 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 7116852 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |