JP4173408B2 - 地震動識別方法及びそのための装置 - Google Patents
地震動識別方法及びそのための装置 Download PDFInfo
- Publication number
- JP4173408B2 JP4173408B2 JP2003175176A JP2003175176A JP4173408B2 JP 4173408 B2 JP4173408 B2 JP 4173408B2 JP 2003175176 A JP2003175176 A JP 2003175176A JP 2003175176 A JP2003175176 A JP 2003175176A JP 4173408 B2 JP4173408 B2 JP 4173408B2
- Authority
- JP
- Japan
- Prior art keywords
- motion
- seismic
- initial
- fitting
- time
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Lifetime
Links
Images
Landscapes
- Geophysics And Detection Of Objects (AREA)
Description
【０００１】
【発明の属する技術分野】
本発明は、地震動識別方法に係り、特に、地震のデータからそれが解析対象となる地震動であるか、ノイズであるかを自動的に識別する方法及びそのための装置に関するものである。
【０００２】
【従来の技術】
本願発明者は、既に、地震波の初動部分の波形形状をパラメータが数個の簡易な関数でフィッティングしてその波形形状を定量化し、得られたパラメータから震央距離とマグニチュードを推定することができる、震央距離及びマグニチュード推定方法とそのための装置（下記特許文献１参照）を提案している。
【０００３】
【特許文献１】
特開２００２−２７７５５７号公報（第３−４頁 図１）
【０００４】
【発明が解決しようとする課題】
本発明は、上記した特許文献１に示される関数フィッティング法に基づいて、解析対象となる地震動であるかノイズであるかを簡便に、かつ自動的に識別することができる地震動識別方法及びそのための装置を提供することを目的とする。
【０００５】
【課題を解決するための手段】
本発明は、上記目的を達成するために、
〔１〕地震動識別方法において、（ａ）地震計から得られる時系列データの絶対値をＶ（ｔ）となし〔ここで、ｔは時間（秒）〕、地震を検知した時刻を時間原点（ｔ＝０）として、そこから、数秒間のデータに簡易な関数を当てはめ、関数形として、Ｖ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）を用い（ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ）、上記式の両辺の常用対数をとって、ｌｏｇ〔Ｖ（ｔ）／ｔ〕＝ｌｏｇＢ−Ａ・ｌｏｇｅ・ｔを得て、未知パラメータを線形化し、その形状パラメータを通常の最小二乗法を用いて算出し、前記初動部分の形状パラメータＡ，Ｂを求め、（ｂ）この初動部分の形状パラメータＡ，Ｂに閾値ＴA ，ＴB を設定し、（ｃ）この設定された閾値ＴA ，ＴB を有するフィッティングしたデータに基づいて、地震動であるか、それ以外のノイズであるかを自動的に識別することを特徴とする。
【０００６】
〔２〕上記〔１〕記載の地震動識別方法において、前記関数のパラメータが係数Ａであり、この係数Ａが閾値ＴA より小さい場合に地震動であると識別することを特徴とする。
【０００７】
〔３〕上記〔１〕記載の地震動識別方法において、前記関数のパラメータが係数Ｂであり、この係数Ｂが閾値ＴB より大きい場合に地震動であると識別することを特徴とする。
【０００８】
〔４〕地震動識別方法において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化し、（ｂ）これによって得られた前記関数Ｙのフィッティング残差Ｚ
【０００９】
【数３】
を設定し、（ｃ）この設定されたフィッティング残差Ｚが閾値ＴZ より小さい場合に地震動であると識別することを特徴とする。
【００１０】
〔５〕地震動識別方法において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化し、（ｂ）前記関数フィッティングによる、前記地震計の検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別することを特徴とする。
【００１１】
〔６〕地震動識別方法において、（ａ）上記〔２〕記載の係数Ａが閾値ＴA より小さく、（ｂ）上記〔３〕記載の係数Ｂが閾値ＴB より大きく、（ｃ）上記〔４〕記載のフィッティング残差Ｚが閾値ＴZ より小さく、（ｄ）上記〔５〕の最大振幅Ａmax が閾値Ｇ1 を超えるときに地震動であると総合的に識別することを特徴とする。
【００１２】
〔７〕地震動識別装置において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、（ｂ）前記初動部分の波形形状パラメータ（係数）Ａ，Ｂを求める手段と、（ｃ）この初動部分の波形形状パラメータ（係数）Ａ，Ｂに閾値ＴA ，ＴB を設定する手段と、（ｄ）この設定された閾値ＴA ，ＴB を有するフィッティングしたデータに基づいて地震動であるか、それ以外のノイズであるかを自動的に識別する手段とを具備することを特徴とする。
【００１３】
〔８〕地震動識別装置において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、（ｂ）これによって得られた前記関数Ｙのフィッティング残差Ｚ
【００１４】
【数４】
を設定する手段と、（ｃ）このフィッティング残差Ｚが、このフィッティング残差Ｚの閾値ＴZ より小さい場合に地震動であると識別する手段とを具備することを特徴とする。
【００１５】
〔９〕地震動識別装置において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、（ｂ）前記関数フィッティングによる、前記地震計の検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別する手段とを具備することを特徴とする。
【００１６】
〔１０〕地震動識別装置において、（ａ）上記〔７〕記載の係数に基づく地震動の識別情報と、（ｂ）上記〔８〕記載のフィッティング残差Ｚに基づく地震動の識別情報と、（ｃ）上記〔９〕記載の加速度エンベロープの最大振幅Ａmax に基づく地震動の識別情報と、（ｄ）上記（ａ）、（ｂ）及び（ｃ）を論理和として地震動を判定する手段とを具備することを特徴とする。
【００１７】
【発明の実施の形態】
以下、本発明の実施の形態について詳細に説明する。
【００１８】
地震波形の時系列データをｖ（ｔ）とする。更に、この時系列データｖ（ｔ）の絶対値をＶ（ｔ）とする。ここで、ｔは時間（秒）を表す。
【００１９】
まず、地震を検知した時刻を時間原点（ｔ＝０）として、そこから、数秒間のデータに簡易な関数を当てはめる。関数形として、ここでは、次式（１）を用いる。なお、絶対値Ｖ（ｔ）の時間的な変動が大きい場合は、適当な方法で絶対値Ｖ（ｔ）を平滑化した後に関数の当てはめを行う。
【００２０】
Ｖ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ） …（１）
ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関するパラメータであり、以下、これらを形状パラメータと仮称する。上記（１）式の両辺の常用対数をとると、次式（２）のように未知パラメータが線形化されるので、その形状パラメータは通常の最小二乗法を用いて算出できる。
【００２１】
ｌｏｇ〔Ｖ（ｔ）／ｔ〕＝ｌｏｇＢ−Ａ・ｌｏｇｅ・ｔ …（２）
上記の方法を用いて多くの地震波形データについて初動部の形状パラメータＡ，Ｂを求める。
【００２２】
図３は実際の地震波に簡易な関数をフィッティングした例を示す図であり、図３（ａ）は実際の初動波（Ｐ波）、図３（ｂ）はその簡易な関数のフィッティング例を示す図である。
【００２３】
このように、多くの地震データに対して、同様のフィッティングを行い、初動部分の形状パラメータを求め、このパラメータを用いて地震動識別処理を行う。
【００２４】
以下、その地震動識別処理の具体例を詳細に説明する。
【００２５】
図１は本発明に係る地震動識別装置のブロック図、図２は本発明に係る地震動識別処理の概略フローチャートである。
【００２６】
これらの図において、１は地震計、２は通信回線、１０は制御処理装置、１１はアンチエイリアシングフィルタ（ローパスフィルタ）、１２はＡ／Ｄ変換器、１３は波形収録部、１４はメモリ、１５はパラメータ演算部、１６は閾値設定部、１７は地震動識別部、１８は条件判定部、１９は情報文送出部、２０はＧＰＳ時計、２１は状態表示装置、２２はモデムである。
【００２７】
この地震動識別方法を、図２に示すフローを参照しながら説明する。
【００２８】
（１）地震計１からの情報を制御処理装置１０に取込み、ディジタル波形データとしてメモリ１４に取得する（ステップＳ１）。
【００２９】
（２）そこで、パラメータ演算部１５において、オフセットの除去を行い（ステップＳ２）、次に、データの絶対値を取得する。このとき絶対値がある基準値以下の場合には、強制的に最小基準値とする（ステップＳ３）。次に、スムージング（例えば、波形のエンベロープを取る）を行う（ステップＳ４）。次に、ｙ＝Ｂｔ・ｅ-At とフィッティングを行い形状パラメータ（係数）Ａ，Ｂを求める（ステップＳ５）。
【００３０】
（３）次いで、閾値設定部１６において、例えば、前記関数の形状パラメータ（係数）Ａ，Ｂに閾値ＴA ，ＴB を設定する（ステップＳ６）。
【００３１】
（４）次に、設定された閾値ＴA ，ＴB に基づいて地震動識別部１７で地震動識別処理を行う（ステップＳ７）。
【００３２】
（５）次に、条件判定部１８において、地震動識別処理された情報に基づき地震動であるか否かの総合判定を行う（ステップＳ８）。
【００３３】
（６）ステップＳ８において、ＹＥＳの場合は、情報文送出部１９において、早期検知情報の送信を行う（ステップＳ９）。
〔具体的な地震動識別方法〕
上記したように、これまでに地震動初動部分（Ｐ波部分）に簡単な関数
Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ） （ｔ：時間）
をフィッティングすることによって、係数（形状パラメータ）Ａ，Ｂを決定し、地震のマグニチュードを地震検知から数秒間で推定する方法を開発している（上記特許文献１参照）。
【００３４】
本発明では、このフィッティングした際に得られる形状パラメータＡ，Ｂやフィッティング残差Ｚに閾値を設定することによって、自動的に地震動の識別を行う。
【００３５】
図４は本発明の第１実施例を示す地震動識別方法の説明図であり、横軸は時間（ｓ）、縦軸は加速度（絶対値）を示す図である。
【００３６】
ここでは、まず、初動部分の最大振幅に関するパラメータである係数Ａによるノイズ識別を行う。
【００３７】
図４において、実線ａとして示すように、係数Ａがある閾値（基準値）ＴA より大きい場合に、ノイズと識別する。
【００３８】
Ａ≧ＴA ⇒ ノイズ
一方、破線ｂとして示すように、係数Ａがある閾値（基準値）ＴA より小さい場合に、地震動と識別する。
【００３９】
Ａ＜ＴA ⇒ 地震動
図５に実際の地震計で観測されたノイズデータに関数をフィッティングして得た係数Ａの頻度分布〔図５（ａ）〕と地震動データに関数をフィッティングして得た係数Ａの頻度分布〔図５（ｂ）〕をそれぞれ示す。
【００４０】
図５（ａ）に示すように、ノイズデータは、係数Ａが１．１をピークに持つ分布になるのに対して、図５（ｂ）に示すように、地震動データは、係数Ａが０．７をピークに持つ分布であり、ノイズデータと地震動データでは係数Ａの分布が異なる。この場合、係数Ａによるノイズ識別の閾値ＴA を１．０とすれば、ノイズデータの７１％をノイズと識別でき、地震動データの９５％を地震動と識別することができる。
【００４１】
図６は本発明の第２実施例を示す地震動識別方法の説明図であり、横軸は時間（ｓ）、縦軸は加速度（絶対値）を示す図である。
【００４２】
ここでは、係数Ｂによるノイズ識別を行う。
【００４３】
図６において、実線ａとして示すように、係数Ｂがある閾値（基準値）ＴB より小さい場合に、ノイズと識別する。
【００４４】
Ｂ≦ＴB ⇒ ノイズ
一方、破線ｂとして示すように、係数Ｂがある閾値（基準値）ＴB より大きい場合に、地震動と識別する。
【００４５】
Ｂ＞ＴB ⇒ 地震動
次いで、本発明の第３実施例を示すフィッティング残差Ｚによるノイズ識別について説明する。
【００４６】
フィッティング残差Ｚがある閾値（基準値）ＴZ より大きい場合に、ノイズと識別する。
【００４７】
Ｚ≧ＴZ ⇒ ノイズ
ただし、Ｚはフィッティングの残差二乗和の平均で、
【００４８】
【数５】
である。
【００４９】
次に、加速度エンベロープの振幅によるノイズの識別について説明する。
【００５０】
図７は本発明の第４実施例を示す地震動識別方法の説明図、図８は本発明の第４実施例を示すノイズ波の識別方法の説明図であり、横軸は時間、縦軸は加速度（絶対値）を示す図である。
【００５１】
図８に示すように、検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax （■）が閾値（基準値）Ｇ1 を超えない場合はノイズと識別する。なお、図７及び図８において、点線（◇）が加速度エンベロープを、実線（■）がその加速度エンベロープの最大振幅Ａmax を示している。すなわち、図７では、一定時間の半ばすぎほどで加速度エンベロープが閾値（基準値）Ｇ1 を超えているので、その最大振幅Ａmax もＧ1 を超えている。つまり、地震波と識別される。
【００５２】
一方、図８に示すような場合には、一定時間内で加速度エンベロープの最大振幅Ａmax が閾値（基準値）Ｇ1 を超えてはいないので、ノイズ波であると識別できる。
【００５３】
このように、地震計による検測時刻から一定時間内に最大振幅Ａmax （■）が閾値Ｇ1 を超える場合には、地震波であると識別する。
【００５４】
図９は本発明の実施例を示す総合的な高精度地震動識別フローチャートである。
【００５５】
（１）まず、パラメータ演算部１５において、関数ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）を地震波初動部分のデータにフィッティングし、係数Ａ，Ｂ、残差Ｚを求める（ステップＳ１１）。
【００５６】
（２）次に、閾値設定部１６において、係数Ａ，Ｂ、残差Ｚに閾値ＴA ，ＴB ，ＴZ を設定する（ステップＳ１２）。
【００５７】
（３）地震動識別部１７で係数Ａによる地震動識別（Ａ＜ＴA ）を行う（ステップＳ１３）。Ａ≧ＴA の場合にはノイズと識別する（ステップＳ１８）。
【００５８】
（４）次に、係数Ｂによる地震動識別（Ｂ＞ＴB ）を行う（ステップＳ１４）。Ｂ≦ＴB の場合にはノイズと識別する（ステップＳ１８）。
【００５９】
（５）次に、フィッティング残差Ｚによる地震動識別（Ｚ＜ＴZ ）を行う（ステップＳ１５）。Ｚ≧ＴZ の場合にはノイズと識別する（ステップＳ１８）。
【００６０】
（６）次に、検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax に基づいて、地震動識別（Ａmax ＞Ｇ1 ）を行う（ステップＳ１６）。Ａmax ≦Ｇ1 の場合には、ノイズと識別する（ステップＳ１８）。
【００６１】
（７）ステップＳ１６においてＹＥＳの場合には、条件判定部１８で地震動であると判断する（ステップＳ１７）。そして、地震動であると判別されると、情報文送出部１９において、早期検知情報の送信を行い、それに対応して措置を講じることができる。
【００６２】
このように、地震動であるかの判定処理を図９に示す処理フローに沿って行うことによって、精度の高い総合的な地震の自動識別処理を行うことができる。
【００６３】
なお、本発明は上記実施例に限定されるものではなく、本発明の趣旨に基づいて種々の変形が可能であり、これらを本発明の範囲から排除するものではない。
【００６４】
【発明の効果】
以上、詳細に説明したように、本発明によれば、以下のような効果を奏することができる。
【００６５】
（Ａ）関数フィッティング法に基づいて、解析対象となる地震動であるかノイズであるかを簡便に、かつ自動的に識別することができる。
【００６６】
（Ｂ）関数フィッティング法による係数Ａが閾値ＴA より小さい場合に地震動であると識別することができる。
【００６７】
（Ｃ）関数フィッティング法による係数Ｂが閾値ＴB より大きい場合に地震動であると識別することができる。
【００６８】
（Ｄ）関数フィッティング法によるフィッティング残差Ｚが閾値ＴZ より小さい場合に地震動であると識別することができる。
【００６９】
（Ｅ）関数フィッティング法による検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別することができる。
【００７０】
（Ｆ）上記（Ｂ）から（Ｅ）をＡＮＤ条件とした場合に、総合的な高精度の地震動の識別を行うことができる。
【図面の簡単な説明】
【図１】 本発明に係る地震動識別装置のブロック図である。
【図２】 本発明に係る地震動識別処理の概略フローチャートである。
【図３】 本発明にかかる簡易な関数のフィッティングの例を示す図である。
【図４】 本発明の第１実施例を示す地震動識別方法の説明図である。
【図５】 （ａ）は実際の地震計で観測されたノイズデータに関数をフィッティングして得た係数Ａの頻度分布を示す図、（ｂ）は地震動データに関数をフィッティングして得た係数Ａの頻度分布を示す図である。
【図６】 本発明の第２実施例を示す地震動識別方法の説明図である。
【図７】 本発明の第４実施例を示す地震動識別方法の説明図である。
【図８】 本発明の第４実施例を示すノイズ波の識別方法の説明図である。
【図９】 本発明の実施例を示す総合的な高精度地震動識別フローチャートである。
【符号の説明】
１ 地震計
２ 通信回線
１０ 制御処理装置
１１ アンチエイリアシングフィルタ（ローパスフィルタ）
１２ Ａ／Ｄ変換器
１３ 波形収録部
１４ メモリ
１５ パラメータ演算部
１６ 閾値設定部
１７ 地震動識別部
１８ 条件判定部
１９ 情報文送出部
２０ ＧＰＳ時計
２１ 状態表示装置
２２ モデム
【発明の属する技術分野】
本発明は、地震動識別方法に係り、特に、地震のデータからそれが解析対象となる地震動であるか、ノイズであるかを自動的に識別する方法及びそのための装置に関するものである。
【０００２】
【従来の技術】
本願発明者は、既に、地震波の初動部分の波形形状をパラメータが数個の簡易な関数でフィッティングしてその波形形状を定量化し、得られたパラメータから震央距離とマグニチュードを推定することができる、震央距離及びマグニチュード推定方法とそのための装置（下記特許文献１参照）を提案している。
【０００３】
【特許文献１】
特開２００２−２７７５５７号公報（第３−４頁 図１）
【０００４】
【発明が解決しようとする課題】
本発明は、上記した特許文献１に示される関数フィッティング法に基づいて、解析対象となる地震動であるかノイズであるかを簡便に、かつ自動的に識別することができる地震動識別方法及びそのための装置を提供することを目的とする。
【０００５】
【課題を解決するための手段】
本発明は、上記目的を達成するために、
〔１〕地震動識別方法において、（ａ）地震計から得られる時系列データの絶対値をＶ（ｔ）となし〔ここで、ｔは時間（秒）〕、地震を検知した時刻を時間原点（ｔ＝０）として、そこから、数秒間のデータに簡易な関数を当てはめ、関数形として、Ｖ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）を用い（ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ）、上記式の両辺の常用対数をとって、ｌｏｇ〔Ｖ（ｔ）／ｔ〕＝ｌｏｇＢ−Ａ・ｌｏｇｅ・ｔを得て、未知パラメータを線形化し、その形状パラメータを通常の最小二乗法を用いて算出し、前記初動部分の形状パラメータＡ，Ｂを求め、（ｂ）この初動部分の形状パラメータＡ，Ｂに閾値ＴA ，ＴB を設定し、（ｃ）この設定された閾値ＴA ，ＴB を有するフィッティングしたデータに基づいて、地震動であるか、それ以外のノイズであるかを自動的に識別することを特徴とする。
【０００６】
〔２〕上記〔１〕記載の地震動識別方法において、前記関数のパラメータが係数Ａであり、この係数Ａが閾値ＴA より小さい場合に地震動であると識別することを特徴とする。
【０００７】
〔３〕上記〔１〕記載の地震動識別方法において、前記関数のパラメータが係数Ｂであり、この係数Ｂが閾値ＴB より大きい場合に地震動であると識別することを特徴とする。
【０００８】
〔４〕地震動識別方法において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化し、（ｂ）これによって得られた前記関数Ｙのフィッティング残差Ｚ
【０００９】
【数３】
【００１０】
〔５〕地震動識別方法において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化し、（ｂ）前記関数フィッティングによる、前記地震計の検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別することを特徴とする。
【００１１】
〔６〕地震動識別方法において、（ａ）上記〔２〕記載の係数Ａが閾値ＴA より小さく、（ｂ）上記〔３〕記載の係数Ｂが閾値ＴB より大きく、（ｃ）上記〔４〕記載のフィッティング残差Ｚが閾値ＴZ より小さく、（ｄ）上記〔５〕の最大振幅Ａmax が閾値Ｇ1 を超えるときに地震動であると総合的に識別することを特徴とする。
【００１２】
〔７〕地震動識別装置において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、（ｂ）前記初動部分の波形形状パラメータ（係数）Ａ，Ｂを求める手段と、（ｃ）この初動部分の波形形状パラメータ（係数）Ａ，Ｂに閾値ＴA ，ＴB を設定する手段と、（ｄ）この設定された閾値ＴA ，ＴB を有するフィッティングしたデータに基づいて地震動であるか、それ以外のノイズであるかを自動的に識別する手段とを具備することを特徴とする。
【００１３】
〔８〕地震動識別装置において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、（ｂ）これによって得られた前記関数Ｙのフィッティング残差Ｚ
【００１４】
【数４】
【００１５】
〔９〕地震動識別装置において、（ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、（ｂ）前記関数フィッティングによる、前記地震計の検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別する手段とを具備することを特徴とする。
【００１６】
〔１０〕地震動識別装置において、（ａ）上記〔７〕記載の係数に基づく地震動の識別情報と、（ｂ）上記〔８〕記載のフィッティング残差Ｚに基づく地震動の識別情報と、（ｃ）上記〔９〕記載の加速度エンベロープの最大振幅Ａmax に基づく地震動の識別情報と、（ｄ）上記（ａ）、（ｂ）及び（ｃ）を論理和として地震動を判定する手段とを具備することを特徴とする。
【００１７】
【発明の実施の形態】
以下、本発明の実施の形態について詳細に説明する。
【００１８】
地震波形の時系列データをｖ（ｔ）とする。更に、この時系列データｖ（ｔ）の絶対値をＶ（ｔ）とする。ここで、ｔは時間（秒）を表す。
【００１９】
まず、地震を検知した時刻を時間原点（ｔ＝０）として、そこから、数秒間のデータに簡易な関数を当てはめる。関数形として、ここでは、次式（１）を用いる。なお、絶対値Ｖ（ｔ）の時間的な変動が大きい場合は、適当な方法で絶対値Ｖ（ｔ）を平滑化した後に関数の当てはめを行う。
【００２０】
Ｖ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ） …（１）
ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関するパラメータであり、以下、これらを形状パラメータと仮称する。上記（１）式の両辺の常用対数をとると、次式（２）のように未知パラメータが線形化されるので、その形状パラメータは通常の最小二乗法を用いて算出できる。
【００２１】
ｌｏｇ〔Ｖ（ｔ）／ｔ〕＝ｌｏｇＢ−Ａ・ｌｏｇｅ・ｔ …（２）
上記の方法を用いて多くの地震波形データについて初動部の形状パラメータＡ，Ｂを求める。
【００２２】
図３は実際の地震波に簡易な関数をフィッティングした例を示す図であり、図３（ａ）は実際の初動波（Ｐ波）、図３（ｂ）はその簡易な関数のフィッティング例を示す図である。
【００２３】
このように、多くの地震データに対して、同様のフィッティングを行い、初動部分の形状パラメータを求め、このパラメータを用いて地震動識別処理を行う。
【００２４】
以下、その地震動識別処理の具体例を詳細に説明する。
【００２５】
図１は本発明に係る地震動識別装置のブロック図、図２は本発明に係る地震動識別処理の概略フローチャートである。
【００２６】
これらの図において、１は地震計、２は通信回線、１０は制御処理装置、１１はアンチエイリアシングフィルタ（ローパスフィルタ）、１２はＡ／Ｄ変換器、１３は波形収録部、１４はメモリ、１５はパラメータ演算部、１６は閾値設定部、１７は地震動識別部、１８は条件判定部、１９は情報文送出部、２０はＧＰＳ時計、２１は状態表示装置、２２はモデムである。
【００２７】
この地震動識別方法を、図２に示すフローを参照しながら説明する。
【００２８】
（１）地震計１からの情報を制御処理装置１０に取込み、ディジタル波形データとしてメモリ１４に取得する（ステップＳ１）。
【００２９】
（２）そこで、パラメータ演算部１５において、オフセットの除去を行い（ステップＳ２）、次に、データの絶対値を取得する。このとき絶対値がある基準値以下の場合には、強制的に最小基準値とする（ステップＳ３）。次に、スムージング（例えば、波形のエンベロープを取る）を行う（ステップＳ４）。次に、ｙ＝Ｂｔ・ｅ-At とフィッティングを行い形状パラメータ（係数）Ａ，Ｂを求める（ステップＳ５）。
【００３０】
（３）次いで、閾値設定部１６において、例えば、前記関数の形状パラメータ（係数）Ａ，Ｂに閾値ＴA ，ＴB を設定する（ステップＳ６）。
【００３１】
（４）次に、設定された閾値ＴA ，ＴB に基づいて地震動識別部１７で地震動識別処理を行う（ステップＳ７）。
【００３２】
（５）次に、条件判定部１８において、地震動識別処理された情報に基づき地震動であるか否かの総合判定を行う（ステップＳ８）。
【００３３】
（６）ステップＳ８において、ＹＥＳの場合は、情報文送出部１９において、早期検知情報の送信を行う（ステップＳ９）。
〔具体的な地震動識別方法〕
上記したように、これまでに地震動初動部分（Ｐ波部分）に簡単な関数
Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ） （ｔ：時間）
をフィッティングすることによって、係数（形状パラメータ）Ａ，Ｂを決定し、地震のマグニチュードを地震検知から数秒間で推定する方法を開発している（上記特許文献１参照）。
【００３４】
本発明では、このフィッティングした際に得られる形状パラメータＡ，Ｂやフィッティング残差Ｚに閾値を設定することによって、自動的に地震動の識別を行う。
【００３５】
図４は本発明の第１実施例を示す地震動識別方法の説明図であり、横軸は時間（ｓ）、縦軸は加速度（絶対値）を示す図である。
【００３６】
ここでは、まず、初動部分の最大振幅に関するパラメータである係数Ａによるノイズ識別を行う。
【００３７】
図４において、実線ａとして示すように、係数Ａがある閾値（基準値）ＴA より大きい場合に、ノイズと識別する。
【００３８】
Ａ≧ＴA ⇒ ノイズ
一方、破線ｂとして示すように、係数Ａがある閾値（基準値）ＴA より小さい場合に、地震動と識別する。
【００３９】
Ａ＜ＴA ⇒ 地震動
図５に実際の地震計で観測されたノイズデータに関数をフィッティングして得た係数Ａの頻度分布〔図５（ａ）〕と地震動データに関数をフィッティングして得た係数Ａの頻度分布〔図５（ｂ）〕をそれぞれ示す。
【００４０】
図５（ａ）に示すように、ノイズデータは、係数Ａが１．１をピークに持つ分布になるのに対して、図５（ｂ）に示すように、地震動データは、係数Ａが０．７をピークに持つ分布であり、ノイズデータと地震動データでは係数Ａの分布が異なる。この場合、係数Ａによるノイズ識別の閾値ＴA を１．０とすれば、ノイズデータの７１％をノイズと識別でき、地震動データの９５％を地震動と識別することができる。
【００４１】
図６は本発明の第２実施例を示す地震動識別方法の説明図であり、横軸は時間（ｓ）、縦軸は加速度（絶対値）を示す図である。
【００４２】
ここでは、係数Ｂによるノイズ識別を行う。
【００４３】
図６において、実線ａとして示すように、係数Ｂがある閾値（基準値）ＴB より小さい場合に、ノイズと識別する。
【００４４】
Ｂ≦ＴB ⇒ ノイズ
一方、破線ｂとして示すように、係数Ｂがある閾値（基準値）ＴB より大きい場合に、地震動と識別する。
【００４５】
Ｂ＞ＴB ⇒ 地震動
次いで、本発明の第３実施例を示すフィッティング残差Ｚによるノイズ識別について説明する。
【００４６】
フィッティング残差Ｚがある閾値（基準値）ＴZ より大きい場合に、ノイズと識別する。
【００４７】
Ｚ≧ＴZ ⇒ ノイズ
ただし、Ｚはフィッティングの残差二乗和の平均で、
【００４８】
【数５】
【００４９】
次に、加速度エンベロープの振幅によるノイズの識別について説明する。
【００５０】
図７は本発明の第４実施例を示す地震動識別方法の説明図、図８は本発明の第４実施例を示すノイズ波の識別方法の説明図であり、横軸は時間、縦軸は加速度（絶対値）を示す図である。
【００５１】
図８に示すように、検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax （■）が閾値（基準値）Ｇ1 を超えない場合はノイズと識別する。なお、図７及び図８において、点線（◇）が加速度エンベロープを、実線（■）がその加速度エンベロープの最大振幅Ａmax を示している。すなわち、図７では、一定時間の半ばすぎほどで加速度エンベロープが閾値（基準値）Ｇ1 を超えているので、その最大振幅Ａmax もＧ1 を超えている。つまり、地震波と識別される。
【００５２】
一方、図８に示すような場合には、一定時間内で加速度エンベロープの最大振幅Ａmax が閾値（基準値）Ｇ1 を超えてはいないので、ノイズ波であると識別できる。
【００５３】
このように、地震計による検測時刻から一定時間内に最大振幅Ａmax （■）が閾値Ｇ1 を超える場合には、地震波であると識別する。
【００５４】
図９は本発明の実施例を示す総合的な高精度地震動識別フローチャートである。
【００５５】
（１）まず、パラメータ演算部１５において、関数ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）を地震波初動部分のデータにフィッティングし、係数Ａ，Ｂ、残差Ｚを求める（ステップＳ１１）。
【００５６】
（２）次に、閾値設定部１６において、係数Ａ，Ｂ、残差Ｚに閾値ＴA ，ＴB ，ＴZ を設定する（ステップＳ１２）。
【００５７】
（３）地震動識別部１７で係数Ａによる地震動識別（Ａ＜ＴA ）を行う（ステップＳ１３）。Ａ≧ＴA の場合にはノイズと識別する（ステップＳ１８）。
【００５８】
（４）次に、係数Ｂによる地震動識別（Ｂ＞ＴB ）を行う（ステップＳ１４）。Ｂ≦ＴB の場合にはノイズと識別する（ステップＳ１８）。
【００５９】
（５）次に、フィッティング残差Ｚによる地震動識別（Ｚ＜ＴZ ）を行う（ステップＳ１５）。Ｚ≧ＴZ の場合にはノイズと識別する（ステップＳ１８）。
【００６０】
（６）次に、検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax に基づいて、地震動識別（Ａmax ＞Ｇ1 ）を行う（ステップＳ１６）。Ａmax ≦Ｇ1 の場合には、ノイズと識別する（ステップＳ１８）。
【００６１】
（７）ステップＳ１６においてＹＥＳの場合には、条件判定部１８で地震動であると判断する（ステップＳ１７）。そして、地震動であると判別されると、情報文送出部１９において、早期検知情報の送信を行い、それに対応して措置を講じることができる。
【００６２】
このように、地震動であるかの判定処理を図９に示す処理フローに沿って行うことによって、精度の高い総合的な地震の自動識別処理を行うことができる。
【００６３】
なお、本発明は上記実施例に限定されるものではなく、本発明の趣旨に基づいて種々の変形が可能であり、これらを本発明の範囲から排除するものではない。
【００６４】
【発明の効果】
以上、詳細に説明したように、本発明によれば、以下のような効果を奏することができる。
【００６５】
（Ａ）関数フィッティング法に基づいて、解析対象となる地震動であるかノイズであるかを簡便に、かつ自動的に識別することができる。
【００６６】
（Ｂ）関数フィッティング法による係数Ａが閾値ＴA より小さい場合に地震動であると識別することができる。
【００６７】
（Ｃ）関数フィッティング法による係数Ｂが閾値ＴB より大きい場合に地震動であると識別することができる。
【００６８】
（Ｄ）関数フィッティング法によるフィッティング残差Ｚが閾値ＴZ より小さい場合に地震動であると識別することができる。
【００６９】
（Ｅ）関数フィッティング法による検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別することができる。
【００７０】
（Ｆ）上記（Ｂ）から（Ｅ）をＡＮＤ条件とした場合に、総合的な高精度の地震動の識別を行うことができる。
【図面の簡単な説明】
【図１】 本発明に係る地震動識別装置のブロック図である。
【図２】 本発明に係る地震動識別処理の概略フローチャートである。
【図３】 本発明にかかる簡易な関数のフィッティングの例を示す図である。
【図４】 本発明の第１実施例を示す地震動識別方法の説明図である。
【図５】 （ａ）は実際の地震計で観測されたノイズデータに関数をフィッティングして得た係数Ａの頻度分布を示す図、（ｂ）は地震動データに関数をフィッティングして得た係数Ａの頻度分布を示す図である。
【図６】 本発明の第２実施例を示す地震動識別方法の説明図である。
【図７】 本発明の第４実施例を示す地震動識別方法の説明図である。
【図８】 本発明の第４実施例を示すノイズ波の識別方法の説明図である。
【図９】 本発明の実施例を示す総合的な高精度地震動識別フローチャートである。
【符号の説明】
１ 地震計
２ 通信回線
１０ 制御処理装置
１１ アンチエイリアシングフィルタ（ローパスフィルタ）
１２ Ａ／Ｄ変換器
１３ 波形収録部
１４ メモリ
１５ パラメータ演算部
１６ 閾値設定部
１７ 地震動識別部
１８ 条件判定部
１９ 情報文送出部
２０ ＧＰＳ時計
２１ 状態表示装置
２２ モデム
Claims (10)
- （ａ）地震計から得られる時系列データの絶対値をＶ（ｔ）となし〔ここで、ｔは時間（秒）〕、地震を検知した時刻を時間原点（ｔ＝０）として、そこから、数秒間のデータに簡易な関数を当てはめ、関数形として、Ｖ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）を用い（ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ）、上記式の両辺の常用対数をとって、ｌｏｇ〔Ｖ（ｔ）／ｔ〕＝ｌｏｇＢ−Ａ・ｌｏｇｅ・ｔを得て、未知パラメータを線形化し、その形状パラメータを通常の最小二乗法を用いて算出し、前記初動部分の形状パラメータＡ，Ｂを求め、
（ｂ）該初動部分の形状パラメータＡ，Ｂに閾値ＴA ，ＴB を設定し、
（ｃ）該設定された閾値ＴA ，ＴB を有するフィッティングしたデータに基づいて、地震動であるか、それ以外のノイズであるかを自動的に識別することを特徴とする地震動識別方法。 - 請求項１記載の地震動識別方法において、前記関数のパラメータが係数Ａであり、該係数Ａが閾値ＴA より小さい場合に地震動であると識別することを特徴とする地震動識別方法。
- 請求項１記載の地震動識別方法において、前記関数のパラメータが係数Ｂであり、該係数Ｂが閾値ＴB より大きい場合に地震動であると識別することを特徴とする地震動識別方法。
- （ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化し、
（ｂ）これによって得られた前記関数Ｙのフィッティング残差Ｚ
（ｃ）該設定されたフィッティング残差Ｚが閾値ＴZ より小さい場合に地震動であると識別することを特徴とする地震動識別方法。 - （ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化し、
（ｂ）前記関数フィッティングによる、前記地震計の検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別することを特徴とする地震動識別方法。 - （ａ）上記請求項２記載の係数Ａが閾値ＴA より小さく、
（ｂ）上記請求項３記載の係数Ｂが閾値ＴB より大きく、
（ｃ）上記請求項４記載のフィッティング残差Ｚが閾値ＴZ より小さく、
（ｄ）上記請求項５記載の最大振幅Ａmax が閾値Ｇ1 を超えるときに地震動であると総合的に識別することを特徴とする地震動識別方法。 - （ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、
（ｂ）前記初動部分の波形形状パラメータ（係数）Ａ，Ｂを求める手段と、
（ｃ）該初動部分の波形形状パラメータ（係数）Ａ，Ｂに閾値ＴA ，ＴB を設定する手段と、
（ｄ）該設定された閾値ＴA ，ＴB を有するフィッティングしたデータに基づいて地震動であるか、それ以外のノイズであるかを自動的に識別する手段とを具備することを特徴とする地震動識別装置。 - （ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、
（ｂ）これによって得られた前記関数Ｙのフィッティング残差Ｚ
（ｃ）該フィッティング残差Ｚが、該フィッティング残差Ｚの閾値ＴZ より小さい場合に地震動であると識別する手段とを具備することを特徴とする地震動識別装置。 - （ａ）地震計から得られる地震波初動部分の波形形状をパラメータが数個の簡易な関数Ｙ（ｔ）＝Ｂ・ｔ×ｅｘｐ（−Ａ・ｔ）〔ここで、Ｂは地震波形の初動振幅の時間変化に関するパラメータ、Ａは初動部分の最大振幅に関係するパラメータ、ｔは時間（秒）〕でフィッティングすることで定量化する手段と、
（ｂ）前記関数フィッティングによる、前記地震計の検測時刻から一定時間内の加速度エンベロープの最大振幅Ａmax が閾値Ｇ1 を超える場合に地震動であると識別する手段とを具備することを特徴とする地震動識別装置。 - （ａ）上記請求項７記載の係数に基づく地震動の識別情報と、
（ｂ）上記請求項８記載のフィッティング残差Ｚに基づく地震動の識別情報と、
（ｃ）上記請求項９記載の加速度エンベロープの最大振幅Ａmax に基づく地震動の識別情報と、
（ｄ）上記（ａ）、（ｂ）及び（ｃ）を論理和として地震動を判定する手段とを具備することを特徴とする地震動識別装置。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003175176A JP4173408B2 (ja) | 2003-06-19 | 2003-06-19 | 地震動識別方法及びそのための装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP2003175176A JP4173408B2 (ja) | 2003-06-19 | 2003-06-19 | 地震動識別方法及びそのための装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JP2005010041A JP2005010041A (ja) | 2005-01-13 |
| JP4173408B2 true JP4173408B2 (ja) | 2008-10-29 |
Family
ID=34098456
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP2003175176A Expired - Lifetime JP4173408B2 (ja) | 2003-06-19 | 2003-06-19 | 地震動識別方法及びそのための装置 |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP4173408B2 (ja) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2012083134A (ja) * | 2010-10-07 | 2012-04-26 | Railway Technical Research Institute | 即時応答性を向上させた震央距離推定方法 |
| JP6189226B2 (ja) * | 2014-02-07 | 2017-08-30 | 戸田建設株式会社 | 地震情報配信システムとノイズ判定方法 |
-
2003
- 2003-06-19 JP JP2003175176A patent/JP4173408B2/ja not_active Expired - Lifetime
Also Published As
| Publication number | Publication date |
|---|---|
| JP2005010041A (ja) | 2005-01-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Bai et al. | Comparison of alternatives to amplitude thresholding for onset detection of acoustic emission signals | |
| KR102027252B1 (ko) | 임의 지진구역 대상 무작위 생성 지진학습데이터를 이용한 지진여부판별 방법 및 이에 의한 진도예측 시스템 | |
| CA2892698A1 (en) | Method and apparatus for estimating the fall risk of a user | |
| US20150308920A1 (en) | Adaptive baseline damage detection system and method | |
| US7609874B2 (en) | System and method for prediction of pitting corrosion growth | |
| JP2020052740A (ja) | 異常検知装置、異常検知方法及びプログラム | |
| JP7052869B2 (ja) | 漏えい調査装置、漏えい調査方法、およびプログラム | |
| US20240053302A1 (en) | Ultrasonic non-destructive test method and system using deep learning, and auto-encoder-based prediction model training method used therefor | |
| CN114878582B (zh) | 一种特钢的缺陷检测分析方法及系统 | |
| JPWO2018008708A1 (ja) | 震央距離推定装置、震央距離推定方法、及びプログラム | |
| JP3695579B2 (ja) | 震央距離及びマグニチュード推定方法とそのための装置 | |
| CN109557583B (zh) | 一种地震属性提取方法及系统 | |
| JP6777224B2 (ja) | 損傷検出装置、方法およびプログラム | |
| CN110455923B (zh) | 一种锚杆锚固质量等级快速评价方法 | |
| JP4173408B2 (ja) | 地震動識別方法及びそのための装置 | |
| JP7375692B2 (ja) | 情報処理装置、情報処理方法、及び情報処理システム | |
| JP4617168B2 (ja) | 軸受損傷評価装置および軸受損傷評価方法 | |
| CN109870404B (zh) | 一种雨棚结构损伤识别方法、装置及终端设备 | |
| JP2011523453A (ja) | 測定信号の相対ノイズ | |
| US6281814B1 (en) | Data conversion method, data converter, and program storage medium | |
| JP5351835B2 (ja) | 音信号区間抽出装置及び音信号区間抽出方法 | |
| CN108985350B (zh) | 一种基于梯度幅值稀疏特征信息识别模糊图像的方法和装置、计算设备及存储介质 | |
| JP2009041978A (ja) | 打音解析による健全性診断方法 | |
| CN111971581A (zh) | 用于验证由雨量传感器提供的数据的设备、方法和计算机程序产品 | |
| JP6832005B2 (ja) | 被験者判断装置、方法、及びプログラム |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050627 |
|
| A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20070202 |
|
| A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080624 |
|
| A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080714 |
|
| TRDD | Decision of grant or rejection written | ||
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 Effective date: 20080812 |
|
| A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20080813 |
|
| R150 | Certificate of patent or registration of utility model |
Ref document number: 4173408 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110822 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110822 Year of fee payment: 3 |
|
| S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110822 Year of fee payment: 3 |
|
| R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20110822 Year of fee payment: 3 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120822 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20120822 Year of fee payment: 4 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20130822 Year of fee payment: 5 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| EXPY | Cancellation because of completion of term |