(57)【要約】【課題】精度よく波の水位を予測することができる。【解決手段】入力部１０は、各観測地点での波の視線方向の流速の入力を受け付ける。予測部３２は、予測対象地点での波の水位を含む状態を予測する。推定部３４は、各観測地点での波の視線方向の流速の入力を受け付けた場合に、入力された各観測地点での波の視線方向の流速と、観測地点について予測された予測対象地点での水位を含む波の状態を、予め定められた観測行列を用いて変換することにより得られる各観測地点での波の視線方向の流速との差分に基づいて、予測対象地点での水位を含む波の状態を推定する。判定部３６は、状態の予測と、状態の推定とを予め定められた条件を満たすまで繰り返させる。状態の予測は、繰り返しの一つ前に推定された状態、又は繰り返しの一つ前に予測された状態に基づいて予測する。【選択図】図１
2020056649.tif 000002
【特許請求の範囲】
【請求項１】
各観測地点での波の視線方向の流速を入力として、予測対象地点での波の水位を予測する津波予測装置であって、
各観測地点での波の視線方向の流速の入力を受け付ける入力部と、
前記予測対象地点での波の水位を含む状態を予測する予測部と、
各観測地点での波の視線方向の流速の入力を受け付けた場合に、
入力された各観測地点での波の視線方向の流速と、観測地点について予測された前記予測対象地点での水位を含む波の状態を、予め定められた観測行列を用いて変換することにより得られる各観測地点での波の視線方向の流速との差分に基づいて、前記予測対象地点での水位を含む波の状態を推定する推定部と、
前記予測部による前記状態の予測と、前記推定部による前記状態の推定とを予め定められた条件を満たすまで繰り返させる判定部と、
を含み、
前記予測部による前記状態の予測は、繰り返しの一つ前に前記推定部で推定された前記状態、又は繰り返しの一つ前に前記予測部で予測された前記状態に基づいて予測する津波予測装置。
【請求項２】
前記状態は、水位と線流量とを含む請求項１に記載の津波予測装置。
【請求項３】
前記状態は、水位と前記視線方向の線流量と前記視線方向と直交する方向の線流量とを含む請求項２記載の津波予測装置。
【請求項４】
前記観測行列は、前記視線方向の流速と前記視線方向の線流量と静水深との関係を線形近似により求めたものとする請求項３に記載の津波予測装置。
【請求項５】
コンピュータを、請求項１〜請求項４のいずれか１項に記載の津波予測装置の各部として機能させるためのプログラム。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、津波予測装置、方法、及びプログラムに係り、特に、波の水位を予測するための津波予測装置、方法、及びプログラムに関する。
【背景技術】
【０００２】
従来より、津波の波高（波の水位の高さ）を予測する技術がある。例えば、海洋レーダで観測される視線方向の津波流速分布から津波の水位分布を推定し、推定結果を初期条件として津波伝播シミュレーションを行うことで、予測対象地点での津波到達時間、及び津波の水位を予測している。
【先行技術文献】
【特許文献】
【０００３】
【特許文献１】特開２０１６−８５２０６号公報
【発明の概要】
【発明が解決しようとする課題】
【０００４】
ただし、特許文献１の技術では、津波の水位分布の推定に際し、運動方程式を簡略化している等の原因により、津波の波高、及び到達時間に関する予測精度があまり高くないという課題があった。原因の内容としては、例えば、円筒座標系におけるθ方向の流速を入力値として使っていないことがある。また、運動方程式の各項のうち複数の項を無視して解いている。予測精度に関しては、各シミュレーション結果の包絡線で最大水位、最小水位を予測するため、波形を正確に表現できるものではなかった。シミュレーション結果の比較に関しては、最新のシミュレーション結果が最も良い結果になるとは限らないため、過去の予測結果と大きさ等を比較する手間が必要であった。
【０００５】
また、経済性や実用性の観点からも、単一の基地局の海洋レーダから観測できる視線方向の流速のみから津波を予測できることが望ましいものである。
【０００６】
本発明は、上記事情を鑑みて成されたものであり、精度よく波の水位を予測することができる津波予測装置、方法、及びプログラムを提供することを目的とする。
【課題を解決するための手段】
【０００７】
上記目的を達成するために、本発明に係る津波予測装置は、各観測地点での波の視線方向の流速を入力として、予測対象地点での波の水位を予測する津波予測装置である。津波予測装置は、各観測地点での波の視線方向の流速の入力を受け付ける入力部を含む。津波予測装置は、前記予測対象地点での波の水位を含む状態を予測する予測部を含む。津波予測装置は、各観測地点での波の視線方向の流速の入力を受け付けた場合に、入力された各観測地点での波の視線方向の流速と、観測地点について予測された前記予測対象地点での水位を含む波の状態を、予め定められた観測行列を用いて変換することにより得られる各観測地点での波の視線方向の流速との差分に基づいて、前記予測対象地点での水位を含む波の状態を推定する推定部を含む。津波予測装置は、前記予測部による前記状態の予測と、前記推定部による前記状態の推定とを予め定められた条件を満たすまで繰り返させる判定部を含む。そして、前記予測部による前記状態の予測は、繰り返しの一つ前に前記推定部で推定された前記状態、又は繰り返しの一つ前に前記予測部で予測された前記状態に基づいて予測するように構成されている。
【０００８】
また、本発明に係る津波予測装置において、前記状態は、水位と線流量とを含むようにしてもよい。
【０００９】
また、本発明に係る津波予測装置において、前記状態は、水位と前記視線方向の線流量と前記視線方向と直交する方向の線流量とを含むようにしてもよい。
【００１０】
また、本発明に係る津波予測装置において、前記観測行列は、前記視線方向の流速と前記視線方向の線流量と静水深との関係を線形近似により求めたものとするようにしてもよい。
【００１１】
本発明に係るプログラムは、コンピュータを、津波予測装置の各部として機能させるためのプログラムである。
【発明の効果】
【００１２】
本発明の津波予測装置、方法、及びプログラムによれば、各観測地点での波の視線方向の流速の入力を受け付ける。予測対象地点での波の水位を含む状態を予測する。各観測地点での波の視線方向の流速の入力を受け付けた場合に、入力された各観測地点での波の視線方向の流速と、予測された予測対象地点での水位を含む波の状態を、予め定められた観測行列を用いて変換することにより得られる各観測地点での波の視線方向の流速との差分に基づいて、予測対象地点での水位を含む波の状態を推定する。状態の予測と、状態の推定とを予め定められた条件を満たすまで繰り返させる。状態の予測は、繰り返しの一つ前に推定された状態、又は繰り返しの一つ前に予測された状態に基づいて予測する。これにより、精度よく波の水位を予測することができる、という効果が得られる。
【図面の簡単な説明】
【００１３】
【図１】本発明の実施の形態に係る津波予測装置の構成を示すブロック図である。
【図２】本発明の実施の形態に係る津波予測装置における津波予測処理ルーチンを示すフローチャートである。
【図３】実験例における海洋レーダによる観測範囲の一例を示す図である。
【図４】実験例におけるデータ同化津波予測の計算領域を示す図である。
【図５】実験例における各波源の位置を示す図である。
【図６】実験例における予測の結果の一例である。
【発明を実施するための形態】
【００１４】
以下、図面を参照して本発明の実施の形態を詳細に説明する。
【００１５】
まず、本発明の実施の形態におけるデータ同化手法、及び予測の原理を説明する。
【００１６】
本実施の形態では、海洋レーダから得られる視線方向の津波の流速分布から水位分布を想定することなく、流速分布を用いて津波伝播のシミュレーションを行うことで、津波の到達時間、及び津波の波高（水位の高さ）を予測する。シミュレーションの際、気象予測などに用いられているデータ同化手法を用いることで解析値と観測値の親和性を高める。
【００１７】
次に、データ同化手法を用いた津波予測の方法について説明する。
【００１８】
データ同化とは、モデル（シミュレーション）に観測値を取り込み、より真値に近い結果を出す手法である。本実施の形態では、解析の計算負荷の観点から、背景場の誤差情報が時間的に変化しないものとする静的な同化手法である最適内挿法を用いることとする。最適内挿法では、最適な推定値ｘａを以下（１）式のように、シミュレーション結果である予報値ｘｂと、観測値ｙとの誤差に重みをかけたものとの和で与える。
【００１９】
【数１】
2020056649.tif 000003
・・・（１）
【００２０】
Ｈは観測行列、Ｗは重み行列である。観測行列Ｈは、計算格子点における予報値ｘｂから、観測点における値への変換行列であり、観測値ｙと予報値ｘｂとが同じ物理量の場合は空間内挿を表す。計算格子点は、観測範囲において所定の格子間隔で設定した格子点である。重み行列Ｗは、推定値ｘａの誤差分散が最小となるように設定する。推定誤差は、以下（２）式で与えられる。
【００２１】
【数２】
2020056649.tif 000004
・・・（２）
【００２２】
シミュレーションの誤差である背景誤差εｂと測定誤差εｏとは無相関とすると、推定誤差共分散は以下（３）式のようになる。
【００２３】
【数３】
2020056649.tif 000005
・・・（３）
【００２４】
ここで、ＢとＲとは以下（４）式、（５）式のように定義する。
【００２５】
【数４】
2020056649.tif 000006
・・・（４）
2020056649.tif 000007
・・・（５）
【００２６】
上記（３）式の推定誤差共分散の対角成分が推定値の誤差分散である。対角成分の和を重み行列Ｗで微分すると以下（６）式のようになる。
【００２７】
【数５】
2020056649.tif 000008
・・・（６）
ここで、ＢとＲとが対称行列であることを用いる。これを０とおくと、最適な重み行列Ｗは以下（７）式を満たす。
【００２８】
【数６】
2020056649.tif 000009
・・・（７）
観測点ｉ，ｊ間の背景誤差共分散行列ＨＢＨＴの成分をｂｉｊ、観測誤差共分散行列Ｒの成分をｒｉｊ、計算格子点ｇと観測点ｉとの間の背景誤差共分散行列ＢＨＴの成分をｂｇｉとすると、計算格子点ｇの対する観測点ｊの観測値の持つ重みｗｇｊは以下（８）式の連立一次方程式から求めることができる。
【００２９】
【数７】
2020056649.tif 000010
・・・（８）
【００３０】
σｉｂを観測点ｉにおける背景誤差の標準偏差、σｇｂを計算格子点ｇにおける背景誤差の標準偏差、σｉｏを観測点ｉにおける観測誤差の標準偏差として、両辺をσｉｂ及びσｇｂで割ると以下（９）式のように変形できる。
【００３１】
【数８】
2020056649.tif 000011
・・・（９）
【００３２】
μｉｊｂは、観測点ｉと観測点ｊとの背景誤差の相関係数であり
2020056649.tif 000012
と表せる。
【００３３】
μｉｊＯは、観測点ｉと観測点ｊとの観測誤差の相関係数であり
2020056649.tif 000013
と表せる。
【００３４】
μｇｉｂは、計算格子点ｇと観測点ｉとの背景誤差の相関係数であり
2020056649.tif 000014
と表せる。
【００３５】
さらに観測点間の観測誤差には相関がないものとして
2020056649.tif 000015
と仮定することで、以下（１０）式のように単純化できる。
【００３６】
【数９】
2020056649.tif 000016
・・・（１０）
【００３７】
以上のデータ同化手法を用いて、水位を含む波の状態を予測する。
【００３８】
次に予測の原理について説明する。特許文献１に記載されているように、津波の挙動について、ｘ軸及びｙ軸を有する２次元直交座標系の以下の質量保存式である（１１）式と運動方程式である（１２）式、及び（１３）式とから構成される長波理論の基礎方程式によって、津波の到達時間、及び津波の水位ηの状態を予測するシミュレーションが導出できる。
【００３９】
【数１０】
2020056649.tif 000017
・・・（１１）
2020056649.tif 000018
・・・（１２）
2020056649.tif 000019
・・・（１３）
【００４０】
（１１）〜（１３）式において、ηは津波の波高であり、Ｍはｘ軸方向の線流量であり、Ｎはｙ軸方向の線流量であり、ｎは海底摩擦係数である。Ｄは全水深であり、静水深ｈ及び波高ηを用いると、Ｄ＝ｈ＋ηとなる。ｔは時間であり、ｇは重力加速度である。
【００４１】
長波理論では、津波の流速は深さ方向（ｚ軸方向）に一定と仮定できるので、津波のｘ軸方向の流速Ｕ及びｙ軸方向の流速ＶはそれぞれＵ＝Ｍ／Ｄ、Ｖ＝Ｎ／Ｄとして算出される。すなわち、計測された海面のｘ軸方向の流速Ｕ及びｙ軸方向の流速Ｖは、ｘｙ平面上の座標により決定付けられる。従って、波高の推定では、ｘ軸方向の流速Ｕ及びｙ軸方向の流速Ｖと津波の波高ηとを関連付けるためのデータベースや経験式を必要とすることなしに、上述した津波の基礎方程式に基づいて、計測された津波のｘ軸方向の流速Ｕ及びｙ軸方向の流速Ｖから波高ηを算出できる。
【００４２】
以上の説明に基づいて、以下、津波予測装置の構成を説明する。
【００４３】
＜本発明の実施の形態に係る津波予測装置の構成＞
【００４４】
次に、本発明の実施の形態に係る津波予測装置の構成について説明する。図１に示すように、本発明の実施の形態に係る津波予測装置１００は、ＣＰＵと、ＲＡＭと、後述する津波予測処理ルーチンを実行するためのプログラムや各種データを記憶したＲＯＭと、を含むコンピュータで構成することが出来る。この津波予測装置１００は、機能的には図１に示すように入力部１０と、演算部２０と、出力部５０とを備えている。
【００４５】
入力部１０は、各観測地点ｉでの波の視線方向の流速ｕの観測値ｙを受け付ける。観測値ｙは、海洋レーダから随時に受け付けるものとする。
【００４６】
演算部２０は、重み算出部３０と、予測部３２と、推定部３４と、判定部３６とを含んで構成されている。
【００４７】
重み算出部３０は、観測値ｙを用いて求められた、背景誤差共分散行列ＨＢＨＴ、観測誤差共分散行列Ｒ、及び背景誤差共分散行列ＢＨＴに基づいて、上記（８）式に従って、重みｗｇｊからなる重み行列Ｗを算出する。重み行列Ｗは、推定値ｘａの誤差分散が最小となるような重み行列である。なお、観測値ｙ以外に、背景誤差共分散行列ＨＢＨＴ、観測誤差共分散行列Ｒ、及び背景誤差共分散行列ＢＨＴを求めるのに必要な値は、実験等により定めればよい。
【００４８】
ここで、観測値ｙの観測結果は流速ｕとして観測される。一方、津波解析のシミュレーションでは、水位ηと、視線方向の線流量Ｍと、視線方向と直交する方向の線流量Ｎとを解析に用いる必要があるため、視線方向の流速ｕから観測行列Ｈを用いて変換する必要がある。流速ｕは以下（１４）式で与えられる。
【００４９】
【数１１】
2020056649.tif 000020
・・・（１４）
【００５０】
Ｄは全水深、ｈは静水深を表す。ただし、上記（１４）式のままでは非線形であるため観測行列Ｈは作成できない。そのため、静水深に対し水位の変化は非常に小さいものとして、以下（１５）式の線形近似で観測行列Ｈを作成し、重み行列Ｗの算出に用いる。
【００５１】
【数１２】
2020056649.tif 000021
・・・（１５）
【００５２】
以下、推定部３４でも観測行列Ｈを用いる際には同様に線形近似したものを用いればよい。
【００５３】
予測部３２は、予測対象地点での波の水位ηと、視線方向の線流量Ｍと、視線方向と直交する方向の線流量Ｎとを含む波の状態の予報値ｘｎｂを予測する。予測対象地点が上述した計算格子点ｇに対応する。具体的には、予測部３２は、一時刻前（ｎ−１）の状態の推定値ｘｎ−１ａに基づき、上述した津波の基礎方程式から導出できるシミュレーションを行って、次の時刻ｎの状態の予報値ｘｎｂを予測する。推定部３４で一時刻前（ｎ−１）の状態の推定を行っていなければ、一時刻前（ｎ−１）の状態の予報値ｘｎ−１ｂを用いる。水位ηは、上記（１１）式の連続式に従って更新できる。線流量Ｍ、及び線流量Ｎは、上記（１２）式、及び（１３）式の運動方程式に従って更新できる。なお、ｎは時刻でなく回数としてもよい。
【００５４】
推定部３４は、各観測地点ｉでの波の視線方向の流速ｕの観測値ｙの入力を受け付けた場合に、上記（１）式に示すデータ同化手法により、重み行列Ｗを係数として、入力された各観測地点ｉでの波の視線方向の流速ｕの観測値ｙｎと、予測された予測対象地点での水位ηを含む波の状態の予報値ｘｎｂを、予め定められた観測行列Ｈを用いて変換することにより得られる各観測地点ｉでの波の視線方向の流速Ｈｘｎｂとの差分に基づいて、予測対象地点での水位ηを含む状態の推定値ｘｎａを推定する。
【００５５】
判定部３６は、予測部３２による波の状態の予報値ｘｎｂの予測と、推定部３４による状態の推定値ｘｎａの推定とを予め定められた条件を満たすまで繰り返させる。条件としては所定の時間や回数を定めておけばよい。予測部３２から繰り返しの度に予報値ｘｎｂを出力部５０に出力させてもよいし、繰り返しの終了後に出力させるようにしてもよい。
【００５６】
＜本発明の実施の形態に係る津波予測装置の作用＞
【００５７】
次に、本発明の実施の形態に係る津波予測装置１００の作用について説明する。津波予測装置１００は、入力部１０で観測値ｙを随時受け付けているときに、図２に示す津波予測処理ルーチンを実行する。例えば、観測値ｙを２分ごとに受け付けているものとする。
【００５８】
まず、ステップＳ１００では、観測値ｙを用いて求められた、背景誤差共分散行列ＨＢＨＴ、観測誤差共分散行列Ｒ、及び背景誤差共分散行列ＢＨＴに基づいて、上記（８）式に従って、重みｗｇｊからなる重み行列Ｗを算出する。
【００５９】
次に、ステップＳ１０２では、繰り返しの単位のカウントであるｎをｎ＝１と設定する。例えば、ｎは、時刻であり、計算時間間隔が１秒である場合には、１秒ごとに１カウントするものとする。
【００６０】
ステップＳ１０４では、予測対象地点での波の水位ηと、視線方向の線流量Ｍと、視線方向と直交する方向の線流量Ｎとを含む波の状態の予報値ｘｎｂを予測する。予測は、一時刻前（ｎ−１）の状態の推定値ｘｎ−１ａに基づき、上述した津波の基礎方程式から導出できるシミュレーションを行って、波の状態の予報値ｘｎｂを予測する。ステップＳ１０８で一時刻前（ｎ−１）の状態の推定を行っていなければ、一時刻前（ｎ−１）の状態の予報値ｘｎ−１ｂを用いる。
【００６１】
ステップＳ１０６では、各観測地点ｉでの波の視線方向の流速ｕの観測値ｙの入力を受け付けたかを判定し、受け付けていればステップＳ１０８へ移行し、受け付けていなければステップＳ１１０へ移行する。
【００６２】
ステップＳ１０８では、上記（１）式に示すデータ同化手法により、重み行列Ｗを係数として、入力された各観測地点ｉでの波の視線方向の流速ｕの観測値ｙｎと、観測地点ｉについて予測された予測対象地点での水位ηを含む波の状態の予報値Ｘｎｂを、予め定められた観測行列Ｈを用いて変換することにより得られる各観測地点ｉでの波の視線方向の流速Ｈｘｎｂとの差分に基づいて、予測対象地点での水位ηを含む状態の推定値ｘｎａを推定する。
【００６３】
ステップＳ１１０では、ｎ＝ｎｅｎｄであるかを判定する。ｎｅｎｄはｎについて予め定められた条件である。ｎｅｎｄであれば条件を満たしたものとして津波予測処理ルーチンを終了し、ｎｅｎｄでなければステップＳ１１２へ移行してｎ＝ｎ＋１とカウントアップしてステップＳ１０４〜Ｓ１１０の処理を繰り返す。
【００６４】
以上説明したように、本発明の実施の形態に係るによれば、各観測地点での波の視線方向の流速の入力を受け付け、予測対象地点での波の水位を含む状態を予測し、各観測地点での波の視線方向の流速の入力を受け付けた場合に、入力された各観測地点での波の視線方向の流速と、予測された予測対象地点での水位を含む波の状態を、予め定められた観測行列を用いて変換することにより得られる各観測地点での波の視線方向の流速との差分に基づいて、予測対象地点での水位を含む波の状態を推定し、状態の予測と、状態の推定とを予め定められた条件を満たすまで繰り返させ、状態の予測は、繰り返しの一つ前に推定された状態、又は繰り返しの一つ前に予測された状態に基づいて予測することにより、精度よく波の水位を予測することができる。
【００６５】
なお、本発明は、上述した実施の形態に限定されるものではなく、この発明の要旨を逸脱しない範囲内で様々な変形や応用が可能である。
【００６６】
例えば、上述した重み算出部３０で重み行列Ｗを一度算出して推定に用いる場合を例に説明したが、これに限定されるものではなく、判定部３６による時刻ｎの繰り返しの度に、背景誤差共分散行列ＨＢＨＴ、観測誤差共分散行列Ｒ、及び背景誤差共分散行列ＢＨＴを更新して、重み行列Ｗを算出するようにしてもよい。
【００６７】
［実験例］
本実施の形態の手法による津波予測の可能性を検討するため、実際に設置されている海洋レーダを対象とした実地形モデルにおいて、海洋レーダの流速観測値として予め実施した津波波源を配置した数値シミュレーションから得られる視線方向流速のみを与え、データ同化による津波予測の実験を行った。実験では、海底露出や陸上への遡上を考慮し、非線形長波理論を用いた検討を行った。データ同化津波予測の計算領域は図４に示すように、海洋レーダによる観測範囲近傍（東西４３．２ｋｍ、南北５７．６ｋｍ）に限定した。格子間隔は、２４０ｍ→８０ｍ→４０ｍ→２０ｍ→１０ｍ→５ｍと順次細分化し、１０ｍ以上の格子領域は陸上完全反射条件、５ｍ格子領域のみ遡上境界条件とした。
【００６８】
津波波源は敷地周辺海域活断層（長さ２９・５５・７２ｋｍ、幅２１・２６・２６ｋｍ、走向０・５５・３０°、上縁深さ２．５ｋｍ、傾斜角４５・３５・３５°、すべり角６２・９６・９０°、すべり量７．７ｍ、Ｍｗ８．０）を想定した。各波源の位置を図５に示す。
【００６９】
計算条件を表１に示す。
【表１】
2020056649.tif 000022
【００７０】
図６に実験における予測の結果を示す。実験では、敷地周辺海域の活断層について、途中までの観測値でデータ同化を行い、どの程度先の時刻まで予測可能かの検討を行った。検討では、背景誤差相関係数μｂは、８ｋｍで固定とした。使用する観測値を地震発生後２分、４分、６分までと変えた場合の水位時刻歴波形を検証している。流速の測定間隔を２分としているため、地震発生後２分までの観測値を使用した場合、データ同化は１回のみとなる。この場合でも、水位の予測は過小評価となるものの、津波の到達時間は概ね予測できている。地震発生後４分までの観測値を用いてデータ同化が２回になると第一波に関しては最後まで同化した結果とほぼ同等となっている。地震発生後６分までの観測値を用いてデータ同化が３回になると地震発生から３４分後まで概ね予測できている。このように、津波の到達時間については、データが入力されれば短時間で予測可能であり、複数回同化が進めば３０分程度先まで水位波形を精度よく予測可能となることが確認できた。
【００７１】
以上のように、データ同化手法を用いることにより短時間で精度よく波の水位を予測できることがわかる。
【符号の説明】
【００７２】
１０   入力部
２０   演算部
３０   重み算出部
３２   予測部
３４   推定部
３６   判定部
５０   出力部
１００ 津波予測装置
【図１】
2020056649.tif 000023
【図２】
2020056649.tif 000024
【図３】
2020056649.tif 000025
【図４】
2020056649.tif 000026
【図５】
2020056649.tif 000027
【図６】
2020056649.tif 000028