(57)【要約】      （修正有）【課題】媒質に伝わる波動を測定することで、その測定した二点間の伝播時間をリアルタイムに推定することができる二点間の波動伝播時間の推定方法を提供する。【解決手段】二点間を波動が伝播することを考える。それぞれの位置をｘ、ｙとして、それぞれの点で観測された波形の相互相関関数を算定する。位置ｘにおけるある時点ｊの波形を基準として、この波形と位置ｙにおけるある時点ｊから時間ｉさかのぼった波形との相互相関関数を時間ｉをサンプル時間間隔を単位として変動させ、相互相関関数の絶対値が最大になる場合の時間ｉが、ある時点ｊにおける波動伝播時間となる。実際に処理を行うにあたっては、前述の位置ｘにおける波形を基準にした場合と、位置ｙにおける波形を基準にした場合の相互相関係数を、それぞれ時間差ｉをサンプル時間を単位として複数さかのぼって算出する。【選択図】図１
2019082353.tif 000002
【特許請求の範囲】
【請求項１】
二点間の波動伝播時間の推定方法において、
空間的に離れた二点間を波動が伝播する場合に、
おのおのの点において該波動を時系列の振動波形として計測し、
一方の点において計測された時系列の振動波形についてサンプリング時間をさかのぼった時系列の振動波形と、
他方の点において計測された時系列の振動波形と、
の相互相関関数を、
各測定時点について該サンプリング時間をさかのぼって連続的に算出したものを第一の相互相関関数群とし、
また、
一方の点において計測された時系列の振動波形と、
他方の点において計測された時系列の振動波形についてサンプリング時間をさかのぼった時系列の振動波形と、
の相互相関関数を、
各測定時点について該サンプリング時間をさかのぼって連続的に算出したものを第二の相互相関関数群とし、
該第一の相互相関関数群と該第二の相互相関関数群との中で、
相互相関係数の絶対値が最大になる該サンプリング時間のさかのぼり量を、
該測定時刻における該二点間の波動伝播時間とする、
ことを特徴とする、
二点間の波動伝播時間の推定方法。
【発明の詳細な説明】
【技術分野】
【０００１】
本発明は波動を伝播する媒質の二点間の波動伝播時間の推定方法に関する。
【背景技術】
【０００２】
波動の伝播時間を推定する方法としては、媒質中の二点で波動を観測し、他方の波動からもう一点で観測した波動を推定し、その差が最小になるサンプリング時間を波動の伝播時間とする方法（例えば特許文献１）や、波動の伝播時間を推定する方法としては、地震波の伝播時間を推定するＮＩＯＭ法（例えば、非特許文献１）や、ダムの波動伝播時間の推定方法（例えば、非特許文献２参照）が知られている。
波動の伝播時間を推定する方法としては、媒質中の二点で波動を観測し、他方の波動からもう一点で観測した波動を推定し、その差が最小になるサンプリング時間を波動の伝播時間とする方法（例えば特許文献１）や、地震波の伝播時間を推定するＮＩＯＭ法（例えば非特許文献１）、ダムの波動伝播時間の推定方法（例えば、非特許文献２参照）が知られている。
【０００３】
まず特許文献１による方法であるが、この方法では媒質内の波動場が大きく乱された場合に誤差が大きくなり、例えば地盤などに適用する場合には液状化の影響を受けやすいという特徴がある。
【０００４】
また、非特許文献１のＮＩＯＭ法は、波動を２点で観測し、それらの観測データに対して一定の時間間隔を取出しフーリエ変換を行い、その結果からその２点間の伝達関数を求めておき、パルス状の入力波形と出力波形の関係をこの伝達関数から求めることで、波動の伝播時間を求めるものである。
【０００５】
しかし、この方法によると一定の測定区間を用いてフーリエ変換を行うためデータ処理に遅れが生じ、また処理そのものも複雑になる。
【０００６】
さらに、非特許文献２の大町ほかによるダムの波動伝播速度の算出方法は、波動に含まれるある振動数が２点間を伝播する場合に位相が遅れることに着目したもので、フーリエ位相スペクトルを用いて振動数と位相差の関係を求めて、そこから伝播速度を算出する。
【０００７】
しかし、この方法の場合でも一定の測定区間を用いてフーリエ変換を行うためデータ処理に遅れが生じ、また処理そのものも複雑になる。さらに、この方法を双方向の伝播波動がほぼ等しい建物に適用した場合は、その建物内の上下２点では位相差が見かけ上なくなり、波動伝播時間の算定ができなくなるという問題もある。
【先行技術文献】
【特許文献】
【０００８】
【特許文献１】特許第６０２４０１２号特許公報
【非特許文献】
【０００９】
【非特許文献１】Kawakami, H. and Haddadi, H. R.: Modeling wave propagation by using Normalized Input-Output Minimization (NIOM), Soil Dyn. Earthq. Engng., 17, pp.117-126, 1998.
【００１０】
【非特許文献２】大町達夫他：直下地震の観測記録に基づくロックフィルダムの非線形地震応答特性，第54回地盤工学シンポジウム 平成21年度論文集，pp.243-250，2009.
【発明の概要】
【発明が解決しようとする課題】
【００１１】
本件発明が解決しようとする課題は、物性の変化や構造物などの劣化を評価するために波動伝播時間を求めようとする際に、用いるデータ個数に制限なくリアルタイム処理ができないことである。
【課題を解決するための手段】
【００１２】
二点間を波動が伝播することを考える。それぞれの位置をｘ、ｙとすると、ある時点ｊにおいてそれぞれの点で観測された波形はｆ（ｘ，ｊ）、ｆ（ｙ，ｊ）とあらわされる。このとき、これら二つの波形の相互相関関数ＴＲ（ｊ）は以下のように算定することができる。
【００１３】
【数１】
2019082353.tif 000003
【００１４】
ここで、演算子［Ｇ（ｘ）］はＧ（ｘ）のある時間にわたる平均演算を意味する。また、
Ｒ（ｉ，ｊ）=［ｆ（ｘ，ｊ）×ｆ（ｙ，ｊ−ｉ）］
Ｒｘ（ｊ）＝［ｆ（ｘ、ｊ）２］
Ｒｙ（ｊ）＝［ｆ（ｙ、ｊ）２］
であり、ｓｉｇｎ（）は符号関数である。
【００１５】
このとき、ｉは二つの波形の時間差を表しており、数１は位置ｘにおけるある時点ｊの波形を基準として、この波形と位置ｙにおけるある時点ｊから時間ｉさかのぼった波形との相互相関関数を算定するものである。この時間ｉをサンプル時間間隔を単位として変動させ、相互相関関数の絶対値が最大になる場合の時間ｉが、ある時点ｊにおける波動伝播時間となる。
【００１６】
実際に処理を行うにあたっては、前述の位置ｘにおける波形を基準にした場合と、位置ｙにおける波形を基準にした場合の相互相関係数を、それぞれ時間ｉをサンプル時間を単位として複数さかのぼって算出する。この処理を行うことでリアルタイムかつ時間遅れなく位置ｘと位置ｙの間の波動伝播時間を算出することができる。
【００１７】
また、時系列の振動波形において、連続するサンプリングデータを補間することで、サンプリング間隔よりも短い時間差の波動伝播時間を推定することができる。なお、補間の方法は、たとえば直線的に二つのサンプリングデータ間を補間するなど、その方法に限定はない。
【発明の効果】
【００１８】
本発明の方法によれば、二点間を伝わる波動を測定することで、その二点間の波動伝播時間をリアルタイムに時間遅れなく推定することができる。このため、従来の方法では二点で観測を行い、限られたデータ長の波形を回収し、オフラインで複雑な処理を行う必要があったものを、非常に長いデータ長の観測波形であっても観測された波動を逐次処理し、連続的に波動伝播時間を算出することができる。また、振動場が乱された場合であっても波動伝播時間を算出することができる。
【００１９】
本発明の方法によれば、媒質に伝わる波動を測定することで、その測定した二点間の伝播時間をリアルタイムに推定することができる。このため非常に長いデータ長の観測波形であっても逐次処理し、連続的に算出することができる。
【図面の簡単な説明】
【００２０】
【図１】本発明の手法により求めた波動伝播時間を、特許文献１で求めた波動伝播時間と比較して示す。
【発明を実施するための形態】
【００２１】
以下、本発明を実施するための形態を示す。
【実施例】
【００２２】
本発明の手法を、東北大学工学部の人間・環境系研究棟で観測された東北地方太平洋沖地震の観測波形に適用した実施例を示す。
【００２３】
本発明の手法によると、図１に示すような波動伝播時間の時系列の変動をリアルタイムに得ることができる。なお、図１は特許文献１の方法により求められた波動伝播時間に、本発明の手法により求めた波動伝播時間を重ね書きしたものである。この図から、本発明の手法による波動伝播時間は特許文献１の方法による波動伝播時間とよく一致しており、本発明の手法の有効性を示している。
【００２４】
なお、本発明の実施例に用いた地震観測波形は独立行政法人建築研究所により提供された強震波形である。
【産業上の利用可能性】
【００２５】
二点間の波動伝播時間をリアルタイムに推定することができるので、例えば波動が伝播する媒質の物性の変化や、具体的には構造物の地震時の被害状況を、リアルタイムに把握することなどに資することができる。
【符号の説明】
【００２６】
１  波動伝播時間を示す軸（単位：１／１００秒）
２  時間軸（単位：秒）
１１  本発明の手法による波形伝播時間
１２  特許文献１の手法による波形伝播時間
【図１】
2019082353.tif 000004
【手続補正書】
【提出日】平成30年2月16日(2018.2.16)
【手続補正１】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
実際に処理を行うにあたっては、前述の位置ｘにおける波形を基準にした場合と、位置ｙにおける波形を基準にした場合の相互相関係数を、それぞれ時間ｉをサンプル時間を単位として複数さかのぼって算出する。波形データをサンプルするごとにこの計算を行うとともに、現在時点を基準とした適当な時間区間で平均化処理を行うことで、リアルタイムかつ時間遅れなく位置ｘと位置ｙの間の波動伝播時間を算出することができる。
【手続補正書】
【提出日】平成30年5月2日(2018.5.2)
【手続補正１】
【補正対象書類名】特許請求の範囲
【補正対象項目名】全文
【補正方法】変更
【補正の内容】
【特許請求の範囲】
【請求項１】
二点間の波動伝播時間の推定方法において、
空間的に離れた二点間を波動が伝播する場合に、
おのおのの点において該波動を時系列の振動波形として計測し、
一方の点において計測された時系列の振動波形についてサンプリング時間をさかのぼった時系列の振動波形と、
他方の点において計測された時系列の振動波形と、
の相互相関関数を、
各測定時点について該サンプリング時間をさかのぼって連続的に算出したものを第一の相互相関関数群とし、
また、
一方の点において計測された時系列の振動波形と、
他方の点において計測された時系列の振動波形についてサンプリング時間をさかのぼった時系列の振動波形と、
の相互相関関数を、
各測定時点について該サンプリング時間をさかのぼって連続的に算出したものを第二の相互相関関数群とし、
該第一の相互相関関数群と該第二の相互相関関数群との中で、
相互相関係数の絶対値が最大になる該サンプリング時間のさかのぼり量を測定時点ごとにリアルタイムに算出し、
算出された相互相関係数の絶対値が最大になるサンプリング時間の該さかのぼり量を、
該測定時点における該二点間の波動伝播時間とする、
ことを特徴とする、
二点間の波動伝播時間の推定方法。
【手続補正２】
【補正対象書類名】明細書
【補正対象項目名】００１６
【補正方法】変更
【補正の内容】
【００１６】
実際に処理を行うにあたっては、前述の位置ｘにおける波形を基準にした場合と、位置ｙにおける波形を基準にした場合の相互相関係数を、それぞれ時間ｉをサンプル時間を単位として複数さかのぼって算出する。この処理を行うことでリアルタイムかつ時間遅れなく位置ｘと位置ｙの間の波動伝播時間を算出することができる。