16:10 〜 16:20
[S17P-09] Kalmanフィルタを用いた津波波高即時把握手法の開発:海底水圧計・加速度計記録の同時活用
日本列島太平洋沖に近年整備されたDONETやS-netといった観測網によって、断層近傍における津波観測が可能になった。一方で、津波発生時に発生する水圧変化には水深変化や海底加速度変化などの非津波成分が重畳することが指摘されており(e.g., Saito and Tsushima, 2016)、水圧計による津波即時把握には課題が残されている。
本研究では、DONET・S-netの各観測点に設置されている水圧計・加速度計の記録をKalmanフィルタを介して組み合わせることで海底変位の逐次推定手法を提案する。これまでの研究で開発した、地震時水圧記録から津波・変位成分のみを抽出する手法(Mizutani et al., in press)により得られた波形に対して、そこに含まれる津波成分をノイズ、変位成分をシグナルとみなしてKalmanフィルタを適用した。また、推定した海底変位とMizutani et al.の手法で抽出した水圧波形を比較することで、津波波高の即時推定を試みた。
カルマンフィルタを用いた地震時変位の推定は、陸上観測点においてGNSS観測と加速度計を組み合わせる形でおこなわれており(e.g., Bock et al., 2011)、フィルタの設計はこれらの先行研究の結果を参考におこなった。GNSS観測と水圧計記録の相違点として、GNSS観測の場合はノイズが定常的に加わっていると考えられるのに対して、水圧記録に含まれるノイズである津波が観測されるのは有限の時間のみであることがあげられる。そのため本研究では、ノイズの非定常性を考慮した適応Kalmanフィルタ(e.g., 佐藤・竹井,1998; Hu et al., 2003)を用いた。さらに、結果の発散を抑制するために、通常のKalmanフィルタを平方根形式に変形した、平方根Kalmanフィルタを用いた(e.g., Haykin, 1996)。
本手法を2016年4月1日に三重県南東沖で発生した地震(Mw 6.0)によるDONET記録に適用した結果、図1のように変位成分(オフセット成分)を低減させることに成功した。しかし、このような結果を得るためには、適応Kalmanフィルタの調節パラメータである忘却係数をあらかじめ適切に設定する必要がある。現段階では、オフセットを十分低減させるためには忘却係数の設定を事後的におこなう必要があるため、実用的な津波即時把握のためには、最適な忘却係数の決定を即時把握に適用できるように改善する必要がある。
本研究では、DONET・S-netの各観測点に設置されている水圧計・加速度計の記録をKalmanフィルタを介して組み合わせることで海底変位の逐次推定手法を提案する。これまでの研究で開発した、地震時水圧記録から津波・変位成分のみを抽出する手法(Mizutani et al., in press)により得られた波形に対して、そこに含まれる津波成分をノイズ、変位成分をシグナルとみなしてKalmanフィルタを適用した。また、推定した海底変位とMizutani et al.の手法で抽出した水圧波形を比較することで、津波波高の即時推定を試みた。
カルマンフィルタを用いた地震時変位の推定は、陸上観測点においてGNSS観測と加速度計を組み合わせる形でおこなわれており(e.g., Bock et al., 2011)、フィルタの設計はこれらの先行研究の結果を参考におこなった。GNSS観測と水圧計記録の相違点として、GNSS観測の場合はノイズが定常的に加わっていると考えられるのに対して、水圧記録に含まれるノイズである津波が観測されるのは有限の時間のみであることがあげられる。そのため本研究では、ノイズの非定常性を考慮した適応Kalmanフィルタ(e.g., 佐藤・竹井,1998; Hu et al., 2003)を用いた。さらに、結果の発散を抑制するために、通常のKalmanフィルタを平方根形式に変形した、平方根Kalmanフィルタを用いた(e.g., Haykin, 1996)。
本手法を2016年4月1日に三重県南東沖で発生した地震(Mw 6.0)によるDONET記録に適用した結果、図1のように変位成分(オフセット成分)を低減させることに成功した。しかし、このような結果を得るためには、適応Kalmanフィルタの調節パラメータである忘却係数をあらかじめ適切に設定する必要がある。現段階では、オフセットを十分低減させるためには忘却係数の設定を事後的におこなう必要があるため、実用的な津波即時把握のためには、最適な忘却係数の決定を即時把握に適用できるように改善する必要がある。