(57)【特許請求の範囲】
【請求項１】
海面に平行で且つ直交する二方向における水平波数の各成分と、浮力周波数およびコリオリパラメータとが外部から入力され、前記水平波数の各成分を用いて複数の合成水平波数を算出し、算出した前記各合成水平波数と、前記浮力周波数および前記コリオリパラメータとを用いて、モード番号に対する複数の固有深度関数および複数の固有周波数を算出する内部波モード演算部と、
前記内部波モード演算部において算出された前記各固有周波数を用いて、前記各固有周波数および前記モード番号に対するスペクトル特性を前記各合成水平波数および前記モード番号に対する内部波スペクトル特性へ変換する際に使用される微分値を前記各固有周波数の差分と前記各合成水平波数の差分を用いて算出する分散特性演算部と、
前記内部波モード演算部において算出された前記各固有周波数と、前記浮力周波数および前記コリオリパラメータとを用いて前記スペクトル特性を算出し、算出した前記スペクトル特性を、前記微分値を用いて前記内部波スペクトル特性に変換するスペクトル特性演算部と、
前記スペクトル特性演算部による変換後の前記内部波スペクトル特性と、前記内部波モード演算部において算出された前記各固有深度関数および前記各固有周波数と、前記水平波数の各成分とを用いて、海中の等密度面の変位である内部波変位を算出する内部波変位算出部と、を有する内部波変位シミュレーション装置。
【請求項２】
前記分散特性演算部は、前記各固有周波数の前進差分に対する前記各合成水平波数の前進差分の比率と、前記各固有周波数の後退差分に対する前記各合成水平波数の後退差分の比率との加重平均値を用いて前記微分値を算出する請求項１に記載の内部波変位シミュレーション装置。
【請求項３】
前記分散特性演算部は、前記各合成水平波数の前進差分に対する前記各固有周波数の前進差分の比率と、前記各合成水平波数の後退差分に対する前記各固有周波数の後退差分の比率との加重平均値を用いて前記微分値を算出する請求項１に記載の内部波変位シミュレーション装置。
【請求項４】
海面に平行で且つ直交する二方向における水平波数の各成分を外部から入力し、入力した前記水平波数の各成分から複数の合成水平波数を算出する水平波数算出工程と、
前記水平波数算出工程において算出した前記各合成水平波数と、外部から入力された浮力周波数およびコリオリパラメータとを用いて、モード番号に対する複数の固有深度関数および複数の固有周波数を求める内部波モード演算工程と、
前記内部波モード演算工程において算出された前記各固有周波数を用いて、前記各固有周波数および前記モード番号に対するスペクトル特性を前記各合成水平波数および前記モード番号に対する内部波スペクトル特性へ変換する際に使用する微分値を前記各固有周波数の差分と前記各合成水平波数の差分を用いて算出する微分値算出工程と、
前記内部波モード演算工程において算出された前記各固有周波数と、前記浮力周波数および前記コリオリパラメータとを用いて前記スペクトル特性を算出するスペクトル特性算出工程と、
前記スペクトル特性算出工程において算出された前記スペクトル特性を、前記微分値を用いて前記内部波スペクトル特性に変換するスペクトル特性変換工程と、
前記スペクトル特性変換工程における変換後の前記内部波スペクトル特性と、前記内部波モード演算工程において算出した前記各固有深度関数および前記各固有周波数と、前記水平波数の各成分とを用いて、海中の等密度面の変位である内部波変位を算出する内部波変位算出工程と、を有する内部波変位シミュレーション方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、海洋内部波による海中の等密度面の変位を予測する内部波変位シミュレーション装置および内部波変位シミュレーション方法に関する。
【背景技術】
【０００２】
海洋における海水は、海面から海底に向かって密度が徐々に大きくなるように成層している。海面上の風および潮汐などの外力が働くと、海面が波立つのと同じように、海中の等密度面には海洋内部波が生じる。海洋内部波は、海中の水温および音速の分布に時間的および空間的な変動をもたらす。そして、海中の水温および音速の分布の変動は、例えば、海中における音波伝搬のシミュレーション結果に影響を与えることが知られており、その影響の大きさを見積もるためには、海洋内部波による等密度面の変位を予測する必要がある。
【０００３】
海洋内部波による等密度面のランダムな変位は、海洋内部波の水平波数およびモード番号に対するスペクトル特性（内部波スペクトル特性）を分散とするランダム変数を要素に含む数式を用いてシミュレーションされる（例えば非特許文献１参照）。そして、内部波スペクトル特性は、海洋内部波の固有周波数およびモード番号に対するスペクトル特性を算出し、算出した結果を所定の数式により変換することで求められる（例えば非特許文献２参照）。
【先行技術文献】
【非特許文献】
【０００４】
【非特許文献１】Colosi and Brown, "Efficient numerical simulation of stochastic internal-wave induced sound-speed perturbation fields," J. Acoust. Soc. Am. 103(4), pp.2232-2235, 1998
【非特許文献２】Flatte, "3.5 Equivalent spectra," in Sound transmission through a fluctuating ocean (Cambridge University Press, 1979)
【発明の概要】
【発明が解決しようとする課題】
【０００５】
しかしながら、非特許文献１および２に記載の技術では、海洋内部波の水平波数と固有周波数との関係を近似的に処理して上記変換を行っているため、近似的な処理の影響により、海洋内部波による等密度面の変位のシミュレーションの最終的な結果に誤差が生じるという課題がある。この課題に鑑みて、従来から、海洋内部波による等密度面の変位のシミュレーション結果に生じる誤差を低減する手法が望まれている。
【課題を解決するための手段】
【０００６】
本発明に係る内部波変位シミュレーション装置は、海面に平行で且つ直交する二方向における水平波数の各成分と、浮力周波数およびコリオリパラメータとが外部から入力され、水平波数の各成分を用いて複数の合成水平波数を算出し、算出した各合成水平波数と、浮力周波数およびコリオリパラメータとを用いて、モード番号に対する複数の固有深度関数および複数の固有周波数を算出する内部波モード演算部と、内部波モード演算部において算出された各固有周波数を用いて、各固有周波数およびモード番号に対するスペクトル特性を各合成水平波数およびモード番号に対する内部波スペクトル特性へ変換する際に使用される微分値を各固有周波数の差分と各合成水平波数の差分を用いて算出する分散特性演算部と、内部波モード演算部において算出された各固有周波数と、浮力周波数およびコリオリパラメータとを用いてスペクトル特性を算出し、算出したスペクトル特性を、微分値を用いて内部波スペクトル特性に変換するスペクトル特性演算部と、スペクトル特性演算部による変換後の内部波スペクトル特性と、内部波モード演算部において算出された各固有深度関数および固有周波数と、水平波数の各成分とを用いて、海中の等密度面の変位である内部波変位を算出する内部波変位算出部と、を有するものである。
【０００７】
また、本発明に係る内部波変位シミュレーション方法は、海面に平行で且つ直交する二方向における水平波数の各成分を外部から入力し、入力した水平波数の各成分から複数の合成水平波数を算出する水平波数算出工程と、水平波数算出工程において算出した各合成水平波数と、外部から入力された浮力周波数およびコリオリパラメータとを用いて、モード番号に対する複数の固有深度関数および複数の固有周波数を求める内部波モード演算工程と、内部波モード演算工程において算出された各固有周波数を用いて、各固有周波数およびモード番号に対するスペクトル特性を各合成水平波数およびモード番号に対する内部波スペクトル特性へ変換する際に使用する微分値を各固有周波数の差分と各合成水平波数の差分を用いて算出する微分値算出工程と、内部波モード演算工程において算出された各固有周波数と、浮力周波数およびコリオリパラメータとを用いてスペクトル特性を算出するスペクトル特性算出工程と、スペクトル特性算出工程において算出されたスペクトル特性を、微分値を用いて内部波スペクトル特性に変換するスペクトル特性変換工程と、スペクトル特性変換工程における変換後の内部波スペクトル特性と、内部波モード演算工程において算出した各固有深度関数および各固有周波数と、水平波数の各成分とを用いて、海中の等密度面の変位である内部波変位を算出する内部波変位算出工程と、を有するものである。
【発明の効果】
【０００８】
本発明は、近似によらず、固有周波数を用いて算出した微分値を用いて内部波スペクトルを算出し、内部波モード演算部において算出された各固有深度関数および各固有周波数と、内部波スペクトルとを用いて内部波変位を算出することができるため、海洋内部波による等密度面の変位のシミュレーション結果に生じる誤差を低減することができる。
【図面の簡単な説明】
【０００９】
【図１】本発明の第１実施形態に係る内部波変位シミュレーション装置の構成を示すブロック図である。
【図２】図１の内部波変位シミュレーション装置の動作を示すフローチャートである。
【図３】従来における固有周波数と水平波数との関係を示すグラフである。
【図４】従来における固有周波数と水平波数の微分値との関係を示すグラフである。
【図５】本発明の第１実施形態における固有周波数と合成水平波数との関係を示すグラフである。
【図６】本発明の第１実施形態における固有周波数と合成水平波数の微分値との関係を示すグラフである。
【図７】本発明の第２実施形態に係る内部波変位シミュレーション装置の構成を示すブロック図である。
【発明を実施するための形態】
【００１０】
［第１実施形態］
図１は、本発明の第１実施形態に係る内部波変位シミュレーション装置１００の構成を示すブロック図である。内部波変位シミュレーション装置１００は、下記式１に基づく演算処理により、海洋内部波（以下「内部波」という。）による海中の等密度面のランダムな変位（内部波変位）ζをシミュレーションするものである。
【００１１】
【数１】
0006547425.tif 000002
【００１２】
ここで、海面に平行で且つ直交する二方向をｘ方向およびｙ方向とし、海面に垂直な方向をｚ方向と定義する。ｘ方向およびｙ方向は、海面に平行な面上の方向（水平方向）において任意に設定される。すなわち、式１において、ｘおよびｙは水平方向における位置であり、ｚは深度である。また、ｔは時間、Ｒｅは複素数の実部をとる関数である。
【００１３】
そして、αは内部波の水平波数である。αｘおよびαｙは、それぞれ、内部波の水平波数α（以下「水平波数α」という。）のｘ方向成分およびｙ方向成分であり、α２＝αｘ２＋αｙ２の関係がある。なお、内部波変位シミュレーション装置１００は、後述する通り、海面に平行で且つ直交する二方向における水平波数の各成分から複数の合成水平波数αｋｌを算出し、水平波数αとして各種演算に用いるように構成されている。
【００１４】
さらに、ｊはモード番号である。モード番号ｊは、予め設定された最大個数Ｊを上限とする任意の自然数をとる（ｊ＝１、２、…、Ｊ）。ｇ（α，ｊ）は、平均が０で、内部波スペクトル特性Ｇ（α，ｊ）を分散とするランダム変数である。Ｗαｊは、水平波数αおよびモード番号ｊに対する固有深度関数であり、ωαｊは、水平波数αおよびモード番号ｊに対する固有周波数である。
【００１５】
より具体的に、内部波変位シミュレーション装置１００は、式１の積分を加算で表した後述する式２７によって、海中の等密度面の変位である内部波変位ζを算出するものである。内部波変位シミュレーション装置１００は、図１に示すように、式２７を構成する各要素を導出するための処理を実行する内部波モード演算部１０、分散特性演算部２０、スペクトル特性演算部３０、および内部波変位算出部４０を有している。また、内部波変位シミュレーション装置１００は、内部波モード演算部１０、分散特性演算部２０、スペクトル特性演算部３０、および内部波変位算出部４０が各種演算に用いる情報を記憶する記憶部５０を有している。
【００１６】
なお、内部波モード演算部１０、分散特性演算部２０、スペクトル特性演算部３０、および内部波変位算出部４０は、これらの機能を実現する回路デバイスなどのハードウェアで実現することもできるし、例えばＤＳＰ等のマイコン又はＣＰＵ等の演算装置上で実行されるソフトウェアとして実現することもできる。また、記憶部５０は、ＨＤＤ（Ｈａｒｄ  Ｄｉｓｋ  Ｄｒｉｖｅ）又はフラッシュメモリ等により構成することができる。
【００１７】
記憶部５０には、水平波数αのｘ方向成分であるｘ系列αｘ，ｋ（ｋ＝１、２、…、Ｋ）、水平波数αのｙ方向成分であるｙ系列αｙ，ｌ（ｌ＝１、２、…、Ｌ）、浮力周波数（ブラント−バイサラ周波数）Ｎ（ｚ）、およびコリオリパラメータ（慣性周波数）ｆが、外部より入力されて記憶されている。コリオリパラメータｆは、Ωを地球の自転角速度とし、φを緯度とした場合に、ｆ＝２Ωｓｉｎ（φ）によって計算されるものである。ｘ系列αｘ，ｋは、予め決めた間隔Δαｘに対してΔαｘ＝αｘ，ｋ＋１−αｘ，ｋ（ｋ＝１、２、…、Ｋ−１）を満たす。ｙ系列αｙ，ｌは、予め決めた間隔Δαｙに対して、Δαｙ＝αｙ，ｌ＋１−αｙ，ｌ（ｌ＝１、２、…、Ｌ）を満たす。
【００１８】
内部波モード演算部１０は、記憶部５０から、ｘ系列αｘ，ｋ、ｙ系列αｙ，ｌ、浮力周波数Ｎ（ｚ）、およびコリオリパラメータｆを取得し、取得したｘ系列αｘ，ｋおよびｙ系列αｙ，ｌを用いて、下記式２により、複数の合成水平波数αｋｌを算出するものである。
【００１９】
【数２】
0006547425.tif 000003
【００２０】
また、内部波モード演算部１０は、式２により算出した合成水平波数αｋｌと、記憶部５０から取得した浮力周波数Ｎ（ｚ）およびコリオリパラメータｆとを用いて、下記式３により、合成水平波数αｋｌおよびモード番号ｊに対するＪ個の固有深度関数Ｗｋｌｊ（ｊ＝１、２、…、Ｊ）と、合成水平波数αｋｌおよびモード番号ｊに対するＪ個の固有周波数ωｋｌｊ（ｊ＝１、２、…、Ｊ）とを算出するものである。すなわち、固有深度関数Ｗｋｌｊおよび固有周波数ωｋｌｊは、合成水平波数αｋｌに対して、式３に示す微分方程式を満足する解である。ここで、Ｊは予め設定された最大個数であり、Ｈは水深である。
【００２１】
【数３】
0006547425.tif 000004
【００２２】
そして、新たな変数λｋｌｊを下記式４のように定義すると、式３は下記式５のように書き換えられる。
【００２３】
【数４】
0006547425.tif 000005
【００２４】
【数５】
0006547425.tif 000006
【００２５】
ここで、深度ｚを予め決めた間隔Δｚで離散化して、ｍΔｚ（ｍ＝１、２、…、Ｍ）により定義したサンプル深度ｚｍ（ｚｍ＝ｍΔｚ）を適用し、固有深度関数Ｗｋｌｊ（ｚｍ) の２階微分を下記式６のように差分で近似する。以下では、固有深度関数Ｗｋｌｊ（ｚｍ)をＷｋｌｊ，ｍとも表記する。そして、式６を境界条件である「Ｗｋｌｊ，０＝Ｗｋｌｊ，Ｍ＝０」と合わせてｍについて連立させると、下記式７の行列表現が得られる。
【００２６】
【数６】
0006547425.tif 000007
【００２７】
【数７】
0006547425.tif 000008
【００２８】
また、行列Ｗは、下記式８に示すＷｋｌｊであり、Ｔは、行列およびベクトルの転置を表す記号である。すなわち、式７の行列Ｗは、行列［Ｗｋｌｊ，１，Ｗｋｌｊ，２，…，Ｗｋｌｊ，Ｍ−１］の転置行列である。さらに、行列Ｂは、下記式９により求めることができる。なお、式９に関して、ｄｉａｇ（ｘ）は、ベクトルｘの要素を対角成分とする対角行列を表す。
【００２９】
【数８】
0006547425.tif 000009
【００３０】
【数９】
0006547425.tif 000010
【００３１】
ここで、式７の両辺にＢ−1をかけると、下記式１０となるため、新たに下記式１１および式１２の通りにＶｋｌｊおよび行列Ｃを定義すると、式７は下記式１３のように書き換えられる。そして、式１３を満足するＶｋｌｊおよびλｋｌｊは、標準的な固有値問題を解くことにより求めることができる。
【００３２】
【数１０】
0006547425.tif 000011
【００３３】
【数１１】
0006547425.tif 000012
【００３４】
【数１２】
0006547425.tif 000013
【００３５】
【数１３】
0006547425.tif 000014
【００３６】
上記の行列表現に基づく数式処理をもとに、内部波モード演算部１０は、あるｋおよびｌに対する固有深度関数Ｗｋｌｊ（ｚｍ）および固有周波数ωｋｌｊを、以下のように算出するものである。すなわち、内部波モード演算部１０は、行列Ｃの第ｍ行第ｎ列の要素Ｃｍ，ｎを、式１２より下記式１４のように計算する。また、内部波モード演算部１０は、行列Ｂを式９により求める。
【００３７】
【数１４】
0006547425.tif 000015
【００３８】
さらに、内部波モード演算部１０は、式１３を満足するＶｋｌｊおよびλｋｌｊを、標準的な固有値問題を解くことにより、λの小さな方からＪ個求める。そして、内部波モード演算部１０は、式１３により求めたＶｋｌｊを用いて、式１１を変形した下記式１５により、固有深度関数Ｗｋｌｊ（ｚｍ）を算出する。
【００３９】
【数１５】
0006547425.tif 000016
【００４０】
また、内部波モード演算部１０は、式４を変形した下記式１６により固有周波数ωｋｌｊを算出する。
【００４１】
【数１６】
0006547425.tif 000017
【００４２】
内部波モード演算部１０は、以上の処理を、全てのｋ＝１、２、…、Ｋおよびｌ＝１、２、…、Ｌに対して行うように構成されている。そして、内部波モード演算部１０は、各合成水平波数αｋｌを、分散特性演算部２０に出力するものである。また、内部波モード演算部１０は、浮力周波数Ｎ（ｚｍ）に対するＪ個の固有周波数ωｋｌｊを、分散特性演算部２０、スペクトル特性演算部３０、内部波変位算出部４０に出力するものである。さらに、内部波モード演算部１０は、浮力周波数Ｎ（ｚｍ）に対するＪ個の固有深度関数Ｗｋｌｊ（ｚｍ）を、内部波変位算出部４０に出力するものである。
【００４３】
分散特性演算部２０は、内部波モード演算部１０から出力されるそれぞれの固有周波数ωｋｌｊに対する微分ｄα／ｄωｋｌｊの値（微分値）Ｄｋｌｊを算出し、算出した微分値Ｄｋｌｊをスペクトル特性演算部３０に出力するものである。なお、微分値Ｄｋｌｊは、後述するように、スペクトル特性演算部３０が、各固有周波数ωｋｌｊおよびモード番号ｊに対するスペクトル特性Ｆ（ωｋｌｊ，ｊ）を、各合成水平波数αｋｌおよびモード番号ｊに対する内部波スペクトル特性Ｇ（αｋｌ，ｊ）へ変換する際に使用される。
【００４４】
より具体的に、分散特性演算部２０は、以下の演算により、合成水平波数αにおける固有周波数ωでの微分ｄα／ｄωを計算するものである。分散特性演算部２０は、固有周波数ωｋｌｊと合成水平波数αｋｌとの組（ωｋｌｊ，αｋｌ）を、モード番号ｊごとに固有周波数ωｋｌｊを基準として昇順に並べる。以下では、昇順に並べたものを(ωｐｊ，αｐ)とする（ｐ＝１、２、…、Ｋ×Ｌ）。また、分散特性演算部２０は、並べ替えの前後におけるインデックスの関係（（ｋ，ｌ）とｐとの関係）を記憶部５０又は内部メモリ等（図示せず）に記録させる。
【００４５】
さらに、分散特性演算部２０は、それぞれのｐ（ｐ＝１、２、…、Ｋ×Ｌ）に対する微分ｄα／ｄωの固有周波数ωｐｊでの微分値Ｄｐｊを算出する。すなわち、分散特性演算部２０は、ｐ（ｐ＝２、３、...、Ｋ×Ｌ−１）に対する微分値Ｄｐｊとして、下記式１７に示す前進差分商Δ＋と下記式１８に示す後退差分商Δ−との加重平均値を下記式１９により算出する。前進差分商Δ＋は、式１７に示すように、固有周波数ωｐｊの前進差分に対する合成水平波数αｐの前進差分の比率であり、後退差分商Δ−は、式１８に示すように、固有周波数ωｐｊの後退差分に対する合成水平波数αｐの後退差分の比率である。また、分散特性演算部２０は、ｐ＝１に対する微分値Ｄｐｊを下記式２０により算出し、Ｐ＝Ｋ×Ｌに対する微分値Ｄｐｊを下記式２１により算出する。
【００４６】
【数１７】
0006547425.tif 000018
【００４７】
【数１８】
0006547425.tif 000019
【００４８】
【数１９】
0006547425.tif 000020
【００４９】
【数２０】
0006547425.tif 000021
【００５０】
【数２１】
0006547425.tif 000022
【００５１】
そして、分散特性演算部２０は、記録しておいたインデックスの関係（（ｋ，ｌ）とｐとの関係）を用いて微分値Ｄｐｊを並べ替えることにより微分値Ｄｋｌｊを算出する。分散特性演算部２０は、以上の処理を全てのモード番号ｊに対して行うように構成されている。また、分散特性演算部２０は、内部波モード演算部１０から出力された固有周波数ωｋｌｊに対する微分値Ｄｋｌｊをスペクトル特性演算部３０に出力するものである。
【００５２】
スペクトル特性演算部３０は、内部波モード演算部１０から出力される固有周波数ωｋｌｊを用いて、固有周波数ωｋｌｊおよびモード番号ｊに対するスペクトル特性Ｆ（ωｋｌｊ，ｊ）を算出するものである。また、スペクトル特性演算部３０は、算出したスペクトル特性Ｆ（ωｋｌｊ，ｊ）を、合成水平波数αｋｌおよびモード番号ｊに対する内部波スペクトル特性Ｇ（αｋｌ，ｊ）に変換し、変換後の内部波スペクトル特性Ｇ（αｋｌ，ｊ）を内部波変位算出部４０に出力するものである。
【００５３】
ところで、スペクトル特性Ｆ（ω，ｊ）は、固有周波数ωおよびモード番号ｊに対して与えられた下記式２２によって求めることができる。
【００５４】
【数２２】
0006547425.tif 000023
【００５５】
ｊ＊は、モード帯域幅を表すパラメータであり、ｊ＊の値を大きくすると、Ｆ（ω，ｊ）のモード番号に対するスペクトル特性が緩やかになる。一般に、深海域では、ｊ＊の値が３に設定される。ｂおよびＮ０は、浮力周波数Ｎ（ｚ）が下記式２３で表されるときの定数であり、ｂはＮ（ｚ）の深度に対する低下度合いを決めるパラメータ、Ｎ０は海面での浮力周波数である。Ｅは、内部波のエネルギーパラメータと呼ばれる定数である。
【００５６】
【数２３】
0006547425.tif 000024
【００５７】
任意のＮ（ｚ）については、式２３の関係から、水深Ｈが十分大きいと仮定すると、下記式２４の関係が得られる。このため、本第１実施形態において、スペクトル特性演算部３０は、式２２を下記式２５のように変形して使用する。
【００５８】
【数２４】
0006547425.tif 000025
【００５９】
【数２５】
0006547425.tif 000026
【００６０】
すなわち、スペクトル特性演算部３０は、記憶部５０から浮力周波数Ｎ（ｚ）およびコリオリパラメータｆを取得し、取得した浮力周波数Ｎ（ｚ）およびコリオリパラメータｆと、内部波モード演算部１０から出力される固有周波数ωｋｌｊとを用いて、式２５により、固有周波数ωｋｌｊおよびモード番号ｊに対するスペクトル特性Ｆ（ωｋｌｊ，ｊ）を算出するものである。ここで、式２５の要素である積ｂＥは、内部波による変位の大きさを調整するパラメータとなり、実際に測定された変位などから決めることができる。
【００６１】
また、スペクトル特性演算部３０は、下記式２６により、スペクトル特性Ｆ（ωｋｌｊ，ｊ）を、内部波による等密度面の変位をシミュレーションする際に必要となる内部波スペクトル特性Ｇ（α，ｊ）に変換するように構成されている。そして、スペクトル特性演算部３０は、式２６による変換で求めた内部波スペクトル特性Ｇ（αｋｌ，ｊ）を内部波変位算出部４０に出力するように構成されている。
【００６２】
【数２６】
0006547425.tif 000027
【００６３】
内部波変位算出部４０は、記憶部５０からｘ系列αｘ，ｋおよびｙ系列αｙ，ｌを取得し、ｘ系列αｘ，ｋおよびｙ系列αｙ，ｌと、内部波モード演算部１０から出力される固有深度関数Ｗｋｌｊ（ｚｍ）と、スペクトル特性演算部３０から出力される内部波スペクトル特性Ｇ（αｋｊ，ｊ）とを用いて、下記式２７により、位置（ｘ，ｙ，ｚｍ）および時間ｔに対する内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ）を算出するものである。また、内部波変位算出部４０は、算出した内部波変位ζ（ｘ，ｙ，ｚｍ，ｔ）を外部に出力するものである。
【００６４】
【数２７】
0006547425.tif 000028
【００６５】
ここで、ランダム変数ｇ（αｋｌ，ｊ）は、平均が０であり、分散が内部波スペクトル特性Ｇ（αｋｌ，ｊ）の正規分布に従うランダムな値をとる。また、ΔαｘおよびΔαｙは、それぞれ、合成水平波数αｋｌのｘ方向成分およびｙ方向成分の間隔である。
【００６６】
続いて、本第１実施形態における内部波変位シミュレーション装置１００の動作を、図２に示すフローチャートに基づいて説明する。
【００６７】
まず、内部波モード演算部１０は、記憶部５０から、ｘ系列αｘ，ｋ、ｙ系列αｙ，ｌ、浮力周波数Ｎ（ｚ）、およびコリオリパラメータｆを取得する。そして、内部波モード演算部１０は、記憶部５０から取得したｘ系列αｘ，ｋおよびｙ系列αｙ，ｌを用いて、式２により、合成水平波数αｋｌを算出する（図２：ステップＳ１０１／水平波数算出工程）。次いで、内部波モード演算部１０は、算出した合成水平波数αｋｌと、記憶部５０から取得した浮力周波数Ｎ（ｚ）およびコリオリパラメータｆとを用いて、式３により、浮力周波数Ｎ（ｚ）に対するＪ個の固有深度関数Ｗｋｌｊおよび固有周波数ωｋｌｊを算出する。すなわち、内部波モード演算部１０は、式１５により固有深度関数Ｗｋｌｊ（ｚｍ）を算出し、式１６により固有周波数ωｋｌｊを算出する（図２：ステップＳ１０２／内部波モード演算工程）。
【００６８】
内部波モード演算部１０は、全てのｋ＝１、２、…、Ｋおよびｌ＝１、２、…、Ｌに対して、上記ステップＳ１０１およびＳ１０２の処理を行う。そして、内部波モード演算部１０は、算出した合成水平波数αｋｌおよびＪ個の固有周波数ωｋｌｊを分散特性演算部２０に出力する。また、内部波モード演算部１０は、算出したＪ個の固有周波数ωｋｌｊをスペクトル特性演算部３０に出力する。さらに、内部波モード演算部１０は、算出したＪ個の固有周波数ωｋｌｊおよびＪ個の固有深度関数Ｗｋｌｊ（ｚｍ）を内部波変位算出部４０に出力する（図２：ステップＳ１０３）。
【００６９】
分散特性演算部２０は、内部波モード演算部１０から出力される合成水平波数αｋｌと固有周波数ωｋｌｊと入力する。そして、分散特性演算部２０は、内部波モード演算部１０から出力される固有周波数ωｋｌｊと合成水平波数αｋｌとの組（ωｋｌｊ，αｋｌ）を、モード番号ｊごとに昇順に並べ替えて（ωｐｊ，αｐ)（ｐ＝１、２、…、Ｋ×Ｌ）とする（図２：ステップＳ１０４）。また、分散特性演算部２０は、並べ替えの前後におけるインデックスの関係（（ｋ，ｌ）とｐとの関係）を記憶部５０又は内部メモリ等（図示せず）に記録させる（図２：ステップＳ１０５）。
【００７０】
次に、分散特性演算部２０は、ｐ（ｐ＝２、３、...、Ｋ×Ｌ−１）に対する微分値Ｄｐｊを式１９により算出する。また、分散特性演算部２０は、ｐ＝１に対する微分値Ｄｐｊを式２０により算出し、Ｐ＝Ｋ×Ｌに対する微分値Ｄｐｊを式２１により算出する（図２：ステップＳ１０６）。次いで、分散特性演算部２０は、記録させておいたインデックスの関係（（ｋ，ｌ）とｐとの関係）を用いて微分値Ｄｐｊを並べ替えることにより微分値Ｄｋｌｊを算出する（図２：ステップＳ１０７／微分値算出工程）。分散特性演算部２０は、上記ステップＳ１０４〜Ｓ１０７までの処理を、全てのモード番号ｊに対して行い、算出した微分値Ｄｋｌｊをスペクトル特性演算部３０に出力する（図２：ステップＳ１０８）。
【００７１】
スペクトル特性演算部３０は、内部波モード演算部１０から出力される固有周波数ωｋｌｊと、分散特性演算部２０から出力される微分値Ｄｋｌｊとを入力し、記憶部５０から、浮力周波数Ｎ（ｚ）およびコリオリパラメータｆを取得する。そして、スペクトル特性演算部３０は、固有周波数ωｋｌｊ、浮力周波数Ｎ（ｚ）、およびコリオリパラメータｆを用いて、式２５により、固有周波数ωｋｌｊおよびモード番号ｊに対するスペクトル特性Ｆ（ωｋｌｊ，ｊ）を算出する（図２：ステップＳ１０９／スペクトル特性算出工程）。
【００７２】
次に、スペクトル特性演算部３０は、算出したスペクトル特性Ｆ（ωｋｌｊ，ｊ）を、微分値Ｄｋｌｊを用いて式２６により内部波スペクトル特性Ｇ（αｋｌ，ｊ）へ変換し（図２：ステップＳ１１０／スペクトル特性変換工程）、変換後の内部波スペクトル特性Ｇ（αｋｌ，ｊ）を内部波変位算出部４０に出力する（図２：ステップＳ１１１）。
【００７３】
内部波変位算出部４０は、内部波モード演算部１０から出力される固有深度関数Ｗｋｌｊ（ｚｍ）および固有周波数ωｋｌｊ（ｊ＝１、２、…、Ｊ）と、スペクトル特性演算部３０から出力される内部波スペクトル特性Ｇ（αｋｊ，ｊ）とを入力し、記憶部５０から、水平波数αのｘ方向成分であるｘ系列αｘ，ｋとｙ方向成分であるｙ系列αｙ，ｌを取得する。そして、内部波変位算出部４０は、内部波モード演算部１０から出力される固有深度関数Ｗｋｌｊ（ｚｍ）および固有周波数ωｋｌｊと、スペクトル特性演算部３０から出力される内部波スペクトル特性Ｇ（αｋｊ，ｊ）と、ｘ系列αｘ，ｋおよびｙ系列αｙ，ｌとを用いて、式２７により、位置（ｘ，ｙ，ｚｍ）および時間ｔに対する内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)を算出する（図２：ステップＳ１１２／内部波変位算出工程）。そして、内部波変位算出部４０は、算出した内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)を外部に出力する（図２：ステップＳ１１３）。
【００７４】
（第１実施形態の効果）
本第１実施形態における内部波変位シミュレーション装置１００は、内部波モード演算部１０において算出された固有周波数ωｋｌｊを用いて分散特性演算部２０が微分値Ｄｋｌｊを算出し、分散特性演算部２０において算出された微分値Ｄｋｌｊを用いてスペクトル特性演算部３０がスペクトル特性Ｆ（ωｋｌｊ，ｊ）を内部波スペクトル特性Ｇ（αｋｌ，ｊ）に変換するという構成を採っている。さらに、内部波変位算出部４０は、内部波モード演算部１０から出力される固有深度関数Ｗｋｌｊ（ｚｍ）および固有周波数ωｋｌｊと、スペクトル特性演算部３０から出力される内部波スペクトル特性Ｇ（αｋｊ，ｊ）とを用いて内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)を算出するように構成されている。よって、内部波変位シミュレーション装置１００によれば、従来のような近似的処理を用いずに、内部波スペクトル特性Ｇ（αｋｌ，ｊ）および内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)を算出することができるため、海洋内部波による等密度面の変位のシミュレーション結果に生じる誤差を低減することができる。
【００７５】
ここで、図３〜図６を参照して、内部波変位シミュレーション装置１００によって得られる効果を詳細に説明する。図３は、従来における固有周波数と水平波数との関係を示すグラフである。図４は、従来における固有周波数と水平波数の微分値との関係を示すグラフである。図５は、本第１実施形態における固有周波数と合成水平波数との関係を示すグラフである。図６は、本第１実施形態における固有周波数と水平波数の微分値との関係を示すグラフである。
【００７６】
従来から、式１を使って内部波による等密度面の変位をシミュレーションする場合は、下記式２８により、スペクトル特性Ｆ（ω，ｊ）から内部波スペクトル特性Ｇ（α，ｊ）への変換を行っている。そして、従来においては、水平波数αと固有周波数ωとの関係を、以下のように近似的に表すことで微分値ｄα／ｄωを計算し、式２８による変換を行っている。
【００７７】
【数２８】
0006547425.tif 000029
【００７８】
式３の左辺のＷαｊ（Ｚ）にかかっている項は、固有深度関数の深度方向の波数βαｊ（以下「鉛直波数βαｊ」という。）により、下記式２９のように表すことができる。
【００７９】
【数２９】
0006547425.tif 000030
【００８０】
また、第ｊモードに対する鉛直波数βαｊについては、下記式３０に示す近似が成り立つため、式２９および式３０により、下記式３１の関係を導出する。さらに、式３１をαについて解くことにより、下記式３２の関係を得る。
【００８１】
【数３０】
0006547425.tif 000031
【００８２】
【数３１】
0006547425.tif 000032
【００８３】
【数３２】
0006547425.tif 000033
【００８４】
ここで、水平波数αは、深度ｚに依存しない量であるため、ｆ＜ωαｊ≪Ｎ（ｚ）と仮定して、下記式３３のように近似すると、微分値ｄα／ｄωαｊは、下記式３４のように表すことができる。
【００８５】
【数３３】
0006547425.tif 000034
【００８６】
【数３４】
0006547425.tif 000035
【００８７】
すなわち、従来は、式３３および式３４を使って、式２８による変換を行っている。なお、式３３および式３４の右辺にあるｂＮ０は、式２４の関係によって置き換えることができる。
【００８８】
しかしながら、式３３および式３４は、水平波数αと固有周波数ωとの近似的な関係を表しているに過ぎない。すなわち、従来は、スペクトル特性Ｆ（ω，ｊ）を内部波スペクトル特性Ｇ（α，ｊ）へ変換する際に介在する近似的な処理の影響により、最終的には、式１でシミュレーションする海中の等密度面の変位に誤差が生じるという課題がある。
【００８９】
ここで、従来構成による内部波変位の算出誤差について具体的に説明する。例えば、簡単のために、浮力周波数Ｎ（ｚ）が深度に依存せず、Ｎ＝Ｎ０で表せる場合を想定する。このとき、式３の微分方程式のｊ番目の解は、下記式３５で表される。∝は比例関係を示す記号である。また、鉛直波数βαｊは下記式３６で表されるため、式２９の関係から、下記式３７および式３８を得ることができる。
【００９０】
【数３５】
0006547425.tif 000036
【００９１】
【数３６】
0006547425.tif 000037
【００９２】
【数３７】
0006547425.tif 000038
【００９３】
【数３８】
0006547425.tif 000039
【００９４】
図３は、式３３による算出値と式３７による算出値とを比較した結果を例示するグラフであり、図３（Ａ）にｊ＝１での関係を示し、図３（Ｂ）にｊ＝５での関係を示す。図４は、式３４による算出値と式３８による算出値とを比較した結果を例示するグラフであり、図４（Ａ）にｊ＝１での関係を示し、図４（Ｂ）にｊ＝５における関係を示す。
【００９５】
また、図３では、式３３による算出値を破線で示し、式３７による算出値を実線で示す。図４では、式３４による算出値を破線で示し、式３８による算出値を実線で示す。図３および図４は、Ｈ＝５０００ｍ、Ｎ０＝０．００５２４ｒａｄ／ｓ（＝３．０ｃｐｈ）、ｆ＝７．２９×１０−５ｒａｄ／ｓの条件の下で計算したものである。図３および図４において、横軸は、固有周波数ωを１時間あたりのサイクル数（ｃｐｈ）で表しており、縦軸は、対数スケールをとっている。
【００９６】
図３に示すように、縦軸に示す水平波数αについては、横軸の固有周波数ωが大きくなりＮ０（３．０ｃｐｈ）に近づくにつれて、式３３による算出結果と式３７による算出結果との差が大きくなることがわかる。また、図４に示すように、縦軸に示す微分値ｄα／ｄωについても、横軸の固有周波数ωが大きくなりＮ０（３．０ｃｐｈ）に近づくにつれて、式３４による算出結果と式３８による算出結果との差が大きくなることがわかる。すなわち、従来においては、固有周波数ωが大きくなるにつれて、内部波の変位に生じる誤差が大きくなっている。
【００９７】
ところで、式３７および式３８を導出する際の想定（Ｎ＝Ｎ０）のように、式３の微分方程式を解析的に解ける形でＮ（ｚ）が与えられた場合には、水平波数αと固有周波数ωとの関係を厳密に表すことができるが、一般にそれは不可能である。したがって、水平波数αと固有周波数ωとの関係を式３３および式３４を使って近似的に表すと、内部波の変位に誤差が生じる。
【００９８】
図５は、内部波モード演算部１０が算出する合成水平波数αと固有周波数ωとの関係と、式３７における水平波数αと固有周波数ωとの関係とを比較した結果を例示するグラフであり、図５（Ａ）にはｊ＝１での関係を示し、図５（Ｂ）にはｊ＝５での関係を示す。図６は、分散特性演算部２０が算出する微分値ｄα／ｄω（微分値Ｄｋｌｊ）と固有周波数ωとの関係と、式３８における微分値ｄα／ｄωと固有周波数ωとの関係とを比較した結果を例示するグラフであり、図６（Ａ）にはｊ＝１での関係を示し、図６（Ｂ）にはｊ＝５における関係を示す。
【００９９】
また、図５では、内部波モード演算部１０による算出値を丸印でプロットしており、式３７による算出値を実線で示している。図６では、内部波モード演算部１０による算出値を丸印でプロットしており、式３８による算出値を実線で示している。図５および図６は、図３および図４の場合と同じ条件下におけるデータを示すものである。
【０１００】
上述した図３および図４では、近似式である式３３および式３４により求めた結果と、式３７および式３８により求めた結果との間には差が生じており、その差は、固有周波数ωが大きくなるにつれて大きくなっている。これに対し、本第１実施形態では、図５および図６に示すように、内部波変位シミュレーション装置１００が求めた結果と、式３７および式３８により求めた結果との間に生じる差が小さくなっている。すなわち、従来の構成と比較して、合成水平波数αｋｌと固有周波数ωｋｌｊとの関係について生じる誤差と、微分値Ｄｋｌｊと固有周波数ωｋｌｊとの関係について生じる誤差が小さくなっていることがわかる。
【０１０１】
以上のように、内部波変位シミュレーション装置１００は、内部波モード演算部１０が合成水平波数αｋｌと固有周波数ωｋｌｊとの関係を演算し、内部波モード演算部１０による演算結果をもとに分散特性演算部２０が微分値Ｄｋｌｊを算出し、分散特性演算部２０において算出された微分値Ｄｋｌｊを用いてスペクトル特性演算部３０が内部波スペクトル特性Ｇ（α，ｊ）を求めるように構成されている。したがって、合成水平波数αｋｌと固有周波数ωｋｌｊとの関係および微分値Ｄｋｌｊを、近似によらず、内部波モードが満足する微分方程式を解くことで得られる固有深度関数Ｗｋｌｊ（ｚｍ）および固有周波数ωｋｌｊとの関係を用いて算出ことができるため、内部波スペクトル特性Ｇ（α，ｊ）および内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)の算出誤差を小さくすることができる。すなわち、内部波変位シミュレーション装置１００によれば、近似式を用いた場合に生じる内部波スペクトル特性Ｇ（α，ｊ）および内部波変位ζの算出誤差を小さくすることができるため、海洋内部波による等密度面の変位のシミュレーション結果に生じる誤差を低減することができる。
【０１０２】
＜変形例＞
前述の第１実施形態は、式２８に相当する式２６を使って、スペクトル特性Ｆ（ω，ｊ）から内部波スペクトル特性Ｇ（α，ｊ）への変換を行うことが前提となっており、かかる前提のもと、固有周波数ωと合成水平波数αとの関係を用いた処理を説明するものである。
【０１０３】
ところで、固有周波数ωと合成水平波数αとの関係は、内部波スペクトル特性Ｇ（α，ｊ）をスペクトル特性Ｆ（ω，ｊ）に変換する場合にも使用することができる。すなわち、内部波スペクトル特性Ｇ（α，ｊ）は、下記式３９を使って、スペクトル特性Ｆ（ω，ｊ）に変換することができる。
【０１０４】
【数３９】
0006547425.tif 000040
【０１０５】
例えば、内部波の変位の空間的な変動が観測されたとき、その変動をスペクトル分析して得られるものは、合成水平波数αに対する内部波スペクトル特性Ｇ（α，ｊ）である。したがって、観測された変動を既存の周波数に対するスペクトル特性と比較したいときは、得られた内部波スペクトル特性Ｇ（α，ｊ）をスペクトル特性Ｆ（ω，ｊ）に変換するとよい。すなわち、内部波の変位の空間的な変動が観測されたときには、式３９による変換処理により、観測された変動を既存の周波数に対するスペクトル特性と比較することができる。
【０１０６】
［第２実施形態］
次に、本発明の第２実施形態に係る内部波変位シミュレーション装置２００について説明する。前述の第１実施形態では、式２８に相当する式２６を使ってスペクトル特性Ｆ（ω，ｊ）を内部波スペクトル特性Ｇ（α，ｊ）に変換する例を示したが、本第２実施形態では、内部波変位シミュレーション装置２００が、式２８と同値な下記式４０を使ってスペクトル特性Ｆ（ω，ｊ）を内部波スペクトル特性Ｇ（α，ｊ）に変換するという点に特徴がある。第１実施形態と同等の構成部材については、同一の符号を用いて説明は省略する。
【０１０７】
【数４０】
0006547425.tif 000041
【０１０８】
分散特性演算部１２０は、上記式１７〜式２１の固有周波数ωと水平波数αとを入れ換えた式により微分値Ｄｐｊを算出するものである。すなわち、分散特性演算部１２０は、固有周波数ωｋｌｊと合成水平波数αｋｌとの組（ωｋｌｊ，αｋｌ）を、モード番号ｊごとに昇順に並べて (ωｐｊ，αｐ)とする（ｐ＝１、２、…、Ｋ×Ｌ）。また、分散特性演算部１２０は、並べ替えの前後におけるインデックスの関係（（ｋ，ｌ）とｐとの関係）を記憶部５０又は内部メモリ等（図示せず）に記録させる。
【０１０９】
さらに、分散特性演算部１２０は、ｐ（ｐ＝２、３、...、Ｋ×Ｌ−１）に対する微分値Ｄｐｊとして、下記式４１に示す前進差分商Δ＋と下記式４２に示す後退差分商Δ−との加重平均値を下記式４３により算出する。前進差分商Δ＋は、式４１に示すように、合成水平波数αｐの前進差分に対する固有周波数ωｐｊの前進差分の比率であり、後退差分商Δ−は、式４２に示すように、合成水平波数αｐの後退差分に対する固有周波数ωｐｊの後退差分の比率である。また、分散特性演算部１２０は、ｐ＝１に対する微分値Ｄｐｊを下記式４４により算出し、Ｐ＝Ｋ×Ｌに対する微分値Ｄｐｊを下記式４５により算出する。
【０１１０】
【数４１】
0006547425.tif 000042
【０１１１】
【数４２】
0006547425.tif 000043
【０１１２】
【数４３】
0006547425.tif 000044
【０１１３】
【数４４】
0006547425.tif 000045
【０１１４】
【数４５】
0006547425.tif 000046
【０１１５】
そして、分散特性演算部１２０は、記録しておいたインデックスの関係（（ｋ，ｌ）とｐとの関係）を用いて微分値Ｄｐｊを並べ替えることにより微分値Ｄｋｌｊを算出する。分散特性演算部１２０は、以上の処理を全てのモード番号ｊに対して行い、算出した微分値Ｄｋｌｊをスペクトル特性演算部１３０に出力するように構成されている。
【０１１６】
また、スペクトル特性演算部１３０は、式２６の代わりに、式４０に相当する下記式４６を用いてスペクトル特性Ｆ（ω，ｊ）を内部波スペクトル特性Ｇ（α，ｊ）に変換するものである。なお、他の構成及び動作内容については、前述した第１実施形態と同様である。
【０１１７】
【数４６】
0006547425.tif 000047
【０１１８】
（第２実施形態の効果）
本第２実施形態における内部波変位シミュレーション装置２００は、内部波モード演算部１０において算出された固有周波数ωｋｌｊを用いて分散特性演算部１２０が微分値Ｄｋｌｊを算出することから、スペクトル特性演算部１３０は、近似によらず、分散特性演算部１２０において算出された微分値Ｄｋｌｊを用いてスペクトル特性Ｆ（ωｋｌｊ，ｊ）を内部波スペクトル特性Ｇ（αｋｌ，ｊ）に変換することができる。また、内部波変位算出部４０は、近似によらず、内部波モード演算部１０から出力される固有深度関数Ｗｋｌｊ（ｚｍ）および固有周波数ωｋｌｊと、スペクトル特性演算部１３０から出力される内部波スペクトル特性Ｇ（αｋｊ，ｊ）とを用いて内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)を算出するように構成されている。すなわち、内部波変位シミュレーション装置２００によれば、従来のような近似的処理を用いずに、内部波スペクトル特性Ｇ（αｋｌ，ｊ）および内部波変位ζ(ｘ，ｙ，ｚｍ，ｔ)を算出することができるため、海洋内部波による等密度面の変位のシミュレーション結果に生じる誤差を低減することができる。
【０１１９】
なお、上述した各実施形態は、内部波変位シミュレーション装置および内部波変位シミュレーション方法における好適な具体例であり、本発明の技術的範囲は、特に本発明を限定する記載がない限り、これらの態様に限定されるものではない。例えば、第１又は第２実施形態では、式１７〜式２１又は式４１〜式４５に示す差分により微分値を求める方法を説明したが、これに限定されず、点（ωｐｊ，αｐ）をスプライン補間したときの点ωｐｊにおける微分値を用いる方法等を採用するようにしてもよい。また、記憶部５０として、内部波モード演算部１０の内部メモリ等を機能させるようにしてもよい。さらに、各実施形態では、内部波モード演算部１０が、ｘ系列αｘ，ｋおよびｙ系列αｙ，ｌを用いて、式２により合成水平波数αｋｌを算出するという構成を採ったが、内部波モード演算部１０の外部に、合成水平波数αｋｌを算出する機能をもつ構成部材を設け、該構成部材にて算出された合成水平波数αｋｌを内部波モード演算部１０が取得するようにしてもよい。
【０１２０】
加えて、各実施形態では、内部波変位算出部４０が、記憶部５０から水平波数αのｘ方向成分αｘ，ｋおよびｙ方向成分αｙ，ｌを取得するという構成を例示しているが、これに限定されるものではない。すなわち、例えば、内部波モード演算部１０が、記憶部５０から取得したｘ系列αｘ，ｋおよびｙ系列αｙ，ｌを内部波変位算出部４０に送信するようにしてもよい。また、第１又は第２実施形態では、スペクトル特性演算部３０又は１３０が、記憶部５０から浮力周波数Ｎ（ｚ）およびコリオリパラメータｆを取得するという構成を例示しているが、これに限定されるものではない。すなわち、例えば、内部波モード演算部１０が、記憶部５０から取得した浮力周波数Ｎ（ｚ）およびコリオリパラメータｆをスペクトル特性演算部３０又は１３０へ出力するようにしてもよい。さらに、各実施形態では、内部波モード演算部１０が、記憶部５０から、ｘ系列αｘ，ｋ、ｙ系列αｙ，ｌ、浮力周波数Ｎ（ｚ）、およびコリオリパラメータｆを取得するという構成を例示しているが、これに限定されるものではない。すなわち、ｘ系列αｘ，ｋ、ｙ系列αｙ，ｌ、浮力周波数Ｎ（ｚ）、およびコリオリパラメータｆが、外部から内部波モード演算部１０へ入力されるようにしてもよい。同様に、浮力周波数Ｎ（ｚ）およびコリオリパラメータｆが、外部からスペクトル特性演算部３０および１３０へ入力されるようにしてもよく、ｘ系列αｘ，ｋおよびｙ系列αｙ，ｌが、外部から内部波変位算出部４０へ入力されるようにしてもよい。
【符号の説明】
【０１２１】
１０  内部波モード演算部、２０、１２０  分散特性演算部、３０、１３０  スペクトル特性演算部、４０  内部波変位算出部、５０  記憶部、１００、２００  内部波変位シミュレーション装置。
【図１】
0006547425.tif 000048
【図２】
0006547425.tif 000049
【図３】
0006547425.tif 000050
【図４】
0006547425.tif 000051
【図５】
0006547425.tif 000052
【図６】
0006547425.tif 000053
【図７】
0006547425.tif 000054