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
-
- G01V1/01—
Description
【０００１】
【発明の属する技術分野】
本発明は、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法及びその予測装置に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】
従来技術による震源時・震源距離・規模の算出は、地震波のＰ波のみ、またはＰ波とＳ波の計測により行われている（例えば、特許文献１参照）。これらの値の算出には、Ｐ波のみの場合、４箇所以上で同じ地震のＰ波を計測する必要があり、Ｐ波とＳ波を用いる場合でも、３箇所以上で同じ地震のＰ波とＳ波を計測する必要がある。
【０００３】
このため、例えば４箇所のＰ波の観測点の震源距離が約１００ｋｍまたはそれ以内にある場合、震源時（地震発生時刻）・震源距離（震源と観測点間の距離）・規模（マグニチュードＭ）の算出には、震源時から少なくともＰ波が観測点に到達するまでの１５秒を要する。しかし、通常はＰ波とＳ波の両方を用いているので、３箇所の観測点の震源距離が約１００ｋｍまたはそれ以内の場合、Ｓ波が観測点に到達するまでの約２５秒を要する。どちらの方法にせよ、地震波が３箇所、または４箇所の観測点に到達するまで、上記の例の場合は、地震発生から１５ないし２５秒間は、震源距離・規模を知ることはできない。
【０００４】
【特許文献１】
特開平６−３２４１６０号公報
【０００５】
【課題を解決するための手段】
本発明は、上記課題を解決するものであって、簡単な方法及び装置により、地震発生とほぼ同時に震源時を検知する共に、短時間に地震波の到来前に、震源距離、規模を予測できるようにするものである。
【０００６】
そのために本発明は、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法であって、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測することを特徴とするものである。
【０００７】
また、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置であって、電界強度を計測する計測手段と、前記計測手段による計測値データを収集・蓄積するデータ蓄積手段と、前記データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測するデータ解析手段とを備えたことを特徴とするものである。
【０００８】
さらに、前記震源時からの経過時間ｔにおける電界強度Ｅ（ｔ）と該経過時間ｔの関係から次式［数３］を基に２個以上の計測値データにより震源距離Ｄと規模Ｍを予測し、
【０００９】
【数３】
Ｅ（ｔ）＝ａ１０M ／｛Ｄ−ｒ（ｔ）｝b
ただし、ａ，ｂ：地質及び電界観測システムに依存する常数
ｒ（ｔ）：時刻ｔにおける震源域とＰ波の面間の距離
＝ｔ×Ｐ波の速度
地震波を計測する地震計と併用して震源時・震源距離・規模を予測することを特徴とするものである。
【００１０】
【発明の実施の形態】
以下、本発明の実施の形態を図面を参照しつつ説明する。図１は本発明に係る電界観測に基づく震源時・震源距離・規模の予測方法の実施の形態を説明するための図であり、１は震源、２は震源から伝播する地震波、３は地震波の波面（Ｗａｖｅ Ｆｒｏｎｔ）であり、Ｐ波の波面であると同時に帯電面でもある。４は帯電面により誘起される誘導電界、５は電界強度を計測する電磁界センサーを示す。
【００１１】
図１において、地震波２は、地震発生と同時に震源１から発生し伝播して行く。この際、Ｐ波の速度が一番速いため、地震波の波面（Ｗａｖｅ Ｆｒｏｎｔ）はＰ波の波面となる。Ｐ波は縦波であり、その波面は帯電面３となる。帯電面３は電界を誘起し誘導電界４が伝播する。
【００１２】
誘導電界４は、後述の観測例に示すように地震と同時に発生し、Ｐ波の伝播と同期して強度が増大し、Ｐ波の到達時に極大となる。そして、Ｐ波の到達時からＳ波の到達まで、電界は減少し、Ｓ波の伝播と同期した電界は観測されない。また、Ｐ波による電界パルスは観測されず、Ｓ波の到達時に電界パルスが発生する。
【００１３】
この誘導電界４の強度を電磁界センサー５により計測し、解析することにより震源時は、電界上昇開始時刻から検知し、前記震源距離及び規模（マグニチュードＭ）は、震源時からの経過時間と電磁界センサー５による計測値データから予測する。
【００１４】
Ｐ波の伝播と同期した誘導電界４の波形は、Ｐ波の波面（Ｗａｖｅ Ｆｒｏｎｔ）に電荷が発生するとすれば、その誘導電界で定量的に説明できる。観測において、振幅は、Ｐ波の波面から電界観測点までの距離の二乗にほぼ反比例している。マグニチュードＭの定義として、地震動の最大速度振幅の常用対数を用いるものが多い。また、流動電位及びピエゾ電荷は最大速度振輻に比例するから、電界が流動電位またはピエゾ電荷によるものならば、電界は１０のＭ乗に比例する。よって、マグニチュードＭは次式［数４］から求められる。
【００１５】
【数４】
Ｅ（ｔ）＝ａ１０M ／｛Ｄ−ｒ（ｔ）｝b
ただし、Ｅ（ｔ）：震源時からの経過時間ｔにおける電界強度
ａ，ｂ：地質及び電界観測システムに依存する常数
Ｄ：震源距離
ｒ（ｔ）：時刻ｔにおける震源域とＰ波の面間の距離
＝ｔ×Ｐ波の速度
上式［数４］は、Ｐ波の到達前でも、震源時（電界上昇開始時刻）からの経過時間ｔと電界の振幅（電界強度）Ｅの関係から、マグニチュードＭと震源距離Ｄは予測できることを示している。即ち、マグニチュードＭと震源距離Ｄを未知数とする上式［数４］は、時間がｔのときの電界強度Ｅ（ｔ）の計測値データが２個以上得られれば解け、マグニチュードＭと震源距離Ｄは求められる。よって、地震発生後、サンプリング周期を０．１秒とすると５秒で１０個の計測値データが得られるので、数秒で、Ｐ波の到来前に、地震動は予測でき、上記従来の例の場合で、約１０秒は短縮するこができ、大部分の地域に対し地震波の到来前の予報が可能になる。
【００１６】
次に、観測装置の構成およびその処理を説明する。図２は本発明に係る電界観測に基づく震源時・震源距離・規模の予測装置の構成例を示す図であり、２１は電極、２２は受信機、２３は地震計、２４はデータ収集部、２５はデータ蓄積部、２６はデータ解析部、２７はペンレコーダを示す。
【００１７】
図２において、電極２１は、電極間の電位差を出力するものである。受信機２２は、電極２１の出力電圧の内の、所要帯域外の出力を減少させると共に、所要帯域内の出力を所要電圧（通常は１０Ｖ以下）まで増幅するものである。地震計２３は、地震波による加速度等の振幅を電圧に変換して出力するものである。データ収集部２４は、所定のサンプリング周期（０，１秒以下）で、受信機２２及び地震計２３の計測値データを収集するものである。データ蓄積部２５は、その収集した計測値データについて有意・無意を判断したうえで、記憶手段に蓄積するものである。
【００１８】
データ解析部２６は、データ蓄積部２５の記憶手段に収集し蓄積した計測値データを解析して地震と同時に発生した電界を検出し、電界上昇開始時刻から震源時を予測し、震源時からの経過時間と計測値データから震源距離と規模を予測して出力する。ペンレコーダ２７は、電極２１、受信機２２及び地震計２３の計測状況を監視するものである。地震計２３の計測値データを用いる場合には、Ｐ波を検知して、Ｐ波の計測値データと電界強度の計測値データから震源距離と規模を予測する。
【００１９】
次に、観測例によって、観測データの処理法を具体的に示す。図３は地震と同時に発生する電界及び該地震のＰ波の観測例を示す図であり、下記の２００２年６月１４日の地震の際の観測例で、該地震の震源時・震源距離・規模を予測するための時刻関係が示されている。３１は地震と同時に発生した電界、３２は南北方向の加速度、３３は東西方向の加速度、３４は上下方向の加速度、３５は震源時、３６はＰ波の到達時刻、３７は電界の特異点、３８はＳ波到達による電界パルスの上昇開始時刻を示す。
【００２０】
上記観測例によれば、地震と同時に発生した電界３１は、地震波が観測されなくなるまで継続した。南北方向の加速度３２、東西方向の加速度３３、上下方向の加速度３４は、何れも電界の観測点（波崎）の東南１４ｋｍの銚子での強震計の記録である。震源時３５（計算値）は、電界上昇開始時刻、１１時４２分４９．０秒と計算誤差の範囲内で一致している。Ｐ波の到達時刻３６は、電界観測点（波崎）での高感度地震計での観測値である。電界の特異点３７は、電界が増加から減少に急変する点である。Ｓ波到達による電界パルスの上昇開始時刻３８は、図２では、Ｓ波の到達時刻と一致しているが、電界の観測点（波崎）との距離差が１１ｋｍある銚子での観測値であり、Ｓ波の波崎への実際の到達時刻は、上昇開始時刻３８の１秒前である。この１秒の差は、Ｓ波による横波の誘起電圧の位相遅れによる。
【００２１】
上記の観測例の場合、震源距離９８ｋｍは、Ｐ波到達時刻３６と震源時３５の差１３．７秒にＰ波の伝播速度７．２ｋｍ／ｓを掛けて得られる。規模は、従来技術により、震源距離と地震波の振幅との関係から得られる他、Ｐ波到達時刻３６と電界の特異点３７との差からも下記の方法で得られる。本観測例の場合は、この差は０．８秒であり、これにＰ波の伝播速度を掛けると震源域の長さ５ｋｍが得られ、この値は規模５（通常、規模の誤差は０．５）の地震の震源域の長さに相当する。上記の他、下記の方法で規模は得られる。
【００２２】
電界の振幅（電界強度）は、規模と共に大きくなり、Ｐ波の波面（Ｗａｖｅ Ｆｒｏｎｔ）から電界観測点までの距離の二乗に反比例するため、Ｐ波の到達前でも、震源時３５（電界上昇開始時刻）からの経過時間と電界の振幅の関係から、規模と震源距離は予測できる。即ち、規模と震源距離を未知数とする上記の関係式は、震源時３５からの経過時間がｔのときの電界の振幅Ｅ（ｔ）の計測値が２個以上得られれば解けて、規模と震源距離は一意的に求められる。よって、地震発生後、数秒で、地震波の到来前に、規模・震源は予測できることになる。
【００２３】
本実施形態の電磁界センサー５による電界強度の計測は、自然雑音及び人工雑音が共に少ない場所で行うようにするのが望ましい。例えば上記観測例によれば、地震に伴う電界が０．４μＶ／ｍ程度であるので、最適観測帯とされる０．０１−３０Ｈｚ帯で、自然雑音及び人工雑音が常時０．１μＶ／ｍ程度以下の観測点とすることが望ましい。
【００２４】
なお、本発明は、上記実施の形態に限定されるものではなく、種々の変形が可能である。例えば上記実施の形態では、１観測地点から震源時・震源距離・規模を予測したが、震源距離から震源位置を求める手法として、従来のＰ波とＳ波から求める手法と同じ手法を採用してもよい。
【００２５】
【発明の効果】
以上の説明から明らかなように、本発明によれば、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法として、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、震源時からの経過時間と計測値データから震源距離と規模を予測し、また、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置として、電界強度を計測する計測手段と、計測手段による計測値データを収集・蓄積するデータ蓄積手段と、データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、震源時からの経過時間と計測値データから震源距離と規模を予測するデータ解析手段とを備えたので、電界強度の計測と解析により短時間に地震波の到来前に、震源・震源距離・規模を予測できる。
【００２６】
このため、本発明は、原子炉、高速鉄道、高速道路、空港などでの地震波到来前の緊急停止処置への利用が可能である。
【図面の簡単な説明】
【図１】 本発明に係る電界観測に基づく震源時・震源距離・規模の予測方法の実施の形態を説明するための図である。
【図２】 本発明に係る電界観測に基づく震源時・震源距離・規模の予測装置の構成例を示す図である。
【図３】 地震と同時に発生する電界及び該地震のＰ波の観測例を示す図である。
【符号の説明】
１…震源、２…地震波、３…Ｐ波面・帯電面、４…誘導電界、５…電磁界センサー、２１…電極、２２…受信機、２３…地震計、２４…データ収集部、２５…データ蓄積部、２６…データ解析部、２７…ペンレコーダ、３１…地霞と同時に発生した電界、３２…南北方向の加速度、３３…束西方向の加速度、３４…上下方向の加速度、３５…震源時、３６…Ｐ波の到達時刻、３７…電界の特異点、３８…Ｓ波到達による電界パルスの上昇開始時刻
【発明の属する技術分野】
本発明は、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法及びその予測装置に関する。
【０００２】
【従来の技術及び発明が解決しようとする課題】
従来技術による震源時・震源距離・規模の算出は、地震波のＰ波のみ、またはＰ波とＳ波の計測により行われている（例えば、特許文献１参照）。これらの値の算出には、Ｐ波のみの場合、４箇所以上で同じ地震のＰ波を計測する必要があり、Ｐ波とＳ波を用いる場合でも、３箇所以上で同じ地震のＰ波とＳ波を計測する必要がある。
【０００３】
このため、例えば４箇所のＰ波の観測点の震源距離が約１００ｋｍまたはそれ以内にある場合、震源時（地震発生時刻）・震源距離（震源と観測点間の距離）・規模（マグニチュードＭ）の算出には、震源時から少なくともＰ波が観測点に到達するまでの１５秒を要する。しかし、通常はＰ波とＳ波の両方を用いているので、３箇所の観測点の震源距離が約１００ｋｍまたはそれ以内の場合、Ｓ波が観測点に到達するまでの約２５秒を要する。どちらの方法にせよ、地震波が３箇所、または４箇所の観測点に到達するまで、上記の例の場合は、地震発生から１５ないし２５秒間は、震源距離・規模を知ることはできない。
【０００４】
【特許文献１】
特開平６−３２４１６０号公報
【０００５】
【課題を解決するための手段】
本発明は、上記課題を解決するものであって、簡単な方法及び装置により、地震発生とほぼ同時に震源時を検知する共に、短時間に地震波の到来前に、震源距離、規模を予測できるようにするものである。
【０００６】
そのために本発明は、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法であって、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測することを特徴とするものである。
【０００７】
また、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置であって、電界強度を計測する計測手段と、前記計測手段による計測値データを収集・蓄積するデータ蓄積手段と、前記データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測するデータ解析手段とを備えたことを特徴とするものである。
【０００８】
さらに、前記震源時からの経過時間ｔにおける電界強度Ｅ（ｔ）と該経過時間ｔの関係から次式［数３］を基に２個以上の計測値データにより震源距離Ｄと規模Ｍを予測し、
【０００９】
【数３】
Ｅ（ｔ）＝ａ１０M ／｛Ｄ−ｒ（ｔ）｝b
ただし、ａ，ｂ：地質及び電界観測システムに依存する常数
ｒ（ｔ）：時刻ｔにおける震源域とＰ波の面間の距離
＝ｔ×Ｐ波の速度
地震波を計測する地震計と併用して震源時・震源距離・規模を予測することを特徴とするものである。
【００１０】
【発明の実施の形態】
以下、本発明の実施の形態を図面を参照しつつ説明する。図１は本発明に係る電界観測に基づく震源時・震源距離・規模の予測方法の実施の形態を説明するための図であり、１は震源、２は震源から伝播する地震波、３は地震波の波面（Ｗａｖｅ Ｆｒｏｎｔ）であり、Ｐ波の波面であると同時に帯電面でもある。４は帯電面により誘起される誘導電界、５は電界強度を計測する電磁界センサーを示す。
【００１１】
図１において、地震波２は、地震発生と同時に震源１から発生し伝播して行く。この際、Ｐ波の速度が一番速いため、地震波の波面（Ｗａｖｅ Ｆｒｏｎｔ）はＰ波の波面となる。Ｐ波は縦波であり、その波面は帯電面３となる。帯電面３は電界を誘起し誘導電界４が伝播する。
【００１２】
誘導電界４は、後述の観測例に示すように地震と同時に発生し、Ｐ波の伝播と同期して強度が増大し、Ｐ波の到達時に極大となる。そして、Ｐ波の到達時からＳ波の到達まで、電界は減少し、Ｓ波の伝播と同期した電界は観測されない。また、Ｐ波による電界パルスは観測されず、Ｓ波の到達時に電界パルスが発生する。
【００１３】
この誘導電界４の強度を電磁界センサー５により計測し、解析することにより震源時は、電界上昇開始時刻から検知し、前記震源距離及び規模（マグニチュードＭ）は、震源時からの経過時間と電磁界センサー５による計測値データから予測する。
【００１４】
Ｐ波の伝播と同期した誘導電界４の波形は、Ｐ波の波面（Ｗａｖｅ Ｆｒｏｎｔ）に電荷が発生するとすれば、その誘導電界で定量的に説明できる。観測において、振幅は、Ｐ波の波面から電界観測点までの距離の二乗にほぼ反比例している。マグニチュードＭの定義として、地震動の最大速度振幅の常用対数を用いるものが多い。また、流動電位及びピエゾ電荷は最大速度振輻に比例するから、電界が流動電位またはピエゾ電荷によるものならば、電界は１０のＭ乗に比例する。よって、マグニチュードＭは次式［数４］から求められる。
【００１５】
【数４】
Ｅ（ｔ）＝ａ１０M ／｛Ｄ−ｒ（ｔ）｝b
ただし、Ｅ（ｔ）：震源時からの経過時間ｔにおける電界強度
ａ，ｂ：地質及び電界観測システムに依存する常数
Ｄ：震源距離
ｒ（ｔ）：時刻ｔにおける震源域とＰ波の面間の距離
＝ｔ×Ｐ波の速度
上式［数４］は、Ｐ波の到達前でも、震源時（電界上昇開始時刻）からの経過時間ｔと電界の振幅（電界強度）Ｅの関係から、マグニチュードＭと震源距離Ｄは予測できることを示している。即ち、マグニチュードＭと震源距離Ｄを未知数とする上式［数４］は、時間がｔのときの電界強度Ｅ（ｔ）の計測値データが２個以上得られれば解け、マグニチュードＭと震源距離Ｄは求められる。よって、地震発生後、サンプリング周期を０．１秒とすると５秒で１０個の計測値データが得られるので、数秒で、Ｐ波の到来前に、地震動は予測でき、上記従来の例の場合で、約１０秒は短縮するこができ、大部分の地域に対し地震波の到来前の予報が可能になる。
【００１６】
次に、観測装置の構成およびその処理を説明する。図２は本発明に係る電界観測に基づく震源時・震源距離・規模の予測装置の構成例を示す図であり、２１は電極、２２は受信機、２３は地震計、２４はデータ収集部、２５はデータ蓄積部、２６はデータ解析部、２７はペンレコーダを示す。
【００１７】
図２において、電極２１は、電極間の電位差を出力するものである。受信機２２は、電極２１の出力電圧の内の、所要帯域外の出力を減少させると共に、所要帯域内の出力を所要電圧（通常は１０Ｖ以下）まで増幅するものである。地震計２３は、地震波による加速度等の振幅を電圧に変換して出力するものである。データ収集部２４は、所定のサンプリング周期（０，１秒以下）で、受信機２２及び地震計２３の計測値データを収集するものである。データ蓄積部２５は、その収集した計測値データについて有意・無意を判断したうえで、記憶手段に蓄積するものである。
【００１８】
データ解析部２６は、データ蓄積部２５の記憶手段に収集し蓄積した計測値データを解析して地震と同時に発生した電界を検出し、電界上昇開始時刻から震源時を予測し、震源時からの経過時間と計測値データから震源距離と規模を予測して出力する。ペンレコーダ２７は、電極２１、受信機２２及び地震計２３の計測状況を監視するものである。地震計２３の計測値データを用いる場合には、Ｐ波を検知して、Ｐ波の計測値データと電界強度の計測値データから震源距離と規模を予測する。
【００１９】
次に、観測例によって、観測データの処理法を具体的に示す。図３は地震と同時に発生する電界及び該地震のＰ波の観測例を示す図であり、下記の２００２年６月１４日の地震の際の観測例で、該地震の震源時・震源距離・規模を予測するための時刻関係が示されている。３１は地震と同時に発生した電界、３２は南北方向の加速度、３３は東西方向の加速度、３４は上下方向の加速度、３５は震源時、３６はＰ波の到達時刻、３７は電界の特異点、３８はＳ波到達による電界パルスの上昇開始時刻を示す。
【００２０】
【００２１】
上記の観測例の場合、震源距離９８ｋｍは、Ｐ波到達時刻３６と震源時３５の差１３．７秒にＰ波の伝播速度７．２ｋｍ／ｓを掛けて得られる。規模は、従来技術により、震源距離と地震波の振幅との関係から得られる他、Ｐ波到達時刻３６と電界の特異点３７との差からも下記の方法で得られる。本観測例の場合は、この差は０．８秒であり、これにＰ波の伝播速度を掛けると震源域の長さ５ｋｍが得られ、この値は規模５（通常、規模の誤差は０．５）の地震の震源域の長さに相当する。上記の他、下記の方法で規模は得られる。
【００２２】
電界の振幅（電界強度）は、規模と共に大きくなり、Ｐ波の波面（Ｗａｖｅ Ｆｒｏｎｔ）から電界観測点までの距離の二乗に反比例するため、Ｐ波の到達前でも、震源時３５（電界上昇開始時刻）からの経過時間と電界の振幅の関係から、規模と震源距離は予測できる。即ち、規模と震源距離を未知数とする上記の関係式は、震源時３５からの経過時間がｔのときの電界の振幅Ｅ（ｔ）の計測値が２個以上得られれば解けて、規模と震源距離は一意的に求められる。よって、地震発生後、数秒で、地震波の到来前に、規模・震源は予測できることになる。
【００２３】
本実施形態の電磁界センサー５による電界強度の計測は、自然雑音及び人工雑音が共に少ない場所で行うようにするのが望ましい。例えば上記観測例によれば、地震に伴う電界が０．４μＶ／ｍ程度であるので、最適観測帯とされる０．０１−３０Ｈｚ帯で、自然雑音及び人工雑音が常時０．１μＶ／ｍ程度以下の観測点とすることが望ましい。
【００２４】
なお、本発明は、上記実施の形態に限定されるものではなく、種々の変形が可能である。例えば上記実施の形態では、１観測地点から震源時・震源距離・規模を予測したが、震源距離から震源位置を求める手法として、従来のＰ波とＳ波から求める手法と同じ手法を採用してもよい。
【００２５】
【発明の効果】
以上の説明から明らかなように、本発明によれば、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法として、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、震源時からの経過時間と計測値データから震源距離と規模を予測し、また、地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置として、電界強度を計測する計測手段と、計測手段による計測値データを収集・蓄積するデータ蓄積手段と、データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、震源時からの経過時間と計測値データから震源距離と規模を予測するデータ解析手段とを備えたので、電界強度の計測と解析により短時間に地震波の到来前に、震源・震源距離・規模を予測できる。
【００２６】
このため、本発明は、原子炉、高速鉄道、高速道路、空港などでの地震波到来前の緊急停止処置への利用が可能である。
【図面の簡単な説明】
【図１】 本発明に係る電界観測に基づく震源時・震源距離・規模の予測方法の実施の形態を説明するための図である。
【図２】 本発明に係る電界観測に基づく震源時・震源距離・規模の予測装置の構成例を示す図である。
【図３】 地震と同時に発生する電界及び該地震のＰ波の観測例を示す図である。
【符号の説明】
１…震源、２…地震波、３…Ｐ波面・帯電面、４…誘導電界、５…電磁界センサー、２１…電極、２２…受信機、２３…地震計、２４…データ収集部、２５…データ蓄積部、２６…データ解析部、２７…ペンレコーダ、３１…地霞と同時に発生した電界、３２…南北方向の加速度、３３…束西方向の加速度、３４…上下方向の加速度、３５…震源時、３６…Ｐ波の到達時刻、３７…電界の特異点、３８…Ｓ波到達による電界パルスの上昇開始時刻
Claims (4)
- 地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する方法であって、電界強度の計測値データを収集・蓄積し、電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間ｔにおける電界強度Ｅ（ｔ）と該経過時間ｔの関係から次式［数１］を基に２個以上の計測値データにより震源距離Ｄと規模Ｍを予測することを特徴とする電界観測に基づく震源時・震源距離・規模の予測方法。
ｒ（ｔ）：時刻ｔにおける震源域とＰ波の面間の距離
＝ｔ×Ｐ波の速度 - 地震波を計測する地震計と併用して震源時・震源距離・規模を予測することを特徴とする請求項１記載の電界観測に基づく震源時・震源距離・規模の予測方法。
- 地震動に伴い観測される電界強度の計測値データの時系列から震源時・震源距離・規模を予測する装置であって、電界強度を計測する計測手段と、前記計測手段による計測値データを収集・蓄積するデータ蓄積手段と、前記データ蓄積手段に蓄積された計測値データを解析して電界上昇開始時刻から震源時を検知し、前記震源時からの経過時間と前記計測値データから震源距離と規模を予測するデータ解析手段とを備え、前記データ解析手段は、前記震源時からの経過時間ｔにおける電界強度Ｅ（ｔ）と該経過時間ｔの関係から次式［数２］を基に２個以上の計測値データにより震源距離Ｄと規模Ｍを予測することを特徴とする電界観測に基づく震源時・震源距離・規模の予測装置。
ｒ（ｔ）：時刻ｔにおける震源域とＰ波の面間の距離
＝ｔ×Ｐ波の速度 - さらに地震波を計測する地震計を備え、前記データ解析手段は、前記地震計による計測データと併用して震源時・震源距離・規模を予測することを特徴とする請求項３記載の電界観測に基づく震源時・震源距離・規模の予測装置。
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 |
|---|---|---|---|---|
| EP1718996A1 (en) * | 2004-02-26 | 2006-11-08 | ExxonMobil Upstream Research Company | 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 |
|---|---|
| US6885945B2 (en) | 2005-04-26 |
| JP2004157094A (ja) | 2004-06-03 |
| US20040093162A1 (en) | 2004-05-13 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Meng et al. | Detecting bridge dynamics with GPS and triaxial accelerometers | |
| 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 | |
| WO2023087890A1 (zh) | 一种基于动应力、振动和oma综合分析判断构架模态共振方法 | |
| KR101951028B1 (ko) | 진도추정과 설비/구조물 피해예측기능을 가진 mems 기반 지진계측 장치 | |
| JP4171800B2 (ja) | 電界観測に基づく震源時・震源距離・規模の予測方法及びその予測装置 | |
| CN109541690B (zh) | 一种浅层介质结构面松散程度评价方法 | |
| WO2012052770A2 (en) | Methods and apparatus for geophysical prospecting to detect bodies of fluids in underground formations | |
| AU655281B2 (en) | Improvements in and relating to particle detection and analysis | |
| Lee et al. | Structural damage detection by power spectral density estimation using output-only measurement | |
| CN104483389B (zh) | 基于组合震源法检测桥梁预应力管道注浆质量的方法 | |
| CN107192763B (zh) | 利用运行列车为激励的钢轨振动衰减率的测试方法 | |
| CN107128329B (zh) | 一种基于应变量测反演变形加速度响应的轨距动态变化监测装置与设计方法 | |
| CN109855771A (zh) | 一种基于钢轨竖向加速度功率谱密度检测温度力的方法 | |
| JPS63186122A (ja) | 構造物の異常診断方式 | |
| CN114537466B (zh) | 高速铁路轨道结构变形损伤监测系统及方法 | |
| EP1769267B1 (en) | Movement detection system and method | |
| CN106053602A (zh) | 一种基于磁致伸缩效应的自比式锚杆无损检测方法 | |
| 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) | 建物の動的耐震性能評価システム |
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 |