特許 6206980 圧力センサの出力周波数算出方法およびそれを用いた気圧観測による津波警報装置、津波警報システム

(19)【発行国】日本国特許庁(JP)

(12)【公報種別】特許公報(B2)

(11)【特許番号】6206980

(24)【登録日】2017年9月15日

(45)【発行日】2017年10月4日

(54)【発明の名称】圧力センサの出力周波数算出方法およびそれを用いた気圧観測による津波警報装置、津波警報システム

(51)【国際特許分類】

   G01L 7/00 20060101AFI20170925BHJP

   G01V 1/00 20060101ALI20170925BHJP

   G01L 9/00 20060101ALI20170925BHJP

【ＦＩ】

   G01L7/00 P

   G01V1/00 D

   G01L9/00 C

   G01L9/00 E

【請求項の数】3

【全頁数】10

(21)【出願番号】特願2015-3183(P2015-3183)

(22)【出願日】2015年1月9日

(65)【公開番号】特開2016-128766(P2016-128766A)

(43)【公開日】2016年7月14日

【審査請求日】2017年6月2日

【早期審査対象出願】

(73)【特許権者】

【識別番号】501138231

【氏名又は名称】国立研究開発法人防災科学技術研究所

(74)【代理人】

【識別番号】100139114

【弁理士】

【氏名又は名称】田中 貞嗣

(74)【代理人】

【識別番号】100139103

【弁理士】

【氏名又は名称】小山 卓志

(74)【代理人】

【識別番号】100088041

【弁理士】

【氏名又は名称】阿部 龍吉

(74)【代理人】

【識別番号】100145920

【弁理士】

【氏名又は名称】森川 聡

(72)【発明者】

【氏名】▲功▼刀 卓

【審査官】 森 雅之

(56)【参考文献】

【文献】 米国特許出願公開第２００７／０１１２７０５（ＵＳ，Ａ１）

【文献】 特開２０１３−９６８０２（ＪＰ，Ａ）

【文献】 特許第２９７８１４０（ＪＰ，Ｂ２）

【文献】 欧州特許出願公開第０７０６０３８（ＥＰ，Ａ１）

【文献】 特許第２５８６８３８（ＪＰ，Ｂ２）

【文献】 特許第２８７６５２８（ＪＰ，Ｂ２）

【文献】 特許第２８７９６７７（ＪＰ，Ｂ２）

(58)【調査した分野】（Int.Cl.，ＤＢ名）

Ｇ０１Ｌ

Ｇ０１Ｖ

(57)【特許請求の範囲】

【請求項1】

圧力センサから出力される矩形波の周波数を算出する圧力センサの出力周波数算出方法であって、
時間フレーム中の矩形波の数をカウントする際には、
時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅を、開始時矩形波の時間幅で除したものによりカウントし、
時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅を、終了時矩形波の時間幅で除したものによりカウントすることを特徴とする圧力センサの出力周波数算出方法。

【請求項2】

地震センサと、前記地震センサからの出力をＡＤ変換することで地震データを生成するＡＤ変換器と、からなる地震計と、
圧力センサと、前記圧力センサから出力される矩形波を算出することで気圧データを生成する周波数算出部と、からなる気圧計と、
前記地震データと、前記気圧データとが入力されると共に、前記地震データと、前記気圧データとに基づいて、津波警報を発する警報判定部と、を有し、
前記周波数算出部は、
時間フレーム中の矩形波の数をカウントする際には、
時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅を、開始時矩形波の時間幅で除したものによりカウントし、
時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅を、終了時矩形波の時間幅で除したものによりカウントすることを特徴とする津波警報装置。

【請求項3】

地震センサと、前記地震センサからの出力をＡＤ変換することで地震データを生成するＡＤ変換器と、からなる地震計と、
圧力センサと、前記圧力センサから出力される矩形波を算出することで気圧データを生成する周波数算出部と、からなる気圧計と、
前記地震データと、前記気圧データとを、通信ネットワークを介して受信する受信部と、
前記受信部で受信した前記地震データと、前記気圧データとに基づいて、津波警報を発する警報判定部と、を有し、
前記周波数算出部は、
時間フレーム中の矩形波の数をカウントする際には、
時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅を、開始時矩形波の時間幅で除したものによりカウントし、
時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅を、終了時矩形波の時間幅で除したものによりカウントすることを特徴とする津波警報システム。

【発明の詳細な説明】

【技術分野】

【0001】

本発明は、気圧変動を計測するための圧力センサの出力周波数算出方法、さらに、地震動および気圧変動を計測し、その計測値から津波襲来の危険度を判定する津波警報装置、津波警報システムに関する。

【背景技術】

【0002】

現行の津波警報は、地震観測データから、地震の位置（緯度、経度、深さ）、および規模（マグニチュード）を推定し、地震が引き起こす津波高を理論的に予想することに寄っている。このため、実際に津波が発生したか否かの情報は用いられていない（非特許文献１）。

【0003】

これを改善するために、沖合の海上、海底に津波計を設置し、津波計のデータを用いることが検討されているが、海上、海底への津波計敷設には多くのコストがかかる。

【0004】

一方、巨大な津波に伴い、気圧変動が観測されることが知られている。これは津波による海面の変化が大気に作用して、その変動が気圧変動として観測されるものである。気圧変動の伝播速度は津波の伝播速度に比べると速いので、気圧変動を監視することで沿岸に津波が到達する前に津波発生の有無を知ることが出来る。ただし、気圧変動を引き起こす要因は火山噴火、隕石の飛来、等多岐に渡るため、気圧観測のみからは津波発生を判定することは難しい（非特許文献２）。

【先行技術文献】

【非特許文献】

【0005】

【非特許文献1】http://www.jma.go.jp/jma/press/1109/12a/torimatome.pdf（東北地方太平洋沖地震による津波被害を踏まえた津波警報の改善の方向性について）

【非特許文献2】「大地震，津波，火山大爆発などから発生した気圧波」地震第２輯、６４巻、４７−６２

【発明の概要】

【発明が解決しようとする課題】

【0006】

上記で問題となった、津波による気圧変動とその他の要因による気圧変動との違いを判定するためには、地震観測と気圧観測を同時に行うことで解決できる。津波発生は大きな地震の発生に伴うため、津波による気圧変動が観測される際には、それに先立ち顕著な地震動が観測されるからである。

【0007】

多くの場合地震観測はＧＰＳ信号等の高精度な時刻基準をもとに観測が行われている。これには、ＡＤ変換器のクロックを同期する方法がとられる。

【0008】

しかしながら、高精度気圧観測に用いられるセンサは周波数変動をカウントする方式のため、ＡＤ変換器のクロックを同期する方法をとることができず、高精度に時刻同期を行うことはされていなかった。気圧変動と地震動の時間関係が不明確では、地震動と気圧変動の双方を用いた処理ができず、津波による気圧変動とそれ以外の気圧変動の区別を行うことが出来ない。

【0009】

本発明は、このような課題を解決するために、高精度な時刻同期を行うことのできる圧力センサの出力周波数算出方法を提供すると共に、この出力周波数算出方法を用いて、地
震動と気圧の同時計測を行い、その計測値から津波襲来の危険度を判定する津波警報装置、津波警報システムを提供することを目的とする。

【課題を解決するための手段】

【0010】

前記課題を解決するために、本発明の請求項１に係る発明は、圧力センサから出力される矩形波の周波数を算出する圧力センサの出力周波数算出方法であって、時間フレーム中の矩形波の数をカウントする際には、時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅を、開始時矩形波の時間幅で除したものによりカウントし、時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅を、終了時矩形波の時間幅で除したものによりカウントすることを特徴とする。

【0011】

また、請求項２に係る発明は、地震センサと、前記地震センサからの出力をＡＤ変換することで地震データを生成するＡＤ変換器と、からなる地震計と、圧力センサと、前記圧力センサから出力される矩形波を算出することで気圧データを生成する周波数算出部と、からなる気圧計と、前記地震データと、前記気圧データとが入力されると共に、前記地震データと、前記気圧データとに基づいて、津波警報を発する警報判定部と、を有し、前記周波数算出部は、時間フレーム中の矩形波の数をカウントする際には、時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅を、開始時矩形波の時間幅で除したものによりカウントし、時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅を、終了時矩形波の時間幅で除したものによりカウントすることを特徴とする津波警報装置である。

【0012】

また、請求項３に係る発明は、地震センサと、前記地震センサからの出力をＡＤ変換することで地震データを生成するＡＤ変換器と、からなる地震計と、圧力センサと、前記圧力センサから出力される矩形波を算出することで気圧データを生成する周波数算出部と、からなる気圧計と、前記地震データと、前記気圧データとを、通信ネットワークを介して受信する受信部と、前記受信部で受信した前記地震データと、前記気圧データとに基づいて、津波警報を発する警報判定部と、を有し、前記周波数算出部は、時間フレーム中の矩形波の数をカウントする際には、時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅を、開始時矩形波の時間幅で除したものによりカウントし、時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅を、終了時矩形波の時間幅で除したものによりカウントすることを特徴とする津波警報システムである。

【発明の効果】

【0015】

本発明に係る圧力センサの出力周波数算出方法よれば、正確な時刻同期を行うことが可能となる。さらに、本発明に係る津波警報装置および津波警報システムによれば、陸上の観測のみで津波発生の有無を判定し、津波警報を配信することが可能となる。

【図面の簡単な説明】

【0016】

【図1】本発明の実施形態に係る津波警報装置１００の主要構成を示す図である。

【図2】従来の気圧センサの周波数算出方法を模式的に説明する図である。

【図3】本発明に係る気圧センサの周波数算出方法を模式的に説明する図である。

【図4】本発明の実施形態に係る津波警報装置１００の地震発生判定処理のフローチャートを示す図である。

【図5】本発明の実施形態に係る津波警報装置１００の津波発生判定処理のフローチャートを示す図である。

【図6】本発明の実施形態に係る地震警報システム３００の主要構成を示す図である。

【発明を実施するための形態】

【0017】

以下、本発明の実施の形態を図面を参照しつつ説明する。図１は本発明の実施形態に係る津波警報装置１００の主要構成を示す図である。津波警報装置１００は、南北、東西、上下方向の地震動を計測する３つの地震計１を備える。それぞれの地震計１は地震センサ２の出力を高精度に時刻同期されたＡＤ変換器３で、デジタル値に変換し、これを地震データとして、処理部１１に送信する。

【0018】

ここで、処理部１１としては、データの送受、データの演算、データの記憶を行うことができる汎用の情報処理装置を用いることができる。

【0019】

また、津波警報装置１００は、一つ以上の気圧を計測する気圧計５を備える。気圧計５は圧力センサ６の出力を高精度に時刻同期された周波数算出部７で、デジタル値に変換し気圧データとして、処理部１１に送信する。処理部１１における警報判定部１２では地震計１からの地震動データ、気圧計５からの気圧データを基に、地震発生判定および津波警報判定を行う。

【0020】

ここで、気圧計５における圧力センサ６としては、パロサイエンティフィック（Paroscientific）社製のものを用いることができる。この圧力センサ６による高精度の気圧観測は圧力に応じてその周波数が変動する矩形波を出力する周波数出力方式のもの（http://www.paroscientific.com/Nano-Resolution.pdf参照）で、圧力センサ６の出力である矩形
波の周波数を算出することで気圧データを算出する。このような算出処理を行うのが周波数算出部７である。

【0021】

警報判定部１２において、津波警報判定がなされた場合は、警報が、通信部２１から通信回線を介して、遠隔地に向けて送信がなされる。また、適宜、表示部２２への表示および警報出力部２３への警報出力がなされる。津波警報装置１００は蓄電池２６でバックアップがなされた電源装置２５から供給される電源２４で動作する。

【0022】

以上のように構成される周波数算出部７における周波数算出方法について、具体的に説明する。ここで、まず、従来の周波数算出部７における周波数算出方法について説明し、次いで、本発明に係る圧力センサ６の周波数算出方法について説明することとする。

【0023】

図２は従来技術（より詳しくは、例えば、http://www.jamstec.go.jp/scdc/top#j.htmlを参照。この例は水圧計測であるが気圧計測と計測方法は同一である）による周波数算出を模式的に表したものである。

【0024】

ここで、時間フレームはＧＰＳ信号等の高精度な時刻標準と同期している。この時間フレームの時間間隔（ΔＴ）は、一般的なＡＤ変換のサンプリング間隔に相当するものである。

【0025】

周波数算出は、時間フレーム（ΔＴ）中に含まれる、矩形波の数の合計を、矩形波の占める時間幅で除して行う。なお、矩形波の時間幅は、矩形波をゲートとして、矩形波の周波数より高速なクロック（例えば１０ＭＨｚ）をカウントすることにより計測することが
できる。ここで時間フレーム（ΔＴ）中に含まれる矩形波とは、矩形波終了の立ち上がりが時間フレーム中に含まれる矩形波であるものとする。すなわち、図２に示す１からＮまでナンバリングされた矩形波が、周波数を算出するためにカウントされる。

【0026】

図２において、Ｎは時間フレーム中に含まれるものとカウントされた矩形波の数、
Ｔ_iはｉ番目の矩形波の時間幅とすれば、圧力センサからの出力周波数Ｆ₁は、下式（１）によって算出される。
Ｆ₁＝Ｎ／（Ｔ₁＋Ｔ₂＋…＋Ｔ_N） ・・・（１）
しかしながら、この従来の周波数算出方法では、矩形波の占める時間幅は、（Ｔ₁＋Ｔ₂＋…＋Ｔ_N）であり、ΔＴとは必ずしも一致せず変動することとなる。このように、従来
方法では、一定間隔のサンプリングを行う事ができないため正確な時刻同期ができない。

【0027】

一方、図３は周波数算出部７で行われる、本発明による周波数算出方法を表したものである。周波数算出部７は、時間フレーム中に含まれる矩形波の数をカウントするが、その際、時間フレームの開始と終了にかかる矩形波の数については、その矩形波が時間フレーム中に占める時間の割合で案分する。

【0028】

より具体的には、図３において、
ΔＴは時間フレームの時間幅、
１とナンバリングされた矩形波は時間フレームの開始時間に重なる開始時矩形波、
Ｎ＋１とナンバリングされた矩形波は時間フレームの終了時間に重なる終了時矩形波、
Ｔ₁は１番目の矩形波（開始時矩形波）の時間幅、
Ｔ_N+1はＮ＋１番目の矩形波（終了時矩形波）の時間幅、
Ｔ_aは１番目の矩形波が着目している時間フレーム中に含まれる時間幅、
Ｔ_bはＮ＋１番目の矩形波が着目している時間フレーム中に含まれる時間幅、
とすれば、圧力センサの出力周波数Ｆ₂を下式（２）によって算出する。
Ｆ₂＝（Ｔ_a／Ｔ₁＋Ｔ_b／Ｔ_N+1＋Ｎ−１）／ΔＴ・・・（２）
すなわち、周波数算出部７における周波数算出方法においては、時間フレーム中の矩形波の数をカウントする際には、時間フレームの開始時間に重なる開始時矩形波は、開始時矩形波が時間フレーム中に含まれる時間幅（Ｔ_a）を、開始時矩形波の時間幅（Ｔ₁）で除したものによりカウントし、時間フレームの終了時間に重なる終了時矩形波は、終了時矩形波が時間フレーム中に含まれる時間幅（Ｔ_b）を、終了時矩形波の時間幅（Ｔ_N+1）で除したものによりカウントするようにしている。

【0029】

上記のような周波数算出部７における周波数算出方法によれば、圧力センサ６の出力をＧＰＳ信号等に同期した一定間隔でサンプリングすることが可能になる。地震動は既存技術によってＧＰＳ信号等に同期した一定間隔でサンプリングすることができるので、高精度に時刻同期した、気圧と地震の同時観測が可能になる。

【0030】

次に、以上のように構成される津波警報装置１００の処理・動作について説明する。図４は本発明の実施形態に係る津波警報装置１００の地震発生判定処理のフローチャートを示す図である。また、このような処理は、処理部１１の警報判定部１２で実行される。

【0031】

図４において、ステップＳ１００で地震発生判定処理が開始されると、続くステップＳ１０１においては、地震計１から地震動データを取得する。次に、ステップＳ１０２においては、地震動データをフィルタ処理（若しくは微分処理）することにより地震波形を生成する。次に、ステップＳ１０３において、前ステップで生成された地震波形の絶対値を演算する。

【0032】

ステップＳ１０４においては、当該絶対値が所定の閾値を超えるか否かが判定される。
ステップＳ１０４の判定がＮＯである場合には、ステップＳ１０１に戻る。

【0033】

一方、ステップＳ１０４の判定がＹＥＳである場合には、地震が発生したものと判定し、ステップＳ１０５に進む。

【0034】

ステップＳ１０５では、所定の一定時間が経過したか否かが判定される。ステップＳ１０５の判定がＮＯである間には、ステップＳ１０５をループし、ステップＳ１０５の判定がＹＥＳとなると、地震が終了したものと判定して、再びステップＳ１０１に戻る。

【0035】

次に、津波警報装置１００の津波発生判定処理について説明する。図５は本発明の実施形態に係る津波警報装置１００の津波発生判定処理のフローチャートを示す図である。また、このような処理は、処理部１１の警報判定部１２で実行される。

【0036】

ステップＳ２００において、津波発生判定処理が開始されると、続くステップＳ２０１においては、気圧計５から気圧データを取得する。次に、ステップＳ２０２においては、気圧データをフィルタ処理（若しくは微分処理）することにより気圧波形を生成する。次に、ステップＳ２０３において、前ステップで生成された気圧波形の絶対値を演算する。

【0037】

ステップＳ２０４においては、当該絶対値が閾値を超えるか否かが判定される。ステップＳ２０４の判定がＮＯである場合には、ステップＳ２０１に戻る。

【0038】

一方、ステップＳ２０４の判定がＹＥＳである場合には、津波発生可能性ありと判定し、続いてステップＳ２０５に進み、先の地震発生判定処理を参照して、地震が発生中であるか否かが判定される。

【0039】

ステップＳ２０５における判定がＮＯである場合には、気圧の変化が地震によるものではなく、従って、津波の発生がないものと判定し、ステップＳ２０１に戻る。

【0040】

一方、ステップＳ２０５における判定がＹＥＳである場合には、気圧の変化が地震によるものであり、津波が発生したものと判定し、ステップＳ２０６に進む。

【0041】

ステップＳ２０６においては、通信部２１、表示部２２や警報出力部２３などを利用して、各種警報を発する。

【0042】

以上のように、本発明に係る圧力センサ６の出力周波数算出方法よれば、正確な時刻同期を行うことが可能となる。さらに、本発明に係る津波警報装置１００によれば、陸上の観測のみで津波発生の有無を判定し、津波警報を配信することが可能となる。

【0043】

次に、本発明に係る圧力センサの出力周波数算出方法を津波警報システムに適用した例について説明する。図６は本発明の実施形態に係る地震警報システム３００の主要構成を示す図である。

【0044】

先の津波警報装置１００のように地震計１、気圧計５や警報判定部１２といった構成を、一つの装置として構成することも可能であるが、地震計１、気圧計５や警報判定部１２といった構成を分散すると共に、各構成同士をデータ通信可能に接続し、図６に示すような津波警報システム３００として構成することも可能である。

【0045】

すなわち、警報判定部１２は、地震計１、気圧計５は同一の地点に存在してもよいし、地震データ、気圧データを通信回線等で遠隔地にあるセンター等に設けられた地震気圧データ処理装置２００に伝送し、そこで地震発生判定と津波警報判定を行ってもよい。

【0046】

図６において、２０１は地震計、２０５は気圧計、１５０は通信ネットワーク、２００は地震気圧データ処理装置である。

【0047】

津波警報システム３００で用いる地震計２０１及び気圧計２０５には、それぞれ地震データ及び気圧データを、通信ネットワーク１５０を介して、地震気圧データ処理装置２００に送信する機能を有している。

【0048】

地震気圧データ処理装置２００は、受信部２１１で地震データ、気圧データを受信する。また、図１で示した津波警報装置１００と、同様の警報判定部２１２をもち、地震発生判定と津波警報判定を行うことができるようになっている。

【0049】

なお、地震気圧データ処理装置２００は、データセンター（通信ネットワーク１５０上に複数あっても一つであってもよい）に設置される情報処理装置であり、データ処理機能、データ記憶機能、データの送受信機能などを備える、例えば、汎用のサーバーを用いることができる。

【0050】

津波警報判定がなされた場合は配信部２１３より警報配信を行う。地震気圧データ処理装置２００における配信部２１３は津波警報を情報端末（携帯電話、スマートフォン、パーソナルコンピュータ、テレビ、ラジオなど）に地震警報を配信する。

【0051】

このシステムでは、地震計２０１と気圧計２０５で観測された、地震データ、気圧データを通信ネットワーク１５０経由でセンターに転送し、離れた場所で地震発生判定と津波警報判定を行うものである。ここで、地震計２０１、気圧計２０５および地震気圧データ処理装置２００はそれぞれ複数設けるようにしてもよい。

【0052】

以上のように構成される、本発明に係る津波警報システム３００によれば、陸上の観測のみで津波発生の有無を判定し、津波警報を配信することが可能となる。

【符号の説明】

【0053】

１・・・地震計
２・・・地震センサ
３・・・ＡＤ変換器
５・・・気圧計
６・・・圧力センサ
７・・・周波数算出部
１１・・・処理部
１２・・・警報判定部
２１・・・通信部
２２・・・表示部
２３・・・警報出力部
２４・・・電源
２５・・・電源装置
２６・・・蓄電池
１００・・・津波警報装置
１５０・・・通信ネットワーク
２００・・・地震気圧データ処理装置
２０１・・・地震計
２０５・・・気圧計
２１１・・・受信部
２１２・・・警報判定部
２１３・・・配信部
３００・・・津波警報システム

【図1】

【図2】

【図3】

【図4】

【図5】

【図6】