津波の検知方法、検知装置および検知プログラム
【課題】従来の津波の検知方法には、正確かつ保守性に優れたものが存在せず、新規な津波の検知方法が望まれていた。
【解決手段】陸地の複数の地点において磁場を測定し、その中の任意の一地点と他の地点の磁場の差信号を求め、該差信号につき、一定の過去から現時点までの平均値を求め、該平均値と現時点の該差信号の値を比較し、その差が所定の値を超過し、所定の時間継続した場合は、巨大津波が発生したと判断する津波の検知方法である。
【解決手段】陸地の複数の地点において磁場を測定し、その中の任意の一地点と他の地点の磁場の差信号を求め、該差信号につき、一定の過去から現時点までの平均値を求め、該平均値と現時点の該差信号の値を比較し、その差が所定の値を超過し、所定の時間継続した場合は、巨大津波が発生したと判断する津波の検知方法である。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は、陸地における磁場変動の測定に基づく津波の検知に関するものである。
【背景技術】
【０００２】
従来、津波の観測方法としては、潮位計や圧力センサーを用いて津波を検出し警報するシステムが知られていたが、例えば巨大な津波が押し寄せてきた場合には、潮位計の針が振り切れていたり、故障したりしていた（非特許文献１参照）。また、沖合に圧力センサーが設けられ、海底ケーブルで情報を地上に送信していたが、センサーが故障した場合には、保守が困難であった。
【０００３】
また、これらの方法ではセンサーや計測部の設置位置が海底や海面に制約されることから、高い設置費用や維持費用が発生すること、船舶との衝突による損傷や太陽光線による劣化が生じること、小さな波や突風による瞬間的な波も雑音として検出してしまい津波と区別する処理が煩雑になることなどの課題があった。
【０００４】
そこで、津波や高波の発生時に生ずるうねり性の波によって渦電流が海面に沿って流れ、渦電流と直角な方向に、正弦波状の周期的な磁界が生ずるのを、スカラ磁力計を海中ブイに収納して計測することにより、装置を海底や海面に設置することなく、錨等を用いて簡単に海に設置できるようにすると共に、津波の識別を簡単に行うことができる津波の検知方法が提案されている（特許文献１参照）。
しかしながら、この検知方法では未だ十分に、リアルタイムに正確な津波情報を得ることができず、津波被害を確実に防止するには至っていない。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９―１９８２０７号公報
【非特許文献】
【０００６】
【非特許文献１】Earth Planets Space、 63、 0-0、 2011 “Significant tsunami observed at the ocean-bottom pressure gauges at2011 Off the Pacific Coast of Tohoku Earthquake”
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明が解決しようとする課題は、巨大津波が押し寄せてきた場合においても針が振り切れたり、故障したりせず、また、保守点検が容易であり、正確に津波の発生を検知することである。
【課題を解決するための手段】
【０００８】
本発明は、基本的には、導電性である海水の巨大な塊の流体運動、すなわち電流により発生する磁場の微少な変動を計測することにより津波の発生を検知するための発明（下記のステップを含む津波を検知するための方法、津波を検知するためのプログラム、このプログラムを実行するための装置）である。
具体的には、本発明は、
津波を検知するために、記憶デバイスと演算処理デバイスとを備えた装置で使用されるプログラムであって、
陸地の測定地点において磁場を経時的に測定して磁場変動を把握するステップ、
現時点から所定時間遡って測定された磁場の平均値を求めるステップ、及び
該平均値と現時点の磁場の値とを比較し現時点の磁場の値と該平均値との差
を求め、当該差が前記磁場変動から算出された通常時の磁場変動と認められる所定の値を、所定の時間継続して超過した場合、津波が発生したと判断するステップを前記装置において実行させる実行可能命令を備えた、前記プログラムである。
前記測定地点は複数箇所である。前記磁場変動は、前記測定地点の中の任意の一地点と他の地点との磁場の差を求めて差信号を得、得られた差信号に基づいて把握される。前記平均値は、前記差信号について求められる。このように、陸地の複数の地点において磁場を測定するのは、地磁気の乱れ等の影響を除くためである。
また、前記プログラムを含む記憶デバイスと演算処理デバイスとを備えた装置は、磁場を測定するための磁場測定器と、前記演算処理デバイスにより演算された結果を表示する出力デバイスと、津波の発生を検知した際に津波の発生を表示する表示デバイスとを具備する。また、この装置は、津波の発生を検知したことを受けて警報を鳴らすことや、ニュース速報をメディアに流すようなシステムを更に具備することもできる。
なお、上述のプログラムを用いた津波検知のフローチャートを図８に示した。
【０００９】
また、上記磁場は、互いに直交する磁場の３成分のうちの1成分、２成分又は３成分に基づいて求めたものである。
【００１０】
また、上記所定の値は、一定の過去から現時点までの前記差の標準偏差（σ）の定数倍（ｋσ：ｋは定数（例えば２から４））である。
【発明の効果】
【００１１】
本発明によれば、陸地において津波の検出ができ、観測装置のメンテナスが簡単になる。また、磁場変化を測定することで津波の発生をリアルタイムで観測できるため、即時に正確な津波予測が可能となる。
【図面の簡単な説明】
【００１２】
【図１】磁場変動状況の観測点を示す概略図である。
【図２】各地点における観測結果を示す波形チャート図である。
【図３】津波発生時の細倉観測点における磁場変動観測例（全磁力Ｆ）である。具体的には、細倉における全磁力（3成分の二乗和の平方根）の変化のチャートを示している。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σ（σは、標準偏差）である。
【図４】津波発生時の細倉における垂直磁力（Ｚ成分）の変化を示すチャートである。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σである。
【図５】津波発生時の細倉と女満別の全磁力の差分を示すチャートである。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σである。
【図６】津波発生時の細倉と女満別の垂直磁力（Ｚ成分）の差分を示すチャートである。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σである。
【図７】津波発生時の細倉と女満別の全磁力差分信号変化に津波の検知と相馬の水位計の相関図である。
【図８】上述のプログラムを用いた津波検知のフローチャートである。
【発明を実施するための形態】
【００１３】
本発明者は、図１に示すように宮城県の細倉（HSK）に観測地点を設け、地球磁場を継続的に観測している。
細倉観測地点では、宮城県栗原市にある細倉鉱山において、2004年3月以降、１対のフラックスゲート磁力計を用いた地磁気３成分観測を行っている。１台は地下70mの坑道内に設置され、もう１台は地表下１ｍに設置されている。地中磁力計は0.5秒サンプリング、地上磁力計は１秒サンプリングで記録し、どちらの観測データもGPS電波によって時刻同期されている。
【００１４】
また、それぞれの記録システムは、太陽光発電やバッテリーを電力の供給源としており、商用電力系とは独立しているので地震に伴う停電等の影響はほとんどない。
【００１５】
2011年3月11日14時46分18秒（日本時間）に2011年東北地方太平洋沖地震（モーメントマグニチュード（Mw）が9.0）が発生した。細倉観測地点と震源域までの距離は約100kmである。
【００１６】
図２に、各観測地点での地震発生時の地球磁場信号を示す。茨城県石岡市の柿岡（KAK）と北海道女満別（MMB）の地磁気データは、気象庁地磁気観測所の提供によるもので、東京都の奥多摩（OKM）は、東海大学地震予知研究センターより提供されたものである。地震発生時刻は、図のＴ０で示してある。地震発生から数分後には、各観測地点に地震が到達しているはずだが、磁場にさほど大きな変動は生じていない。しかしながら、細倉観測点の磁場信号を参照すると、地震発生から10分ないし18分にかけて、巨大地震、それに引き続いて発生した津波による注目すべき磁場変化が見られた。
【００１７】
地磁気の測定には、３成分観測を行っている。したがって、３成分のうちの１成分のみ（例えば、垂直磁力成分：Ｚ成分）または２成分（例えば、垂直成分Ｚと水平成分のうちの一つの成分）に注目しても良いし、３成分すべてに注目しても良い。
【００１８】
観測されたデータを観測しているだけでも、磁場変動の異常を観測可能であるが、観測時点から所定時間前からの磁場の平均値（Ｍ）を求め、その平均値と現時点での磁場を比較することにより、より客観的に津波を検知することができる。
【００１９】
図３は、巨大地震発生日の細倉観測点における全磁力（３成分の二乗和の平方根）のデータである。その図において、実線は、全磁力の観測波形であり、破線は、２０分間の移動平均値（Ｍで表す。）であり、二点鎖線は、Ｍ±３σである。図３より、一観測所のみの信号において津波と判断する閾値として、３σは大きすぎ、２σ程度とすべきであることがわかる。
【００２０】
図４は、巨大地震発生日の細倉観測点における垂直磁力強度（Ｚ成分）のデータである。
この場合も、全磁力と同様に、一観測所のみの信号において津波と判断する閾値として、３σは大きすぎ、２σ程度とすべきであることがわかる。
【００２１】
以上のように、一観測所のデータだけからでも、津波の発生を検知できるが、さらに良い方法がある。すなわち、複数の観測点に共通な、津波以外の原因による雑音、例えば、地磁気の乱れを消去するために、津波に近い観測点と津波から遠い観測点において観測された磁場変化の差分信号を利用することで、津波に基づく磁場変化を感度よく検出することが可能となった。
【００２２】
図５は、細倉と女満別の観測データの全磁力の差分信号のグラフである。実線は観測波形であり、破線は、３０分間の移動平均値、二点鎖線は、Ｍ±３σである。普段は、平均値からほとんどずれずに差分信号が推移するが、巨大津波発生時には、大きくずれることが確認された。これを検出することで、マルチ観測点における磁場変化の差信号を利用した津波検知が可能である。
【００２３】
図６は、細倉と女満別の観測データの垂直磁場Ｚの差分信号のグラフである。実線は観測波形であり、破線は、３０分間の移動平均値、二点鎖線は、Ｍ±３σである。図５と同じく、普段は、平均値からほとんどずれずに差分信号が推移するが、巨大津波発生時には、大きくずれることが確認された。このように、３成分のうちの１成分であるＺ成分を検出することで、マルチ観測点における磁場変化の差信号を利用した津波検知が可能である。
【００２４】
これを数式で現すと、時刻ｔにおける上記誤差信号をＤ（ｔ）とすると、所定の時刻から現時点までの誤差信号の平均値Ｍは、
Ｍ＝｛Ｄ（ｔ−ｎＴ）＋Ｄ（ｔ−（ｎ−１）Ｔ）＋・・・＋Ｄ（ｔ−Ｔ）｝／ｎ
となる（サンプリング間隔をＴとする。）。また、標準偏差σは、
σ２＝｛（Ｄ（ｔ−ｎＴ）―Ｍ）２＋・・・＋（Ｄ（ｔ−Ｔ）―Ｍ）２｝／ｎ
から求められる。
したがって、本発明は、
絶対値（Ｄ（ｔ）―Ｍ）＞３σ
または
絶対値（Ｄ（ｔ）―Ｍ）＞ｋσ（ｋは、任意の定数（例えば２から４））
となった時（ｔ）に、津波が発生したと判断する。
なお、上述のｎは、10分から30分となる離散時間となる。換言すれば、サンプリング間隔Ｔが1秒ならｎ＝600から1800となり、Ｔが2秒ならｎ＝300から900となる。
【００２５】
サンプリング間隔Ｔは、０．５秒から２秒程度であればよい。比較判断する大きさは、平常時の標準偏差が１ｎＴ程度であるので、具体的に、例えば、２ｎＴとしても良いが、標準偏差を規準にして定めた方が好ましい。上記の例では、３σとした。
【００２６】
さらに、津波と検知するときは、ｋσと決めた基準の値（閾値）を所定の時間（例えば３０秒程度）連続して超えた場合とするのが好ましい。なぜなら、地磁気観測を行っている場合に、非常に鋭いピークを持ったパルスが観測されることがしばしばあるからである。地磁気観測では有名な話であり、原因は定かではないが、車とか人間とか鳥とか言われている。
【００２７】
図７は、細倉と女満別の観測データの全磁力の差分信号のグラフと相馬における水位計のグラフを時刻同期して表示したものである。相馬は、震源に比較的近く、水位が１０ｍにも達している。図７より、津波は、１５時頃に発生し、１６時前に相馬に到達したことがわかる。したがって、津波到達の５０分前には津波警報を発することが可能である。
【００２８】
以上から明らかなように、マルチ観測点における磁場変化の差信号を利用して津波検知が可能となる。しかも、観測点間で所定の差以上の信号が出た地域の近くで津波が発生したことが分かる。このように本発明は、磁場の差信号を利用する点、リアルタイムの観測が可能である。
【産業上の利用可能性】
【００２９】
本発明者は、細倉鉱山に作られた観測点で2011年3月11日に発生した東北地方太平洋沖地震（東日本大震災）において、巨大津波に伴う磁場変動を観測することに成功した。これは地震電磁気観測にとって画期的なことであり、この成果をもとに、例えば巨大な津波への緊急警報が可能となるため、本分野の飛躍的な発展が大いに期待できる。
【符号の説明】
【００３０】
MMB 女満別（北海道）
HSK 細倉（宮城県）
KAK 柿岡（茨城県）
OKM 奥多摩（東京都）
ESS 江刺（岩手県）
Aftershock Area 余震の領域
Tohoku EQ 東北地方地震の震源地
Tsunami 津波
【技術分野】
【０００１】
本発明は、陸地における磁場変動の測定に基づく津波の検知に関するものである。
【背景技術】
【０００２】
従来、津波の観測方法としては、潮位計や圧力センサーを用いて津波を検出し警報するシステムが知られていたが、例えば巨大な津波が押し寄せてきた場合には、潮位計の針が振り切れていたり、故障したりしていた（非特許文献１参照）。また、沖合に圧力センサーが設けられ、海底ケーブルで情報を地上に送信していたが、センサーが故障した場合には、保守が困難であった。
【０００３】
また、これらの方法ではセンサーや計測部の設置位置が海底や海面に制約されることから、高い設置費用や維持費用が発生すること、船舶との衝突による損傷や太陽光線による劣化が生じること、小さな波や突風による瞬間的な波も雑音として検出してしまい津波と区別する処理が煩雑になることなどの課題があった。
【０００４】
そこで、津波や高波の発生時に生ずるうねり性の波によって渦電流が海面に沿って流れ、渦電流と直角な方向に、正弦波状の周期的な磁界が生ずるのを、スカラ磁力計を海中ブイに収納して計測することにより、装置を海底や海面に設置することなく、錨等を用いて簡単に海に設置できるようにすると共に、津波の識別を簡単に行うことができる津波の検知方法が提案されている（特許文献１参照）。
しかしながら、この検知方法では未だ十分に、リアルタイムに正確な津波情報を得ることができず、津波被害を確実に防止するには至っていない。
【先行技術文献】
【特許文献】
【０００５】
【特許文献１】特開２００９―１９８２０７号公報
【非特許文献】
【０００６】
【非特許文献１】Earth Planets Space、 63、 0-0、 2011 “Significant tsunami observed at the ocean-bottom pressure gauges at2011 Off the Pacific Coast of Tohoku Earthquake”
【発明の概要】
【発明が解決しようとする課題】
【０００７】
本発明が解決しようとする課題は、巨大津波が押し寄せてきた場合においても針が振り切れたり、故障したりせず、また、保守点検が容易であり、正確に津波の発生を検知することである。
【課題を解決するための手段】
【０００８】
本発明は、基本的には、導電性である海水の巨大な塊の流体運動、すなわち電流により発生する磁場の微少な変動を計測することにより津波の発生を検知するための発明（下記のステップを含む津波を検知するための方法、津波を検知するためのプログラム、このプログラムを実行するための装置）である。
具体的には、本発明は、
津波を検知するために、記憶デバイスと演算処理デバイスとを備えた装置で使用されるプログラムであって、
陸地の測定地点において磁場を経時的に測定して磁場変動を把握するステップ、
現時点から所定時間遡って測定された磁場の平均値を求めるステップ、及び
該平均値と現時点の磁場の値とを比較し現時点の磁場の値と該平均値との差
を求め、当該差が前記磁場変動から算出された通常時の磁場変動と認められる所定の値を、所定の時間継続して超過した場合、津波が発生したと判断するステップを前記装置において実行させる実行可能命令を備えた、前記プログラムである。
前記測定地点は複数箇所である。前記磁場変動は、前記測定地点の中の任意の一地点と他の地点との磁場の差を求めて差信号を得、得られた差信号に基づいて把握される。前記平均値は、前記差信号について求められる。このように、陸地の複数の地点において磁場を測定するのは、地磁気の乱れ等の影響を除くためである。
また、前記プログラムを含む記憶デバイスと演算処理デバイスとを備えた装置は、磁場を測定するための磁場測定器と、前記演算処理デバイスにより演算された結果を表示する出力デバイスと、津波の発生を検知した際に津波の発生を表示する表示デバイスとを具備する。また、この装置は、津波の発生を検知したことを受けて警報を鳴らすことや、ニュース速報をメディアに流すようなシステムを更に具備することもできる。
なお、上述のプログラムを用いた津波検知のフローチャートを図８に示した。
【０００９】
また、上記磁場は、互いに直交する磁場の３成分のうちの1成分、２成分又は３成分に基づいて求めたものである。
【００１０】
また、上記所定の値は、一定の過去から現時点までの前記差の標準偏差（σ）の定数倍（ｋσ：ｋは定数（例えば２から４））である。
【発明の効果】
【００１１】
本発明によれば、陸地において津波の検出ができ、観測装置のメンテナスが簡単になる。また、磁場変化を測定することで津波の発生をリアルタイムで観測できるため、即時に正確な津波予測が可能となる。
【図面の簡単な説明】
【００１２】
【図１】磁場変動状況の観測点を示す概略図である。
【図２】各地点における観測結果を示す波形チャート図である。
【図３】津波発生時の細倉観測点における磁場変動観測例（全磁力Ｆ）である。具体的には、細倉における全磁力（3成分の二乗和の平方根）の変化のチャートを示している。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σ（σは、標準偏差）である。
【図４】津波発生時の細倉における垂直磁力（Ｚ成分）の変化を示すチャートである。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σである。
【図５】津波発生時の細倉と女満別の全磁力の差分を示すチャートである。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σである。
【図６】津波発生時の細倉と女満別の垂直磁力（Ｚ成分）の差分を示すチャートである。観測日時は、地震（津波）が発生した２０１１年３月１１日の１４時〜１６時の２時間である。実線は観測波形、破線は移動平均値（Ｍ）、二点鎖線はＭ±３σである。
【図７】津波発生時の細倉と女満別の全磁力差分信号変化に津波の検知と相馬の水位計の相関図である。
【図８】上述のプログラムを用いた津波検知のフローチャートである。
【発明を実施するための形態】
【００１３】
本発明者は、図１に示すように宮城県の細倉（HSK）に観測地点を設け、地球磁場を継続的に観測している。
細倉観測地点では、宮城県栗原市にある細倉鉱山において、2004年3月以降、１対のフラックスゲート磁力計を用いた地磁気３成分観測を行っている。１台は地下70mの坑道内に設置され、もう１台は地表下１ｍに設置されている。地中磁力計は0.5秒サンプリング、地上磁力計は１秒サンプリングで記録し、どちらの観測データもGPS電波によって時刻同期されている。
【００１４】
また、それぞれの記録システムは、太陽光発電やバッテリーを電力の供給源としており、商用電力系とは独立しているので地震に伴う停電等の影響はほとんどない。
【００１５】
2011年3月11日14時46分18秒（日本時間）に2011年東北地方太平洋沖地震（モーメントマグニチュード（Mw）が9.0）が発生した。細倉観測地点と震源域までの距離は約100kmである。
【００１６】
図２に、各観測地点での地震発生時の地球磁場信号を示す。茨城県石岡市の柿岡（KAK）と北海道女満別（MMB）の地磁気データは、気象庁地磁気観測所の提供によるもので、東京都の奥多摩（OKM）は、東海大学地震予知研究センターより提供されたものである。地震発生時刻は、図のＴ０で示してある。地震発生から数分後には、各観測地点に地震が到達しているはずだが、磁場にさほど大きな変動は生じていない。しかしながら、細倉観測点の磁場信号を参照すると、地震発生から10分ないし18分にかけて、巨大地震、それに引き続いて発生した津波による注目すべき磁場変化が見られた。
【００１７】
地磁気の測定には、３成分観測を行っている。したがって、３成分のうちの１成分のみ（例えば、垂直磁力成分：Ｚ成分）または２成分（例えば、垂直成分Ｚと水平成分のうちの一つの成分）に注目しても良いし、３成分すべてに注目しても良い。
【００１８】
観測されたデータを観測しているだけでも、磁場変動の異常を観測可能であるが、観測時点から所定時間前からの磁場の平均値（Ｍ）を求め、その平均値と現時点での磁場を比較することにより、より客観的に津波を検知することができる。
【００１９】
図３は、巨大地震発生日の細倉観測点における全磁力（３成分の二乗和の平方根）のデータである。その図において、実線は、全磁力の観測波形であり、破線は、２０分間の移動平均値（Ｍで表す。）であり、二点鎖線は、Ｍ±３σである。図３より、一観測所のみの信号において津波と判断する閾値として、３σは大きすぎ、２σ程度とすべきであることがわかる。
【００２０】
図４は、巨大地震発生日の細倉観測点における垂直磁力強度（Ｚ成分）のデータである。
この場合も、全磁力と同様に、一観測所のみの信号において津波と判断する閾値として、３σは大きすぎ、２σ程度とすべきであることがわかる。
【００２１】
以上のように、一観測所のデータだけからでも、津波の発生を検知できるが、さらに良い方法がある。すなわち、複数の観測点に共通な、津波以外の原因による雑音、例えば、地磁気の乱れを消去するために、津波に近い観測点と津波から遠い観測点において観測された磁場変化の差分信号を利用することで、津波に基づく磁場変化を感度よく検出することが可能となった。
【００２２】
図５は、細倉と女満別の観測データの全磁力の差分信号のグラフである。実線は観測波形であり、破線は、３０分間の移動平均値、二点鎖線は、Ｍ±３σである。普段は、平均値からほとんどずれずに差分信号が推移するが、巨大津波発生時には、大きくずれることが確認された。これを検出することで、マルチ観測点における磁場変化の差信号を利用した津波検知が可能である。
【００２３】
図６は、細倉と女満別の観測データの垂直磁場Ｚの差分信号のグラフである。実線は観測波形であり、破線は、３０分間の移動平均値、二点鎖線は、Ｍ±３σである。図５と同じく、普段は、平均値からほとんどずれずに差分信号が推移するが、巨大津波発生時には、大きくずれることが確認された。このように、３成分のうちの１成分であるＺ成分を検出することで、マルチ観測点における磁場変化の差信号を利用した津波検知が可能である。
【００２４】
これを数式で現すと、時刻ｔにおける上記誤差信号をＤ（ｔ）とすると、所定の時刻から現時点までの誤差信号の平均値Ｍは、
Ｍ＝｛Ｄ（ｔ−ｎＴ）＋Ｄ（ｔ−（ｎ−１）Ｔ）＋・・・＋Ｄ（ｔ−Ｔ）｝／ｎ
となる（サンプリング間隔をＴとする。）。また、標準偏差σは、
σ２＝｛（Ｄ（ｔ−ｎＴ）―Ｍ）２＋・・・＋（Ｄ（ｔ−Ｔ）―Ｍ）２｝／ｎ
から求められる。
したがって、本発明は、
絶対値（Ｄ（ｔ）―Ｍ）＞３σ
または
絶対値（Ｄ（ｔ）―Ｍ）＞ｋσ（ｋは、任意の定数（例えば２から４））
となった時（ｔ）に、津波が発生したと判断する。
なお、上述のｎは、10分から30分となる離散時間となる。換言すれば、サンプリング間隔Ｔが1秒ならｎ＝600から1800となり、Ｔが2秒ならｎ＝300から900となる。
【００２５】
サンプリング間隔Ｔは、０．５秒から２秒程度であればよい。比較判断する大きさは、平常時の標準偏差が１ｎＴ程度であるので、具体的に、例えば、２ｎＴとしても良いが、標準偏差を規準にして定めた方が好ましい。上記の例では、３σとした。
【００２６】
さらに、津波と検知するときは、ｋσと決めた基準の値（閾値）を所定の時間（例えば３０秒程度）連続して超えた場合とするのが好ましい。なぜなら、地磁気観測を行っている場合に、非常に鋭いピークを持ったパルスが観測されることがしばしばあるからである。地磁気観測では有名な話であり、原因は定かではないが、車とか人間とか鳥とか言われている。
【００２７】
図７は、細倉と女満別の観測データの全磁力の差分信号のグラフと相馬における水位計のグラフを時刻同期して表示したものである。相馬は、震源に比較的近く、水位が１０ｍにも達している。図７より、津波は、１５時頃に発生し、１６時前に相馬に到達したことがわかる。したがって、津波到達の５０分前には津波警報を発することが可能である。
【００２８】
以上から明らかなように、マルチ観測点における磁場変化の差信号を利用して津波検知が可能となる。しかも、観測点間で所定の差以上の信号が出た地域の近くで津波が発生したことが分かる。このように本発明は、磁場の差信号を利用する点、リアルタイムの観測が可能である。
【産業上の利用可能性】
【００２９】
本発明者は、細倉鉱山に作られた観測点で2011年3月11日に発生した東北地方太平洋沖地震（東日本大震災）において、巨大津波に伴う磁場変動を観測することに成功した。これは地震電磁気観測にとって画期的なことであり、この成果をもとに、例えば巨大な津波への緊急警報が可能となるため、本分野の飛躍的な発展が大いに期待できる。
【符号の説明】
【００３０】
MMB 女満別（北海道）
HSK 細倉（宮城県）
KAK 柿岡（茨城県）
OKM 奥多摩（東京都）
ESS 江刺（岩手県）
Aftershock Area 余震の領域
Tohoku EQ 東北地方地震の震源地
Tsunami 津波
【特許請求の範囲】
【請求項１】
津波を、当該津波が陸地に到達する前に検知するための方法であって、
陸地の測定地点において磁場を経時的に測定して磁場変動を把握するステップ、
現時点から所定時間遡って測定された磁場の平均値を求めるステップ、及び
該平均値と現時点の磁場の値とを比較し現時点の磁場の値と該平均値との差を求め、当該差が前記磁場変動から算出された通常時の磁場変動と認められる所定の値を、所定の時間継続して超過した場合、津波が発生したと判断するステップ
を備え、
前記測定地点が複数箇所であり、
前記磁場変動は、前記測定地点の中の任意の一地点と他の地点との磁場の差を求めて差信号を得、得られた差信号に基づいて把握され、
前記平均値は前記差信号について求められる、
前記方法。
【請求項２】
津波を検知するために、記憶デバイスと演算処理デバイスとを備えた装置で使用されるプログラムであって、
陸地の測定地点において磁場を経時的に測定して磁場変動を把握するステップ、
現時点から所定時間遡って測定された磁場の平均値を求めるステップ、及び
該平均値と現時点の磁場の値とを比較し現時点の磁場の値と該平均値との差を求め、当該差が前記磁場変動から算出された通常時の磁場変動と認められる所定の値を、所定の時間継続して超過した場合、津波が発生したと判断するステップ
を前記装置において実行させる実行可能命令を備え、
前記測定地点が複数箇所であり、
前記磁場変動は、前記測定地点の中の任意の一地点と他の地点との磁場の差を求めて差信号を得、得られた差信号に基づいて把握され、
前記平均値は前記差信号について求められる、
前記プログラム。
【請求項３】
記憶デバイスと演算処理デバイスとを備え、請求項２に記載のプログラムを実行するための装置あって、
磁場を測定するための磁場測定器と、
前記演算処理デバイスにより演算された結果を表示する出力デバイスと、
津波の発生を検知した際に津波の発生を表示する表示デバイスと
を具備する、津波を検知するための装置。
【請求項１】
津波を、当該津波が陸地に到達する前に検知するための方法であって、
陸地の測定地点において磁場を経時的に測定して磁場変動を把握するステップ、
現時点から所定時間遡って測定された磁場の平均値を求めるステップ、及び
該平均値と現時点の磁場の値とを比較し現時点の磁場の値と該平均値との差を求め、当該差が前記磁場変動から算出された通常時の磁場変動と認められる所定の値を、所定の時間継続して超過した場合、津波が発生したと判断するステップ
を備え、
前記測定地点が複数箇所であり、
前記磁場変動は、前記測定地点の中の任意の一地点と他の地点との磁場の差を求めて差信号を得、得られた差信号に基づいて把握され、
前記平均値は前記差信号について求められる、
前記方法。
【請求項２】
津波を検知するために、記憶デバイスと演算処理デバイスとを備えた装置で使用されるプログラムであって、
陸地の測定地点において磁場を経時的に測定して磁場変動を把握するステップ、
現時点から所定時間遡って測定された磁場の平均値を求めるステップ、及び
該平均値と現時点の磁場の値とを比較し現時点の磁場の値と該平均値との差を求め、当該差が前記磁場変動から算出された通常時の磁場変動と認められる所定の値を、所定の時間継続して超過した場合、津波が発生したと判断するステップ
を前記装置において実行させる実行可能命令を備え、
前記測定地点が複数箇所であり、
前記磁場変動は、前記測定地点の中の任意の一地点と他の地点との磁場の差を求めて差信号を得、得られた差信号に基づいて把握され、
前記平均値は前記差信号について求められる、
前記プログラム。
【請求項３】
記憶デバイスと演算処理デバイスとを備え、請求項２に記載のプログラムを実行するための装置あって、
磁場を測定するための磁場測定器と、
前記演算処理デバイスにより演算された結果を表示する出力デバイスと、
津波の発生を検知した際に津波の発生を表示する表示デバイスと
を具備する、津波を検知するための装置。
【図１】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【図２】
【図３】
【図４】
【図５】
【図６】
【図７】
【図８】
【公開番号】特開２０１３−８３６１８（Ｐ２０１３−８３６１８Ａ）
【公開日】平成２５年５月９日（２０１３．５．９）
【国際特許分類】
【出願番号】特願２０１１−２４１６８８（Ｐ２０１１−２４１６８８）
【出願日】平成２３年１１月２日（２０１１．１１．２）
【出願人】（３０５０２７４０１）公立大学法人首都大学東京 (385)
【公開日】平成２５年５月９日（２０１３．５．９）
【国際特許分類】
【出願日】平成２３年１１月２日（２０１１．１１．２）
【出願人】（３０５０２７４０１）公立大学法人首都大学東京 (385)
[ Back to top ]