JP4171800B2 - 電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 - Google Patents
電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 Download PDFInfo
- Publication number
- JP4171800B2 JP4171800B2 JP2002325538A JP2002325538A JP4171800B2 JP 4171800 B2 JP4171800 B2 JP 4171800B2 JP 2002325538 A JP2002325538 A JP 2002325538A JP 2002325538 A JP2002325538 A JP 2002325538A JP 4171800 B2 JP4171800 B2 JP 4171800B2
- Authority
- JP
- Japan
- Prior art keywords
- epicenter
- time
- electric field
- distance
- scale
- 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
- 230000005684 electric field Effects 0.000 title claims description 82
- 238000000034 method Methods 0.000 title claims description 14
- 238000005259 measurement Methods 0.000 claims description 19
- 238000007405 data analysis Methods 0.000 claims description 8
- 230000001133 acceleration Effects 0.000 description 10
- 238000013500 data storage Methods 0.000 description 6
- 238000010586 diagram Methods 0.000 description 6
- 230000005672 electromagnetic field Effects 0.000 description 5
- 238000004364 calculation method Methods 0.000 description 3
- 238000013480 data collection Methods 0.000 description 3
- 238000007796 conventional method Methods 0.000 description 2
- 239000012212 insulator Substances 0.000 description 2
- 230000000630 rising effect Effects 0.000 description 2
- 238000005070 sampling Methods 0.000 description 2
- 230000001360 synchronised effect Effects 0.000 description 2
- 238000004458 analytical method Methods 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000006698 induction Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 230000000644 propagated effect Effects 0.000 description 1
- 230000001902 propagating effect Effects 0.000 description 1
- 230000035945 sensitivity Effects 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V1/00—Seismology; Seismic or acoustic prospecting or detecting
- G01V1/01—Measuring or predicting earthquakes
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Acoustics & Sound (AREA)
- Environmental & Geological Engineering (AREA)
- Geology (AREA)
- Remote Sensing (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Geophysics And Detection Of Objects (AREA)
Description
【発明の属する技術分野】
本発明は、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法及びその予測装置に関する。
【0002】
【従来の技術及び発明が解決しようとする課題】
従来技術による震源時・震源距離・規模の算出は、地震波のP波のみ、またはP波とS波の計測により行われている(例えば、特許文献1参照)。これらの値の算出には、P波のみの場合、4箇所以上で同じ地震のP波を計測する必要があり、P波とS波を用いる場合でも、3箇所以上で同じ地震のP波とS波を計測する必要がある。
【0003】
このため、例えば4箇所のP波の観測点の震源距離が約100kmまたはそれ以内にある場合、震源時(地震発生時刻)・震源距離(震源と観測点間の距離)・規模(マグニチュードM)の算出には、震源時から少なくともP波が観測点に到達するまでの15秒を要する。しかし、通常はP波とS波の両方を用いているので、3箇所の観測点の震源距離が約100kmまたはそれ以内の場合、S波が観測点に到達するまでの約25秒を要する。どちらの方法にせよ、地震波が3箇所、または4箇所の観測点に到達するまで、上記の例の場合は、地震発生から15ないし25秒間は、震源距離・規模を知ることはできない。
【0004】
【特許文献1】
特開平6−324160号公報
【0005】
【課題を解決するための手段】
本発明は、上記課題を解決するものであって、簡単な方法及び装置により、地震発生とほぼ同時に震源時を検知する共に、短時間に地震波の到来前に、震源距離、規模を予測できるようにするものである。
【0006】
そのために本発明は、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法であって、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測することを特徴とするものである。
【0007】
また、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置であって、電界強度を計測する計測手段と、前記計測手段による計測値データを収集・蓄積するデータ蓄積手段と、前記データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測するデータ解析手段とを備えたことを特徴とするものである。
【0008】
さらに、前記震源時からの経過時間tにおける電界強度E(t)と該経過時間tの関係から次式[数3]を基に2個以上の計測値データにより震源距離Dと規模Mを予測し、
【0009】
【数3】
E(t)=a10M /{D−r(t)}b
ただし、a,b:地質及び電界観測システムに依存する常数
r(t):時刻tにおける震源域とP波の面間の距離
=t×P波の速度
地震波を計測する地震計と併用して震源時・震源距離・規模を予測することを特徴とするものである。
【0010】
【発明の実施の形態】
以下、本発明の実施の形態を図面を参照しつつ説明する。図1は本発明に係る電界観測に基づく震源時・震源距離・規模の予測方法の実施の形態を説明するための図であり、1は震源、2は震源から伝播する地震波、3は地震波の波面(Wave Front)であり、P波の波面であると同時に帯電面でもある。4は帯電面により誘起される誘導電界、5は電界強度を計測する電磁界センサーを示す。
【0011】
図1において、地震波2は、地震発生と同時に震源1から発生し伝播して行く。この際、P波の速度が一番速いため、地震波の波面(Wave Front)はP波の波面となる。P波は縦波であり、その波面は帯電面3となる。帯電面3は電界を誘起し誘導電界4が伝播する。
【0012】
誘導電界4は、後述の観測例に示すように地震と同時に発生し、P波の伝播と同期して強度が増大し、P波の到達時に極大となる。そして、P波の到達時からS波の到達まで、電界は減少し、S波の伝播と同期した電界は観測されない。また、P波による電界パルスは観測されず、S波の到達時に電界パルスが発生する。
【0013】
この誘導電界4の強度を電磁界センサー5により計測し、解析することにより震源時は、電界上昇開始時刻から検知し、前記震源距離及び規模(マグニチュードM)は、震源時からの経過時間と電磁界センサー5による計測値データから予測する。
【0014】
P波の伝播と同期した誘導電界4の波形は、P波の波面(Wave Front)に電荷が発生するとすれば、その誘導電界で定量的に説明できる。観測において、振幅は、P波の波面から電界観測点までの距離の二乗にほぼ反比例している。マグニチュードMの定義として、地震動の最大速度振幅の常用対数を用いるものが多い。また、流動電位及びピエゾ電荷は最大速度振輻に比例するから、電界が流動電位またはピエゾ電荷によるものならば、電界は10のM乗に比例する。よって、マグニチュードMは次式[数4]から求められる。
【0015】
【数4】
E(t)=a10M /{D−r(t)}b
ただし、E(t):震源時からの経過時間tにおける電界強度
a,b:地質及び電界観測システムに依存する常数
D:震源距離
r(t):時刻tにおける震源域とP波の面間の距離
=t×P波の速度
上式[数4]は、P波の到達前でも、震源時(電界上昇開始時刻)からの経過時間tと電界の振幅(電界強度)Eの関係から、マグニチュードMと震源距離Dは予測できることを示している。即ち、マグニチュードMと震源距離Dを未知数とする上式[数4]は、時間がtのときの電界強度E(t)の計測値データが2個以上得られれば解け、マグニチュードMと震源距離Dは求められる。よって、地震発生後、サンプリング周期を0.1秒とすると5秒で10個の計測値データが得られるので、数秒で、P波の到来前に、地震動は予測でき、上記従来の例の場合で、約10秒は短縮するこができ、大部分の地域に対し地震波の到来前の予報が可能になる。
【0016】
次に、観測装置の構成およびその処理を説明する。図2は本発明に係る電界観測に基づく震源時・震源距離・規模の予測装置の構成例を示す図であり、21は電極、22は受信機、23は地震計、24はデータ収集部、25はデータ蓄積部、26はデータ解析部、27はペンレコーダを示す。
【0017】
図2において、電極21は、電極間の電位差を出力するものである。受信機22は、電極21の出力電圧の内の、所要帯域外の出力を減少させると共に、所要帯域内の出力を所要電圧(通常は10V以下)まで増幅するものである。地震計23は、地震波による加速度等の振幅を電圧に変換して出力するものである。データ収集部24は、所定のサンプリング周期(0,1秒以下)で、受信機22及び地震計23の計測値データを収集するものである。データ蓄積部25は、その収集した計測値データについて有意・無意を判断したうえで、記憶手段に蓄積するものである。
【0018】
データ解析部26は、データ蓄積部25の記憶手段に収集し蓄積した計測値データを解析して地震と同時に発生した電界を検出し、電界上昇開始時刻から震源時を予測し、震源時からの経過時間と計測値データから震源距離と規模を予測して出力する。ペンレコーダ27は、電極21、受信機22及び地震計23の計測状況を監視するものである。地震計23の計測値データを用いる場合には、P波を検知して、P波の計測値データと電界強度の計測値データから震源距離と規模を予測する。
【0019】
次に、観測例によって、観測データの処理法を具体的に示す。図3は地震と同時に発生する電界及び該地震のP波の観測例を示す図であり、下記の2002年6月14日の地震の際の観測例で、該地震の震源時・震源距離・規模を予測するための時刻関係が示されている。31は地震と同時に発生した電界、32は南北方向の加速度、33は東西方向の加速度、34は上下方向の加速度、35は震源時、36はP波の到達時刻、37は電界の特異点、38はS波到達による電界パルスの上昇開始時刻を示す。
【0020】
上記観測例によれば、地震と同時に発生した電界31は、地震波が観測されなくなるまで継続した。南北方向の加速度32、東西方向の加速度33、上下方向の加速度34は、何れも電界の観測点(波崎)の東南14kmの銚子での強震計の記録である。震源時35(計算値)は、電界上昇開始時刻、11時42分49.0秒と計算誤差の範囲内で一致している。P波の到達時刻36は、電界観測点(波崎)での高感度地震計での観測値である。電界の特異点37は、電界が増加から減少に急変する点である。S波到達による電界パルスの上昇開始時刻38は、図2では、S波の到達時刻と一致しているが、電界の観測点(波崎)との距離差が11kmある銚子での観測値であり、S波の波崎への実際の到達時刻は、上昇開始時刻38の1秒前である。この1秒の差は、S波による横波の誘起電圧の位相遅れによる。
【0021】
上記の観測例の場合、震源距離98kmは、P波到達時刻36と震源時35の差13.7秒にP波の伝播速度7.2km/sを掛けて得られる。規模は、従来技術により、震源距離と地震波の振幅との関係から得られる他、P波到達時刻36と電界の特異点37との差からも下記の方法で得られる。本観測例の場合は、この差は0.8秒であり、これにP波の伝播速度を掛けると震源域の長さ5kmが得られ、この値は規模5(通常、規模の誤差は0.5)の地震の震源域の長さに相当する。上記の他、下記の方法で規模は得られる。
【0022】
電界の振幅(電界強度)は、規模と共に大きくなり、P波の波面(Wave Front)から電界観測点までの距離の二乗に反比例するため、P波の到達前でも、震源時35(電界上昇開始時刻)からの経過時間と電界の振幅の関係から、規模と震源距離は予測できる。即ち、規模と震源距離を未知数とする上記の関係式は、震源時35からの経過時間がtのときの電界の振幅E(t)の計測値が2個以上得られれば解けて、規模と震源距離は一意的に求められる。よって、地震発生後、数秒で、地震波の到来前に、規模・震源は予測できることになる。
【0023】
本実施形態の電磁界センサー5による電界強度の計測は、自然雑音及び人工雑音が共に少ない場所で行うようにするのが望ましい。例えば上記観測例によれば、地震に伴う電界が0.4μV/m程度であるので、最適観測帯とされる0.01−30Hz帯で、自然雑音及び人工雑音が常時0.1μV/m程度以下の観測点とすることが望ましい。
【0024】
なお、本発明は、上記実施の形態に限定されるものではなく、種々の変形が可能である。例えば上記実施の形態では、1観測地点から震源時・震源距離・規模を予測したが、震源距離から震源位置を求める手法として、従来のP波とS波から求める手法と同じ手法を採用してもよい。
【0025】
【発明の効果】
以上の説明から明らかなように、本発明によれば、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法として、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、震源時からの経過時間と計測値データから震源距離と規模を予測し、また、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置として、電界強度を計測する計測手段と、計測手段による計測値データを収集・蓄積するデータ蓄積手段と、データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、震源時からの経過時間と計測値データから震源距離と規模を予測するデータ解析手段とを備えたので、電界強度の計測と解析により短時間に地震波の到来前に、震源・震源距離・規模を予測できる。
【0026】
このため、本発明は、原子炉、高速鉄道、高速道路、空港などでの地震波到来前の緊急停止処置への利用が可能である。
【図面の簡単な説明】
【図1】 本発明に係る電界観測に基づく震源時・震源距離・規模の予測方法の実施の形態を説明するための図である。
【図2】 本発明に係る電界観測に基づく震源時・震源距離・規模の予測装置の構成例を示す図である。
【図3】 地震と同時に発生する電界及び該地震のP波の観測例を示す図である。
【符号の説明】
1…震源、2…地震波、3…P波面・帯電面、4…誘導電界、5…電磁界センサー、21…電極、22…受信機、23…地震計、24…データ収集部、25…データ蓄積部、26…データ解析部、27…ペンレコーダ、31…地霞と同時に発生した電界、32…南北方向の加速度、33…束西方向の加速度、34…上下方向の加速度、35…震源時、36…P波の到達時刻、37…電界の特異点、38…S波到達による電界パルスの上昇開始時刻
Claims (4)
- 地震波を計測する地震計と併用して震源時・震源距離・規模を予測することを特徴とする請求項1記載の電界観測に基づく震源時・震源距離・規模の予測方法。
- 地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置であって、電界強度を計測する計測手段と、前記計測手段による計測値データを収集・蓄積するデータ蓄積手段と、前記データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測するデータ解析手段とを備え、前記データ解析手段は、前記震源時からの経過時間tにおける電界強度E(t)と該経過時間tの関係から次式[数2]を基に2個以上の計測値データにより震源距離Dと規模Mを予測することを特徴とする電界観測に基づく震源時・震源距離・規模の予測装置。
r(t):時刻tにおける震源域とP波の面間の距離
=t×P波の速度 - さらに地震波を計測する地震計を備え、前記データ解析手段は、前記地震計による計測データと併用して震源時・震源距離・規模を予測することを特徴とする請求項3記載の電界観測に基づく震源時・震源距離・規模の予測装置。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002325538A JP4171800B2 (ja) | 2002-11-08 | 2002-11-08 | 電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 |
US10/634,837 US6885945B2 (en) | 2002-11-08 | 2003-08-06 | Method for estimating origin time, hypocentral distance, and scale based on electric field observation, and apparatus for prediction |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2002325538A JP4171800B2 (ja) | 2002-11-08 | 2002-11-08 | 電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2004157094A JP2004157094A (ja) | 2004-06-03 |
JP4171800B2 true JP4171800B2 (ja) | 2008-10-29 |
Family
ID=32211948
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2002325538A Expired - Lifetime JP4171800B2 (ja) | 2002-11-08 | 2002-11-08 | 電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 |
Country Status (2)
Country | Link |
---|---|
US (1) | US6885945B2 (ja) |
JP (1) | JP4171800B2 (ja) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7573780B2 (en) * | 2004-02-26 | 2009-08-11 | Exxonmobil Upstream Research Co. | Electrode configurations for suppression of electroseismic source noise |
US7392139B2 (en) * | 2006-03-27 | 2008-06-24 | Swiss Reinsurance Company | System and method for providing earthquake data |
JP4385233B2 (ja) * | 2007-07-31 | 2009-12-16 | 日本電気株式会社 | 地震予測情報通知サーバ、地震予測情報通知システム、地震予測情報の通知方法およびプログラム |
JP5470632B2 (ja) * | 2008-10-01 | 2014-04-16 | 公立大学法人首都大学東京 | 緊急地震警報を発するための方法および装置 |
CN105890916B (zh) * | 2016-04-08 | 2018-07-31 | 中国铁道科学研究院铁道科学技术研究发展中心 | 一种高速铁路车载地震紧急处置装置测试***及方法 |
CN109683195B (zh) * | 2019-01-09 | 2020-04-24 | 武汉理工大学 | 一种地震动摇摆分量获取方法 |
Family Cites Families (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JPS63184088A (ja) * | 1987-01-27 | 1988-07-29 | Radio Res Lab | 地震前兆の電界変動の観測法 |
JPH0194286A (ja) * | 1987-10-07 | 1989-04-12 | Yuuseishiyou Tsushin Sogo Kenkyusho | 地震前兆の長波・地電流の発生領域のトモグラフィ法 |
US5694129A (en) * | 1995-08-29 | 1997-12-02 | Science And Technology Agency National Research Institute For Earth Science And Disaster Prevention | Method of imminent earthquake prediction by observation of electromagnetic field and system for carrying out the same |
-
2002
- 2002-11-08 JP JP2002325538A patent/JP4171800B2/ja not_active Expired - Lifetime
-
2003
- 2003-08-06 US US10/634,837 patent/US6885945B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2004157094A (ja) | 2004-06-03 |
US20040093162A1 (en) | 2004-05-13 |
US6885945B2 (en) | 2005-04-26 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
US6962082B2 (en) | Device and method for acoustic diagnosis and measurement by pulse electromagnetic force | |
US11543286B2 (en) | Distributed acoustic sensing based acoustic wave speed scanning and mapping of civil infrastructures | |
Zhu et al. | Non-contact detection of surface waves in concrete using an air-coupled sensor | |
KR101951028B1 (ko) | 진도추정과 설비/구조물 피해예측기능을 가진 mems 기반 지진계측 장치 | |
WO2023087890A1 (zh) | 一种基于动应力、振动和oma综合分析判断构架模态共振方法 | |
JP4171800B2 (ja) | 電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 | |
AU655281B2 (en) | Improvements in and relating to particle detection and analysis | |
CN107128329B (zh) | 一种基于应变量测反演变形加速度响应的轨距动态变化监测装置与设计方法 | |
CN104483389B (zh) | 基于组合震源法检测桥梁预应力管道注浆质量的方法 | |
CN109855771A (zh) | 一种基于钢轨竖向加速度功率谱密度检测温度力的方法 | |
CN106053602A (zh) | 一种基于磁致伸缩效应的自比式锚杆无损检测方法 | |
JPS63186122A (ja) | 構造物の異常診断方式 | |
CN114537466B (zh) | 高速铁路轨道结构变形损伤监测***及方法 | |
EP1769267B1 (en) | Movement detection system and method | |
CN104570145B (zh) | 一种光纤传感生命探测设备 | |
Ryden et al. | Portable Seismic Acquistion System (PSAS) for Pavement MASW | |
Deraemaeker | Vibration based structural health monitoring using large sensor arrays: overview of instrumentation and feature extraction based on modal filters | |
CN112577752B (zh) | 在城市复杂噪声环境中用地震检波器测量汽车转速的方法 | |
Epple et al. | Ultrasonic monitoring of large-scale structures-input to engineering assessment | |
JP2002348949A (ja) | 建物の動的耐震性能評価システム | |
Ng | Characterization of the Subsurface from Evaluating Low Amplitude Seismic Recordings | |
JP3106183B2 (ja) | 地盤探査方法及び装置 | |
Mehrkens et al. | Investigation of a pre-stressed concrete structure using passive image interferometry | |
WANG et al. | Stability assessment of road rock slopes using seismic motions | |
Fiandaca et al. | A Vehicle Bridge Interaction based approach for the monitoring of bridges through an electric mobile platform |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20050304 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20080409 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20080606 |
|
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: 20080709 |
|
A01 | Written decision to grant a patent or to grant a registration (utility model) |
Free format text: JAPANESE INTERMEDIATE CODE: A01 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 4171800 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
S533 | Written request for registration of change of name |
Free format text: JAPANESE INTERMEDIATE CODE: R313533 |
|
R350 | Written notification of registration of transfer |
Free format text: JAPANESE INTERMEDIATE CODE: R350 |
|
EXPY | Cancellation because of completion of term |