JP2022165187A - 損傷領域推定システム、推定装置及び損傷領域推定方法 - Google Patents
損傷領域推定システム、推定装置及び損傷領域推定方法 Download PDFInfo
- Publication number
- JP2022165187A JP2022165187A JP2021070431A JP2021070431A JP2022165187A JP 2022165187 A JP2022165187 A JP 2022165187A JP 2021070431 A JP2021070431 A JP 2021070431A JP 2021070431 A JP2021070431 A JP 2021070431A JP 2022165187 A JP2022165187 A JP 2022165187A
- Authority
- JP
- Japan
- Prior art keywords
- unit
- elastic wave
- estimation system
- damaged area
- rail
- 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.)
- Pending
Links
- 238000000034 method Methods 0.000 title claims abstract description 28
- 238000012545 processing Methods 0.000 claims description 68
- 238000004891 communication Methods 0.000 claims description 53
- 230000007797 corrosion Effects 0.000 claims description 48
- 238000005260 corrosion Methods 0.000 claims description 48
- 230000006378 damage Effects 0.000 claims description 47
- 238000001514 detection method Methods 0.000 claims description 23
- 239000000284 extract Substances 0.000 claims description 6
- 230000001960 triggered effect Effects 0.000 claims 1
- 230000005540 biological transmission Effects 0.000 description 45
- 238000000605 extraction Methods 0.000 description 25
- 238000010586 diagram Methods 0.000 description 23
- 230000006870 function Effects 0.000 description 13
- 230000008030 elimination Effects 0.000 description 9
- 238000003379 elimination reaction Methods 0.000 description 9
- 238000007493 shaping process Methods 0.000 description 9
- 238000012360 testing method Methods 0.000 description 8
- 230000001133 acceleration Effects 0.000 description 6
- 238000009434 installation Methods 0.000 description 6
- 238000007689 inspection Methods 0.000 description 5
- 230000035945 sensitivity Effects 0.000 description 5
- 238000012795 verification Methods 0.000 description 5
- 230000005856 abnormality Effects 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 238000007639 printing Methods 0.000 description 4
- 238000007796 conventional method Methods 0.000 description 3
- 230000006866 deterioration Effects 0.000 description 3
- 238000005401 electroluminescence Methods 0.000 description 2
- 230000003628 erosive effect Effects 0.000 description 2
- 238000011156 evaluation Methods 0.000 description 2
- 238000004088 simulation Methods 0.000 description 2
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 2
- 241001669679 Eleotris Species 0.000 description 1
- 239000000853 adhesive Substances 0.000 description 1
- 230000001070 adhesive effect Effects 0.000 description 1
- 238000005452 bending Methods 0.000 description 1
- 239000013078 crystal Substances 0.000 description 1
- 230000001186 cumulative effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 230000009977 dual effect Effects 0.000 description 1
- 238000003487 electrochemical reaction Methods 0.000 description 1
- 238000001914 filtration Methods 0.000 description 1
- 230000010354 integration Effects 0.000 description 1
- 230000001678 irradiating effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000004807 localization Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 238000010248 power generation Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 150000003839 salts Chemical class 0.000 description 1
- 239000000523 sample Substances 0.000 description 1
- 239000004065 semiconductor Substances 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/043—Analysing solids in the interior, e.g. by shear waves
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61K—AUXILIARY EQUIPMENT SPECIALLY ADAPTED FOR RAILWAYS, NOT OTHERWISE PROVIDED FOR
- B61K9/00—Railway vehicle profile gauges; Detecting or indicating overheating of components; Apparatus on locomotives or cars to indicate bad track sections; General design of track recording vehicles
- B61K9/08—Measuring installations for surveying permanent way
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
- B61L23/041—Obstacle detection
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
- B61L23/042—Track changes detection
- B61L23/044—Broken rails
-
- B—PERFORMING OPERATIONS; TRANSPORTING
- B61—RAILWAYS
- B61L—GUIDING RAILWAY TRAFFIC; ENSURING THE SAFETY OF RAILWAY TRAFFIC
- B61L23/00—Control, warning or like safety means along the route or between vehicles or trains
- B61L23/04—Control, warning or like safety means along the route or between vehicles or trains for monitoring the mechanical state of the route
- B61L23/042—Track changes detection
- B61L23/045—Rail wear
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/04—Analysing solids
- G01N29/12—Analysing solids by measuring frequency or resonance of acoustic waves
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/14—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object using acoustic emission techniques
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N29/00—Investigating or analysing materials by the use of ultrasonic, sonic or infrasonic waves; Visualisation of the interior of objects by transmitting ultrasonic or sonic waves through the object
- G01N29/44—Processing the detected response signal, e.g. electronic circuits specially adapted therefor
- G01N29/4409—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison
- G01N29/4436—Processing the detected response signal, e.g. electronic circuits specially adapted therefor by comparison with a reference signal
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/02—Indexing codes associated with the analysed material
- G01N2291/028—Material parameters
- G01N2291/0289—Internal structure, e.g. defects, grain size, texture
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01N—INVESTIGATING OR ANALYSING MATERIALS BY DETERMINING THEIR CHEMICAL OR PHYSICAL PROPERTIES
- G01N2291/00—Indexing codes associated with group G01N29/00
- G01N2291/26—Scanned objects
- G01N2291/262—Linear objects
- G01N2291/2623—Rails; Railroads
Landscapes
- Engineering & Computer Science (AREA)
- Physics & Mathematics (AREA)
- Mechanical Engineering (AREA)
- General Health & Medical Sciences (AREA)
- Immunology (AREA)
- Chemical & Material Sciences (AREA)
- Analytical Chemistry (AREA)
- Biochemistry (AREA)
- Health & Medical Sciences (AREA)
- General Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Pathology (AREA)
- Acoustics & Sound (AREA)
- Signal Processing (AREA)
- Investigating Or Analyzing Materials By The Use Of Ultrasonic Waves (AREA)
- Machines For Laying And Maintaining Railways (AREA)
- Train Traffic Observation, Control, And Security (AREA)
Abstract
【課題】破断が発生する前段階で鉄道に関する対象物の損傷領域を簡便に推定することができる損傷領域推定システム、推定装置及び損傷領域推定方法を提供することである。【解決手段】実施形態の損傷領域推定システムは、複数のセンサと、位置標定部と、推定部とを持つ。複数のセンサは、損傷領域の推定対象となる鉄道に関する対象物で発生する弾性波を検出する。位置標定部は、前記複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する。推定部は、前記複数の弾性波の発生源の位置に基づいて、前記対象物における損傷領域を推定する。【選択図】図2
Description
本発明の実施形態は、損傷領域推定システム、推定装置及び損傷領域推定方法に関する。
従来、レール異常に伴う鉄道の輸送トラブルが発生している。レール異常の中でも、レール破断は脱線等の重大事故を引き起こす危険性が高いため、早急な対応が求められている。レール破断を検知する技術も提案されているが、従来の手法では、レールが破断したことを検知するため、既にレールの損傷が大きくなっていることが想定される。
損傷の初期段階で異常を捉えることができればより効果的な対策が期待でき、重大事故防止の観点からも望ましい。損傷の初期段階で異常を捉えることができる従来技術も提案されている。しかしながら、従来の手法では、異常を検知できない場所があったり、異常を検知するために数多くのセンサを設置する必要があり、損傷領域を簡便に推定することができない場合があった。なお、このような問題は、レールに限らず、架線等においても生じる。
本発明が解決しようとする課題は、破断が発生する前段階で鉄道に関する対象物の損傷領域を簡便に推定することができる損傷領域推定システム、推定装置及び損傷領域推定方法を提供することである。
実施形態の損傷領域推定システムは、複数のセンサと、位置標定部と、推定部とを持つ。複数のセンサは、損傷領域の推定対象となる鉄道に関する対象物で発生する弾性波を検出する。位置標定部は、前記複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する。推定部は、前記複数の弾性波の発生源の位置に基づいて、前記対象物における損傷領域を推定する。
以下、実施形態の損傷領域推定システム、推定装置及び損傷領域推定方法を、図面を参照して説明する。
(概要)
実施形態における損傷領域推定システムは、軌条で発生した弾性波に基づいて、軌条における損傷領域を推定する。軌条とは、例えば鉄道車両が走行するレールである。軌条において弾性波が発生する現象は、例えば図1に示すように、軌条を流れる電流が漏出して電食が発生することによって生じる。
実施形態における損傷領域推定システムは、軌条で発生した弾性波に基づいて、軌条における損傷領域を推定する。軌条とは、例えば鉄道車両が走行するレールである。軌条において弾性波が発生する現象は、例えば図1に示すように、軌条を流れる電流が漏出して電食が発生することによって生じる。
図1は、実施形態における損傷領域推定システムの概要を説明するための図である。図1には、鉄道車両RSがレールR上を走行している様子が示されている。鉄道車両RSは、変電所1から供給される電流を利用して走行する。より具体的には、変電所1から出力される電流は、電車線2及び鉄道車両RSが備えるパンタグラフ3を介して鉄道車両RSに供給される。鉄道車両RSで利用されなかった電流は、レールRを介して変電所1に戻る。電流がレールRを介して変電所1に戻る際、一部の電流は漏れ電流として地面4に流れる。このような現象により、レールRにおいて電食が発生する。図1では、電食が発生した箇所を電食部5として示している。
電食部5では、弾性波AEが発生する。電食部5で発生した弾性波AEは、レールR内部を伝搬していく。そこで、実施形態における損傷領域推定システムでは、レールR又はレールRに接触している物体に複数のセンサ10を設置して、レールR内部を伝搬した弾性波AEを複数のセンサ10で検出する。レールRに接触している物体は、例えばレールRをまくらぎに固定するためのレールクリップや鉄道車両RSの車輪などである。損傷領域推定システムでは、複数のセンサ10それぞれによって検出された弾性波AEに基づいて位置標定を行う。これにより、弾性波AEの発生源(以下「弾性波源」という)が標定される。そして、損傷領域推定システムでは、位置標定の結果を用いて、レールRにおける損傷領域を推定する。損傷領域とは、電食によりレールRにおいて損傷が発生している領域である。
上記のような手法で損傷領域を推定することにより、従来のようにレールが破断する前段階でレールの異常を発見することができる。
以下、実施形態における損傷領域推定システムの詳細について説明する。
以下、実施形態における損傷領域推定システムの詳細について説明する。
(第1の実施形態)
図2は、第1の実施形態における損傷領域推定システム100の構成を表す図である。損傷領域推定システム100は、複数のセンサ10-1~10-n(nは1以上の整数)、信号処理部20及び推定装置30を備える。第1の実施形態において、センサ10と信号処理部20との間、及び、信号処理部20と推定装置30との間は、有線により接続される。なお、以下の説明では、センサ10-1~10-nを区別しない場合にはセンサ10と記載する。第1の実施形態において、センサ10及び信号処理部20は、軌条で発生する弾性波を検出する検出装置として構成される。
図2は、第1の実施形態における損傷領域推定システム100の構成を表す図である。損傷領域推定システム100は、複数のセンサ10-1~10-n(nは1以上の整数)、信号処理部20及び推定装置30を備える。第1の実施形態において、センサ10と信号処理部20との間、及び、信号処理部20と推定装置30との間は、有線により接続される。なお、以下の説明では、センサ10-1~10-nを区別しない場合にはセンサ10と記載する。第1の実施形態において、センサ10及び信号処理部20は、軌条で発生する弾性波を検出する検出装置として構成される。
センサ10は、レールR又はレールRに接触している物体に設置される。以下の説明では、センサ10が、レールRに設置されているものとする。信号処理部20及び推定装置30は、センサ10と同じ場所に設けられてもよいし、センサ10が設置されている場所と異なる場所(例えば、管理者がいる管理所等)に設けられてもよい。
センサ10は、鉄道車両RSが走行する軌条で発生する弾性波を検出する。センサ10は、検出した弾性波を電気信号に変換する。センサ10には、例えば10kHz~1MHzの範囲に感度を有する圧電素子が用いられる。センサ10としてより好適なものは、100kHz~200kHzに感度を有する圧電素子である。センサ10は、周波数範囲内に共振ピークをもつ共振型、共振を抑えた広帯域型等の種類があるが、センサ10の種類はいずれでもよい。センサ10が弾性波を検出する方法は、電圧出力型、抵抗変化型及び静電容量型等があるが、いずれの検出方法でもよい。センサ10は、増幅器を内蔵していてもよい。
センサ10に代えて加速度センサ又はマイクロホンが用いられてもよい。この場合、加速度センサは、レールR内部で発生する弾性波を検出する。そして、加速度センサは、センサ10と同様の処理を行うことによって、検出した弾性波を電気信号に変換する。
信号処理部20は、センサ10から出力された電気信号を入力とする。信号処理部20は、入力した電気信号に対して信号処理を行う。信号処理部20が行う信号処理は、例えば、ノイズ除去、到達時刻の決定、パラメータ抽出等である。信号処理部20は、信号処理により得られた弾性波の特徴量のデータを送信データとして、有線ケーブル40を介して推定装置30に出力する。
信号処理部20は、アナログ回路又はデジタル回路を用いて構成される。デジタル回路は、例えばFPGA(Field Programmable Gate Array)やマイクロコンピュータにより実現される。不揮発型のFPGAを用いることで、待機時の消費電力を抑えることができる。デジタル回路は、専用のLSI(Large-Scale Integration)により実現されてもよい。信号処理部20は、フラッシュメモリ等の不揮発メモリや、取り外し可能なメモリを搭載してもよい。
図3は、第1の実施形態における信号処理部20の機能を表す概略ブロック図である。信号処理部20は、フィルタ21、A/D変換器22、波形整形フィルタ23、ゲート生成回路24、到達時刻決定部25、特徴量抽出部26、送信データ生成部27、メモリ28及び出力部29を備える。
フィルタ21では、センサ10から出力された電気信号において信号帯域以外のノイズ成分が除去される。フィルタ21は、例えばバンドパスフィルタ(BPF:Band pass filter)である。
A/D変換器22は、ノイズ成分が除去された電気信号を量子化してデジタル信号に変換する。A/D変換器22は、デジタル信号を波形整形フィルタ23に出力する。
波形整形フィルタ23は、入力された時系列データのデジタル信号から所定の信号帯域外のノイズ成分を除去する。波形整形フィルタ23は、例えばバンドパスフィルタ(BPF)である。波形整形フィルタ23は、例えばフィルタ21と同じ周波数帯域を通過させるように設定されているものとする。波形整形フィルタ23は、ノイズ成分除去後の信号(以下「ノイズ除去信号」という。)をゲート生成回路24及び特徴量抽出部26に出力する。
ゲート生成回路24は、波形整形フィルタ23から出力されたノイズ除去信号を入力とする。ゲート生成回路24は、入力したノイズ除去信号に基づいてゲート信号を生成する。ゲート信号は、ノイズ除去信号の波形が持続しているか否かを示す信号である。
ゲート生成回路24は、例えばエンベロープ検出器及びコンパレータにより実現される。エンベロープ検出器は、ノイズ除去信号のエンベロープを検出する。エンベロープは、例えばノイズ除去信号を二乗し、二乗した出力値に対して所定の処理(例えばローパスフィルタを用いた処理やヒルベルト変換)を行うことで抽出される。コンパレータは、ノイズ除去信号のエンベロープが所定の閾値以上であるか否かを判定する。
ゲート生成回路24は、ノイズ除去信号のエンベロープが所定の閾値以上となった場合、ノイズ除去信号の波形が持続していることを示す第1のゲート信号を到達時刻決定部25及び特徴量抽出部26に出力する。一方、ゲート生成回路24は、ノイズ除去信号のエンベロープが所定の閾値未満になった場合、ノイズ除去信号の波形が持続していないことを示す第2のゲート信号を到達時刻決定部25及び特徴量抽出部26に出力する。
到達時刻決定部25は、不図示の水晶発振器などのクロック源から出力されるクロックと、ゲート生成回路24から出力されたゲート信号とを入力とする。到達時刻決定部25は、第1のゲート信号が入力されている間に入力されたクロックを用いて、弾性波到達時刻を決定する。到達時刻決定部25は、決定した弾性波到達時刻を時刻情報として送信データ生成部27に出力する。到達時刻決定部25は、第2のゲート信号が入力されている間に処理を行わない。到達時刻決定部25は、クロック源からの信号をもとに、電源投入時からの累積の時刻情報を生成する。具体的には、到達時刻決定部25は、クロックのエッジをカウントするカウンタとし、カウンタのレジスタの値を時刻情報とすればよい。カウンタのレジスタは所定のビット長を有するように決定される。
特徴量抽出部26は、波形整形フィルタ23から出力されたノイズ除去信号と、ゲート生成回路24から出力されたゲート信号とを入力とする。特徴量抽出部26は、第1のゲート信号が入力されている間に入力されたノイズ除去信号を用いて、ノイズ除去信号の特徴量を抽出する。特徴量抽出部26は、第2のゲート信号が入力されている間に処理を行わない。特徴量は、ノイズ除去信号の特徴を示す情報である。
特徴量は、例えば波形の振幅[mV]、波形の立ち上がり時間[usec]、ゲート信号の持続時間[usec]、ゼロクロスカウント数[times]、波形のエネルギー[arb.]、周波数[Hz]及びRMS(Root Mean Square:二乗平均平方根)値等である。特徴量抽出部26は、抽出した特徴量に関するパラメータを送信データ生成部27に出力する。特徴量抽出部26は、特徴量に関するパラメータを出力する際に、特徴量に関するパラメータにセンサIDを対応付ける。センサIDは、レールRに設置されているセンサ10を識別するための識別情報を表す。これにより、特徴量に関するパラメータが、どのセンサ10により検出された弾性波の特徴量であるのかを特定することができる。
波形の振幅は、例えばノイズ除去信号の中で最大振幅の値である。波形の立ち上がり時間は、例えばゲート信号の立ち上がり開始からノイズ除去信号が最大値に達するまでの時間T1である。ゲート信号の持続時間は、例えばゲート信号の立ち上がり開始から振幅が予め設定される値よりも小さくなるまでの時間である。ゼロクロスカウント数は、例えばゼロ値を通る基準線をノイズ除去信号が横切る回数である。
波形のエネルギーは、例えば各時点においてノイズ除去信号の振幅を二乗したものを時間積分した値である。なお、エネルギーの定義は、上記例に限定されず、例えば波形の包絡線を用いて近似されたものでもよい。周波数は、ノイズ除去信号の周波数である。RMS値は、例えば各時点においてノイズ除去信号の振幅を二乗して平方根により求めた値である。
送信データ生成部27は、センサIDと、時刻情報と、特徴量に関するパラメータとを入力とする。送信データ生成部27は、入力したセンサIDと、時刻情報と、特徴量に関するパラメータとを含む送信データを生成する。送信データ生成部27は、生成した送信データをメモリ28に記録してもよいし、メモリ28に記録せずに出力部29に出力してもよい。
メモリ28は、送信データを記憶する。メモリ28は、例えばデュアルポートRAM(Random Access Memory)である。
出力部29は、メモリ28に記憶されている送信データ、又は、送信データ生成部27から出力された送信データを推定装置30に逐次出力する。
図2に戻って説明を続ける。
推定装置30は、有線通信部31、制御部32、記憶部33及び表示部34を備える。
推定装置30は、有線通信部31、制御部32、記憶部33及び表示部34を備える。
有線通信部31は、信号処理部20から出力された送信データを、有線ケーブル40を介して受信する。有線通信部31は、受信した送信データを制御部32に出力する。
制御部32は、推定装置30全体を制御する。制御部32は、CPU(Central Processing Unit)等のプロセッサやメモリを用いて構成される。制御部32は、プログラムを実行することによって、取得部321、イベント抽出部322、位置標定部323、分布生成部324及び推定部325として機能する。
取得部321、イベント抽出部322、位置標定部323、分布生成部324及び推定部325の機能部のうち一部または全部は、ASIC(Application Specific Integrated Circuit)やPLD(Programmable Logic Device)、FPGAなどのハードウェア(回路部;circuitryを含む)によって実現されてもよいし、ソフトウェアとハードウェアとの協働によって実現されてもよい。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM、CD-ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置などの非一時的な記憶媒体である。プログラムは、電気通信回線を介して送信されてもよい。
取得部321、イベント抽出部322、位置標定部323、分布生成部324及び推定部325の機能の一部は、予め推定装置30に搭載されている必要はなく、追加のアプリケーションプログラムが推定装置30にインストールされることで実現されてもよい。
取得部321は、各種情報を取得する。例えば、取得部321は、有線通信部31によって受信された送信データを取得する。取得部321は、取得した送信データを記憶部33に保存する。
イベント抽出部322は、記憶部33に記憶されている送信データの中から1イベントにおける送信データを抽出する。イベントとは、レールRで発生した弾性波発生事象を表す。本実施形態における弾性波発生事象は、例えばレールRで発生した電食である。1回のイベントが発生した場合、複数のセンサ10で略同時刻に弾性波が検出されることになる。すなわち、記憶部33には、略同時刻に検出された弾性波に関する送信データが記憶されていることになる。そこで、イベント抽出部322は、所定の時間窓を設け、到達時刻が時間窓の範囲内に存在する全ての送信データを1イベントにおける送信データとして抽出する。イベント抽出部322は、抽出した1イベントにおける送信データを位置標定部323に出力する。
時間窓の範囲Twは、対象とするレールRにおける弾性波伝搬速度vと、最大のセンサ間隔dmaxを用いて、Tw≧dmax/vの範囲になるように決定してもよい。誤検出を避けるためには、Twをできるだけ小さい値に設定することが望ましいため、実質的にはTw=dmax/vとすることができる。弾性波伝搬速度vは、予め求められていてもよい。
位置標定部323は、センサ位置情報と、イベント抽出部322によって抽出された複数の送信データそれぞれに含まれるセンサID及び時刻情報とに基づいて弾性波源の位置標定を行う。なお、弾性波は、例えば振幅やエネルギーが所定値以下のデータを除去するなどしてもよい。
センサ位置情報には、センサIDに対応付けてセンサ10の設置位置に関する情報が含まれる。センサ位置情報は、例えばレールRの特定位置からの水平方向および垂直方向の距離など等の情報を含む。位置標定部323は、センサ位置情報を予め保持している。センサ位置情報は、位置標定部323が弾性波源の位置標定を行う前であればどのタイミングで位置標定部323に記憶されてもよい。
センサ位置情報は、記憶部33に記憶されていてもよい。この場合、位置標定部323は、位置標定を行うタイミングで記憶部33からセンサ位置情報を取得する。弾性波源の位置の標定には、カルマンフィルタ、最小二乗法などが用いられてもよい。位置標定部323は、計測期間中に得られた弾性波源の位置情報を分布生成部324に出力する。
分布生成部324は、位置標定部323から出力された複数の弾性波源の位置情報を入力とする。分布生成部324は、入力した複数の弾性波源の位置情報を用いて、弾性波源分布を生成する。弾性波源分布は、弾性波源の位置が示された分布を表す。より具体的には、弾性波源分布は、損傷領域の推定対象となるレールRを表した仮想的なデータ上において弾性波源の位置を示す点が示された分布である。なお、分布生成部324は、弾性波源分布を用いて、弾性波源の位置をコンター図で表すことによって弾性波源密度分布を生成してもよい。
推定部325は、分布生成部324によって生成された弾性波源分布に基づいてレールRにおける損傷領域を推定する。具体的には、推定部325は、弾性波源の密度が閾値以上の領域を損傷領域として推定する。
記憶部33は、取得部321によって取得された送信データを記憶する。記憶部33は、磁気ハードディスク装置や半導体記憶装置などの記憶装置を用いて構成される。
表示部34は、推定部325の制御に従って評価結果を表示する。表示部34は、液晶ディスプレイ、有機EL(Electro Luminescence)ディスプレイ等の画像表示装置である。表示部34は、画像表示装置を推定装置30に接続するためのインタフェースであってもよい。この場合、表示部34は、評価結果を表示するための映像信号を生成し、自身に接続されている画像表示装置に映像信号を出力する。
図4は、第1の実施形態における損傷領域推定システム100の処理の流れを示すシーケンス図である。なお、図4の処理開始時には、レールRに電流が流れているものとする。
各センサ10は、弾性波を検出する(ステップS101)。例えば、各センサ10は、電食により発生した弾性波を検出する。各センサ10は、検出した弾性波を電気信号に変換して信号処理部20に出力する。信号処理部20は、各センサ10から出力された各電気信号に対して信号処理を行う(ステップS102)。具体的には、信号処理部20は、各電気信号に対してノイズ除去、到達時刻の決定、パラメータ抽出等の信号処理を行う。信号処理部20は、上記の信号処理をセンサ10から電気信号が得られる度に実行する。
各センサ10は、弾性波を検出する(ステップS101)。例えば、各センサ10は、電食により発生した弾性波を検出する。各センサ10は、検出した弾性波を電気信号に変換して信号処理部20に出力する。信号処理部20は、各センサ10から出力された各電気信号に対して信号処理を行う(ステップS102)。具体的には、信号処理部20は、各電気信号に対してノイズ除去、到達時刻の決定、パラメータ抽出等の信号処理を行う。信号処理部20は、上記の信号処理をセンサ10から電気信号が得られる度に実行する。
信号処理部20は、信号処理後のデータを用いてセンサ10毎の送信データを生成する(ステップS103)。信号処理部20は、生成したセンサ10毎の送信データを、有線ケーブル40を介して推定装置30に送信する(ステップS104)。ここで、信号処理部20は、送信データを生成する度に推定装置30に送信してもよいし、ある期間分の送信データを生成したタイミングでまとめて推定装置30に送信してもよい。
推定装置30の有線通信部31は、信号処理部20から送信された送信データを受信する。取得部321は、有線通信部31によって受信された送信データを取得する。取得部321は、取得した送信データを記憶部33に記録する(ステップS105)。イベント抽出部322は、記憶部33に記憶されている1イベントにおける複数の送信データを抽出する(ステップS106)。イベント抽出部322は、抽出した1イベントにおける送信データを位置標定部323に出力する。
位置標定部323は、イベント抽出部322によって抽出された1イベントにおける複数の送信データそれぞれに含まれるセンサID、時刻情報と、センサ位置情報とに基づいて、弾性波源の位置を標定する(ステップS107)。位置標定部323は、弾性波源の標定結果を分布生成部324に出力する。イベント抽出部322及び位置標定部323は、ステップS106及びステップS107の処理を所定期間繰り返し実行する。
分布生成部324は、所定期間分の弾性波源の位置情報を用いて、弾性波源分布を生成する(ステップS108)。分布生成部324は、生成した弾性波源分布を推定部325に出力する。推定部325は、分布生成部324によって生成された弾性波源分布に基づいてレールRにおける損傷領域を推定する(ステップS109)。具体的には、推定部325は、弾性波源の密度が閾値以上の領域を損傷領域として推定する。一例として、推定部325は、弾性波源分布上をメッシュで区切り、区切られた領域毎に弾性波源の密度を算出する。そして、推定部325は、弾性波源の密度が閾値以上の領域を損傷領域として推定してもよい。なお、弾性波源分布がコンター図で表されている場合には、ある閾値以上の領域を損傷領域として推定してもよい。
推定部325は、推定結果を出力する(ステップS110)。例えば、推定部325は、推定結果を表示部34に表示させてもよい。
次に、図5から図8を用いて、レール模擬試験体による原理検証試験の結果をもとに、レールRの健全性の評価方法を詳しく説明する。ここで、レール模擬試験体とは図5に示す実レールを模擬した実験モデルである。図5に示す実験モデルは、レールRと、コンクリートブロックCBと、レールクリップRCと、ゴムシートGSとで構成される。コンクリートブロックCBは、まくらぎに見立てたブロックである。レールクリップRCは、レールRをコンクリートブロックCBに固定するための部材であり、ボルトBと、ナットNとにより固定される。ゴムシートGSは、レールR及びレールクリップRCと、コンクリートブロックCBとを機械的又は電気的に絶縁するための部材である。コンクリートブロックCBと、ボルトBとの間には接着剤Aが流入されている。
実施形態における検査手順は以下の通りである。第一に、検査対象となるレールRに複数のセンサ10を設置する。第二に、レールRに所定の電流を印加する。第三に、損傷(例えば、電食部)から発生した弾性波をセンサ10で検出する。第四に、センサ10で検出した弾性波を分析して弾性波源位置を標定する。以降、第二から第四までの工程を繰り返し、損傷領域を推定する。
以下、図6を参照しながら、実験モデルにて実施した原理検証試験の概要を説明する。本試験では、レールRに4個のセンサ10-1~10-4を配置した。図6では、手前に設置したセンサ10-2及び10-4を示している。センサ10-1及び10-3は、センサ10-2及び10-4と対面する位置に設置される。ただし、センサの設置箇所及び個数はこれに限らず、レールクリップRCやまくらぎなどの電食により発生した弾性波を検出可能な場所及び個数(複数)であればよい。
続いて、定電流電源にてレールR(正極)と、レールクリップRC(負極)との間に直流電流を一定時間印加した。ここで、直流電流とした理由は、直流電気鉄道(検測車含む)から変電所間に流れる帰線(一般的にはレールR)電流の漏れ電流を模擬するためである。なお、交流電気鉄道を想定して交流電流を印加してもよい。さらに、現場の湿潤環境(トンネル内の漏水箇所など)を模擬するため、塩水(濃度:3%)を含ませたスポンジSPをレールRとレールクリップRCで挟んだ。
図7は、上記の原理検証試験により得られた1つの弾性波の波形の一例を示す図である。図7の左に示す図は、横軸を時間、縦軸を振幅とした弾性波の波形を表し、図7の右に示す図は、横軸を周波数、縦軸を振幅とした弾性波の波形を表す。図7の右に示すように、電食により発生した弾性波は、低域(例えば、100kHz未満)に主共振周波数(図7の破線6)を持っていることがわかる。
各センサ10によって検出された弾性波を分析(例えば、2次元位置標定、信号振幅分類)した結果を図8に示す。図8では、振幅が40[dB]以上の弾性波源を示している。図8に示す例では、複数の電食部5の位置が標定されており、円7で示す領域が損傷領域として推定される。このように、レールRに設置したセンサ10-1~10-4を利用して電食部5の位置及び損傷領域を推定出来ていることがわかる。
従来、レールR破断を検知する技術として、列車検知用の軌道回路を利用した方法が広く適用されてきた。一方で、軌道回路の利用廃止の動向に合わせて、レールR破断時の衝撃振動を起点に伝搬した超音波を検知するといった軌道回路を利用しない方法も提案さている。この方法では、き裂など比較的程度の大きな損傷に対しては有効に機能する。しかしながら、例えば腐食などのき裂進展以前の損傷については、検出感度などの理由から適用は困難と考えられる。
き裂進展以前の損傷を検知する手法としては、例えば、特許文献2で開示されている超音波を利用した方法が提案されていた。特許文献2の手法では、レールRの頂面である踏面に当てた超音波探触子からレールR内に向けて超音波を照射し、損傷部からの反射エコーを検知することで損傷の有無を知ることができる。特許文献3では、レールR底部に設置した磁気センサを利用して検査箇所の漏洩電流を計測し、レールRの電食状態を判定する方法が提案されていた。
しかしながら、特許文献2に開示されている方法では、レールR底端部などの超音波が到達しにくい箇所への適用は容易ではない。具体的には、レールRとレールクリップ間で発生する電食検知などには適用が容易ではない。一方、特許文献3に開示されている方法であれば、レールRとレールクリップ間で発生する電食検知は可能と推察される。しかしながら、電食損傷位置を特定するためにはレールクリップ毎に磁気センサを配置しなければならず、作業が容易ではない。
それに対して、第1の実施形態における損傷領域推定システム100では、鉄道車両RSが走行するレールRで発生する弾性波を検出する複数のセンサ10と、複数のセンサ10それぞれによって検出された複数の弾性波に基づいて、複数の弾性波の発生源の位置を標定する位置標定部323と、複数の弾性波の発生源の位置に基づいて、レールRにおける損傷領域を推定する推定部325とを備える。そのため、破断が発生する前段階でレールRにおける損傷領域を簡便に推定することが可能になる。
さらに、損傷領域推定システム100では、複数のセンサ10それぞれによって検出された複数の弾性波の特徴量を抽出する信号処理部20をさらに備える。位置標定部323は、弾性波の特徴量に基づいて複数の弾性波の発生源の位置を標定した後に発生源分布を導出する。推定部325は、発生源分布を用いてレールRにおける損傷領域を推定する。このように、弾性波源の位置が示された分布を利用するため、弾性波源の数に応じて損傷領域を推定することができる。
さらに、損傷領域推定システム100では、推定部325が、発生源分布において密度が閾値以上の領域を損傷領域として推定する。弾性波源の密度が閾値以上である領域は、弾性波がより多く発生している領域である。弾性波は、材料の変形、破壊又は電気化学反応時などに損傷部分から放出される波である。弾性波源が多い領域ほど、損傷が発生している可能性がある。このように、損傷領域推定システム100では、発生源分布を用いて、容易に損傷領域を推定することができる。
(第2の実施形態)
第1の実施形態では、信号処理部と推定装置とが有線により接続されていた。第2の実施形態では、信号処理部と推定装置とを無線により接続する構成について説明する。
第1の実施形態では、信号処理部と推定装置とが有線により接続されていた。第2の実施形態では、信号処理部と推定装置とを無線により接続する構成について説明する。
図9は、第2の実施形態における損傷領域推定システム100aの構成を表す図である。損傷領域推定システム100aは、複数のセンサ10-1~10-n、信号処理部20a、無線通信部50、電源供給部60及び推定装置30aを備える。第2の実施形態において、センサ10、信号処理部20a、無線通信部50及び電源供給部60は、レールRで発生する弾性波を検出する検出装置として構成される。以下、第1の実施形態との相違点について説明する。
信号処理部20aは、信号処理により得られたデータを送信データとして無線通信部50に出力する。
無線通信部50は、信号処理部20aから出力された送信データを、所定のタイミングで推定装置30aに送信する。無線通信部50は、第1の無線通信部の一態様である。
電源供給部60は、信号処理部20a及び無線通信部50に電源を供給する。電源供給部60は、自立電源であることが好ましい。電源供給部60は、例えば1次電池、帰線電流による二次電池の充電又は振動発電機に代表されるエナジーハーベスタを利用する装置である。
(振動発電機)
電源供給部60は、信号処理部20a及び無線通信部50に電源を供給可能な位置であればどこに設置されてもよい。電源供給部60がレールR又はレールRに接触している物体に設置されている場合、電源供給部60をレールRの上下振動を利用して発電する振動発電機とすることで、電池交換が不要でかつ配線を台車上で完結することができ、設置コストを低減することができ、好適である。振動発電機における可動子の固有振動数が、設置する場所における固有振動周波数±10%の範囲に含まれるように構成することで、より大きな発電量が得られる。特に振動発電機における可動子の固有振動数が設置場所の1次曲げ固有振動数近傍に設定されることが望ましい。
電源供給部60は、信号処理部20a及び無線通信部50に電源を供給可能な位置であればどこに設置されてもよい。電源供給部60がレールR又はレールRに接触している物体に設置されている場合、電源供給部60をレールRの上下振動を利用して発電する振動発電機とすることで、電池交換が不要でかつ配線を台車上で完結することができ、設置コストを低減することができ、好適である。振動発電機における可動子の固有振動数が、設置する場所における固有振動周波数±10%の範囲に含まれるように構成することで、より大きな発電量が得られる。特に振動発電機における可動子の固有振動数が設置場所の1次曲げ固有振動数近傍に設定されることが望ましい。
推定装置30aは、無線通信部35、制御部32、記憶部33及び表示部34を備える。推定装置30aは、有線通信部31に代えて無線通信部35を備える点で第1の実施形態と構成が異なる。推定装置30aは、それ以外の構成については、第1の実施形態における推定装置30と構成は同じである。
無線通信部35は、無線通信部50から送信された送信データを無線により受信する。無線通信部35は、第2の無線通信部の一態様である。
第2の実施形態における損傷領域推定システム100aの処理は、検出装置から推定装置30aへの送信データの送信が有線から無線になった点を除けば第1の実施形態と同様である。
以上のように構成された損傷領域推定システム100aによれば、弾性波の特徴量のデータを含む送信データを無線により送信する無線通信部50と、少なくとも無線通信部50に給電する電源供給部60と、無線通信部50から送信された送信データを受信する無線通信部35を備えている。これにより、無線通信部50を備える検出装置は、電源供給部60から給電される電力で動作可能になる。さらに、無線化することにより、ケーブルの敷設や管理の作業コストを削減することができる。
さらに、損傷領域推定システム100aでは、少なくとも無線通信部35、位置標定部323及び推定部325を備える推定装置30aが鉄道車両RSに設置される。これにより、推定装置30aを載せた鉄道車両RS(例えば、営業車や検測車など)を少なくとも1台以上走行させて、検出装置からデータを吸い上げればよく、データ収集効率を向上させることができる。
さらに、損傷領域推定システム100aでは、電源供給部60として、1次電池、帰線電流による二次電池の充電、振動発電機に代表されるエナジーハーベスタ等の自立電源を用いる。これにより、電源が近くにない場所であっても容易に送信データを送信することが可能になる。
(第2の実施形態における損傷領域推定システム100aの変形例)
損傷領域推定システム100aにおいて、間欠動作を可能とするウェイクアップ機能を備えるように構成されてもよい。間欠動作を可能とするウェイクアップ機能を備えた構成を図10に示す。図10は、変形例における損傷領域推定システム100bの構成を表す図である。損傷領域推定システム100bは、複数のセンサ10-1~10-n、信号処理部20a、無線通信部50、電源供給部60、検出部70、制御部80及び推定装置30bを備える。センサ10、信号処理部20b、無線通信部50、電源供給部60、検出部70及び制御部80は、レールRで発生する弾性波を検出する検出装置として構成される。検出装置において、信号処理部20b及び無線通信部50が処理開始前はスリープ状態である点、検出部70及び制御部80を新たに備える点が第2の実施形態との相違点である。
損傷領域推定システム100aにおいて、間欠動作を可能とするウェイクアップ機能を備えるように構成されてもよい。間欠動作を可能とするウェイクアップ機能を備えた構成を図10に示す。図10は、変形例における損傷領域推定システム100bの構成を表す図である。損傷領域推定システム100bは、複数のセンサ10-1~10-n、信号処理部20a、無線通信部50、電源供給部60、検出部70、制御部80及び推定装置30bを備える。センサ10、信号処理部20b、無線通信部50、電源供給部60、検出部70及び制御部80は、レールRで発生する弾性波を検出する検出装置として構成される。検出装置において、信号処理部20b及び無線通信部50が処理開始前はスリープ状態である点、検出部70及び制御部80を新たに備える点が第2の実施形態との相違点である。
検出部70は、環境の変化を検出する。例えば、検出部70は、加速度センサである。環境の変化とは、加速度が所定値以上変化したことである。検出部70は、鉄道車両RSの走行振動を監視し、所定値以上の加速度を検出した場合に、環境が変化したことを制御部80に通知する。
制御部80は、検出部70により環境の変化が検出された場合に、スリープ状態である機能部を起動させる。例えば、制御部80は、電源供給部60に対して信号処理部20b及び無線通信部50への給電を行うように制御することによってスリープ状態である信号処理部20b及び無線通信部50を起動させる。なお、信号処理部20b又は無線通信部50の一方がスリープ状態であってもよい。
このように構成されることによって、信号処理部20b及び無線通信部50は常に起動している必要はなく、消費電力を抑制することができる。
なお、電源供給部60として振動発電機を利用する場合には、制御部80は振動発電機からの電源の出力信号が得られたことを契機にスリープ状態である機能部を起動させる。この場合、損傷領域推定システム100bは、検出部70を備えなくてよい。
このように構成されることによって、上記の検出部70を備えなくて済む。そのため、検出部70を備える構成よりもコストを削減することができる。
(第3の実施形態)
第1の実施形態及び第2の実施形態では、センサで検出される弾性波を全て利用する構成を説明した。しかしながら、実際に現場で運用するためには、各種外乱への対策が必要である。例えば、センサで検出される弾性波には、電食により発生する弾性波(以下「電食弾性波」という。)の他に、自然腐食に起因する弾性波(以下「自然腐食弾性波」という。)等がある。自然腐食が発生する現場で作業する場合には、電食弾性波と、自然腐食弾性波とが混在する可能性が高い。電食による損傷領域の推定精度の劣化を抑制するためには、これらの弾性波を分離する必要がある。そこで、第3の実施形態では、センサで検出される弾性波のうち、自然腐食弾性波を除くための構成について説明する。
第1の実施形態及び第2の実施形態では、センサで検出される弾性波を全て利用する構成を説明した。しかしながら、実際に現場で運用するためには、各種外乱への対策が必要である。例えば、センサで検出される弾性波には、電食により発生する弾性波(以下「電食弾性波」という。)の他に、自然腐食に起因する弾性波(以下「自然腐食弾性波」という。)等がある。自然腐食が発生する現場で作業する場合には、電食弾性波と、自然腐食弾性波とが混在する可能性が高い。電食による損傷領域の推定精度の劣化を抑制するためには、これらの弾性波を分離する必要がある。そこで、第3の実施形態では、センサで検出される弾性波のうち、自然腐食弾性波を除くための構成について説明する。
図11は、図6に示す環境下で電流印加の無い場合に得られた弾性波の波形の一例を示す図である。図11の左に示す図は、横軸を時間、縦軸を振幅とした弾性波の波形を表し、図11の右に示す図は、横軸を周波数、縦軸を振幅とした弾性波の波形を表す。図11に示す弾性波は、湿潤環境下で生じた自然腐食弾性波であり電食弾性波とは異なる。図11の右に示すように、自然腐食弾性波は、高域(例えば、100kHz以上)に主共振周波数(図11の破線8)を持っていることがわかる。このように、電食弾性波と、自然腐食弾性波とでは、主共振周波数が異なる。そこで、センサ10により得られた弾性波の周波数成分を確認することで、電食弾性波と、自然腐食弾性波とを分離することができる。
図12は、第3の実施形態における損傷領域推定システム100cの構成を表す図である。損傷領域推定システム100cは、複数のセンサ10-1~10-n、信号処理部20a、無線通信部50、電源供給部60及び推定装置30cを備える。推定装置30c以外の構成は、第2の実施形態と同様である。
推定装置30cは、無線通信部35、制御部32c、記憶部33及び表示部34を備える。制御部32cは、推定装置30c全体を制御する。制御部32cは、CPU等のプロセッサやメモリを用いて構成される。制御部32cは、プログラムを実行することによって、取得部321、イベント抽出部322、信号判定部326、位置標定部323c、分布生成部324及び推定部325として機能する。
取得部321、イベント抽出部322、信号判定部326、位置標定部323c、分布生成部324及び推定部325の機能部のうち一部または全部は、ASICやPLD、FPGAなどのハードウェアによって実現されてもよいし、ソフトウェアとハードウェアとの協働によって実現されてもよい。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM、CD-ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置などの非一時的な記憶媒体である。プログラムは、電気通信回線を介して送信されてもよい。
取得部321、イベント抽出部322、信号判定部326、位置標定部323c、分布生成部324及び推定部325の機能の一部は、予め推定装置30cに搭載されている必要はなく、追加のアプリケーションプログラムが推定装置30cにインストールされることで実現されてもよい。
制御部32cは、位置標定部323に代えて位置標定部323cを備える点、信号判定部326を新たに備える点で制御部32と構成が異なる。以下、相違点について説明する。
信号判定部326は、複数のセンサ10それぞれによって検出された弾性波毎に主共振周波数が所定の範囲内であるか否かを判定する。例えば、信号判定部326は、センサ10によって検出された弾性波毎に、主共振周波数が低域(例えば、100kHz)未満であるか否かを判定する。
位置標定部323cは、信号判定部326によって主共振周波数が所定の範囲内であると判定された弾性波を用いて弾性波源の位置を標定する。主共振周波数が所定の範囲内であると判定された弾性波は、上述したように、電食弾性波である。このように、位置標定部323cは、電食弾性波を用いて弾性波源の位置を標定する。
以上のように構成された損傷領域推定システム100cによれば、複数のセンサ10それぞれによって検出された弾性波毎に、主共振周波数が所定の範囲内であるか否かを判定する信号判定部326を備える。そして、位置標定部323cは、信号判定部326によって主共振周波数が所定の範囲内であると判定された弾性波を用いて発生源の位置を標定する。これにより、損傷領域推定システム100cでは、電食弾性波と、自然腐食弾性波とを分離して用いることができる。したがって、電食に起因する損傷領域を精度よく推定することが可能になる。
(第3の実施形態における損傷領域推定システム100cの変形例)
上記の実施形態では、電食弾性波と、自然腐食弾性波とを分離して、電食に起因する損傷領域を推定する構成を示した。推定装置30cは、自然腐食に起因する損傷領域を推定するように構成されてもよい。このように構成される場合、位置標定部323cは、信号判定部326によって主共振周波数が所定の範囲内ではないと判定された弾性波を用いて弾性波源の位置を標定すればよい。そして、推定部325は、位置標定部323cによって標定された弾性波源の位置に基づいて損傷領域を推定する。このように、位置標定部323cによって標定された弾性波源の位置は、自然腐食に起因する弾性波の発生源である。そのため、自然腐食に起因する損傷領域を精度よく推定することが可能になる。
上記の実施形態では、電食弾性波と、自然腐食弾性波とを分離して、電食に起因する損傷領域を推定する構成を示した。推定装置30cは、自然腐食に起因する損傷領域を推定するように構成されてもよい。このように構成される場合、位置標定部323cは、信号判定部326によって主共振周波数が所定の範囲内ではないと判定された弾性波を用いて弾性波源の位置を標定すればよい。そして、推定部325は、位置標定部323cによって標定された弾性波源の位置に基づいて損傷領域を推定する。このように、位置標定部323cによって標定された弾性波源の位置は、自然腐食に起因する弾性波の発生源である。そのため、自然腐食に起因する損傷領域を精度よく推定することが可能になる。
推定装置30cは、電食に起因する損傷領域と、自然腐食に起因する損傷領域との両方を推定するように構成されてもよい。このように構成される場合、位置標定部323cは、第1段階として信号判定部326によって主共振周波数が所定の範囲内であると判定された弾性波を用いて弾性波源の位置を標定する。さらに、位置標定部323cは、第2段階として信号判定部326によって主共振周波数が所定の範囲内ではないと判定された弾性波を用いて弾性波源の位置を標定する。そして、推定部325は、位置標定部323cによって別々に標定された弾性波源の位置に基づいて、電食に起因する損傷領域と、自然腐食に起因する損傷領域との両方を推定する。推定部325は、推定結果を別々に表示部34に表示する。
(第4の実施形態)
鉄道車両の帰線電流による電食誘起を想定した場合、鉄道走行位置と電食弾性波の発生源とが所定の距離以内では走行外乱をセンサにより同時に検出して電食弾性波のセンサ到達時刻が正確に読み取れず、位置標定の精度劣化を招く可能性がある。一方で、走行外乱の周波数帯域が電食弾性波の主要帯域から離れている場合には、周波数フィルタリングによる外乱除去は十分可能と考えられる。しかし、両者が接近している場合は難しい。このような場合に有効な手段の一つとして、時間領域での外乱除去が挙げられる。具体的には、ある時刻で取得した鉄道車両の位置から所定の距離以上離れたセンサを選別することで、外乱の混入回避が期待できる。そこで、第4の実施形態では、鉄道車両の位置に応じて、位置標定に利用するセンサを選択する構成について説明する。
鉄道車両の帰線電流による電食誘起を想定した場合、鉄道走行位置と電食弾性波の発生源とが所定の距離以内では走行外乱をセンサにより同時に検出して電食弾性波のセンサ到達時刻が正確に読み取れず、位置標定の精度劣化を招く可能性がある。一方で、走行外乱の周波数帯域が電食弾性波の主要帯域から離れている場合には、周波数フィルタリングによる外乱除去は十分可能と考えられる。しかし、両者が接近している場合は難しい。このような場合に有効な手段の一つとして、時間領域での外乱除去が挙げられる。具体的には、ある時刻で取得した鉄道車両の位置から所定の距離以上離れたセンサを選別することで、外乱の混入回避が期待できる。そこで、第4の実施形態では、鉄道車両の位置に応じて、位置標定に利用するセンサを選択する構成について説明する。
図13は、第4の実施形態における損傷領域推定システム100dの構成を表す図である。損傷領域推定システム100dは、複数のセンサ10-1~10-n、信号処理部20a、無線通信部50、電源供給部60、推定装置30d及び位置情報提供装置90を備える。推定装置30d及び位置情報提供装置90以外の構成は、第2の実施形態と同様である。
位置情報提供装置90は、鉄道車両RSの位置情報を取得して、取得した位置情報を推定装置30dに提供する。位置情報提供装置90は、例えばGPS(Global Positioning System)により位置情報を取得してもよい。なお、位置情報提供装置90は、鉄道車両RSの位置情報を取得できればその他の方法で鉄道車両RSの位置情報を取得してもよい。例えば、位置情報提供装置90は、軌道回路やCBTC(Communication Based Train Control)の列車検知情報を参照してもよい。
推定装置30dは、無線通信部35、制御部32d、記憶部33及び表示部34を備える。制御部32dは、推定装置30d全体を制御する。制御部32dは、CPU等のプロセッサやメモリを用いて構成される。制御部32dは、プログラムを実行することによって、取得部321、イベント抽出部322、信号選別部327、位置標定部323d、分布生成部324及び推定部325として機能する。
取得部321、イベント抽出部322、信号選別部327、位置標定部323d、分布生成部324及び推定部325の機能部のうち一部または全部は、ASICやPLD、FPGAなどのハードウェアによって実現されてもよいし、ソフトウェアとハードウェアとの協働によって実現されてもよい。プログラムは、コンピュータ読み取り可能な記録媒体に記録されてもよい。コンピュータ読み取り可能な記録媒体とは、例えばフレキシブルディスク、光磁気ディスク、ROM、CD-ROM等の可搬媒体、コンピュータシステムに内蔵されるハードディスク等の記憶装置などの非一時的な記憶媒体である。プログラムは、電気通信回線を介して送信されてもよい。
取得部321、イベント抽出部322、信号選別部327、位置標定部323d、分布生成部324及び推定部325の機能の一部は、予め推定装置30dに搭載されている必要はなく、追加のアプリケーションプログラムが推定装置30dにインストールされることで実現されてもよい。
制御部32dは、位置標定部323に代えて位置標定部323dを備える点、信号選別部327を新たに備える点で制御部32と構成が異なる。以下、相違点について説明する。
信号選別部327は、位置情報提供装置90から提供された鉄道車両RSの位置情報に基づいて、センサ10と鉄道車両RSとの距離が所定値以上離れているセンサ10を複数のセンサ10の中から選択する。
信号選別部327によるセンサ10の選択方法について説明する。図14及び図15は、信号選別部327の具体的な処理を説明するための図である。図14及び図15において、ある地点の変電所1と他の地点の変電所1との間の距離をLとし、鉄道車両RSの位置をXTとし、あるセンサ10の位置をXs0とし、あるセンサ10と他のセンサ10との間の距離をΔXsmaxとする。鉄道車両RSの位置は、鉄道車両RSの中心の位置であってもよい。この場合、鉄道車両RSの位置と、先頭の車輪との間の距離をΔXTfとし、鉄道車両RSの位置と、後尾の車輪との間の距離をΔXTrとする。
図14では、鉄道車両RSの位置XTが地点L/2の位置に達していない状態を示している。この場合、信号選別部327は、以下の式(1)を満たすセンサ10を複数のセンサ10の中から選択する。
Xs0>XT+ΔXTf+ΔXsmax・・・式(1)
図15では、鉄道車両RSの位置XTが地点L/2の位置を通り過ぎている状態を示している。この場合、信号選別部327は、以下の式(2)を満たすセンサ10を複数のセンサ10の中から選択する。
Xs0<XT-ΔXTr-ΔXsmax・・・式(2)
位置標定部323dは、信号選別部327によって選択されたセンサ10によって検出された弾性波に基づいて弾性波源の位置を標定する。
以上のように構成された損傷領域推定システム100dによれば、鉄道車両RSの位置情報を提供する位置情報提供装置90と、位置情報提供装置90から提供された鉄道車両RSの位置情報に基づいて、センサ10と鉄道車両RSとの距離が所定値以上離れているセンサ10を複数のセンサ10の中から選択する信号選別部327と、を備える。位置標定部323dは、選択されたセンサ10によって検出された弾性波に基づいて弾性波源の位置を標定する。このように、センサ10と鉄道車両RSとの距離が所定値以上離れているセンサ10を複数のセンサ10の中から選択することで、走行外乱の影響を抑制して位置標定ができる。そのため、位置標定の精度劣化を抑制することができる。
(第1の実施形態から第4の実施形態のいずれかに共通する変形例)
上記の各実施形態では、損傷領域の推定対象が、軌条である場合を例に説明した。損傷領域の推定対象は、鉄道に関する対象物であればよく、軌条に限定される必要はない。例えば、損傷領域の推定対象となる鉄道に関する対象物は、架線であってもよい。第2の実施形態から第4の実施形態のように、無線方式を用いる場合には、軌条の他に、架線の損傷領域を推定することも可能である。架線においても、軌条と同様に腐食により弾性波が発生する。センサ10の設置対象が架線に代わる点を除けば、処理は第2の実施形態から第4の実施形態と同様である。
上記の各実施形態では、損傷領域の推定対象が、軌条である場合を例に説明した。損傷領域の推定対象は、鉄道に関する対象物であればよく、軌条に限定される必要はない。例えば、損傷領域の推定対象となる鉄道に関する対象物は、架線であってもよい。第2の実施形態から第4の実施形態のように、無線方式を用いる場合には、軌条の他に、架線の損傷領域を推定することも可能である。架線においても、軌条と同様に腐食により弾性波が発生する。センサ10の設置対象が架線に代わる点を除けば、処理は第2の実施形態から第4の実施形態と同様である。
第3の実施形態及び第4の実施形態では、検出装置と推定装置とが有線により接続されるように構成されてもよい。第3の実施形態及び第4の実施形態は、第2の実施形態と同様に変形されてもよい。
第1の実施形態から第4の実施形態では、レールRを流れる電流として、鉄道車両RSと変電所との間を流れる帰線電流を例に説明したが、レールRに電流を流すことができればそのほかの手法で与えられる電流であってもよい。例えば、図6に示したように、定電流電源を用いて、任意のタイミングでレールRに電流を流してもよい。この方法であれば、検査したいタイミング検査が可能になる。
推定装置30,30a,30b,30cが備える一部の機能は、別の筺体に備えられていてもよい。第1の実施形態における推定装置30を例にすると、例えば、推定装置30の制御部32には取得部321及び推定部325の機能が備えられ、イベント抽出部322、位置標定部323及び分布生成部324の機能が別の筺体に備えられてもよい。このように構成される場合、推定部325は、弾性波源分布を別の筺体から取得し、取得した弾性波源分布を用いて損傷領域を推定する。
第1の実施形態から第4の実施形態において、センサ10として低周波のセンサ10を用いるように構成されてもよい。低周波のセンサ10とは、例えば共振周波数が70kHzのセンサである。上述したように、電食弾性波は、低域に主共振周波数をもつという特徴がある。したがって、センサ10の中心周波数を電食弾性波の主共振周波数に合致させることで自然腐食弾性波の分離及び電食弾性波を高感度に検出することができる。具体的には、センサ10として低周波数のセンサ10(共振周波数:70kHz)を用いればよい。図16は、低周波数のセンサ10と、高周波数のセンサ10とを用いて電食弾性波を検出した結果を示す図である。図16(A)は低周波数のセンサ10を用いて電食弾性波を検出した結果を表し、図16(B)は高周波数のセンサ10を用いて電食弾性波を検出した結果を表す。図16(A)及び図16(B)に示すように、低周波数のセンサ10を適用することで、高周波数のセンサ10(共振周波数:150kHz)よりも電食弾性波を高感度に検出できていることがわかる。
信号処理部20は、推定装置30,30c,30dに備えられてもよい。信号処理部20aは、推定装置30aに備えられてもよい。信号処理部20bは、推定装置30bに備えられてもよい。このように構成される場合、信号処理部20,20a,20bは、センサ10から出力された電気信号を、センサ10から直接、又は、不図示の中継装置を介して取得する。
推定部325は、出力制御部として動作してもよい。出力制御部は、出力部を制御して、推定結果を出力する。ここで、出力部には、表示部34、通信部および印刷部が含まれる。出力部が通信部である場合、出力制御部は通信部を制御して、推定結果を他の装置に送信する。出力部が印刷部である場合、出力制御部は印刷部を制御して、推定結果を印刷する。なお、推定装置30,30a,30b,30c,30dは、表示部34、通信部および印刷部の一部又は全てを備えて上記の動作を実行してもよい。
以上説明した少なくともひとつの実施形態によれば、損傷領域の推定対象となる鉄道に関する対象物(例えば、鉄道車両RSが走行する軌条)で発生する弾性波を検出する複数のセンサと、複数のセンサそれぞれによって検出された複数の弾性波に基づいて、複数の弾性波の発生源の位置を標定する位置標定部と、複数の弾性波の発生源の位置に基づいて、対象物における損傷領域を推定する推定部とを持つことにより、破断が発生する前段階で鉄道に関する対象物の損傷領域を簡便に推定することができる。
本発明のいくつかの実施形態を説明したが、これらの実施形態は、例として提示したものであり、発明の範囲を限定することは意図していない。これら実施形態は、その他の様々な形態で実施されることが可能であり、発明の要旨を逸脱しない範囲で、種々の省略、置き換え、変更を行うことができる。これら実施形態やその変形は、発明の範囲や要旨に含まれると同様に、特許請求の範囲に記載された発明とその均等の範囲に含まれるものである。
10、10-1~10-n…センサ,20、20a、20b…信号処理部,21…フィルタ,22…A/D変換器,23…波形整形フィルタ,24…ゲート生成回路,25…到達時刻決定部,26…特徴量抽出部,27…送信データ生成部,28…メモリ,29…出力部,30、30a、30b、30c、30d…推定装置,50…無線通信部,60…電源供給部,70…検出部,80…制御部,90…位置情報提供装置,31…有線通信部,32、32c、32d…制御部,33…記憶部,34…表示部,321…取得部,322…イベント抽出部,323、323c、323d…位置標定部,324…分布生成部,325…推定部,326…信号判定部,327…信号選別部
Claims (18)
- 損傷領域の推定対象となる鉄道に関する対象物で発生する弾性波を検出する複数のセンサと、
前記複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する位置標定部と、
前記複数の弾性波の発生源の位置に基づいて、前記対象物における損傷領域を推定する推定部と、
を備える損傷領域推定システム。 - 前記複数のセンサそれぞれによって検出された複数の弾性波の特徴量を抽出する信号処理部をさらに備え、
前記位置標定部は、前記弾性波の特徴量に基づいて前記複数の弾性波の発生源の位置を標定し、標定した前記複数の弾性波の発生源の位置を表す発生源分布を導出し、
前記推定部は、前記発生源分布を用いて前記対象物における損傷領域を推定する、
請求項1に記載の損傷領域推定システム。 - 前記推定部は、前記発生源分布において密度が閾値以上の領域を前記損傷領域として推定する、
請求項2に記載の損傷領域推定システム。 - 前記弾性波は、前記対象物を流れる電流が漏出して電食が発生する際に放出される波である、
請求項1から3のいずれか一項に記載の損傷領域推定システム。 - 前記対象物を流れる電流は、鉄道車両と変電所との間を流れる帰線電流である、
請求項4に記載の損傷領域推定システム。 - 前記複数のセンサは、少なくとも前記対象物、又は、前記対象物に接触している物体のいずれかに設置される、
請求項1から5のいずれか一項に記載の損傷領域推定システム。 - 前記複数のセンサそれぞれによって検出された複数の弾性波により得られる弾性波の特徴量のデータを無線により送信する第1の無線通信部と、
少なくとも前記第1の無線通信部に給電する電源供給部と、
前記第1の無線通信部から送信された前記弾性波の特徴量のデータを受信する第2の無線通信部と、
をさらに備え、
前記位置標定部は、前記第2の無線通信部により受信された前記弾性波の特徴量のデータを用いて前記複数の弾性波の発生源の位置を標定する、
請求項1から6のいずれか一項に記載の損傷領域推定システム。 - 前記第2の無線通信部、前記位置標定部及び前記推定部は、鉄道車両に設置される、
請求項7に記載の損傷領域推定システム。 - 前記電源供給部は、自立電源である、
請求項7又は8に記載の損傷領域推定システム。 - 環境の変化を検出する検出部と、
前記検出部により環境の変化が検出された場合に、少なくとも前記第1の無線通信部を起動させる制御部と、
をさらに備える、
請求項7から9のいずれか一項に記載の損傷領域推定システム。 - 前記電源供給部が振動発電機である場合、前記振動発電機が動作したことを契機に少なくとも前記第1の無線通信部を起動させる制御部と、
をさらに備える、
請求項7から9のいずれか一項に記載の損傷領域推定システム。 - 前記複数のセンサそれぞれによって検出された弾性波毎に主共振周波数が所定の範囲内であるか否かを判定する信号判定部、
をさらに備え、
前記位置標定部は、前記信号判定部によって主共振周波数が所定の範囲内であると判定された弾性波を用いて発生源の位置を標定する、
請求項1から11のいずれか一項に記載の損傷領域推定システム。 - 鉄道車両の位置情報を提供する位置情報提供装置と、
前記位置情報提供装置から提供された前記鉄道車両の位置情報に基づいて、センサと前記鉄道車両との距離が所定値以上離れているセンサを前記複数のセンサの中から選択する信号選別部と、
をさらに備え、
前記位置標定部は、選択されたセンサによって検出された弾性波に基づいて発生源の位置を標定する、
請求項1から12のいずれか一項に記載の損傷領域推定システム。 - 前記複数のセンサの中心周波数は、前記弾性波の主共振周波数に略等しい、
請求項1から13のいずれか一項に記載の損傷領域推定システム。 - 前記対象物は、鉄道車両が走行する軌条である、
請求項1から14のいずれか一項に記載の損傷領域推定システム。 - 損傷領域の推定対象となる鉄道に関する対象物で発生する弾性波を検出する複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定する位置標定部と、
前記複数の弾性波の発生源の位置に基づいて、前記対象物における損傷領域を推定する推定部と、
を備える推定装置。 - 損傷領域の推定対象となる鉄道に関する対象物で発生する弾性波を検出する複数のセンサそれぞれによって検出された複数の弾性波に基づいて、前記複数の弾性波の発生源の位置を標定し、
標定により得られた前記複数の弾性波の発生源の位置に基づいて前記対象物における損傷領域を推定する、
損傷領域推定方法。 - 前記対象物は軌条であり、
前記軌条に電流を印加し、
前記軌条を流れる電流が漏出して電食が発生する際に放出される弾性波の発生源の位置を標定し、
標定した前記複数の弾性波の発生源の位置を表す発生源分布を導出し、
前記発生源分布を用いて前記軌条における損傷領域を推定する、
請求項17に記載の損傷領域推定方法。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021070431A JP2022165187A (ja) | 2021-04-19 | 2021-04-19 | 損傷領域推定システム、推定装置及び損傷領域推定方法 |
US17/653,493 US11946906B2 (en) | 2021-04-19 | 2022-03-04 | Damaged region determination system, determination apparatus and damaged region determination method |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2021070431A JP2022165187A (ja) | 2021-04-19 | 2021-04-19 | 損傷領域推定システム、推定装置及び損傷領域推定方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
JP2022165187A true JP2022165187A (ja) | 2022-10-31 |
Family
ID=83601272
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2021070431A Pending JP2022165187A (ja) | 2021-04-19 | 2021-04-19 | 損傷領域推定システム、推定装置及び損傷領域推定方法 |
Country Status (2)
Country | Link |
---|---|
US (1) | US11946906B2 (ja) |
JP (1) | JP2022165187A (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102023125478A1 (de) | 2022-09-30 | 2024-04-04 | Semiconductor Energy Laboratory Co., Ltd. | Halbleitervorrichtung, Herstellungsverfahren der Halbleitervorrichtung und eines elektronischen Geräts |
Families Citing this family (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP7447045B2 (ja) * | 2021-03-22 | 2024-03-11 | 株式会社東芝 | 検査システム、検査装置及び検査方法 |
CN116039698B (zh) * | 2023-03-31 | 2023-07-07 | 成都盛锴科技有限公司 | 一种利用声音特征进行轨道线路健康检测的方法 |
Family Cites Families (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2002286700A (ja) | 2001-03-28 | 2002-10-03 | Railway Technical Res Inst | Ae音による基礎構造物の破壊探知システム |
JP4718857B2 (ja) | 2005-02-10 | 2011-07-06 | 東京計器株式会社 | レール底部腐食検知装置及びレール底部腐食検知方法 |
JP4367710B2 (ja) | 2005-03-15 | 2009-11-18 | 財団法人鉄道総合技術研究所 | レールの電食状態判定システム |
JP5681512B2 (ja) | 2011-02-01 | 2015-03-11 | 日本信号株式会社 | レール破断検知装置 |
JP2014080133A (ja) | 2012-10-17 | 2014-05-08 | Tokyo Keiki Inc | 鉄道レール破断検出装置及び方法 |
NO341972B1 (en) * | 2016-09-07 | 2018-03-05 | Stiftelsen Norsar | A railway track condition monitoring system for detecting a partial or complete disruption of a rail of the railway track |
WO2019136321A1 (en) * | 2018-01-05 | 2019-07-11 | The Board Of Regents Of The Nevada System Of Higher Education On Behalf Of The University Of Nevada, Las Vegas | Acoustic detection of defects in rail at high speed |
JP6857777B2 (ja) | 2018-09-13 | 2021-04-14 | 株式会社東芝 | 構造物評価システム、構造物評価装置及び構造物評価方法 |
-
2021
- 2021-04-19 JP JP2021070431A patent/JP2022165187A/ja active Pending
-
2022
- 2022-03-04 US US17/653,493 patent/US11946906B2/en active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
DE102023125478A1 (de) | 2022-09-30 | 2024-04-04 | Semiconductor Energy Laboratory Co., Ltd. | Halbleitervorrichtung, Herstellungsverfahren der Halbleitervorrichtung und eines elektronischen Geräts |
Also Published As
Publication number | Publication date |
---|---|
US20220334082A1 (en) | 2022-10-20 |
US11946906B2 (en) | 2024-04-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
JP2022165187A (ja) | 損傷領域推定システム、推定装置及び損傷領域推定方法 | |
US10330646B2 (en) | Structure evaluation system, structure evaluation apparatus, and structure evaluation method | |
US6962082B2 (en) | Device and method for acoustic diagnosis and measurement by pulse electromagnetic force | |
US20220135094A1 (en) | Inspection of rail health | |
CN112543725A (zh) | 轨道车声学监测系统及其使用方法 | |
CN106627662B (zh) | 一种列车定位装置及方法 | |
CN110687208A (zh) | 一种基于双曲线定位的无基准Lamb波损伤监测方法 | |
US20230024577A1 (en) | System and method for detecting a defect in a railway track rail | |
US20180217103A1 (en) | Deterioration diagnosis method, a deterioration diagnosis system, and a sensor | |
EP3321673B1 (en) | Structure evaluation system and structure evaluation method | |
CN113281401B (zh) | 无砟轨道隐蔽病害的检测方法、系统与装置 | |
US11519883B2 (en) | Structure evaluation system, structure evaluation apparatus, and structure evaluation method | |
JP7447045B2 (ja) | 検査システム、検査装置及び検査方法 | |
CN109307715B (zh) | 储罐底板主被动声融合检测方法 | |
JP2021017063A (ja) | レール破断の検知装置及びレール破断の検知方法 | |
CN105424803B (zh) | 一种列车轮对探伤处理方法及装置 | |
CN109307709B (zh) | 储罐底板主被动声融合检测系统 | |
Hong et al. | A guided wave approach for real-time health monitoring of high-speed train bogie frames | |
US20230349788A1 (en) | System and method for detecting faults in extended waveguides | |
CN204495852U (zh) | 一种机车测速设备的检测装置 | |
CN115184252B (zh) | 基于脉冲回波法的垂直接地极腐蚀损伤检测方法及系统 | |
JP7383654B2 (ja) | レール破断の検知装置及びレール破断の検知方法 | |
US20230194382A1 (en) | Structure evaluation system and structure evaluation method | |
CN103147740B (zh) | 一种测试液面回波时间的方法 | |
JP6146761B2 (ja) | 地震計のノイズ識別方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20230313 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20231101 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20240109 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20240307 |