JP5798545B2 - 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 - Google Patents
大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 Download PDFInfo
- Publication number
- JP5798545B2 JP5798545B2 JP2012243002A JP2012243002A JP5798545B2 JP 5798545 B2 JP5798545 B2 JP 5798545B2 JP 2012243002 A JP2012243002 A JP 2012243002A JP 2012243002 A JP2012243002 A JP 2012243002A JP 5798545 B2 JP5798545 B2 JP 5798545B2
- Authority
- JP
- Japan
- Prior art keywords
- time
- displacement
- acceleration
- dep
- amplitude
- 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.)
- Active
Links
- 238000000034 method Methods 0.000 title claims description 108
- 238000006073 displacement reaction Methods 0.000 claims description 724
- 230000001133 acceleration Effects 0.000 claims description 543
- 238000001514 detection method Methods 0.000 claims description 434
- 230000033001 locomotion Effects 0.000 claims description 316
- 230000008859 change Effects 0.000 claims description 195
- 238000006243 chemical reaction Methods 0.000 claims description 131
- 230000006870 function Effects 0.000 claims description 121
- 230000002159 abnormal effect Effects 0.000 claims description 56
- 230000008569 process Effects 0.000 claims description 46
- 230000001186 cumulative effect Effects 0.000 claims description 41
- 238000009499 grossing Methods 0.000 claims description 29
- 230000007246 mechanism Effects 0.000 claims description 22
- 238000010606 normalization Methods 0.000 claims description 18
- 230000011664 signaling Effects 0.000 claims description 15
- 230000002441 reversible effect Effects 0.000 claims description 13
- 230000009467 reduction Effects 0.000 claims description 10
- 230000006378 damage Effects 0.000 claims description 4
- 238000009826 distribution Methods 0.000 description 64
- 230000007704 transition Effects 0.000 description 58
- 239000002243 precursor Substances 0.000 description 49
- 230000000694 effects Effects 0.000 description 46
- 238000010586 diagram Methods 0.000 description 38
- 238000012795 verification Methods 0.000 description 35
- 238000009825 accumulation Methods 0.000 description 33
- 230000000704 physical effect Effects 0.000 description 32
- 239000002245 particle Substances 0.000 description 31
- 230000014509 gene expression Effects 0.000 description 28
- 230000000737 periodic effect Effects 0.000 description 28
- 230000007423 decrease Effects 0.000 description 22
- 238000005516 engineering process Methods 0.000 description 21
- 238000007792 addition Methods 0.000 description 19
- 238000005293 physical law Methods 0.000 description 17
- 238000004458 analytical method Methods 0.000 description 16
- 238000012545 processing Methods 0.000 description 15
- 230000004069 differentiation Effects 0.000 description 13
- 230000006872 improvement Effects 0.000 description 11
- 230000003111 delayed effect Effects 0.000 description 10
- 238000001914 filtration Methods 0.000 description 10
- 238000004873 anchoring Methods 0.000 description 9
- 239000000284 extract Substances 0.000 description 9
- 230000003068 static effect Effects 0.000 description 9
- 230000002123 temporal effect Effects 0.000 description 9
- 239000000463 material Substances 0.000 description 8
- 238000012544 monitoring process Methods 0.000 description 8
- 238000012935 Averaging Methods 0.000 description 7
- 238000000605 extraction Methods 0.000 description 6
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 description 6
- 229910052753 mercury Inorganic materials 0.000 description 6
- 230000002547 anomalous effect Effects 0.000 description 5
- 230000002829 reductive effect Effects 0.000 description 5
- 230000003595 spectral effect Effects 0.000 description 5
- 208000012661 Dyskinesia Diseases 0.000 description 4
- 230000007774 longterm Effects 0.000 description 4
- 230000002265 prevention Effects 0.000 description 4
- 239000013598 vector Substances 0.000 description 4
- 102100027731 Endogenous retrovirus group K member 16 Rec protein Human genes 0.000 description 3
- 101000580913 Homo sapiens Endogenous retrovirus group K member 16 Rec protein Proteins 0.000 description 3
- 230000036772 blood pressure Effects 0.000 description 3
- 230000003247 decreasing effect Effects 0.000 description 3
- 230000010354 integration Effects 0.000 description 3
- 238000011002 quantification Methods 0.000 description 3
- 238000010183 spectrum analysis Methods 0.000 description 3
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 3
- 241000219307 Atriplex rosea Species 0.000 description 2
- 230000005856 abnormality Effects 0.000 description 2
- 238000013459 approach Methods 0.000 description 2
- 238000012790 confirmation Methods 0.000 description 2
- 238000011161 development Methods 0.000 description 2
- 230000035487 diastolic blood pressure Effects 0.000 description 2
- 230000007613 environmental effect Effects 0.000 description 2
- 230000010355 oscillation Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 238000001228 spectrum Methods 0.000 description 2
- 238000007619 statistical method Methods 0.000 description 2
- 230000005653 Brownian motion process Effects 0.000 description 1
- 206010044565 Tremor Diseases 0.000 description 1
- 230000002238 attenuated effect Effects 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 238000009530 blood pressure measurement Methods 0.000 description 1
- 238000005537 brownian motion Methods 0.000 description 1
- 244000309464 bull Species 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 230000008878 coupling Effects 0.000 description 1
- 238000010168 coupling process Methods 0.000 description 1
- 238000005859 coupling reaction Methods 0.000 description 1
- 125000004122 cyclic group Chemical group 0.000 description 1
- 238000013016 damping Methods 0.000 description 1
- 238000007405 data analysis Methods 0.000 description 1
- 238000009795 derivation Methods 0.000 description 1
- 230000001066 destructive effect Effects 0.000 description 1
- 239000010432 diamond Substances 0.000 description 1
- 238000005315 distribution function Methods 0.000 description 1
- 238000007667 floating Methods 0.000 description 1
- 239000012530 fluid Substances 0.000 description 1
- 230000005484 gravity Effects 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 230000005923 long-lasting effect Effects 0.000 description 1
- KFOPKOFKGJJEBW-ZSSYTAEJSA-N methyl 2-[(1s,7r,8s,9s,10r,13r,14s,17r)-1,7-dihydroxy-10,13-dimethyl-3-oxo-1,2,6,7,8,9,11,12,14,15,16,17-dodecahydrocyclopenta[a]phenanthren-17-yl]acetate Chemical compound C([C@H]1O)C2=CC(=O)C[C@H](O)[C@]2(C)[C@@H]2[C@@H]1[C@@H]1CC[C@H](CC(=O)OC)[C@@]1(C)CC2 KFOPKOFKGJJEBW-ZSSYTAEJSA-N 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000002093 peripheral effect Effects 0.000 description 1
- 238000003672 processing method Methods 0.000 description 1
- 238000004445 quantitative analysis Methods 0.000 description 1
- 230000008707 rearrangement Effects 0.000 description 1
- 230000000630 rising effect Effects 0.000 description 1
- 239000011435 rock Substances 0.000 description 1
- 230000001932 seasonal effect Effects 0.000 description 1
- 230000035939 shock Effects 0.000 description 1
- 238000004611 spectroscopical analysis Methods 0.000 description 1
- 230000008961 swelling Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
- 238000012360 testing method Methods 0.000 description 1
- 230000009466 transformation Effects 0.000 description 1
- 238000000844 transformation Methods 0.000 description 1
- 230000005428 wave function Effects 0.000 description 1
Images
Description
0番目(最初)の地震が、
(震源時=1983/01/06 09:35:47.50、LAT=北緯34.453度、LON=東経132.607度、DEP=27.0km、MAG=3.5)、
1番目の地震が、
(震源時=1983/01/08 10:36:54.60、LAT=北緯33.817度、LON=東経131.843度、DEP=83.0km、MAG=4.0)、
2番目の地震が、
(震源時=1983/01/17 12:19:49.10、LAT=北緯33.468度、LON=東経132.448度、DEP=48.0km、MAG=4.1)
であったとする。
0番目=(LAT=34.453、LON=132.607、DEP=27.0、INT=X、MAG=3.5)、
1番目=(LAT=33.817、LON=131.843、DEP=83.0、INT=49.0186、MAG=4.0)、
2番目=(LAT=33.468、LON=132.448、DEP=48.0、INT=217.7151、MAG=4.1)
となる。0番目のINT=Xは、その前の地震情報がないので震源時(発生 時刻)の差である間隔時間INTが算出できない事を意味する。従って、1番目が時系列データの始まりとなる。この様に作成する時系列の表記方法を次に述べる。
LATとすると、
[LAT]=[LAT(1)=33.817、LAT(2)=33.468、・・・]、
LONとすると
[LON]=[LON(1)=131.843、LON(2)=132.448、・・・]、
DEPとすると、
[DEP]=[DEP(1)=83.0、DEP(2)=48.0、・・・]、
INTとすると、
[INT]=[INT(1)=49.0186、INT(2)=217.7151、・・・]、
MAGとすると、
[MAG]=[MAG(1)=4.0、MAG(2)=4.1、・・・]、
と表記される。地震の発生順序を示すインデックスjの1、2、・・・が各時系列の時間tとなる。従って、5次元空間(LAT、LON、DEP、INT、MAG)に描かれた地震発生の軌跡の時刻tの各座標軸成分は、LAT(t)、LON(t)、INT(t)、DEP(t)、MAG(t)となる。これら各成分で定められた位置ベクトルは、抽出する地震を選択するために設定した各しきい値、もしくは所定の各基準値からの変位ベクトルともみなす事ができる。従って、仮想粒子の運動が描く軌跡の時刻t(t=j)における各位置成分は、変位成分と同等に取り扱えるので、成分cがLATであれば、その変位D(c,t)を、D(LAT,t)と表記するが、位置ベクトルのLAT成分のLAT(t)と同等である。これら震源要素cの時系列[c]を、次の[数1]で表記する。
(1)大地震発生の時刻(震源時)の予知
(2)マグニチュードの予知
(3)加速度A(DEP,t)とA(INT,t)の波形の反転が、[図3]の様に、略、完全な反転とならない場合、CQKとCQT予兆の推移を如何に確定し予知に利用するか
(4)巨大地震の予知
(5)大地震、巨大地震につながる地殻の蓄積応力の臨界状態を如何にモニターし、大地震、巨大地震の予知に利用するか
これら項目に関する理論及び技術的な背景を順に説明する。
[特許文献1]は、地震発生の2次の変化率である加速度A(c,t)を、[数1]の時系列[c]から、[図6]の差分間隔を、n=k=sと等しくした2次差分検出ウエーブレットを用いて検出した。この時、加速度A(c,t)は、各震源要素間に共通な略2sの短い固有な周期(時間的な尺度)を持って変動するが、大地震発生直前には、通常の地震発生とは異なる2種類のA(INT,t)とA(DEP,t)とA(MAG,t)との振幅及び位相関係が成立する。この2つの特殊な関係が、[図3]に表示した大地震発生のCQK及びCQT予兆であった。この図で、大地震が発生した時刻は、それら予兆を検出した時刻ta1のA(INT,ta1)の振幅のピークが、半周期後に、反転し、逆のピークに到達した時刻ta2である。[特許文献1]は、加速度A(INT,t)が、略2sの周期を持った変動をしているので、予兆検出時刻ta1の半周期前のta0に2sを加算するか、予兆検出時刻ta1に半周期sを加算して大地震発生時刻のta2を予知した。しかし、加速度A(INT,t)の周期は、略2sであるが、振動毎に、変動する。従って、その様な場合も考量し、予兆検出時刻までの前半の半周期を、ta1-ta0とし、この前半の半周期が、後半の半周期と等しいと仮定して、大地震の発生時刻ta2を予知した。しかし、この後半の半周期が、前半と異なる場合、大地震発生の時刻(震源時)の予知には改善が求められる。この改善に、[数44]の変位D(c,t)と[数47]の加速度A(c,t)の間に成立する[数43]のニュートンの運動の第二法則を利用する。
(1)時刻tx2(tx2 > ta2)のD(c,tx2)が、最新の変位となる。
(2)時刻tx1(=tx2-s)のA(c,tx1)が最新の加速度となる。
(3)[数43]の運動の第二法則により最新の加速 度A(c,tx1)が作用している変位は、時刻tx1の変位D(c,tx1)であるが、加速度A(c,tx1)の算出には、変位D(c,tx2)も含まれる。
(4)[数1]の時系列[c]の最後(最新)の変位d(c,m)を与える最新の時刻mは、m=tx2+wとなる。
(5)大地震が発生する時刻は、震源要素cをINTとした加速度A(INT,t)が谷(正の方向のピーク)となる時刻ta2となる。
この大地震発生時刻ta2は、[図3]のCQT予兆の大地震が発生する時刻に対応する。 ニュートンの運動の第二法則から、時刻ta2の加速度A(INT,ta2)は、変位D(INT,ta2)に作用し、変位D(INT,ta2)は、その振動の山(負の方向のピーク)に在る。[図3]のCQK予兆の大地震発生の場合は、大地震発生時刻ta2で、[図8]のA(c,ta2)とD(c,ta2)の位相関係が、それぞれ反転する。
[特許文献1]の発明における、マグニチュードの程度の予知に関しては、[非特許文献11]にもあるように、地震発生の変化に、地殻の応力場の変化が反映されているというよく知られた自然法則を利用している。具体例としては、北緯32度〜36度、東経131.5度〜136.5度で囲まれた中国・四国・近畿地方に発生したマグニチュードが3.5以上の地震の発生時刻の時間間隔INTの時系列d(INT,m)を利用する。INTのランダム的な変動を平滑するために、例えば、 任意な個数、2s個(2s=60〜70)の移動平均、もしくは累積加算を取る。時間に関する固有な尺度(60〜70イベント)は、大地震の予兆を検出する加速度A(c,t)の変動周期に相当する。時刻mのINTの2s個の累積値を、CI(m,2s)と表記する。CI(m,2s)は、2s個の地震発生に、費やした時間なので、CI(m,2s)が、増加すれば、単位時間当たりの地震の数が、減少し、地震活動が静穏になり、応力の蓄積が進む。従って、CI(m,2s)を利用すれば、大地震に至る応力の 蓄積が如何に推移するかを観察する事が出来る。例えば、上記中国・四国・近畿地方の場合、次の3つの段階に分類できる。まず1段階目は、地震の数年前から観測される静穏の始まり、次の2段階目は、1年から半年程前に観測される蓄積された 応力がピークに到達し応力の蓄積が減少し始めたと判断される臨界静穏状態、そして、3段階目は、応力の蓄積が急激に減少し、大地震の発生数日〜数時間前まで急減少が継続する段階に分類できる。すると、臨界静穏状態に対応する形状の振幅値、大地震の発生数日までに観測される形状全体の幅やその面積の大きさが、蓄積応力が作用する断層形状(長さや幅)や断層面積に比例していると仮定できる。又、後述するように、CI(m,2s)は、地殻に蓄積され開放される歪エネルギーの密度に比例する量となるので、その形状の総面積が、大地震により解放される総歪エネルギーと仮定できる。従って、過去に発生した大地震のM6クラスとかM7クラスとかのマグニチュードの程度とこれら形状の大きさとの比例関係から、今にも発生しそうな大地震のマグニチュードの程度の予知をする事ができる。一方、マグニチュードは、[非特許文献1]によると、MKS単位で、次の[数53]で与えられる。
[図3]のCQKとCQT予兆において、加速度A(DEP,t)とA(INT,t)とが、完全な逆位相とならない波形の反転も観測されている。例えば、[特許文献1]の改善の実施例として後述する2011年3月10日に発生した東北地方のCQKのM6.4地震や2011年7月5日に発生した和歌山県のCQKのM5.5地震等がある。これらの地震では、地震発生日は、[図3]のA(INT,t)が負のピーク値をとる時刻ta2ではなく、A(DEP,t)が、時刻ta2の正のピーク値から更に半周期進んだ負のピーク値を取る時刻となる。その時刻は、[図8]のA(c,t)とD(c,t)で震源要素cをDEPとした関係図を用いると、時刻ta3となる。この発生時刻のta2からta3への遅延は、上記CQK型の物理モデルによると、時刻ta1からta2の変化量ΔD(DEP,t)で生じるΔF(DEP,t)は、時刻ta2で今にも発生しそうな大地震の断層に作用するせん断応力を誘発する。しかし、そのせん断応力は、変化量ΔD(DEP,t)で生じるとされたΔF(DEP,t)よりも少ない量となり、変化量ΔD(DEP,t)に相当する断層幅を持った断層の静止摩擦応力より小さくなる。従って、そのCQK大地震は、時刻ta2で発生しない。時刻ta2からta3までの変化量ΔD(DEP,t)で生じたΔF(DEP,t)が誘発するせん断応力が更に加わると、断層に作用するせん断応力は、断層面の静止摩擦応力より大きくなり、そのCQK大地震が時刻ta3で発生すると推論される。ただし、この場合でも、時刻ta1からta2までの変化量ΔD(DEP,t)が、断層幅Wに略等しくなる。従って、断層幅Wを持った断層に作用するせん断応力を生み出す源となるΔF(DEP,t)を発生させる物理モデルは、CQKとCQT予兆の加速度A(DEP,t)とA(INT,t)の波形の反転が、[図3]に図示した理想的な反転と異なる場合、[特許文献1]の大地震発生時刻と震源の予知技術を改善する。
マグニチュードが8.5より大きな巨大地震にも[数55]や[数56]が成立すると仮定すると、巨大地震の断層長Lの予測が、巨大地震のマグニチュードの予知となる。従って、断層長Lの予測をするために、[特許文献1]に用いた変位、1次差分、2次差分検出ウエーブレットを上記2つの物理的性質を満足する各検出ウエーブレットに変更して、巨大地震発生の自然現象を定量化し、巨大地震の予知に利用する。
太平洋に面する東海岸側の下には、太平洋プレートと固着している領域がある。その固着域は、[図11]に茶色で示され、その下の矢印は、太平洋プレートの西方向への運動方向を示し、固着域もその方向へ動いていた。固着域の西方向への運動が、太平洋に面する東海岸側を引きずり込んでいた(沈下させていた)。一方、大陸プレートの東方向へ矢印で示した 運動は、日本海に面する西海岸側を押し上げていた(隆起させていた)。このスローな変形は、長年にわたり生じた変形で、実線からの通常な変形として破線で示されている。
(2)[図12]に描いた予兆的膨らみの発生と発生後の太平洋プレートの西方向への急加速と急停止とからなる異常加速運動
2009年12月8日頃から太平洋に面する東海岸側(大陸プレートの東端)と日本海に面する西海岸側に地殻の膨らみが観察された。大陸プレートの東方向への「押し」に対し、太平洋プレートの西方向への「押し」が強まったプレートの相対運動による膨らみは、「予兆的な膨らみ」として、赤の破線で表示されている。その膨らみは、東海岸側で1mm程、西海岸側で3mm程あった。又、この膨らみは、本州の東海岸側(四国の太平洋側も含む)と日本海岸側でも観察されている箇所がある。又、その「膨らみ」を観測した後、太平洋プレートの赤色の破線矢印で示されている西方向への異常運動が観測された。2010年7月8日頃から、その西方向への運動が加速され始め、2010年12月22日頃に通常の約3倍の西方向への移動速度に達した。その直後から、急減速し、2011年1月27日頃までに急停止した。そして、太平洋プレートの西方向への運動は、3月11日の巨大地震と津波発生まで停止し続けていた。
(3)[図13]に描いたプレート境界の固着領域の壊れ、誘起されたメガスラストが巨大地震と津波を発生させた。
大陸プレートの東端に位置する東日本の地殻の膨らみを支えていた固着域は、太平洋プレートの西方向の運動に耐え切れなくなり、[実施例]で後述する様に、先ず浅い箇所が破壊され、深い箇所が一挙に壊れ、プレート境界で赤色の破線矢印で示す方向にメガスラストを誘起し、東方向の断層運動をし、巨大地震と津波とを発生させた。壊れた固着域は、丸印の破線で表示されている。
[特許文献1]によると、地震発生の変化に、地殻の応力場の変化が反映されているという自然法則を定量化した具体例は、北緯32度〜36度、東経131.5度〜136.5度で囲まれた中国・四国・近畿地方に発生したマグニチュードが、3.5以上の地震の発生時刻の時間間隔d(INT,m)を、時刻mまで2s個、累積したCI(m,2s)であった。CI(m,2s)は、2s個の地震発生に、費やした時間なので、単位時間当たりの地震の数が減少すると、CI(m,2s)は増加し、領域の地殻に応力が蓄積されている状況を示す。逆に、大地震が発生し余震の多発等で地震発生数が増加すれば、蓄積された応力が大地震の発生により開放されている状況を示す。従って、CI(m,2s)は、大地震に至るまでの応力の蓄積と大地震発生後の応力開放のサイクルの推移を表す。このサイクルの形状が、[特許文献1]の発明では、マグニチュードの程度の予知に利用された。しかし、CI(m,2s)は、地殻に蓄積され開放される歪エネルギーの密度に比例する量である事も、本願の発明者が、発見したので、その自然法則を次に述べる。
ステップ1:
任意選択した時刻jの値<CI(j,2s),1> (平均個数は1個なので、CI(j,2s)と同一値)と最も近接した値を、Δ<CI(k-j,2s),1>〜0の条件で探す。
ステップ2:
その値が、時刻kの値<CI(k,2s),1> = CI(k,2s)であったとする。
ステップ3:
所定のしきい値THを、例えば、TH=10とかTH=15とかに設定する。設定に関して、時系列[CI(m,2s)]の標準偏差等の統計量を考慮した値でも良い。
ステップ4:
時刻kを時刻j以外の全ての時刻kxに変え、Δ<CI(kx-j,2s),1>がTHより大きくなる場合の数を数える。その総数をN1とし、N1をN1/N1で規格化する。従って、r=1の場合、規格化した総数は1となる。
なお、時系列<CI(j,2s),1>は、変動しているので、変動が激しくなると、N1の値は、増加し、時系列がランダムであれば、N1は最大となる。
ステップ5:
時刻j+nとk+nに変位検出ウエーブレットを、それぞれ追加配置し、r=2とした時刻jとkの差Δ<CI(k-j,2s),2>を、所定のしきい値THと比較する。Δ<CI(k-j,2s),2>がTHより大であれば、時刻j+nと時刻k+nの変位成分が大きく寄与している事になる。この様にTHより大きくなる総数N2を、時刻kを時刻j以外の全ての時刻kxに変え、ペアーの差の絶対値Δ<CI(kx-j,2s),2>がTHより大となる回数を数え、総回数をN2と、その規格値をN2/N1とする。その値が、例えば、N2/N1〜0.5であったとする。N2は、2個の移動平均後の<CI(j,2s),2>から得るΔ<CI(k-j,2s),2>がTHより大となる総数なので、その平均操作により、N1より少ない数値を取る。N2/N1〜0.5の場合、時系列[CI(m,2s)]の変動を、50%支配している事になり、力学変数の数は、2となる。
ステップ6:
時刻j+2nとk+2nに変位検出ウエーブレットを、それぞれ追加配置したr=3の場合でも([図28]のr=4をr=3とした)Δ<CI(k-j,2s),3>が、THよりも大となると、時刻j+2nとk+2nの変位成分が、依然、寄与している事になる。上記ステップ5と同様のプロセスで、 時刻kを時刻j以外の全ての時刻kxに変えると、Δ<CI(k-j,2s),3>がTHより大となる条件に適合する総数N3は減少する。規格値がN3/N1〜0.2だったとすると、3番目の独立変数の時系列[CI(m,2s)]の変動への寄与は、20%程度で、力学変数の数は3となる。
ステップ7:
時刻j+3nとk+3nに変位検出ウエーブレットを追加配置する。上記ステップ6と同様なプロセスをr=4で実施する事になる。[図28]のΔ<CI(k-j,2s),4>が、THより大となる条件を満足する総数を数える。総数N4は激減する。その規格化値が、N4/N1〜0.05だったとする。4番目の力学変数の時系列[CI(m,2s)]への変動への寄与は、微小だが、依然存在し、力学変数の総数は4となる。
ステップ8:
時刻j+4nとk+4nに変位検出ウエーブレットを追加配置する。上記ステップ6、7と同様なプロセスをr=5で実施する。Δ<CI(k-j,2s),5>が、THより大となる条件を満足する総数がN5で、規格値は、r=4の場合と変わらず、N5/N1 〜0.05だったとする。5番目の独立変数の時系列[CI(m,2s)]への変動への寄与は、r=4の場合と同程度で存在し、力学変数の総数は5となる。
ステップ9:
時刻j+5nとk+5nに変位検出ウエーブレットを追加配置する。上記ステップ8と同様なプロセスをr=6で実施する。Δ<CI(k-j,2s),6>が、THより大となる条件を満足する総数をN6とする。その規格化値が、N6/N1〜0.05となり、r=4、r=5 の場合と同程度だったとする。6番目の独立変数の時系列[CI(m,2s)]への変動への寄与は、r=4、r=5 の場合と同程度で、力学変数の総数は6となる。
ステップ10:
時刻j+6nとk+6nに変位検出ウエーブレットを追加配置し、上記ステップ9と同様なプロセスをr=7で実施し、ステップ9と同様な規格値N7/N1〜0.05を得たとする。ステップ7でr=4の時と、規格値は変化してなく、この様な規格値を作り出す時系列[CI(m,2s)]の変動は、ランダム変動が力学変数に重複している状況を示す。最終的に力学変数の数は、3と断定できる。また、ランダム変動の寄与が無ければ、r=4以降の、規格値はN4/N1、N5/N1、N6/N1は、全てゼロとなり時系列[CI(m,2s)]に変動を与えている力学変数の数は3となる。時系列[CI(m,2s)]には、2sを適当な大きさ以上に取ればランダム変動の寄与は無い。
m = 555 (1994/04/06)の時、CI(m,70)max=70 x 283.0961時間、
m = 551 (1994/01/31)の時、CD(m,70)max=70 x 43.1942km
となる。なお、70倍した283.0961時間と43.1942kmは、d(INT,j)とd(DEP,j)の70個の移動平均値の最大値である。[図32]の左側の縦軸目盛は、[数60]をNCI(m,70)、[数62]をNCD(m,70)とラベル付けしたグラフの目盛である。又、右側の縦軸の6以上の目盛は、時 系列[d(MAG,m)]をMAGとラベル付けした矢印の高さのグラフの目盛である。従って、MAGのグラフは、マグニチュードである矢印の高さが6以上の大地震になると表示される。更に、その上にある 右側の縦軸目盛にLONとラベル付けされている範囲は、震源パラメータの時系列[d(LON,m)]を点グラフにしたLONの目盛で、経度が131.5度〜136.5度に相当する。大地震の余震や群発地震は、震源が略同じなので、それら余震や群発地震の震源の経度は殆んど変化しない。従って、[d(LON,m)]の点グラフの変動は無くなる。又、[図29]-[図31]に使用されたマグニチュードが6以上の大地震の発生年代ラベルは、そのまま[図32]の時系列[d(MAG,m)]に表記されている。又、2つの群発地震S1984とS2001は、時系列[d(LON,m)]の点グラフに矢印と共にラベルが付けられている。
(2)マグニチュードの予知
(3)加速度A(DEP,t)とA(INT,t)の波形の反転が、[図3]の様に、略、完全な反転とならない場合、CQKとCQT予兆の推移を如何に確定し予知に利用するか
(4)巨大地震の予知
(5)大地震、巨大地震につながる地殻の蓄積応力の臨界状態を如何にモニターし大地震と巨大地震の予知に利用するか
[背景技術]の[0109]-[0112]で説明した課題内容をまとめる。[特許文献1]の震源時(大地震発生時刻)の予知は、次に述べる自然法則に基づいた。地殻の応力変化を反映する地震発生の変化を、[数1]の時系列[c]から、[図6]の差分間隔を、n=k=sと等しくした2次差分検出ウエーブレットを用いて、2次の変化率である加速度A(c,t)として抽出した。すると、加速度A(c,t)は、各震源要素間に共通な略2sの短い固有な周期(時間的な尺度)を持って変動するが、大地震発生直前には、通常の地震発生の変化とは異なるA(INT,t)とA(DEP,t)とA(MAG,t)との振幅及び位相関係が2種類成立した。この2つの特殊な関係が、[図3]に表示した大地震発生のCQK及びCQT予兆であった。この図で、大地震が発生する時刻は、それら予兆を検出したA(INT,ta1)の振幅のピークが、時刻ta1から半周期後に、反転し、逆のピークに到達した時刻ta2である。抽出した加速度A(INT,t)の周期は、略2sであるが、振動毎に、変動する。従って、[特許文献1]は、予兆検出時刻ta1までの前半の半周期を、ta1-ta0とし、この前半の半周期が、後半の半周期と等しいと仮定して、大地震の発生時刻ta2を予知するものであった。しかし、この後半の半周期が、前半と異なる場合、大地震発生の時刻(震源時)の予知には改善が求められる。この予知の方法、装置、及びプログラムに改善が求められる状況下、本願の発明者は、地震発生の変化を記述するニュートンの運動の第二法則を次の理論により確立した。
(1)先ず、時刻tx2(tx2>ta2)のD(c,tx2)を、最新の変位とする。
(2)すると、時刻tx1(=tx2-s)のA(c,tx1)が最新の加速度となる。
(3)この時、[数43]の運動の第二法則により最新の加速度A(c,tx1)が作用している変位は、時刻tx1の変位D(c,tx1)であるが、加速度A(c,tx1)の算出には、変位D(c,tx2)も含まれる。
(4)尚、[数1]の時系列[c]の最後(最新)の変位d(c,m)を与える最新の時刻mは、m=tx2+wとなる。
(5)予兆(CQT)を抽出直後、大地震が発生する時刻は、震源要素cをINTとした加速度A(INT,t)が谷(正の方向のピーク)となる時刻ta2となる。
この大地震が発生する[図8]の時刻ta2は、[図3]のCQT予兆の大地震が発生する時刻ta2に対応する。[図3]のCQK予兆の場合は、[図8]のA(c,t)とD(c,t)の位相関係を、それぞれ反転した関係となる。
[背景技術]の[0113]-[0126]で説明した課題内容をまとめる。[特許文献1]の大地震のマグニチュードの予知は、マグニチュードの程度を 予知するもので、次に述べる自然法則に基づいた。地震の発生時刻の時間間隔INTの時系列d(INT,m)の2s個の累積値を、時刻mのCI(m,2s)とすると、CI(m,2s)は、地殻に蓄積された歪エネルギー密度に比例する量となり、その形状の総面積が、大地震により解放される総歪エネルギーとなる。従って、過去に発生した大地震のマグニチュードの程度(例えば、M6クラスとかM7クラス)とこれら形状の大きさとの比例関係から、今にも発生しそうな大地震のマグニチュードの程度の予知をするものであった。この予知の方法、装置およびプログラムを、[数55]か[数56]を用いたより定量的な予知の方法、装置、及びプログラムに改善するには、大地震の予兆の推移をニュートンの運動の第二法則を用いて記述した大地震発生の自然法則を、確立し、大地震の予知に利用する事が望まれる。この状況下、本願の発明者は、先ず、微分操作が持つ物理的性質を正しく反映する[図7]の差分操作を確立した。その差分操作は、時刻τの変位D(c,τ)を[数45]、速度V(c,τ)を[数46]、加速度A(c,τ)を、[数47]もしくは、[数48]で与える。この時、[図7]の各検出ウエーブレットのΔt(=2w+1)、及び、差分を取る間隔n、k、sは、任意な大きさに設定できるので、様々な変動周波数を持つ粒子の運動(地震発生の変化)から、変位検出ウエーブレットの幅と差分を取る間隔で選択的に抽出した周波数領域にある運動を、変位(位置)D(c,τ)に作用する加速度(力)A(c,τ)として検出する。その周波数領域の選択機能は、変位、速度、 加速度を検出する物理的ウエーブレットの周波数特性式、[数48]-[数51]で与えられる。従って、本願の発明者は、大地震発生直前のCQK及びCQT予兆の推移に関する自然法則を、所定の周期で変動している震源要素c(c=DEP,INT,MAG)の変位D(c,τ)に作用する力A(c,τ)として、[数44]で与えられるニュートンの運動の第二法則で記述するに至った。更に、[図9]と[図10]とに図示する大地震の発生の物理モデルを確立させた。そのモデルを使用すると、CQK及びCQT予兆を検出後、[図10]を用い、今にも発生しそうな大地震の断層(長さLと幅Wの形状)に作用するせん断応力は、変化量ΔD(DEP,t)を起こした力ΔF(DEP,t)から誘発され、ΔD(DEP,t)が、その大地震の断層幅W kmと略等しくなる自然法則を、発見した。この法則を利用し、今にも発生しそうな大地震の断層の幅Wを予測し、[数54]に基づいたM6.8とかM7.1等のより定量的なマグニチュードの予知の方法、システム、プログラム、及び記録媒体を開発するに至った。
[背景技術]の[0127]で説明した課題内容をまとめる。この課題の場合、CQKとCQT予兆の大地震が発生する時刻は、[図8]のA(c,t)とD(c,t)で震源要素cをDEPとした関係図を用いると、時刻ta3となる。この発生時刻のta2からta3への遅延は、[図10]の物理モデルによると、時刻ta1からta2のΔD(DEP,t)で生じたΔF(DEP,t)は、時刻ta2で、今にも発生しそうな大地震の断層に作用するせん断応力を誘発するが、変化量ΔD(DEP,t)で生じるとされたΔF(DEP,t)よりも少ない量となり、変化量ΔD(DEP,t)に相当する断層幅を持った断層の静止摩擦応力より小さくなる。従って、そのCQK大地震は、時刻ta2で発生せず、時刻ta2からta3までの変化量ΔD(DEP,t)によるΔF(DEP,t)が更に加わり、断層に作用するせん断応力は、断層面の静止摩擦応力より大きくなり、そのCQK大地震が時刻ta3で発生すると推論される。ただし、この場合でも、時刻ta1からta2までの変化量ΔD(DEP,t)が、断層幅Wに略等しくなる。従って、[特許文献1]で利用した加速度A(DEP,t)とA(INT,t)の波形の反転が、[図3]の様に、略、完全な反転とならない場合、[図10]の大地震の物理モデルを利用すると、その発生時刻はA(INT,ta2)のta2ではなく、A(DEP,ta3)の時刻ta3となる、大地震の予知の方法、システム、プログラム、及び記録媒体を開発するに至った。尚、A(DEP,ta3)の時刻ta3を予知する時、第1の課題と同様に、差分間隔を取った時間だけ先に進んだ変位D(DEP,t)が、負か正の方向の振幅値のピークに到達した観測時刻を大地震発生の予知時刻としてもよい。
[背景技術]の[0128]-[0148]で説明した課題内容をまとめる。この課題は、新たに解決すべきものである。[特許文献1]の大地震の予知方法に関して、任意選択する 領域の広さは、約4〜5度の緯度幅と経度幅からなるメッシュであったが、2011年3月11日に発生した東日本巨大地震の断層長Lは、約500kmで、メッシュを構成する緯度差に換算すると約4.5度となる。巨大地震の発生場所(領域)は、過去に繰り返し発生したほぼ直線的な広範囲のプレート境界領域に限られる。従って、広範囲にわたるプレート境界で発生する巨大地震の予知には、「プレート境界の地殻表面の変化に出現する異常が、巨大地震の予兆となり、その予兆の広がりを、今にも発生しそうな巨大地震の断層長(L)とする」自然現象を利用した巨大地震の予知が求められる。地殻 表面の変化に出現する異常を、巨大地震の予兆として検出するために、その領域に設置した多数のGPSステーションの位置を地殻表面を 構成する質点とみなし、その質点変動(地殻変動)を、それらGPSステーションの所定の位置座標からの変位として定量化する。各GPSステーションの東西(E)、南北(N)、上下方向(h)からなる直交座標系(E,N,h)の各軸の成分をcとすると、時刻jの変位成分cは、D(c,j)となり、各成分の変位時系列は、[数57]で与えられる。変位時系列には、地殻表面の質点運動とは無関係なGPS観測に付随するランダムノイズや、周期変動が多く含まれる。従って、地殻変動に出現する異常を検出する時、第1の課題と同様なランダム変動する箇所に微分を用いた物理法則を適用する事が不可能となる障害が生じる。この様な状況下、本願の発明者は、[図7]に例示した任意な周波数選択機能を持つ物理的ウエーブレットを用いて、地殻表面の質点の変位時系列[c]から得た変位、速度、加速度成分を利用し、プレート境界を含む地殻表面の変化に出現する異常とその推移を、(1)巨大地震の予兆の発生、(2)その予兆の確定、(3)巨大地震が何時発生しても不思議でない状態、の3段階に確定できた。それらは、次のものである。
(1)東日本巨大地震発生(2011年3月11日に発生)の約1年4ヶ月前から生じ始めた大陸プレートの東の端の上に位置する東日本の東西方向の断面からみた、東海岸側で約1mm、西海岸側で約3mmの膨らみを、その巨大地震の予兆の発生とする。
(2)巨大地震発生約8カ月前から生じ始めた太平洋プレートの西方向への異常な急加速後、77日程前に、東方向へのスピードが通常の約3倍に達し、その直後の急減速から停止を予兆の確定とする。
(3)巨大地震発生44日程前に太平洋プレートの西方向への移動が停止し、その停止状態の巨大地震発生まで継続を、巨大地震発生直前の、何時、巨大地震が、発生してもおかしくない予兆の最終状態とする。
従って、巨大地震の上記予兆発生と推移とを確定する物理法則と、巨 大地震発生の上記自然法則とを利用し、広範囲なプレート境界領域に、次に述べる第5の課題である[数60]や[数61]を利用したAMR現象の検出を確立し、巨大地震の予知の方法、システム、プログラム、及び記録媒体を開発するに至った。
[背景技術]の[0149]-[0163]で説明した課題内容をまとめる。[特許文献1]で利用したせん断応力の蓄積の推移をモニターする手段は、地震の発生時刻の時間間隔INTのd(INT,j)を時刻mまで2s個、累積加算したCI(m,2s)をグラフ表示する事であった。[特許文献1]では、その時系列[CI(m,2s)]のグラフを、マグニチュードの大きさの程度の 予知に利用した。本願の発明者は、[図28]を用いて時系列[CI(m,2s)]の独立変数の数を特定する理論を確立し、時系列[CI(m,2s)]のモニター手段を、物理法則から得た歪エネルギー密度を規格化した[数60]や[数61]を用いたモニター手段に改善し、地殻に蓄積された歪エネルギーが 大地震や巨大地震の発生前に急激に解放され始めるAMR現象の検出技術を確立した。更に、領域の広さと、[数61]や[数62]で使用するローパスフィルターのカットオフ周波数とを変えた複数の観測によるAMR現象の検出技術を利用し、大地震と巨大地震のある程度の震源領域と発生時期の予知の方法、システム、プログラム、及び記録媒体を開発するに至った。
対象とする大陸プレート上の領域の地殻状態の変化を連続観測するために、その領域に稠密に設置されたGPSの各ステーションと、対象領域を含む大陸プレート下に沈み込んでいる海洋プレート上に設置されたGPSの各ステーションで、ステーションの東西方向をE、南北方向をN、上下方向をhとした直交座標系(E,N,h)の座標軸成分(要素)をパラメータcとし、j番目に観測したステーションの位置座標値をD(c,j)として数値化する地殻変動成分信号化手段と、
前記、震源要素信号化手段により抽出されたD(c,j)を、地震の発生を観測した順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは、所定の選択処理または平滑化処理を施し、[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]を得る震源要素の時系列化手段と、
前記、地殻変動成分信号化手段により抽出されたD(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは所定の選択処理または平滑化処理を施し、[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]を得る地殻変動成分の時系列化手段と、
前記、所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段と
前記、震源要素と地殻変動成分の時系列化手段により得られた時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差に比例する量を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段か、前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、1次差分値ΔD(c,m,n)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、 [C]を濾過したF1(c,m)を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段かのいずれかの1次差分出力手段と
前記、1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差に比例する量を、時刻mの2次差分値、Δ2D(c,m,n,k)として求めるD(c,m)の2次差分出力手段か前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、2次差分値Δ2D(c,m,n,k)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、 [C]を濾過したF2(c,m)を、時刻mの2次差分値Δ2D(c,m,n,k)として出力するD(c,m)の2次差分出力手段かのいずれかの2次差分出力手段と
からなる、地震発生と地殻状態の変化を定量化する方法において、
前記、2次差分出力手段から得た、時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする、2次差分値の時刻シフト手段と、
前記、1次差分出力手段から得た、時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする、1次差分値の時刻シフト手段と、
前記、震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として出力する変位変換手段と
前記、2次差分値の時刻シフト手段から得たΔ2D(c,t,n,k)を、前記変位変換手段から得たD(c,t)に作用している加速度A(c,t)に変換出力する加速度変換手段と
前記、1次差分値の時刻シフト手段から得たΔD(c,t,n)を、前記、時系列の変位変換手段から得たD(c,t)に作用している速度V(c,t)に変換出力する速度変換手段と
前記、変位変換手段で得たD(c,t)と、前記、速度変換手段で得た速度V(c,t)の各振幅と各位相の関係を出力する変位と速度の振幅位相比較検出手段と
前記、速度変換手段で得た速度V(c,t)と、前記、加速度変換手段で得た加速度A(c,t)との各振幅と各位相の関係を出力する速度と加速度の振幅位相比較検出手段と
前記、変位変換手段で得たD(c,t)と、前記、加速度変換手段で得た加速度A(c,t)の各振幅と各位相の関係を出力する変位と加速度の振幅位相比較検出手段と
前記、変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(INT,t)が、次の半周期後に反転する時刻と予測する 時刻ta2か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(DEP,t)が、次の半周期後に反転する時刻と予測する 時刻ta2か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(DEP,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(DEP,t)が反転したピークを取ると予測する時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(DEP,t)が半周期後に反転したピークを取った 時刻をta2とし、前記2次差分検出手段と加速度変換手段で 差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(DEP,t)が、次の半周期後に反転する時刻と予測する 時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(INT,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(INT,t)が反転したピークを取ると予測する時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(INT,t)が半周期後に反転したピークを取った 時刻をta2とし、前記2次差分検出手段と加速度変換手段で 差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(INT,t)が、次の半周期後に反転する時刻と予測する 時刻ta3か、を今にも発生しそうな大地震の発生時刻とする大地震発生時刻の予知手段と
前記、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記、変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする、今にも発生しそうな大地震の震源の予知手段と、
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値か、を今にも発生しそうな大地震の断層幅Wとする予測断層幅Wの検出手段と
予測断層幅Wの検出手段から出力されたWを用いて今にも発生しそうな大地震のマグニチュードの予知手段と
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の 時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とから、変位D(h,t)の膨らみを検出する手段か、
もしくは前記変位と速度の振幅位相比較検出手段と、前記、変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とから、変位D(h,t)の膨らみを検出する手段と、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記、変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする、前記、変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記、変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、大陸プレートの運動変化を検出する手段か、もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する手段と、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段と
海洋プレート上にある領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段か、
もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係からか、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段と
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する手段と
を備えてなる大地震と巨大地震の予知方法である。
m番目の地震とm−a番目の地震の発生時刻の差をCI(m,a)とする時間間隔累積検出手段と、
前記、CI(m,a)値がNCI(m,a)となるように規格化するNCI(m,a)の規格化手段と
前記、NCI(m,a)から大地震や巨大地震発生前の自然現象の一つである地震発生の変化が加速されるAMR現象(臨界状態)を検出し大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知方法において
対象とする領域の広さを、順次変える、観測領域の縮小と拡大手段と
前記、震源要素信号化手段と時系列化手段により抽出した地震の震源の深さ(DEP)の信号d(DEP,j)を地震の任意な発生順序mまで、a個ずつ移動平均するか、
a個ずつ、累積加算しCD(m,a)とする震源の深さDEPの累積検出手段と
前記、累積値CD(m,a)がNCD(m,a)となるように規格化するNCD(m,a)の規格化手段と
観測領域の縮小と拡大手段から出力された各領域のCI(m,a)とCD(m,a)のそれぞれ規格化したNCI(m,a)とNCD(m,a)とを、領域に蓄積と開放される規格化された歪エネルギー密度とする歪エネルギー変換手段と
領域の歪エネルギー変換手段から出力した地殻の蓄積歪エネルギー密度が巨大地震、大地震発生直前の臨界状態に到達したかを検出する地殻の臨界状態検出手段と、
前記、地殻の臨界状態検出手段で臨界状態を検出後大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知方法である。
対象とする大陸プレート上の領域の地殻状態の変化を連続観測するために、その領域に稠密に設置されたGPSの各ステーションと、対象領域を含む大陸プレート下に沈み込んでいる海洋プレート上に設置されたGPSの各ステーションで、ステーションの東西方向をE、南北方向をN、上下方向をhとした直交座標系(E,N,h)の座標軸成分(要素)をパラメータcとし、j番目に観測したステーションの位置座標値をD(c,j)として数値化する地殻変動成分信号化手段と、
前記、震源要素信号化手段により抽出されたD(c,j)を、地震の発生を観測した順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは、所定の選択処理または平滑化処理を施し、[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]を得る震源要素の時系列化手段と、
前記、地殻変動成分信号化手段により抽出されたD(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは所定の選択処理または平滑化処理を施し、[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]を得る地殻変動成分の時系列化手段と、
前記、所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段と
前記、震源要素と地殻変動成分の時系列化手段により得られた時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差に比例する量を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段か、前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、1次差分値ΔD(c,m,n)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、 [C]を濾過したF1(c,m)を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段かのいずれかの1次差分出力手段と
前記、1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差に比例する量を、時刻mの2次差分値、Δ2D(c,m,n,k)として求めるD(c,m)の2次差分出力手段か前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、2次差分値Δ2D(c,m,n,k)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、 [C]を濾過したF2(c,m)を、時刻mの2次差分値Δ2D(c,m,n,k)として出力するD(c,m)の2次差分出力手段かのいずれかの2次差分出力手段と
からなる、地震発生と地殻状態の変化を定量化する装置において、
前記、2次差分出力手段から得た、時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする、2次差分値の時刻シフト手段と、
前記、1次差分出力手段から得た、時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする、1次差分値の時刻シフト手段と、
前記、震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として出力する変位変換手段と
前記、2次差分値の時刻シフト手段から得たΔ2D(c,t,n,k)を、前記変位変換手段から得たD(c,t)に作用している加速度A(c,t)に変換出力する加速度変換手段と
前記、1次差分値の時刻シフト手段から得たΔD(c,t,n)を、前記、時系列の変位変換手段から得たD(c,t)に作用している速度V(c,t)に変換出力する速度変換手段と
前記、変位変換手段で得たD(c,t)と、前記、速度変換手段で得た速度V(c,t)の各振幅と各位相の関係を出力する変位と速度の振幅位相比較検出手段と
前記、速度変換手段で得た速度V(c,t)と、前記、加速度変換手段で得た加速度A(c,t)との各振幅と各位相の関係を出力する速度と加速度の振幅位相比較検出手段と
前記、変位変換手段で得たD(c,t)と、前記、加速度変換手段で得た加速度A(c,t)の各振幅と各位相の関係を出力する変位と加速度の振幅位相比較検出手段と
前記、変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(INT,t)が、次の半周期後に反転する時刻と予測する 時刻ta2か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(DEP,t)が、次の半周期後に反転する時刻と予測する 時刻ta2か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(DEP,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(DEP,t)が反転したピークを取ると予測する時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(DEP,t)が半周期後に反転したピークを取った 時刻をta2とし、前記2次差分検出手段と加速度変換手段で 差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(DEP,t)が、次の半周期後に反転する時刻と予測する 時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(INT,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(INT,t)が反転したピークを取ると予測する時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(INT,t)が半周期後に反転したピークを取った 時刻をta2とし、前記2次差分検出手段と加速度変換手段で 差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(INT,t)が、次の半周期後に反転する時刻と予測する 時刻ta3か、を今にも発生しそうな大地震の発生時刻とする大地震発生時刻の予知手段と
前記、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記、変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする、今にも発生しそうな大地震の震源の予知手段と、
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値かを今にも発生しそうな大地震の断層幅Wとする予測断層幅Wの検出手段と予測断層幅Wの検出手段から出力されたWを用いて今にも発生しそうな大地震のマグニチュードの予知手段と
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の 時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とから、変位D(h,t)の膨らみを検出する手段か、
もしくは前記変位と速度の振幅位相比較検出手段と、前記、変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とから、変位D(h,t)の膨らみを検出する手段と、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記、変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする、前記、変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記、変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、大陸プレートの運動変化を検出する手段か、
もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する手段と、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段と
海洋プレート上にある領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段か、もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係からか、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段と
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する手段と
を備えてなる大地震と巨大地震の予知装置である。
m番目の地震とm−a番目の地震の発生時刻の差をCI(m,a)とする時間間隔累積検出手段と、
前記、CI(m,a)値がNCI(m,a)となるように規格化するNCI(m,a)の規格化手段と
前記、NCI(m,a)から大地震や巨大地震発生前の自然現象の一つである地震発生の変化が加速されるAMR現象(臨界状態)を検出し大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知装置において
対象とする領域の広さを、順次変える、観測領域の縮小と拡大手段と
前記、震源要素信号化手段と時系列化手段により抽出した地震の震源の深さ(DEP)の信号d(DEP,j)を地震の任意な発生順序mまで、a個ずつ移動平均するか、
a個ずつ、累積加算しCD(m,a)とする震源の深さDEPの累積検出手段と
前記、累積値CD(m,a)がNCD(m,a)となるように規格化するNCD(m,a)の規格化手段と
観測領域の縮小と拡大手段から出力された各領域のCI(m,a)とCD(m,a)のそれぞれ規格化したNCI(m,a)とNCD(m,a)とを、領域に蓄積と開放される規格化された歪エネルギー密度とする歪エネルギー変換手段と
領域の歪エネルギー変換手段から出力した地殻の蓄積歪エネルギー密度が巨大地震、大地震発生直前の臨界状態に到達したかを検出する地殻の臨界状態検出手段と、
前記、地殻の臨界状態検出手段で臨界状態を検出後大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知装置である。
対象とする大陸プレート上の領域の地殻状態の変化を連続観測するために、その領域に稠密に設置されたGPSの各ステーションと、対象領域を含む大陸プレート下に沈み込んでいる海洋プレート上に設置されたGPSの各ステーションで、ステーションの東西方向をE、南北方向をN、上下方向をhとした直交座標系(E,N,h)の座標軸成分(要素)をパラメータcとし、j番目に観測したステーションの位置座標値をD(c,j)として数値化する地殻変動成分信号化手段と、
前記、震源要素信号化手段により抽出されたD(c,j)を、地震の発生を観測した順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは、所定の選択処理または平滑化処理を施し、[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]を得る震源要素の時系列化手段と、
前記、地殻変動成分信号化手段により抽出されたD(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは所定の選択処理または平滑化処理を施し、[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]を得る地殻変動成分の時系列化手段と、
前記、所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段と
前記、震源要素と地殻変動成分の時系列化手段により得られた時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差に比例する量を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段か、前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、1次差分値ΔD(c,m,n)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、 [C]を濾過したF1(c,m)を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段かのいずれかの1次差分出力手段と
前記、1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差に比例する量を、時刻mの2次差分値、Δ2D(c,m,n,k)として求めるD(c,m)の2次差分出力手段か前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、2次差分値Δ2D(c,m,n,k)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、 [C]を濾過したF2(c,m)を、時刻mの2次差分値Δ2D(c,m,n,k)として出力するD(c,m)の2次差分出力手段かのいずれかの2次差分出力手段とからなる、地震発生と地殻状態の変化を定量化するプログラムにおいて、
前記、2次差分出力手段から得た、時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする、2次差分値の時刻シフト手段と、
前記、1次差分出力手段から得た、時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする、1次差分値の時刻シフト手段と、
前記、震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として出力する変位変換手段と
前記、2次差分値の時刻シフト手段から得たΔ2D(c,t,n,k)を、前記変位変換手段から得たD(c,t)に作用している加速度A(c,t)に変換出力する加速度変換手段と
前記、1次差分値の時刻シフト手段から得たΔD(c,t,n)を、前記、時系列の変位変換手段から得たD(c,t)に作用している速度V(c,t)に変換出力する速度変換手段と
前記、変位変換手段で得たD(c,t)と、前記、速度変換手段で得た速度V(c,t)の各振幅と各位相の関係を出力する変位と速度の振幅位相比較検出手段と
前記、速度変換手段で得た速度V(c,t)と、前記、加速度変換手段で得た加速度A(c,t)との各振幅と各位相の関係を出力する速度と加速度の振幅位相比較検出手段と
前記、変位変換手段で得たD(c,t)と、前記、加速度変換手段で得た加速度A(c,t)の各振幅と各位相の関係を出力する変位と加速度の振幅位相比較検出手段と
前記、変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(INT,t)が、次の半周期後に反転する時刻と予測する 時刻ta2か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(DEP,t)が、次の半周期後に反転する時刻と予測する 時刻ta2か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(DEP,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(DEP,t)が反転したピークを取ると予測する時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(DEP,t)が半周期後に反転したピークを取った 時刻をta2とし、前記2次差分検出手段と加速度変換手段で 差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(DEP,t)が、次の半周期後に反転する時刻と予測する 時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(INT,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(INT,t)が反転したピークを取ると予測する時刻ta3か、加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(INT,t)が半周期後に反転したピークを取った 時刻をta2とし、前記2次差分検出手段と加速度変換手段で 差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(INT,t)が、次の半周期後に反転する時刻と予測する 時刻ta3か、を今にも発生しそうな大地震の発生時刻とする大地震発生時刻の予知手段と
前記、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記、変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする、今にも発生しそうな大地震の震源の予知手段と、
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値かを今にも発生しそうな大地震の断層幅Wとする予測断層幅Wの検出手段と
予測断層幅Wの検出手段から出力されたWを用いて今にも発生しそうな大地震のマグニチュードの予知手段と
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の 時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とから、変位D(h,t)の膨らみを検出する手段か、もしくは前記変位と速度の振幅位相比較検出手段と、
前記、変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とから、変位D(h,t)の膨らみを検出する手段と、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記、変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする、前記、変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記、変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、大陸プレートの運動変化を検出する手段か、
もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する手段と、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段と
海洋プレート上にある領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段か、
もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係からか、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段と
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する手段と
を備えてなる大地震と巨大地震の予知プログラムである。
m番目の地震とm−a番目の地震の発生時刻の差をCI(m,a)とする時間間隔累積検出手段と、
前記、CI(m,a)値がNCI(m,a)となるように規格化するNCI(m,a)の規格化手段と
前記、NCI(m,a)から大地震や巨大地震発生前の自然現象の一つである地震発生の変化が加速されるAMR現象(臨界状態)を検出し大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知プログラムにおいて
対象とする領域の広さを、順次変える、観測領域の縮小と拡大手段と
前記、震源要素信号化手段と時系列化手段により抽出した地震の震源の深さ(DEP)の信号d(DEP,j)を地震の任意な発生順序mまで、a個ずつ移動平均するか、
a個ずつ、累積加算しCD(m,a)とする震源の深さDEPの累積検出手段と
前記、累積値CD(m,a)がNCD(m,a)となるように規格化するNCD(m,a)の規格化手段と
観測領域の縮小と拡大手段から出力された各領域のCI(m,a)とCD(m,a)のそれぞれ規格化したNCI(m,a)とNCD(m,a)とを、領域に蓄積と開放される規格化された歪エネルギー密度とする歪エネルギー変換手段と
領域の歪エネルギー変換手段から出力した地殻の蓄積歪エネルギー密度が巨大地震、大地震発生直前の臨界状態に到達したかを検出する地殻の臨界状態検出手段と、
前記、地殻の臨界状態検出手段で臨界状態を検出後大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知プログラムである。
(1)震源時(大地震発生時刻)の予知、
(2)マグニチュードの予知、
(3)加速度A(DEP,t)とA(INT,t)の波形の反転が、[図3]の様に、略、完全な反転とならない場合のCQKとCQT予兆の確定と予知、
(4)巨大地震の予知
であった。しかし、地震発生の変化と地殻表面(GPS)の変化とを連続観測すると、それぞれの観測値の変化は、連続だが、その時間微分が不可能となる。従って、微分を用いる変化率の導出に運動の物理法則を適用できないので、予知に利用可能な大地震発生と巨大地震発生に関する自然法則を確立できなかった。次に述べる各手段は、これらの課題を手段毎に解決する効果がある。
上記発明の震源要素信号化手段は、対象とする領域の地震の発生を、領域で観測される仮想粒子の出現とみなし、j番目に出現した仮想粒子の要素cに関する観測値をD(c,j)として数値化する。
従って、上記、震源要素の時系列化手段は、
仮想粒子が観測空間に描くジグザグな運動(軌跡)の各座標軸成分を、直接、もしくは所定の選択処理または平滑化処理を施し、各震源要素cの時系列、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
で与える。
仮想粒子が観測空間に描くジグザグな軌跡の各座標軸成分を、直接、得る場合、上記時系列[C]は[数1]で与えられる時系列
[c]=[d(c,1)、d(c,2)、d(c,3)、・・、d(c,m)、・・]
と同一な時系列となる。
上記発明の地殻変動成分信号化手段は、
広範囲なプレート境界に設置されたGPSの各ステーションの位置を地殻の質点の位置とみなし、その 質点の位置座標の変動を、時刻jの東西(E)、南北(N)、上下(h)方向からなる直行座標系(E,N,h)における各座標軸成分パラメータcとする変位D(c,j)とする。
上記発明の地殻変動成分の時系列化手段は、
前記地殻変動成分信号化手段により抽出されたD(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、得るか、もしくは、所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る。
直接、得る場合は、上記時系列[C]は、[数57]で与えられる時系列
[c]=[d(c,1)、d(c,2)、d(c,3)、・・、d(c,m)、・・]
と同一な時系列となる。
所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段は、
例えば、平滑処理に[数44]に示される2w+1個の移動平均を使用して時刻τの変位D(c,τ)を得る場合、その処理の遅れ時間Δは、wとなる。
上記発明のD(c,m)の1次差分出力手段は、
震源要素と地殻変動成分の時系列化手段から出力された時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差を、D(c,m)の1次差分値ΔD(c,m,n)とする。更に、移動平均を用いた平滑化処理してD(c,m)が得られたとすると、その平滑機能は、[数48]で与えられる。この時、ΔD(c,m,n)を得る1次差分機能は、[数49]で与えられるバンドパスフィルター機能を持つ。従って、[数49]と、略、等価なフィルターを用いて、[C]を濾過したF1(c,m)を、D(c,m)の1次差分値ΔD(c,m,n)とする事もできる。従って、微分不可能なD(c,m)の時間間隔nの1次の変化率を求める。なお、1次差分値ΔD(c,m,n)は、n=1の場合、D(c,m)が微分可能であれば、その1次微分によって得る事のできる速度に相当するD(c,m)の1次の変化率となる。
上記発明のD(c,m)の2次差分出力手段は、
1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差を、時刻mの2次差分値、Δ2D(c,m,n,k)とする。更に、D(c,m)は、平滑機能が、[数48]で与えられる移動平均を用いて得られたとする。この時、[数50]で与えられるバンドパスフィルター機能と略等価なフィルターを用いて、[C]を濾過したF2(c,m)を、D(c,m)の2次差分値Δ2D(c,m,n,k)とすることもできる。従って、微分不可能なΔD(c,m,n)の時間間隔kの1次の変化率を求めることができるので、微分不可能なD(c,m)の2次の変化率をもとめる事ができる。なお、D(c,m)とΔD(c,m,n)とが微分可能であれば、n=k=1の場合、その1次微分によって得る事のできるD(c,m)の加速度に相当する2次の変化率を求める事ができる。
しかし、D(c,m)のこれら1次、2次の変化率は、物理で定義される速度、加速度とはならない。従って、D(c,m)の1次、2次の変化率を、変位D(c,m)の速度、加速度に変換するために、次の時刻シフトの手段を、変位D(c,m)、その1次、2次差分で与えられた変化率に導入する。
先ず、上記発明の2次差分値の時刻シフト手段は、
時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする効果がある。なお、Δは、D(c,m)がd(c,m)から直接得られていない場合、その処理遅れの時間は、遅れ時間検出手段の出力Δとなる。[数48]で与えられるΔtが2w+1個の移動平均の場合は、Δはwとなる。又、直接得られる場合は、遅れは無いので、Δ=0となる。
上記発明の1次差分値の時刻シフト手段は、
時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする効果がある。
上記発明の変位変換手段は、
前記震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として変換出力する効果がある。
上記発明の加速度変換手段は、
上記Δ2D(c,t,n,k)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)に作用している加速度A(c,t)に変換出力する効果がある。
上記発明の速度変換手段は、
ΔD(c,t,n)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)に作用している速度V(c,t)に変換出力する効果がある。
従って、これら変換出力された変位D(c,t)、速度V(c,t)、加速度A(c,t)は、平滑処理として移動平均を用いると、[数45]-[数47]で記述される変位D(c,t)、速度V(c,t)、加速度A(c,t)に変換するので、上記各手段における各時刻のシフトと変換操作は、微分不可能であったD(c,m)を、微分可能なD(c,t)に変換し、その1次微分を速度、2次微分を加速度と定義を満足する速度V(c,t)と加速度A(c,t)とにそれぞれ変換する効果がある。
更に、変位D(c,t)の出力に使用された平滑を(Δt=2w+1)個の移動平均とすると、3つの変換操作は、任意パラメータw、n、kで周波数領域を選択する次の3つの機能を有する。
(1)変位D(c,t)への変換操作は、[数48]のローパスフィルター。
(2)速度V(c,t)への変換操作は、[数49]のバンドパスフィルター。
(3)加速度A(c,t)への変換操作は、[数50]のバンドパスフィルター。
従って、任意パラメータw、n、kで選択された周波数領域において、時系列要素が震源要素D(c,j)の場合、地震の発生の変化を、大地震の予知に利用できる周波数領域において、変位D(c,t)、速度V(c,t)、加速度A(c,t)で記述する効果がある。即ち、ノイズに埋もれた大地震発生の自然現象を抽出可能とし、更に、物理法則で記述できる効果がある。
また時系列要素が地殻変動成分D(c,j)場合、震源要素D(c,j)の場合と同様に、地殻状態の変化を、その変化の異常検出に利用できる周波数領域において、変位D(c,t)、速度V(c,t)、加速度A(c,t)で記述する効果がある。即ち、ノイズに埋もれた巨大地震を発生させる地殻変動に関する自然法則を抽出可能とし、物理法則で記述できる効果がある。
上記発明の変位と速度の振幅位相比較検出手段は、
要素や成分のパラメータcを変えて変位D(c,t)と速度V(c,t)との振幅と位相の比較をするのみならず、変位D(c,t)と速度V(c,t)の各振幅と各位相の関係を、位相平面図(D(c,t),V(c,t))にも出力する。
上記発明の速度と加速度の振幅位相比較検出手段は、
要素や成分のパラメータcを変えて速度V(c,t)と加速度A(c,t)との振幅と位相の比較をするのみならず、速度V(c,t)と加速度A(c,t)との各振幅と各位相の関係を平面図(V(c,t),A(c,t))にも出力する。
上記発明の変位と加速度の振幅位相比較検出手段は、
要素や成分のパラメータcを変えて変位D(c,t)と加速度A(c,t)との振幅と位相の比較をするのみならず、変位D(c,t)と加速度A(c,t)との各振幅と各位相の関係を平面図(D(c,t),A(c,t))にも出力し、[数43]のニュートンの運動の第二法則を出力する効果がある。
従って、これら振幅位相比較比較の手段は、大地震と巨大地震の発生に関する自然法則を、次に、説明するように予知に利用できる効果がある。
上記、発明の大地震発生時刻の予知手段は、
前記変位と加速度の振幅位相比較検出手段の出力から得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取るような特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)もしくはA(DEP,t)の振幅値が、ピークを取った時刻から、その半周期後に反転したピークに大地震が発生する自然法則に、[数43]のニュートンの運動の第二法則を導入して、その大地震発生の時刻を予知できる効果がある。
又、
前記変位と加速度の振幅位相比較検出手段の出力から得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取るような特別な位相振幅関係が、明確でない場合、上記発明の予測断層幅Wの検出手段で、利用する、[数43]のニュートンの運動の第二法則から、断層面をスリップさせるせん断応力と静止摩擦応力のバランス状態を把握して、大地震発生の時刻を予知できる効果がある。
大地震の震源の予知手段は、
前記、A(DEP,t)とA(INT,t)とA(MAG,t)との振幅値が、負の値を取るような特別な位相振幅関係を検出から、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を、線形延長した(LAT、LON、DEP)情報を予測震源とする効果がある。上記各手段における時刻は、イベント数であるが、震源時の間隔時間INTの平均値を用いて実時間への変換ができるので、今にも発生しそうな大地震の発生時刻、震源を予知する効果がある。
上記発明の予測断層幅Wの検出手段は、
[数43]のニュートンの運動の第二法則を利用すると、力の変化量ΔF(DEP,τ)が、D(DEP,τ)の変化量ΔD(DEP,τ)に比例するので、このΔF(DEP,τ)が、断層面をスリップさせるせん断応力の変化量に比例する。従って、D(DEP,τ)の変化量ΔD(DEP,τ)を、今にも発生しそうな大地震の断層幅Wとする効果がある。
具体的にΔD(DEP,τ)を算出するには、前記変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、
変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か
変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、
変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値か
をD(DEP,τ)の変化量ΔD(DEP,τ)とする。
上記発明のマグニチュードの予知手段は、
予測断層幅Wの検出手段から出力されたWと[数56]とを利用して、今にも発生しそうな大地震のマグニチュードMを算出できる効果がある。
上記発明の変位D(h,t)の膨らみを検出する手段は、
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の 時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とからの位相平面上(D(h,t)、V(h,t))に描かれる軌跡解析からか、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とからか、位相空間(D(h,t)、V(h,t)、A(h,t))の軌跡解析から変位D(h,t)の膨らみを検出する効果がある。
上記発明の大陸プレートの運動変化を検出する手段は、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係からか、もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する効果がある。
上記発明の大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段は、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する効果がある。
上記発明の海洋プレートの異常加速運動検出手段は、
海洋プレート上で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係からか、もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係から、海洋プレートの大陸プレート方向への 速度成分が、通常の速度の倍以上に急加速し、その直後、急減速、急停止する海洋プレートの異常加速運動を検出する効果がある。
上記発明の巨大地震の発生時刻、マグニチュードを予知する手段は、
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する効果がある。
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
とする。遅れ時間検出手段は、時系列化手段の処理に遅れが生じる場合、その遅れ時間をΔとする。例えば、平滑処理に、[数44]で、与えられるd(c,j)の2w+1個の移動平均や、[数48]のΔtを2w+1としたフィルターを利用すると、その遅れ時間Δは、Δ=Δt/2=wとなる。d(c,j)信号を、直接に、時系列[C]にする場合は、遅れ時間は無く、Δ=0で、[数1]や[数57]で与えられる。
ΔD(c,m,n)= D(c,m)-D(c,m-n)、
に比例する量を、D(c,m)の1次差分値ΔD(c,m,n)として出力するか、
[数49]のΔtを2w+1としたバンドパスフィルターと、略等価なフィルターを用いて、[C]を濾過したF1(c,m)を、D(c,m)の1次差分値ΔD(c,m,n)として出力するかのいずれかである。
Δ2D(c,m,n,k)=ΔD(c,m,n)- ΔD(c,m-k,n)、
に比例する量を、D(c,m)の2次差分値Δ2D(c,m,n,k)として出力するか、
[数50]のΔtを2w+1としたバンドパスフィルターと略等価なフィルターを用いて[C]を濾過したF2(c,m)を、D(c,m)の2次差分値Δ2D(c,m,n,k)として出力するかのいずれかである。
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tmin(>tmax)の変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、
変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か
変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、
変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値かを
│ΔD(DEP,t)│とし、このDEPの変化量に、ニュートンの運動の第二法則を適用すると、│ΔD(DEP,t)│を、次に述べる様に、今にも発生しそうな大地震の予測断層幅Wとする事が出来る。
大地震の発生時刻の予知手段は、
(1)加速度の振幅位相比較検出手段の出力から得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取るような特別な位相振幅関係([図3]の関係)が成立後、逆位相にあったA(INT,t)もしくはA(DEP,t)の振幅値のピークが次の半周期後に 反転する時点が、大地震の発生時刻となる自然法則か
(2)加速度の振幅位相比較検出手段の出力から得た加速度A(DEP,t)と、加速度A(INT,t)とが略逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取るような特別な位相振幅関係([図3]の関係が不明瞭)が成立後、略逆位相にあったA(DEP,t)の振幅値のピークが、次の1周期後に、戻る時点が、大地震の発生時刻となる自然法則か
に、ブロック8から出力されるニュートンの運動の第二法則を適用する。上記自然法則(1)に、運動の第二法則を適用すると、震源要素cをINTもしくはDEPしたA(c,t)の振幅値のピークが、次の半周期後に反転する時刻を、差分を取る時間間隔Δa、Δa = (n+k)/2、だけ進んだ時刻の変位D(c,t+Δa)が、A(c,t)と反対方向のピークに到達した時刻を、A(c,t)が反転ピークに到達する時刻とする。上記自然法則(2)に、運動の第二法則を適用すると、震源要素cをDEPとしたA(c,t)の振幅値のピークが、次の1周期後に戻る時刻を、差分を取る 時間間隔Δaだけ進んだ先の時刻の変位D(c,t+Δa)が、A(c,t)と反対方向のピークに到達した時刻を、A(c,t)が1周期後に到達する時刻とする。
大地震の震源を予知する手段は、
上記A(DEP,t)と、A(INT,t)とA(MAG,t)とが特別な位相振幅関係が成立した時刻から、大地震の発生時刻の予知手段から出力された予知時刻まで、前記、 変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする。
(1)中・四国・近畿地方に出現している予兆の進行状況に2月7日以降変化はありません。
(2)他の地域にも予兆が出現していますが、解析が間に合っていません。
なお、(2)の予知結果は、公開していない。
(1)発生時刻が、1994年10月24日のm=583から19イベント後の132日後のm=602(1995年4月25日頃)
(2)マグニチュードが、6.8
(3)震源が、(LAT=34.53度、LON=135.18度、DEP=20.8 km)
となる。なお上記(2)のイベント数時間の実時間への変更は、1994年10月24日頃のマグニチュード3.5以上の地震の発生率が、[図32]のNCI(m,70)から、m=583付近で9.66日毎に1回発生すると仮定できるので、19イベントは、1994年10月24日から183.5日後となる事による。
m YYYY/MM/DD-hh:mm:ss.ss LAT(度) LON(度) DEP(Km) INT(h) MAG
597 1994/12/15 05:58:39.60 34.475 135.88 8.1 133.1432 3.6
598 1995/01/03 13:09:49.29 34.093 135.194 12.67 463.186 4
599 1995/01/12 06:35:32.14 33.979 135.198 12.6 209.4286 3.5
600 1995/01/16 01:19:00.81 33.662 132.195 59.65 90.7246 3.5
601 1995/01/16 18:28:02.68 34.598 135.039 18.4 17.1505 3.7
602 1995/01/17 05:46:51.86 34.595 135.038 16.06 11.3137 7.2
603 1995/01/17 05:49:14.63 34.617 135.082 13.85 0.0397 4.4
604 1995/01/17 05:49:35.17 34.621 135.076 13.31 0.0057 4.2
605 1995/01/17 05:49:48.81 34.658 135.118 11.5 0.0038 4.5
606 1995/01/17 05:50:23.95 34.653 135.133 12.96 0.0098 5.2
607 1995/01/17 05:51:05.23 34.589 135.002 12.16 0.0115 3.8
上記、予知情報と、テーブルのm=602の兵庫県南部地震の震源情報は、精度良く一致しているが、イベント数の時間から、実時間への変換は、1994年10月24日の時点では悪い。しかし、イベント数の予知時刻に接近すると、実時間への変換精度を上げることができる。兵庫県南部地震の予知検証を、[図35]を用いて行う。そのため、[図35]は、[図34]の対象とした期間を、1991年1月7日(m=440)から、1995年1月17日(m=640)まで拡張した。従って、[図34]に示した予知情報に、兵庫県南部地震の発生が、[図35]に加わる。この拡張した期間に、加わったグラフ表示機能は、次の機能である。最下段のMAGのd(MAG,m)の細線グラフは、MAGの値が6以上になると、その直前の値から太線に変換されて連結される。M7.2の兵庫県南部地震は、m=602(実時間で1995年1月17日)で、発生した。従って、その大地震発生直前のm=601(実時間で1995年1月16日)の値d(MAG,601)=3.7からd(MAG,602)=7.2まで、太い線グラフとなる。その太線には、横矢印で、1995/1/17 M7.2と兵庫県南部地震の発生日とマグニチュードを表示した。又、各震源要素cのc行の変位グラフもd(c,601)−d(c,602)が細線から太線表示となる。ただし、大地震発生直前(m=601)のM3.7地震は前震で、その震源は、本震の大地震と略同じであったので、d(MAG,m)の様にM3.7からM7.2へと大きな数値の変化がないので、グラフの細線は、それらの個所で、太線でなく太い点表示になっている。更に、各震源要素cのc行の加速度A(c,t)も、t= m-w-sの関係から、A(c,554)-A(c,555)の細線グラフが太い点となる。
NDEP-25(t) = 47.16km、
NINT-25(t) = 331.33 h (時間)、
NMAG-25(t) = 4.4
である。ただし、NMAG-25(t)は、3からオフセット値を規格化しているので、最大値は1.4が1となる規格化である。左側の縦軸は、規格化した変位の目盛りで、右側の縦軸は、時系列[d(LON,m)]と[d(MAG,m)]グラフの目盛である。これら表示グラフに用いたシンボルとグラフの名前は、左上の枠内に、上から、○印が、NDEP-25(t)、△印が、NINT-25(t)、□印が、NMAG-25(t)、点が、d(LON,m)、矢印が、d(MAG,m)、として表示してある。なお、右側の縦軸目盛のラベルMAGは5以上なので、地震のマグニチュードが5以上になると矢印で表示される。右側の縦軸目盛の上のラベルLONは、震源パラメータの時系列[d(LON,m)]の相対目盛で、経度目盛りの131.5度〜136.5度に相当する。横軸は、時間tとmで、t=m-12の関係がある。時間mの範囲は、m=400-700で、1990/4/28から1995/2/2の期間に相当する。兵庫県南部地震のCQKD予兆出現個所には、NINT-25(t)が、t=530で下方向のピークを取る個所に上向きの二重線矢印とCQKDのラベルが付けられている。時系列[d(MAG,m)]の兵庫県南部地震の発生個所は、時刻m = 602なので、d(MAG,601)からd(MAG,602)の矢印の横に1995/1/17 M7.2のラベルを付けた。規格化したNDEP-25(t)、NINT-25(t)、NMAG-25(t)のグラフの兵庫県南部地震の発生時刻は、時刻t=590なので、その個所に、下向きの破線矢印を付けてある。なお、NINT-25(t)とNDEP-25(t)は、[数60]と[数61]の時間mを時間tに変換した歪エネルギー密度でもある。
(1)発生時刻が、24イベント後のm=1708で、148日後の2011年7月26日頃
(2)マグニチュードが、5.8
(3)震源が、(LAT=33.8度、LON=135.6度、DEP=6km)
となる。
m YYYY/MM/DD-hh:mm:ss.ss LAT(度) LON(度) DEP(Km) INT(h) MAG
1705 2011/06/29 23:36:11.12 33.414 131.912 64.5 290.451 3.8
1706 2011/07/01 14:22:09.82 33.391 132.228 46.4 38.7663 4.1
1707 2011/07/04 11:07:11.82 33.931 135.547 46.5 68.7506 4.1
1708 2011/07/05 19:18:43.44 33.99 135.234 7.3 32.1921 5.5
1709 2011/07/05 19:34:55.64 33.996 135.242 7.1 0.2701 4.7
1710 2011/07/05 22:31:39.59 32.665 132.303 38.5 2.9455 3.5
1711 2011/07/24 23:32:12.40 33.92 136.139 42.3 457.0091 4.8
従って、その予知地震は、和歌山県で2011年7月5日に発生した震源の深さが、7.3km、マグニチュードが5.5の地震であった事が判明し、時刻のイベント数予知、震源とマグニチュードの予知は、正確であった。イベント数の時間から、実時間への変換は、2011年2月28日の時点では悪いが、予知時刻に接近すると、精度を上げることができる。
年代 m 年/月/日 緯度 経度 深さ(km) M
S1997 51 1997/03/03 34.956 139.155 5.45 3.9
S1998a 457 1998/04/21 34.959 139.161 7.42 3.9
S1998b 610 1998/08/07 36.237 137.661 4.55 4.2
S2000 1184 2000/06/27 34.095 139.47 10.26 3.9
2003 3994 2003/11/12 33.17 137.057 397.8 6.5
2004a 4140 2004/05/30 34.105 141.862 23.3 6.7
2004b 4256 2004/09/05 33.03 136.8 37.6 6.9
2004c 4276 2004/09/05 33.144 137.142 43.5 7.4
2004d 4517 2004/10/23 37.289 138.87 13.1 6.8
2005 4844 2005/01/19 33.953 141.993 30 6.8
2007a 5547 2007/03/25 37.221 136.686 10.7 6.9
2007b 5749 2007/07/16 37.557 138.609 16.8 6.8
2008a 6068 2008/05/08 36.275 141.979 18 6.3
2008b 6073 2008/05/08 36.228 141.608 50.6 7
2009 6427 2009/08/09 33.128 138.404 332.9 6.9
2010a 6649 2010/03/14 37.724 141.818 39.8 6.7
2010b 6708 2010/06/13 37.396 141.796 40.3 6.4
2011 6919 2011/03/11 36.108 141.265 43.2 7.4
である。又、[図45]は、[図44]の震源分布に、2011年3月7日から2012年3月26日までに発生したマグニチュード3.0以上の地震を、追加した。従って、[図45]は、2011年3月11日に発生した東北巨大地震の余震が追加されている。
(1)発生時刻が、7か10イベント後の2011年3月7日か3月10日頃
(2)マグニチュードが、7.1
(3)震源が、(LAT=35.83度、LON=141.62度、DEP=29km)
であった。
m YYYY/MM/DD-hh:mm:ss.ss LAT(度) LON(度) DEP(Km) INT(h) MAG
6905 2011/02/28 11:32:59.30 36.164 137.479 4.9 7.4924 3.6
6906 2011/03/01 23:12:54.39 37.16 141.466 47.2 35.6653 4.1
6907 2011/03/02 08:25:26.17 33.044 141.552 60 9.2088 4.7
6908 2011/03/03 16:37:26.62 32.925 138.386 318 32.2001 3.5
6909 2011/03/06 16:35:32.38 35.527 141.168 27.9 71.9683 3.5
6910 2011/03/07 05:10:58.39 36.523 140.571 58 12.5906 3.9
6911 2011/03/07 15:13:04.02 35.683 139.277 139.9 10.0349 3.8
6912 2011/03/09 06:31:36.42 33.764 141.168 42 39.309 3.5
6913 2011/03/11 14:54:42.22 36.717 140.576 9.1 56.3849 5.7
6914 2011/03/11 14:57:17.71 36.249 137.582 2.3 0.0432 4.7
6915 2011/03/11 14:58:05.81 37.676 141.911 23 0.0134 6.4
6916 2011/03/11 15:05:06.28 37.521 141.626 17 0.1168 5.9
6917 2011/03/11 15:11:18.95 36.864 140.62 7.9 0.1035 5.6
6918 2011/03/11 15:12:58.41 37.205 141.66 27.1 0.0276 6.7
6919 2011/03/11 15:15:34.46 36.108 141.265 43.2 0.0433 7.4
6920 2011/03/11 15:17:17.48 35.959 141.07 32.6 0.0286 5.7
上記予知の検証を、[図51]に図示する。先ず、発生時刻の予知に関して、2月28日のm=6905から8イベント後の3月11日に東北巨大地震の余震(M5.7)が発生し始めた。そのM5.7余震発生から2イベント後にM6.4地震、5イベント後にM6.7地震、その直後の6イベント後にM7.4地震が巨大地震の大きな余震として発生した。従って、2月28日から、7イベントか10イベント後の2011年3月7日頃か3月10日頃にM7.1の大地震が発生すると予知した予知精度は、大地震が3月11日に発生した東北巨大地震の余震であったにも関わらず正確である。最下段のMAGのd(MAG,m)の細線グラフは、MAGの値が6以上になると、その直前の値から太線に変換されて連結される。従って、m=6914の値d(MAG,6914)=4.7からd(MAG,6915)=6.4まで、更に、m=6917の値d(MAG,6917)=5.6からd(MAG,6918)=6.7、d(MAG,6919)=7.4まで太い線グラフとなる。太線には、横矢印で、2011/3/11 M6.4と2011/3/11 M7.4とそれら大地震(余震)の発生日とマグニチュードが表示してある。なお、マグニチュードの表示範囲は、その値が3.5から6.5までである。又、各震源要素cの各行の変位グラフも細線から太線表示となる。各震源要素cのc行の加速度A(c,t)も、t= m-w-sの関係から、d(c,m)の太線グラフから27イベント遅れた個所で太線表示となっている。従って、LAT、LON、DEP行の破線矢印は、d(c,m)の太線の位置付近を、指しているので震源の予知は、略正確である。なお、その予知した震源は、破線矢印の先の個所の読み取り値で、緯度が(LAT=35.83)、経度が(LON=141.62)、震源の深さが(DEP=29 km)であり、M7.4地震の震源は、(LAT=36.108度、LON=141.265度、DEP=43.2km)である。発生日の予知に関して、INT行の下向き破線矢印の起点がA(INT,t)の太線の始まりで、その予知も正確であった事を図示している。
NDEP-15(t) = 222.86km、
NINT-15(t) = 2.79 h (時間)、
NMAG-15(t) = 5.55
である。左側の縦軸は、規格化した変位の目盛りで、右側の縦軸は、時系列[d(LON,m)]と[d(MAG,m)]グラフの目盛である。これら表示グラフに用いたシンボルとグラフの名前は、左上の枠内に、上から、○印が、NDEP-15(t)、△印が、NINT-15(t)、□印が、NMAG-15(t)、点が、d(LON,m)、矢印が、d(MAG,m)、として表示してある。なお、右側の縦軸目盛のラベルMAGは6以上なので、地震のマグニチュードが6以上になると矢印の高さに変換表示される。右側の縦軸目盛の上のラベルLONは、震源パラメータの時系列[d(LON,m)]を点グラフにした相対目盛で、その経度範囲は、136.5度〜142度に相当する。これら表示グラフに用いたシンボルとグラフの名前は、左上の枠内に、上から、○印が、NDEP-15(t)、△印が、NINT-15(t)、□印が、NMAG-15(t)、点が、d(LON,m)、矢印が、d(MAG,m)、として表示してある。これらグラフの横軸は、時間tとmで、t=m-7の関係がある。時間mの範囲は、m=6600-7000で、2010/01/04から2011/03/11の期間に相当する。大余震M7.4のCQTDの個所には、NDEP-15(t)が、t=6870で下方向のピークを取る個所に上向きの二重線矢印とCQKDのラベルが、付けられている。そのCQKD予兆のM7.4大地震の発生個所は、時刻m = 6919なので、下向き破線矢印で示した時刻t=6912の個所である。又、その大地震は、2011のラベルを付けた太線矢印のd(MAG,m=6919)である。時刻t=6625と6670付近に示した下向き二重線矢印で示した2つのCQTDは、それぞれ2010aと2010bとラベル付けされたM6.8とM6.4の大地震の予兆であった事を検証している。又、これら予兆時に検出された今にも発生しそうだった大地震の予測断層幅Wは、その幅Wに相当するNDEP-15(t)の個所に破線の両矢印とWとでラベル表示してある。
1997年1月1日から2000年6月30日を[図53]。
2000年6月30日から2000年8月16日までを[図54]。
2000年8月16日から2004年10月13日までを[図55]。
2003年11月17日から2008年3月15日までを[図56]。
2006年8月14日から2011年3月14日までを [図57]。
これらの図は、前記、中国・四国・近畿地方の歪エネルギー密度の推移から大地震発生直前の地殻の臨界状態の抽出を示した[図32]と同様な図である。
m = 3940 (2003/09/06)の時、CI(m,40)max = 40 x 49.0541時間、
m = 6460 (2009/08/20)の時、CD(m,40)max = 40 x 126.815 km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,40)、NCD(m,40)とを、[図53]-[図57]にグラフ表示した。左側の縦軸目盛は、NCI(m,40)とNCD(m,40)グラフの目盛である。又、右側のMAG(マグニチュード)とラベル付けした縦軸目盛は、時系列[d(MAG,m)]の6以上の目盛である。従って、[d(MAG,m)]のグラフは、その6以上の値が矢印の高さに変換され、オフセット表示される。右側の縦軸目盛にLON(経度)とラベル付けした目盛りは、震源パラメータの時系列[d(LON,m)]目盛で、その点グラフの範囲は、経度範囲の136.5度〜142度に相当する。横軸は、各グラフに共通な時間mで、地震の発生順番を示すインデックスである。大地震の余震や群発地震は、震源が略同じなので、それら余震や群発地震の震源の経度は殆んど変化しない。従って、[d(LON,m)]の点グラフの変動は無くなる。マグニチュードが6以上の主な大地震には、その発生年代が、[図55]-[図57]の時系列[d(MAG,m)]にラベル表記されている。しかし、群発地震は、マグニチュードが6以下の地震からなる事が多い。従って、それら群発地震を、発生年代毎に、群発地震(Swarm)のSを年代に付け、[図53]-[図54]の時系列[d(LON,m)]の点グラフで変動の無い個所に矢印で示した。なお、2000年6月27日から始まり2カ月余りに及んだ神津島沖のマグニチュードが6を超える大きな群発地震にも、[図53]-[図54]にS2000とラベル付けしている。群発地震発生中の歪エネルギー密度の変化は、発生前の変化と比較すると少ないので、[図54]のNCI(m,40)は、10倍に拡大、NCD(m,40)は、5倍に拡大してある。この様に、拡大した歪エネルギー密度から、群発地震発生期間中でも大地震が発生する直前の地殻の臨界状態(AMR現象)を抽出している実施例が、[図54]である。同様に、大地震や巨大地震の余震 発生中の地殻の歪エネルギー密度の変化を拡大すると、大きな余震が発生する直前のAMR現象を抽出できるので、大きな余震発生時刻の予知に利用できる。この様に、[図53]-[図57]にグラフ表示したNCD(m,40)とNCI(m,40)は、選択した関東・甲信越・中越地方(32.5度−38度、136.5度−142度)に今にも発生しそうな群発地震、大地震、巨大地震発生直前の地殻の臨界状態を、地殻に蓄積、解放される歪エネルギーの推移から抽出した実施例である。各グラフの最上段に記した選択地震、領域、期間は、[図32]と同様に、選択した地震のマグニチュードMと領域(緯度、経度)の範囲、横軸のイベント時間mの範囲を実時間に変換した期間である。
年代 m 年/月/日 緯度 経度 深さ(km) M
S1998 119 1998/08/07 36.237 137.661 4.55 4.2
2004 452 2004/10/23 37.289 138.87 13.1 6.8
2007a 895 2007/03/25 37.221 136.686 10.7 6.9
2007b 1085 2007/07/16 37.557 138.609 16.8 6.8
2011 1396 2011/03/12 36.986 138.598 8.4 6.7
である。
(1)発生時刻が、14イベント後の77日後の5月16日頃
(2)マグニチュードが、6.7
(3)震源が(LAT=36.14度、LON=138.33度、DEP=27km)
であった。
m YYYY/MM/DD-hh:mm:ss.ss LAT(度) LON(度) DEP(Km) INT(h) MAG
1380 2011/02/28 11:32:59.30 36.164 137.479 4.9 7.4924 3.6
1381 2011/03/01 16:50:19.93 36.157 137.445 4.6 29.2891 3.4
1382 2011/03/11 14:57:17.71 36.249 137.582 2.3 238.116 4.7
1383 2011/03/11 15:06:34.00 36.258 137.582 2.3 0.1545 3.3
1384 2011/03/11 16:33:40.95 37.948 138.096 14.9 1.4519 4.3
1385 2011/03/11 16:34:34.85 38.013 139.427 16.1 0.015 4.1
1386 2011/03/11 17:58:07.70 36.813 139.469 3.1 1.3925 4.8
1387 2011/03/11 19:09:39.05 36.122 136.512 8.6 1.192 4
1388 2011/03/11 19:17:29.63 36.124 136.502 8.9 0.1307 3.4
1389 2011/03/11 19:38:06.20 36.259 137.564 2 0.3435 3.6
1390 2011/03/11 19:40:53.78 36.817 139.461 4.7 0.0465 3.4
1391 2011/03/11 20:10:54.85 36.864 139.327 5.6 0.5003 3.6
1392 2011/03/11 22:09:28.59 36.249 137.56 2.6 1.976 3.4
1393 2011/03/12 00:11:42.85 36.248 137.556 2.3 2.0373 3.8
1394 2011/03/12 00:24:40.01 36.804 139.312 6.4 0.2159 4.5
1395 2011/03/12 00:47:47.56 36.766 139.475 4.9 0.3854 3.3
1396 2011/03/12 03:59:15.62 36.986 138.598 8.4 3.1911 6.7
1397 2011/03/12 04:01:14.19 36.917 138.598 10.5 0.0329 3.6
上記予知の検証を、[図61]に図示する。先ず、発生時刻の予知に関して、2月28日のm=1380から2イベント後の3月11日から東北巨大地震の余震が発生し始めた。その 余震発生から14イベント後に、2月28日から16イベント後に、M6.7の長野県北部の地震が発生した。従って、2月28日から、14イベント後M6.7の大地震が発生すると予知した予知精度は、イベント数では、2イベント早く予知しただけで、その予知精度は良い。しかし、イベントの 時刻から実時間への変換に東北巨大地震の余震の影響が考慮されていないので、実時間の予知精度は悪い。この悪い精度は、3月11日に発生していた東北巨大地震の余震をモニターしていれば、M6.7 地震の発生日の予知を、3月12日とする事が可能であった。最下段のMAGのd(MAG,m)の細線グラフは、MAGの値が6以上になると、その直前の値から太線に変換されて連結される。従って、m=1395の値d(MAG,1395)=3.3からd(MAG,1396)=6.7まで太い線グラフとなるが、表示範囲は、その値が3.5から6.5までである。太線には、横矢印で、2011/3/12 M6.7と大地震(余震)の発生日とマグニチュードが表示してある。従って、マグニチュードの予知は、正確であった。又、他の震源要素cの各行の変位d(c,m)グラフも細線から太線表示となる。各震源要素cのc行の加速度A(c,t)も、t= m-w-sの関係から、d(c,m)の太線グラフから40イベント遅れた個所で太線表示となっている。従って、LAT、LON、DEP行の破線矢印の先は、d(c,m)の太線の位置付近を、指しているので震源の予知も、略正確である。なお、その予知した震源は、緯度が(LAT=36.14)、経度が(LON=138.33)、震源の深さが(DEP=27 km)であり、長野県北部に発生したM6.7地震の震源は、(LAT= 36.986度、LON= 138.598度、DEP= 8.4km)である。発生日の予知に関して、INT行の下向き破線矢印がA(INT,t)の太線の始まりで、その予知も、2イベント早かったが、略正確であった事を図が検証している。
NDEP-21(t) = 243.75km、
NINT-21(t) = 137.47 h (時間)、
NMAG-21(t) = 4.53
である。ただし、NMAG-21(t)は、3からオフセット値を規格化しているので、最大値は1.53が1となる規格化である。左側の縦軸は、規格化した変位の目盛りで、右側の縦軸は、時系列[d(LON,m)]と[d(MAG,m)]グラフの目盛である。又、右側の縦軸目盛のラベルMAGは5以上なので、地震のマグニチュードが5以上になると矢印の高さに変換表示される。右側の縦軸目盛の上のラベルLONは、震源パラメータの時系列[d(LON,m)]を点グラフにした相対目盛で、経度範囲の136度〜140度に相当する。これら表示グラフに用いたシンボルとグラフの名前は、左上の枠内に、上から、○印が、NDEP-21(t)、△印が、NINT-21(t)、□印が、NMAG-21(t)、点が、d(LON,m)、矢印が、d(MAG,m)、として表示してある。これらグラフの横軸は、時間tとmで、t=m-10の関係がある。時間mの範囲は、m=1200-1500で、2008/03/17から2011/04/29の期間に相当する。大余震M6.7のCQTD予兆の個所には、NDEP-21(t)が、t=1325付近で下方向のピークを取る個所に上向きの二重線矢印とCQKDのラベルが、付けられている。そのCQKD予兆のM6.7大地震の発生個所は、時刻m = 1396なので、下向き破線矢印で示した時刻t=1386の個所である。又、その大地震は、2011のラベルを付けた太線矢印のd(MAG,m=1396)である。予兆時に検出された今にも発生しそうだった大地震の予測断層幅Wは、その幅Wに相当するNDEP-21(t)の個所に破線の両矢印とWとでラベル表示してある。
m = 319 (2001/12/08)の時、CI(m,60)max = 60 x 221.016時間
m = 1269 (2009/03/30)の時、CD(m,60)max = 60 x 90.379km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,60)、NCD(m,60)とを、[図65]にグラフ表示した。左側の縦軸目盛は、NCI(m,60)とNCD(m,60)グラフの目盛である。又、右側の縦軸のMAGとラベル付けした5以上の目盛は、時系列[d(MAG,m)]を矢印の高さに変換したグラフ目盛である。従って、矢印の高さがマグニチュードのオフセット値となり、MAGが5以上の大地震になると表示される。更に、その上にある右側の縦軸目盛にLONとラベル付けされている目盛りの範囲は、震源パラメータの時系列[d(LON,m)]を点グラフにした経度範囲の136度〜140度に相当する。横軸は、各グラフに共通な時間mで、地震の発生順番を示すインデックスである。大地震の余震や群発地震は、震源が略同じなので、それら余震や群発地震の震源の経度は殆んど変化しない。従って、[d(LON,m)]の点グラフの変動が無くなる群発地震と、マグニチュードが6以上の主な大地震には、その発生年代が、[図65]の時系列[d(LON,m)]か[d(MAG,m)]にラベル表記されている。[図65]の時系列[d(LON,m)]の点グラフの変化の無い個所に矢印で示したSと発生年代ラベルは、その時発生した群発地震を表す。[図65]にグラフ表示したNCD(m,60)とNCI(m,60)は、選択した新潟・中越地方(32度−36度、136度−140度)に、今にも発生しそうな群発地震、大地震発生直前の地殻の臨界状態(AMR現象)を地殻に蓄積、解放される歪エネルギー密度の推移から抽出した実施例である。ただし、大地震2007bは、2007aと連動しているので、2007bのAMR現象の抽出は、できていない。
年代 m 年/月/日 緯度 経度 深さ(km) M
1998 224 1998/09/03 39.803 140.905 7.86 6.2
2001a 584 2001/08/14 40.993 142.44 37.69 6.4
2001b 613 2001/12/02 39.396 141.267 121.5 6.4
2003a 791 2003/05/26 38.818 141.654 72 7.1
2003b 926 2003/07/26 38.402 141.174 11.9 6.4
2005a 1281 2005/08/16 38.151 142.28 41.6 7.2
2005b 1347 2005/12/02 38.073 142.353 40.3 6.6
2008 1669 2008/06/14 39.03 140.881 7.8 7.2
2010 2110 2010/07/05 39.657 142.653 34.5 6.4
2011a 2197 2011/03/10 38.271 142.879 28.9 6.4
2011b 2201 2011/03/11 38.103 142.861 23.7 9
2011c 2714 2011/04/07 38.204 141.92 65.9 7.5
である。
(1)発生時刻が、3から6イベント後の2011年3月9日から3月21日頃
(2)マグニチュードが、6.7
(3)震源が(LAT=38.24度、LON=142.82度、DEP=19km)
であった。
m YYYY/MM/DD-hh:mm:ss.ss LAT(度) LON(度) DEP(Km) INT(h) MAG
2191 2011/02/24 11:11:32.47 40.897 139.121 28.6 102.6253 4.1
2192 2011/03/02 03:11:42.06 40.161 139.347 20.1 136.0027 3.6
2193 2011/03/05 09:25:56.57 38.537 142.105 42.3 78.2374 3.8
2194 2011/03/09 12:05:48.22 38.283 142.927 21.3 98.6643 5
2195 2011/03/09 13:45:50.59 38.438 142.941 22.1 1.6673 5.7
2196 2011/03/10 03:16:14.18 38.271 142.879 28.9 13.5066 6.4
2197 2011/03/10 06:00:54.64 38.259 142.913 28 2.7446 5
2198 2011/03/10 07:42:45.32 38.174 142.957 21.2 1.6974 4.2
2199 2011/03/11 13:12:30.23 38.263 142.952 18.2 29.4958 4.2
2200 2011/03/11 14:46:18.12 38.103 142.861 23.7 1.5633 9
2201 2011/03/11 14:57:25.81 39.197 142.382 23.8 0.1855 5.7
2202 2011/03/11 15:06:10.75 39.042 142.398 26.6 0.1458 6.7
2203 2011/03/11 15:08:53.50 39.839 142.781 31.7 0.0452 7.4
2204 2011/03/11 15:12:54.40 39.442 142.095 33.4 0.0669 4.8
2205 2011/03/11 15:23:05.40 39.003 142.447 29.8 0.1697 6.1
上記表から発生時刻の予知に関して、2月24日のm=2191から3イベント後の3月9日に東北巨大地震の前震のM5が発生した。更に5イベント後の3月10日にM6.4地震、そして、9イベント後の3月11日に、M9東北巨大地震が発生した。従って、震源を(LAT=38.24度、LON=142.82度、DEP=19km)としたM6.7の大地震が、2月24日から、3イベントから6イベント後に、発生するとした予知は、正確であったと言える。
NDEP-21(t) = 73.95km、
NINT-21(t) = 138.67 h (時間)、
NMAG-21(t) = 5.74
である。ただし、NMAG-21(t)は、3.3からオフセット値を規格化しているので、最大値は2.44が1となる規格化である。左側の縦軸は、規格化した変位の目盛りで、右側の縦軸は、時系列[d(LON,m)]と[d(MAG,m)]グラフの目盛である。又、右側の縦軸目盛のラベルMAGは6以上なので、地震のマグニチュードが6以上になると矢印の高さに変換表示される。右側の縦軸目盛の上のラベルLONは、震源パラメータの時系列[d(LON,m)]を点グラフにした相対目盛で、その経度範囲は、138度〜143度となる。これら表示グラフに用いたシンボルとグラフの名前は、左上の枠内に、上から、○印が、NDEP-21(t)、△印が、NINT-21(t)、□印が、NMAG-21(t)、点が、d(LON,m)、縦矢印が、d(MAG,m)、として表示してある。これらグラフの横軸は、時間tとmで、t=m-10の関係がある。時間mの範囲は、m=1900-2300で、2008/12/18から2011/03/14の期間に相当する。2つのCQTD予兆が、上下方向の二重線矢印で 表示されている。一番目のCQTD予兆は、2010とラベル表示されたM6.4地震の予兆で、2番目のCQTD予兆が、2011とd(MAG,2200)の矢印にラベル表示された東北巨大地震M9もしくはその前震のM6.4地震の予兆である。又、これら予兆時に 検出された今にも発生しそうだった大地震の予測断層幅Wは、その幅Wに相当するNDEP-21(t)の個所に破線の両矢印とWとでラベル表示してある。CQKD予兆の東北巨大地震M9の発生個所は、時刻m = 2200なので、下向き破線矢印で示した時刻t=2190の個所である。2011の東北巨大地震M9発生直後のNDEP-21(t)は、発生直前の加算平均した地震の影響が無くなると、浅い処から深い所へと移動している。この観測は、太平洋プレートと東北地方の東海岸の下の大陸プレートとの固着域が、浅い箇所から深い箇所に向けて破壊した事を示している。又、NDEP-21(t)は、[数61]の時間mを10イベント遅れの時間tに変換した歪エネルギー密度の関数でもある。従って、NDEP-21(t)グラフの変化は、東北巨大地震が発生している現象のグラフで、断層の浅い箇所から深い箇所へと向かって、蓄積されていた歪エネルギーが、一挙に開放された事を表す。巨大地震は、上向きの破線矢印で示した時刻t=2190で発生しているが、その直前の歪エネルギー密度NDEP-21(t)は、地殻の深い所でピークに到達し、深い所から、浅い方向へと、一挙に開放され始めているAMR現象を、検出している。このNDEP-21(t)が、AMR現象の抽出に使用した観測窓の時間幅は、次に述べる[図71]-[図73]のNCD(m,60)の時間幅60イベントより狭い21イベントである。
m = 582 (2001/08/01)の時、CI(m,60)max = 60 x 107.808時間、
m = 856 (2003/05/28)の時、CD(m,60)max = 60 x 69.75667km
となる。
これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,60)、NCD(m,60)とを、[図71]に、歪エネルギー密度としてグラフ表示した。左側の縦軸目盛は、歪 エネルギー密度、NCI(m,60)とNCD(m,60)の相対目盛である。又、右側の縦軸のMAGとラベル付けした6以上の目盛は、時系列[d(MAG,m)]のMAGが6以上になると、その値を矢印の高さに変換表示したグラフ目盛である。更に、MAGが6以上のd(MAG,m)を、破線矢印で歪エネルギー密度のグラフ(NCD(m,60)とNCI(m,60))付近まで延長し、AMR現象が大地震発生前に、抽出されたかを、明確にした。右側の縦軸目盛にLONとラベル付けされている目盛りの範囲は、震源パラメータの時系列[d(LON,m)]を点グラフにした、経度の目盛範囲、138度〜143度に相当する。横軸は、各グラフに共通な時間mで、地震の発生順番を示すインデックスである。大地震の余震や群発地震は、震源が略同じなので、それら余震や群発地震の震源の経度は殆んど変化しない。従って、[d(LON,m)]の点グラフの変動は無くなる。ただし、巨大地震の場合は、断層形状が大きく余震発生領域が拡がるので、点グラフは変動する。マグニチュードが6以上の主な大地震には、その発生年代を、[図71]の時系列[d(MAG,m)]の矢印に表記した。発生年代が同一の場合は、アルファベットa、bを付加した。地殻に蓄積された歪エネルギー、[図71]のNCD(m,60)とNCI(m,60)、のピークが、大地震発生の直前(矢印の少し前)から、減少し始めている。M7.1の大地震2003a(発生日=2003/05/26、LAT=38.818度、LON=141.654度、DEP=72km)の震源の深さDEPは深い。従って、NCD(m,60)で表した歪エネルギー密度は、M7.1 地震発生直後から増加し最大ピークに到達する。しかし、そのピークは、震源の浅いM6.4地震2003b(発生日=2003/07/26、LAT=38.402度、LON=141.174度、DEP=11.9km)とその前震の発生前から、歪エネルギー密度は減少し始めている。従って、今にも発生しそうな大地震発生直前の地殻の臨界状態(AMR現象)を歪エネルギーの推移から抽出した例である。
m = 570 (2001/06/25)の時、CI(m,40)max = 40 x 114.48時間
m = 855 (2003/05/28)の時、CD(m,40)max = 40 x 70.77km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,40)、NCD(m,40)とを、[図72]と[図73]に、歪エネルギー密度の推移グラフとして表示した。[図72]は、期間が1997年1月2日(m=0)から2006年12月27日(m=1500)までで、[図73]は、2006年12月27日(m=1500)から、2011年5月10日(m=3000)までである。[図72]と[図73]のグラフ表示、大地震のラベル、各目盛の表示方法は、2s=60の場合の[図71]と同一である。唯一の違いは、[図73]に表示した東北巨大地震2011b発生直後から、NCI(m,40)を10倍拡大表示した事である。この拡大したNCI(m,40)の歪エネルギー密度の推移から、[図71]のNCI(m,60)から明確に抽出できなかった大きなM7.5余震(2011cとラベル表示したm=2713で発生した余震、もしくは、大地震)の発生直前の地殻の臨界状態(AMR現象)を、抽出できた例を示した。このM7.5地震の発生日と震源は、(2011/04/07、LAT=38.204度、LON=141.92度、DEP=65.9km)であった。発生後、時系列[d(LON,m)]の変動が少なくなっているので、このM7.5 の大地震は、巨大地震2011bの余震ではなく、余震発生中に発生した、別の大地震である。大地震2011cの発生直後から、NCD(m,40)とNCI(m,40)の振幅変動が反転しているのは、大地震2011cとその余震の震源が深かったからである。
関東・甲信越・中越地方(32.5度−38度、136.5度−142度)のCQT予兆、
新潟・中越地方(36度−40度、136度−140度) のCQT予兆、
東北地方(38度−42度、138度−143度)のCQK予兆
から、それら予兆に該当する大地震を予知できた。発生日時の予知に関して、これら複数の大地震は、最終的に、2011年3月9日の時点で、3月11日頃発生すると予知でき、連動した大地震もしくは巨大地震として、3月11日に発生したM9の東北巨大地震の発生を推測できた。しかし、そのマグニチュードは、直接予知できなかった。それは、次の理由による。2011年3月11日のM9東日本巨大地震の断層長Lは、約500kmとなり、メッシュを構成する緯度差に換算すると約4.5度となる。一方、上記大地震の予知のため選択した領域の広さは、約4〜5度の緯度、経度幅からなるメッシュであった。従って、選択した領域(メッシュ)内の地震発生の変化から、いまにも発生しそうであった巨大地震の断層幅Wを直接抽出できなかった。
日本列島及び島々に張り巡らされた、国土地理院の約1200のGPSステーションの地殻変動解析に多くの時間を必要とするために、上記2月28日の時点で、その解析は一切実施されていなかった。その地殻変動の異常抽出は、3月11日の東北巨大地震発生後に、[図1]のブロック12−15で実施した。上記大地震と巨大地震を区別するマグニチュードの予知は、「プレート境界の地殻表面の変化に出現する異常が、巨大地震の予兆となり、その予兆の広がりが、今にも発生しそうな巨大地震の断層長(L)となる」自然現象を利用した。このマグニチュードと、発生時期の予知は、[図1]のブロック16とブロック21で実施し、[背景技術]の[0128]-[0148]で、既に、詳細したので、省略する。
中国・四国・近畿地方(32度-36度、131.5度-136.5度)(深発地震は含まない)、
近畿・福井・東南海地方 (31度-37度、134度-138度)、
関東・甲信越・中越地方(32.5度−38度、136.5度−142度)、
新潟・中越地方(36度−40度、136度−140度)、
東北地方(38度−42度、138度−143度)、
東北・北海道地方(37度−46度、138度−150度)、
日本列島を含む広範囲な領域(16度−52度、116度−156度)
とする。これら各領域の歪エネルギー密度から2011年3月11日に発生したM9の東北巨大地震の予知が可能であったとする検証を、実施例として述べる。
(1)中国・四国・近畿地方の[図33]の歪エネルギー密度の推移からは、2011年7月5日に発生した和歌山県で発生したM5.5地震の小規模なAMRの現象の検出(AMR-22g)は存在するが、2011年3月11日に発生しそうな大地震に関連した大規模なAMRの現象の検出は無い。ただし、この地方の歪エネルギー密度には、近畿地方に存在する深発地震は含まれていない。
(2)近畿・福井・東南海地方の深発地震も含めた歪エネルギー密度NCD(m,60)とNCI(m,60)には、2011年1月12日頃、共にピークに到達し、その直後、急減少を始め、18イベント後に3月11日を迎えた大規模なAMR現象の検出がある。
(3)関東・甲信越・中越地方の[図57]の2011とd(MAG,m)にラベル表示したM7.4の発生直前の大規模なAMR現象の検出がある。
(4)新潟・中越地方の[図65]の2011とd(MAG,m)にラベル表示したM6.7の発生直前の大規模なAMR現象の検出がある。
(5)東北地方の[図71]と[図74]の2011aと2011bとd(MAG,m)にラベル表示したM9巨大地震発生直前の大規模なAMR現象の検出がある。
(6)東北・北海道地方にも後述の[図80]にM9巨大地震発生直前、その前震を含めた大規模なAMR現象の検出がある。
(7)日本列島を含む広範囲な領域で、後述の[図82]-[図85]にM9巨大地震発生直前、その前震を含めた大規模なAMR現象の検出がある。
従って、東日本の太平洋プレートと大陸プレートの境界で大規模なAMR現象を検出している。このプレート境界の直線的な長さは、[背景技術]で説明した様に、500km余りとなる。従って、断層の長さ、L=500kmを[数55]に代入しマグニチュード9を得る。また、東北・北海道地方と日本列島を含む広範囲な領域のAMR現象抽出後から数日から数時間以内に、その巨大地震が発生すると予知できる。
年代 m 年/月/日 緯度 経度 深さ(km) M
1997a 12 1997/03/16 34.925 137.528 39.12 6
1997b 28 1997/05/24 34.498 137.503 23.1 6
2000 295 2000/06/07 36.824 135.566 21.3 6.2
2003 591 2003/11/12 33.17 137.057 397.8 6.5
2004a 648 2004/09/05 33.03 136.8 37.6 6.9
2004b 668 2004/09/05 33.144 137.142 43.5 7.4
2004c 783 2004/09/07 33.206 137.296 41 6.6
2004d 819 2004/09/08 33.114 137.291 36.1 6.5
2007 1111 2007/07/16 36.866 135.105 373.8 6.7
である。
m = 1084 (2007/03/30)の時、CI(m,60)max = 60 x 171.98時間、
m = 650 (2004/09/05)の時、CD(m,60)max = 60 x 239.55km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,60)、NCD(m,60)とを、[図80c]に、歪エネルギー密度としてグラフ表示した。左側の縦軸目盛は、歪 エネルギー密度、NCI(m,60)とNCD(m,60)の相対目盛である。又、右側の縦軸の6以上の目盛は、時系列[d(MAG,m)]をMAGとラベル付けした矢印の高さのグラフの目盛である。従って、MAGのグラフは、マグニチュードのオフセット値が 矢印の高さに変換表示され、MAGが6以上になると表示される。右側の 縦軸目盛にLONとラベル付けされている目盛り範囲は、震源パラメータの時系列[d(LON,m)]を点グラフにした、LONの目盛で、経度の134度〜138度に相当する。横軸は、各グラフに共通な時間mで、地震の発生順番を示すインデックスである。対象期間は、m=9000〜11000の2009年1月10日から2011年3月18日である。マグニチュードが6以上で2000年以降の大地震には、その発生年代を、[図77]の時系列[d(MAG,m)]の矢印に表記した。地殻に蓄積された歪エネルギー密度、[図77]のNCD(m,60)とNCI(m,60)とから、抽出できたM6以上のAMR現象は、
2000、2003、2004a、2004b、2007の各大地震と、
AMR-22iとラベル表示した近畿・福井・東南海地方に出現した2011年3月11日のM9東北巨大地震の大規模なAMR現象である。
このAMR現象は、2011年1月12日(m=1338)頃、共にピークに到達し、その直後、急減少を始め、18イベント後に3月11日(m=1356)を迎え、50イベント後の7月11日まで、継続した。更に、NCD(m,60)の急減少は、76イベント後の10月26日まで継続している。この大規模なAMR現象は、上記中国・四国・近畿地方の予知の実施例で取り上げた2011年7月5日和歌山県のM5.5地震をも含めたマグニチュードが5以上で5.5以下の地震7個程含み、2011年3月11日のM9東北巨大地震発生前後に生じた大規模なAMR現象である。この様に、各地域の地震の系は、閉じた系でなく開いた系で、太平洋プレート下端の深発地震を介し、隣接する地域の影響を受ける。従って、巨大地震の予測断層長は、太平洋プレートと大陸プレートの境界、約500kmであると断定できる。
年代 m 年/月/日 緯度 経度 深さ(km) M
2000 2095 2000/01/28 43.006 146.749 58.51 7
2003a 4203 2003/05/26 38.818 141.654 72 7.1
2003b 4540 2003/09/26 41.78 144.078 42 8
2003c 4564 2003/09/26 41.707 143.695 21.4 7.3
2004 5987 2004/11/29 42.944 145.28 48.2 7.1
2005a 6516 2005/08/16 38.151 142.28 41.6 7.2
2005b 6726 2005/11/15 38.031 144.889 83 7.1
2008a 8414 2008/06/14 39.03 140.881 7.8 7.2
2008b 8729 2008/09/11 41.776 144.151 30.9 7.1
2010 9831 2010/06/18 44.597 149.196 30 7
2011a 10307 2011/03/09 38.328 143.28 8.3 7.5
2011b 10411 2011/03/11 38.103 142.861 23.7 9
である。
m = 9795 (2010/05/30)の時、CI(m,30)max = 30 x 26.06時間、
m = 7231 (2006/06/14)の時、CD(m,30)max = 30 x 97.02km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,30)、NCD(m,30)とを、[図80]に、歪エネルギー密度としてグラフ表示した。左側の縦軸目盛は、歪 エネルギー密度、NCI(m,30)とNCD(m,30)の相対目盛である。又、右側のMAG(マグニチュード)とラベル付けした縦軸目盛は、時系列[d(MAG,m)]の6以上の目盛である。従って、[d(MAG,m)]のグラフは、その6以上の値が矢印の高さに変換表示される。右側の縦軸目盛にLON(経度)とラベル付けした目盛りは、震源パラメータの時系列[d(LON,m)]目盛で、その点グラフの範囲は、経度の138度〜150度に相当する。横軸は、各グラフに共通な時間mで、地震の発生順番を示すインデックスである。対象期間は、m=9000〜11000の2009年1月10日から2011年3月18日である。マグニチュードが7以上の大地震と巨大地震には、その発生年代を、[図80]の時系列[d(MAG,m)]の矢印に表記した。地殻に蓄積された歪エネルギー、[図80]のNCD(m,30)とNCI(m,30)の、これらM7以上のAMR現象は、
2010のM7大地震(m=9831、発生日=2010/06/18)の場合、NCD(m,30)とNCI(m,30)のピーク値は、
NCD(m=9817,30) = 0.66 (2010/06/08)
NCI(m=9795,30) = 1 (2010/05/30)
となる。
2011aのM7.5大地震(m=10307、発生日=2011/03/09)の場合、NCD(m,30)とNCI(m,30)のピーク値は、
NCD(m=10238,30) = 0.90 (2011/01/19)
NCI(m=10251,30) = 0.84 (2011/02/01)
となる。
2011aのM7.5大地震は、2011bの東北巨大地震の前震である。前震である事をを検証するために、[図80]の期間をm=10000〜11000とし、更に、2011a地震発生後NCI(m,30)が最小値に到達直後のm=10338からm=11000の間のNCI(m,30)を25倍に拡大し、[図81]に表示する。25倍したNCI(m,30)のピーク値はm=10413で、2011bの東北巨大地震が発生した時刻である。従って、巨大地震2011bのAMR現象を検出できていない。しかし、NCI(m,30)は、INTを累積加算しているので、前震2011aの余震の影響を少なくするために、 加算数を30より少なくすると、そのピーク値を取る時刻をm=10413より早める事が出来る。例えば、累積加算数を20と異なる観測窓を使用する場合、[数60]の2sを20としたNCI(m,20)のピーク値は、東北巨大地震発生時刻のm=10413より3イベント早いm=10410となり巨大地震のAMR現象が検出される。
m = 11244 (2006/08/10)の時、CI(m,100)max = 100 x 13.008時間、
m = 1119 (1998/02/23)の時、CD(m,100)max = 100 x 120.30km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,100)、NCD(m,100)とを、[図82]に、上記日本列島を含む広範囲な領域(緯度16度-52度、経度116度-156度)の歪エネルギー密度としてグラフ表示した。表示期間はm=15300〜17800で、2009/09/19から2011/03/13までである。マグニチュードが7以上の主要な 地震には、d(MAG,m)に発生した年代をラベル表示した。同じ年代に発生した地震には、アルファベットのa、bで区別している。それら年代を表示した地震の震源情報は、
年代 m 年/月/日 緯度 経度 深さ(km) M
2010a 15755 2010/02/27 25.919 128.68 37 7.4
2010b 16115 2010/06/18 44.597 149.196 30 7
2010c 16645 2010/11/30 28.358 139.589 493.9 7.1
2010d 16695 2010/12/22 27.052 143.935 8 7.9
2011a 17287 2011/03/09 38.328 143.28 8.3 7.5
2011b 17361 2011/03/11 38.103 142.861 23.7 9
である。
m = 11241 (2006/08/07)の時、CI(m,30)max = 30 x 16.60時間、
m = 366 (1997/05/10)の時、CD(m,30)max = 30 x 172.77km
となる。これら最大値が1となるように規格化した時系列、[数60]と[数61]のNCI(m,30)、NCD(m,30)とを、[図83]に、歪エネルギー密度の推移グラフとして表示した。[図83]の表示期間は、[図82]と同様にm=15300〜17800で実時間に変換すると2009年9月19日から2011年3月13までである。歪エネルギー密度の変動からのAMR検出は、NCI(m,30)とNCD(m,30)の方が、NCI(m,100)とNCD(m,100)より、高周波の検出成分を多く含む。従って、深発地震2010cと父島近海地震2010dのAMR現象が分離されて検出されている。その検出をより細かく見るために、[図83]の表示期間からm=16300〜17500で2010年8月12日から2011年3月11の期間を[図84]に拡大表示した。この[図84]の各大地震のAMR現象は、
2010c(m=16645、2010/11/30)のM7.1深発地震の場合、
NCI(m,30)のピークは、m=16635 (2010/11/26)
NCD(m,30)のピークは、m=16636 (2010/11/26)
2010d(m=16695、2010/12/22)のM7.9地震の場合、
NCI(m,30)のピークは、m=16695 (2010/12/22)
NCD(m,30)のピークは、m=16690 (2010/12/17)
2011a(m=17287、2011/03/09)のM7.5地震の場合、
NCI(m,30)のピークは、m=17283 (2011/03/07)
NCD(m,30)のピークは、m=17271 (2011/03/03)
となり、各大地震のAMR現象が、それら大地震発生直前に検出されている。
20倍したNCI(m,30)がm=17359、
5倍したNCD(m,30)がm=17358
となる。従って、巨大地震2011bのAMR現象を検出できているので、2011a地震は、巨大地震の前震と断定できる。観測窓の2s=30を、更に縮小すると、上記AMR現象検出を早める事が可能となる。
2 d(c,j)信号を直接、所定の選択、もしくは平滑処理(処理遅れ時間はΔ)し、[c]=[D(c,1),D(c,2),・・D(c,m),・・]とする時系列化手段、
3 時系列[c]からD(c,m)-D(c,m-n)を1次差分値ΔD(c,m,n)として出力し、ΔD(c,m,n)-ΔD(c,m-k,n)を2次差分値Δ2D(c,m,n,k)として出力するか、
ΔD(c,m,n)検出と等価な濾過機能を設定し、[c]から濾過したF1(c,m)を1次差分値として出力し、Δ2D(c,m,n,k)検出と等価な濾過機能を設定し、[c]から濾過したF2(c,m)を2次差分値として出力するかの
1次、2次差分出力手段
4 時刻mを時刻t、t=m-Δ-(n+k)/2、(2で処理遅れが無ければΔ=0)にシフトする時刻シフト手段
5 D(c,m)を変位D(c,t)に、1次差分値を速度V(c,t)に、2次差分値を加速度A(c,t)に変換する変位、速度、加速度変換手段
6 D(c,t)とV(c,t)の振幅と位相を比較検出する手段
7 V(c,t)とA(c,t)の振幅と位相を比較検出する手段
8 D(c,t)とA(c,t)の振幅と位相を比較検出する手段
9 大地震のA(c,t)のCQKとCQT予兆検出手段と、運動の第二法則を用いて、それら予兆を、D(c,t)のCQKDとCQTDに変換し、|ΔD(DEP,t)|= Wとする予測断層幅Wを検出する手段
10 Mの予知手段
11 発生時刻と震源の予知手段
12 大陸プレートの膨らみ検出手段
13 大陸プレートの運動変化検出手段
14 運動変化に連動した膨らみの広がり(長さL)の検出手段
15 海洋プレートの異常加速運動検出手段
16 巨大地震のMと発生時刻の予知手段
17 観測領域の縮小と拡大手段
18 d(INT,j)とd(DEP,j)のa個の累積加算からCI(m,a)とCD(m,a)を検出する手段
19 規格化したNCI(m,a)とNCD(m,a)とを領域の歪エネルギー密度に変換する手段
20 領域の歪エネルギー密度から地殻状態の臨界到達時刻の予知手段
21 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体
22a-22i AMR現象検出個所
Claims (7)
- 観測対象とする領域に発生した地震の震源情報である緯度(LAT)、経度(LON)と深さ(DEP)、震源時(地震の発生時刻)の間隔(INT)、マグニチュード(MAG)を震源要素とし、更に、発震機構解の情報から得る地震の破壊メカニズムに関する要素を、節面1(N1)と節面2(N2)の走行(R)と傾斜角(DEC)、P軸(PA)とT軸(TA)の方位(DIR)と鉛直角(VER)とし、
断層運動の要素として、断層のすべり角をλとし、更に各節面のN1-R、N1-DEC、N2−R2、N2−DEC2、N1-λ、N2−λ2とするか、
または断層面をNとし、N-R、N-DEC、N-λとし、
それら震源情報の要素を表すパラメータをcとして、その選択領域の観測対象とする地震の発生の変化を、対象とする要素cの時系列データで定量化するために、抽出したj番目の地震の要素cに関する地震情報をd(c,j)として数値化する震源要素信号化手段と、
対象とする大陸プレート上の領域の地殻状態の変化を連続観測するために、その領域に稠密に設置されたGPSの各ステーションと、対象領域を含む大陸プレート下に沈み込んでいる海洋プレート上に設置されたGPSの各ステーションで、ステーションの東西方向をE、南北方向をN、上下方向をhとした直交座標系(E,N,h)の座標軸成分(要素)をパラメータcとし、j番目に観測したステーションの位置座標値をd(c,j)として数値化する地殻変動成分信号化手段と、
前記、震源要素信号化手段により抽出されたd(c,j)を、地震の発生を観測した順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは、所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る震源要素の時系列化手段と、
前記、地殻変動成分信号化手段により抽出されたd(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る地殻変動成分の時系列化手段と、
前記、所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段と
前記、震源要素と地殻変動成分の時系列化手段により得られた時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差に比例する量を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段か、
前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、1次差分値ΔD(c,m,n)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、
[C]を濾過したF1(c,m)を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段かの
いずれかの1次差分出力手段と
前記、1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差に比例する量を、時刻mの2次差分値、Δ2D(c,m,n,k)として求めるD(c,m)の2次差分出力手段か
前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、2次差分値Δ2D(c,m,n,k)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、
[C]を濾過したF2(c,m)を、時刻mの2次差分値Δ2D(c,m,n,k)として出力するD(c,m)の2次差分出力手段かの
いずれかの2次差分出力手段と
からなる、地震発生と地殻状態の変化を定量化する方法において、
前記、2次差分出力手段から得た、時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする、2次差分値の時刻シフト手段と、
前記、1次差分出力手段から得た、時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする、1次差分値の時刻シフト手段と、
前記、震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として出力する変位変換手段と
前記、2次差分値の時刻シフト手段から得たΔ2D(c,t,n,k)を、前記変位変換手段から得たD(c,t)に作用している加速度A(c,t)に変換出力する加速度変換手段と
前記、1次差分値の時刻シフト手段から得たΔD(c,t,n)を、前記、時系列の変位変換手段から得たD(c,t)に作用している速度V(c,t)に変換出力する速度変換手段と
前記、変位変換手段で得たD(c,t)と、前記、速度変換手段で得た速度V(c,t)の各振幅と各位相の関係を出力する変位と速度の振幅位相比較検出手段と
前記、速度変換手段で得た速度V(c,t)と、前記、加速度変換手段で得た加速度A(c,t)との各振幅と各位相の関係を出力する速度と加速度の振幅位相比較検出手段と
前記、変位変換手段で得たD(c,t)と、前記、加速度変換手段で得た加速度A(c,t)の各振幅と各位相の関係を出力する変位と加速度の振幅位相比較検出手段と
前記、変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(INT,t)が、次の半周期後に反転する時刻と予測する時刻ta2か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(DEP,t)が、次の半周期後に反転する時刻と予測する時刻ta2か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(DEP,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(DEP,t)が反転したピークを取ると予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(DEP,t)が半周期後に反転したピークを取った時刻をta2とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(DEP,t)が、次の半周期後に反転する時刻と予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(INT,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(INT,t)が反転したピークを取ると予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(INT,t)が半周期後に反転したピークを取った時刻をta2とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(INT,t)が、次の半周期後に反転する時刻と予測する時刻ta3か、
を今にも発生しそうな大地震の発生時刻とする大地震発生時刻の予知手段と
前記、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記、変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする、今にも発生しそうな大地震の震源の予知手段と、
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、
変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か
変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、
変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値かを
今にも発生しそうな大地震の断層幅Wとする予測断層幅Wの検出手段と
予測断層幅Wの検出手段から出力されたWを用いて今にも発生しそうな大地震のマグニチュードの予知手段と
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とから、変位D(h,t)の膨らみを検出する手段か、
もしくは前記変位と速度の振幅位相比較検出手段と、前記、変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とから、変位D(h,t)の膨らみを検出する手段と、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記、変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする、前記、変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記、変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、大陸プレートの運動変化を検出する手段か、
もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する手段と、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段と
海洋プレート上にある領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段か、
もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係からか、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段と
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する手段と
を備えてなる大地震と巨大地震の予知方法。 - 請求項1に記載の大地震と巨大地震の予知方法において、前記、震源要素信号化手段と時系列化手段により抽出した地震発生時刻の時間間隔(INT)の信号d(INT,j)を地震の任意な発生順序mまで、地震の任意な発生個数をa(a=1、2、・・)とし、a個ずつ、移動平均するか、a個ずつ、累積加算しCI(m,a)とする時間間隔累積検出手段か、もしくは、
m番目の地震とm−a番目の地震の発生時刻の差をCI(m,a)とする時間間隔累積検出手段と、
前記、CI(m,a)値がNCI(m,a)となるように規格化するNCI(m,a)の規格化手段と
前記、NCI(m,a)から大地震や巨大地震発生前の自然現象の一つである地震発生の変化が加速されるAMR現象(臨界状態)を検出し大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知方法において
対象とする領域の広さを、順次変える、観測領域の縮小と拡大手段と
前記、震源要素信号化手段と時系列化手段により抽出した地震の震源の深さ(DEP)の信号d(DEP,j)を地震の任意な発生順序mまで、a個ずつ移動平均するか、
a個ずつ、累積加算しCD(m,a)とする震源の深さDEPの累積検出手段と
前記、累積値CD(m,a)がNCD(m,a)となるように規格化するNCD(m,a)の規格化手段と
観測領域の縮小と拡大手段から出力された各領域のCI(m,a)とCD(m,a)のそれぞれ規格化したNCI(m,a)とNCD(m,a)とを、領域に蓄積と開放される規格化された歪エネルギー密度とする歪エネルギー変換手段と
領域の歪エネルギー変換手段から出力した地殻の蓄積歪エネルギー密度が巨大地震、大地震発生直前の臨界状態に到達したかを検出する地殻の臨界状態検出手段と、
前記、地殻の臨界状態検出手段で臨界状態を検出後大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知方法。 - 観測対象とする領域に発生した地震の震源情報である緯度(LAT)、経度(LON)と深さ(DEP)、震源時(地震の発生時刻)の間隔(INT)、マグニチュード(MAG)を震源要素とし、更に、発震機構解の情報から得る地震の破壊メカニズムに関する要素を、節面1(N1)と節面2(N2)の走行(R)と傾斜角(DEC)、P軸(PA)とT軸(TA)の方位(DIR)と鉛直角(VER)とし、
断層運動の要素として、断層のすべり角をλとし、更に各節面のN1-R、N1-DEC、N2−R2、N2−DEC2、N1-λ、N2−λ2とするか、
または断層面をNとし、N-R、N-DEC、N-λとし、
それら震源情報の要素を表すパラメータをcとして、その選択領域の観測対象とする地震の発生の変化を、対象とする要素cの時系列データで定量化するために、抽出したj番目の地震の要素cに関する地震情報をd(c,j)として数値化する震源要素信号化手段と、
対象とする大陸プレート上の領域の地殻状態の変化を連続観測するために、その領域に稠密に設置されたGPSの各ステーションと、対象領域を含む大陸プレート下に沈み込んでいる海洋プレート上に設置されたGPSの各ステーションで、ステーションの東西方向をE、南北方向をN、上下方向をhとした直交座標系(E,N,h)の座標軸成分(要素)をパラメータcとし、j番目に観測したステーションの位置座標値をd(c,j)として数値化する地殻変動成分信号化手段と、
前記、震源要素信号化手段により抽出されたd(c,j)を、地震の発生を観測した順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは、所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る震源要素の時系列化手段と、
前記、地殻変動成分信号化手段により抽出されたd(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る地殻変動成分の時系列化手段と、
前記、所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段と
前記、震源要素と地殻変動成分の時系列化手段により得られた時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差に比例する量を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段か、
前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、1次差分値ΔD(c,m,n)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、
[C]を濾過したF1(c,m)を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段かの
いずれかの1次差分出力手段と
前記、1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差に比例する量を、時刻mの2次差分値、Δ2D(c,m,n,k)として求めるD(c,m)の2次差分出力手段か
前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、2次差分値Δ2D(c,m,n,k)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、
[C]を濾過したF2(c,m)を、時刻mの2次差分値Δ2D(c,m,n,k)として出力するD(c,m)の2次差分出力手段かの
いずれかの2次差分出力手段と
からなる、
地震発生と地殻状態の変化を定量化する装置において、
前記、2次差分出力手段から得た、時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする、2次差分値の時刻シフト手段と、
前記、1次差分出力手段から得た、時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする、1次差分値の時刻シフト手段と、
前記、震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として出力する変位変換手段と
前記、2次差分値の時刻シフト手段から得たΔ2D(c,t,n,k)を、前記変位変換手段から得たD(c,t)に作用している加速度A(c,t)に変換出力する加速度変換手段と
前記、1次差分値の時刻シフト手段から得たΔD(c,t,n)を、前記、時系列の変位変換手段から得たD(c,t)に作用している速度V(c,t)に変換出力する速度変換手段と
前記、変位変換手段で得たD(c,t)と、前記、速度変換手段で得た速度V(c,t)の各振幅と各位相の関係を出力する変位と速度の振幅位相比較検出手段と
前記、速度変換手段で得た速度V(c,t)と、前記、加速度変換手段で得た加速度A(c,t)との各振幅と各位相の関係を出力する速度と加速度の振幅位相比較検出手段と
前記、変位変換手段で得たD(c,t)と、前記、加速度変換手段で得た加速度A(c,t)の各振幅と各位相の関係を出力する変位と加速度の振幅位相比較検出手段と
前記、変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(INT,t)が、次の半周期後に反転する時刻と予測する時刻ta2か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(DEP,t)が、次の半周期後に反転する時刻と予測する時刻ta2か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(DEP,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(DEP,t)が反転したピークを取ると予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(DEP,t)が半周期後に反転したピークを取った時刻をta2とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(DEP,t)が、次の半周期後に反転する時刻と予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(INT,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(INT,t)が反転したピークを取ると予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(INT,t)が半周期後に反転したピークを取った時刻をta2とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(INT,t)が、次の半周期後に反転する時刻と予測する時刻ta3か、
を今にも発生しそうな大地震の発生時刻とする大地震発生時刻の予知手段と
前記、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記、変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする、今にも発生しそうな大地震の震源の予知手段と、
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、
変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か
変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、
変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値かを
今にも発生しそうな大地震の断層幅Wとする予測断層幅Wの検出手段と
予測断層幅Wの検出手段から出力されたWを用いて今にも発生しそうな大地震のマグニチュードの予知手段と
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とから、変位D(h,t)の膨らみを検出する手段か、
もしくは前記変位と速度の振幅位相比較検出手段と、前記、変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とから、変位D(h,t)の膨らみを検出する手段と、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記、変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする、前記、変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記、変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、大陸プレートの運動変化を検出する手段か、
もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する手段と、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段と
海洋プレート上にある領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段か、
もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係からか、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段と
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する手段と
を備えてなる大地震と巨大地震の予知装置。 - 請求項3に記載の大地震と巨大地震の予知装置において、前記、震源要素信号化手段と時系列化手段により抽出した地震発生時刻の時間間隔(INT)の信号d(INT,j)を地震の任意な発生順序mまで、地震の任意な発生個数をa(a=1、2、・・)とし、a個ずつ、移動平均するか、a個ずつ、累積加算しCI(m,a)とする時間間隔累積検出手段か、もしくは、
m番目の地震とm−a番目の地震の発生時刻の差をCI(m,a)とする時間間隔累積検出手段と、
前記、CI(m,a)値がNCI(m,a)となるように規格化するNCI(m,a)の規格化手段と
前記、NCI(m,a)から大地震や巨大地震発生前の自然現象の一つである地震発生の変化が加速されるAMR現象(臨界状態)を検出し大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知装置において
対象とする領域の広さを、順次変える、観測領域の縮小と拡大手段と
前記、震源要素信号化手段と時系列化手段により抽出した地震の震源の深さ(DEP)の信号d(DEP,j)を地震の任意な発生順序mまで、a個ずつ移動平均するか、
a個ずつ、累積加算しCD(m,a)とする震源の深さDEPの累積検出手段と
前記、累積値CD(m,a)がNCD(m,a)となるように規格化するNCD(m,a)の規格化手段と
観測領域の縮小と拡大手段から出力された各領域のCI(m,a)とCD(m,a)のそれぞれ規格化したNCI(m,a)とNCD(m,a)とを、領域に蓄積と開放される規格化された歪エネルギー密度とする歪エネルギー変換手段と
領域の歪エネルギー変換手段から出力した地殻の蓄積歪エネルギー密度が巨大地震、大地震発生直前の臨界状態に到達したかを検出する地殻の臨界状態検出手段と、
前記、地殻の臨界状態検出手段で臨界状態を検出後大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知装置。 - 請求項5の発明は、観測対象とする領域に発生した地震の震源情報である緯度(LAT)、経度(LON)と深さ(DEP)、震源時(地震の発生時刻)の間隔(INT)、マグニチュード(MAG)を震源要素とし、更に、発震機構解の情報から得る地震の破壊メカニズムに関する要素を、節面1(N1)と節面2(N2)の走行(R)と傾斜角(DEC)、P軸(PA)とT軸(TA)の方位(DIR)と鉛直角(VER)とし、
断層運動の要素として、断層のすべり角をλとし、更に各節面のN1-R、N1-DEC、N2−R2、N2−DEC2、N1-λ、N2−λ2とするか、
または断層面をNとし、N-R、N-DEC、N-λとし、
それら震源情報の要素を表すパラメータをcとして、その選択領域の観測対象とする地震の発生の変化を、対象とする要素cの時系列データで定量化するために、抽出したj番目の地震の要素cに関する地震情報をd(c,j)として数値化する震源要素信号化手段と、
対象とする大陸プレート上の領域の地殻状態の変化を連続観測するために、その領域に稠密に設置されたGPSの各ステーションと、対象領域を含む大陸プレート下に沈み込んでいる海洋プレート上に設置されたGPSの各ステーションで、ステーションの東西方向をE、南北方向をN、上下方向をhとした直交座標系(E,N,h)の座標軸成分(要素)をパラメータcとし、j番目に観測したステーションの位置座標値をd(c,j)として数値化する地殻変動成分信号化手段と、
前記、震源要素信号化手段により抽出されたd(c,j)を、地震の発生を観測した順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは、所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る震源要素の時系列化手段と、
前記、地殻変動成分信号化手段により抽出されたd(c,j)を、観測時間の順序j(j=1、2、3、・・、m、・・)に基づいて、直接、もしくは所定の選択処理または平滑化処理を施し、
[C]=[D(c,1)、D(c,2)、D(c,3)、・・、D(c,m)、・・]
を得る地殻変動成分の時系列化手段と、
前記、所定の選択または平滑処理をした場合、その処理に費やす遅れ時間をΔとする遅れ時間検出手段と
前記、震源要素と地殻変動成分の時系列化手段により得られた時系列[C]の時刻mのD(c,m)と、時刻mから任意間隔n(n=1,2,3、・・)だけ離れたD(c,m-n)との差に比例する量を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段か、
前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、1次差分値ΔD(c,m,n)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、
[C]を濾過したF1(c,m)を、時刻mの1次差分値ΔD(c,m,n)として出力するD(c,m)の1次差分出力手段かの
いずれかの1次差分出力手段と
前記、1次差分出力手段から出力された、ΔD(c,m,n)と、時刻mから任意間隔k(k=1、2、・・)だけ離れたΔD(c,m-k,n)との差に比例する量を、時刻mの2次差分値、Δ2D(c,m,n,k)として求めるD(c,m)の2次差分出力手段か
前記、震源要素の時系列化手段か地殻変動成分の時系列化手段により得られた時系列[C]を、前記、2次差分値Δ2D(c,m,n,k)を得るバンドパスフィルター機能と略等価な機能を設定したフィルターを用いて、
[C]を濾過したF2(c,m)を、時刻mの2次差分値Δ2D(c,m,n,k)として出力するD(c,m)の2次差分出力手段かの
いずれかの2次差分出力手段と
からなる、
地震発生と地殻状態の変化を定量化するプログラムにおいて、
前記、2次差分出力手段から得た、時刻mの2次差分値Δ2D(c,m,n,k)を、時刻t= m-Δ-(n+k)/2の2次差分値Δ2D(c,t,n,k)とする、2次差分値の時刻シフト手段と、
前記、1次差分出力手段から得た、時刻mの1次差分値ΔD(c,m,n)を、時刻t= m-Δ-(n+k)/2の1次差分値ΔD(c,t,n)とする、1次差分値の時刻シフト手段と、
前記、震源要素と地殻変動成分の時系列化手段から得た時刻mのD(c,m)を、時刻t= m-Δ-(n+k)/2の変位D(c,t)として出力する変位変換手段と
前記、2次差分値の時刻シフト手段から得たΔ2D(c,t,n,k)を、前記変位変換手段から得たD(c,t)に作用している加速度A(c,t)に変換出力する加速度変換手段と
前記、1次差分値の時刻シフト手段から得たΔD(c,t,n)を、前記、時系列の変位変換手段から得たD(c,t)に作用している速度V(c,t)に変換出力する速度変換手段と
前記、変位変換手段で得たD(c,t)と、前記、速度変換手段で得た速度V(c,t)の各振幅と各位相の関係を出力する変位と速度の振幅位相比較検出手段と
前記、速度変換手段で得た速度V(c,t)と、前記、加速度変換手段で得た加速度A(c,t)との各振幅と各位相の関係を出力する速度と加速度の振幅位相比較検出手段と
前記、変位変換手段で得たD(c,t)と、前記、加速度変換手段で得た加速度A(c,t)の各振幅と各位相の関係を出力する変位と加速度の振幅位相比較検出手段と
前記、変位と加速度の振幅位相比較検出手段で得た加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となり、加速度A(MAG,t)の振幅値が、負の値を取る位相振幅関係が作用している各変位D(DEP,t)、D(INT,t)、D(MAG,t)の振幅関係と位相関係において、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(INT,t)が、次の半周期後に反転する時刻と予測する時刻ta2か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその特別な関係を検出した時刻ta1とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta1から次の半周期後に既に反転した時刻ta2を、時刻ta1のA(DEP,t)が、次の半周期後に反転する時刻と予測する時刻ta2か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(DEP,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(DEP,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(DEP,t)が反転したピークを取ると予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(DEP,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(DEP,t)が半周期後に反転したピークを取った時刻をta2とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(DEP,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(DEP,t)が、次の半周期後に反転する時刻と予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出した後、逆位相にあったA(INT,t)の振幅値がピークを取った時刻をその検出時刻ta1とし、その時刻ta1から、半周期後にA(INT,t)が反転したピークを取った時刻をta2とし、その時刻ta2から、更に、次の半周期後にA(INT,t)が反転したピークを取ると予測する時刻ta3か、
加速度A(DEP,t)と、加速度A(INT,t)とが逆位相となる特別な位相振幅関係を検出直後、逆位相にあったA(INT,t)が、振幅値のピークを取った時刻をその特別な関係を検出した時刻ta1とし、その時刻ta1のA(INT,t)が半周期後に反転したピークを取った時刻をta2とし、前記2次差分検出手段と加速度変換手段で差分を取った間隔時間Δa、Δa=(n+k)/2、だけ進んだ変位D(INT,t+Δa)が、時刻ta2から次の半周期後に既に反転した時刻ta3を、時刻ta2のA(INT,t)が、次の半周期後に反転する時刻と予測する時刻ta3か、
を今にも発生しそうな大地震の発生時刻とする大地震発生時刻の予知手段と
前記、大地震発生時刻の予知手段で得た大地震の発生予測時刻まで、前記、変位変換手段で得たD(LAT,t)、D(LON,t)、D(DEP,t)を線形延長した(LAT、LON、DEP)情報を予測震源とする、今にも発生しそうな大地震の震源の予知手段と、
変位D(MAG,t)が最大値となる時刻tmaxと最小値となる時刻tminの変位D(DEP,tmax)とD(DEP,tmin)との差の絶対値か、
変位D(MAG,t)が最大値となる時刻tmaxと時刻tmaxからD(DEP,t)の振幅値が反転し始める時刻をtminとし、D(DEP,tmax)とD(DEP,tmin)との差の絶対値か
変位D(MAG,t)の最大値付近から最小値付近に変化する時刻に対応して変化するD(DEP,t)の振幅値の差の絶対値か、
変位D(MAG,t)の最大値付近の、D(DEP,t)の振幅値とその半周期後の振幅値との差の絶対値かを
今にも発生しそうな大地震の断層幅Wとする予測断層幅Wの検出手段と
予測断層幅Wの検出手段から出力されたWを用いて今にも発生しそうな大地震のマグニチュードの予知手段と
大陸プレート上にある広範囲の対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と前記1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(h,t)と、速度V(h,t)とから、変位D(h,t)の膨らみを検出する手段か、
もしくは前記変位と速度の振幅位相比較検出手段と、前記、変位と加速度の振幅位相比較検出手段で得た変位D(h,t)と、速度V(h,t)と、加速度A(h,t)とから、変位D(h,t)の膨らみを検出する手段と、
膨らみが検出された大陸プレート上にある対象領域で、時系列の要素パラメータcを(E,N,h)とし、前記、変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする、前記、変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記、変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、大陸プレートの運動変化を検出する手段か、
もしくは、膨らみが検出された大陸プレート上にある対象領域で、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)と変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とV(N,t)、A(E,t)とA(N,t)との関係から、大陸プレートの運動変化を検出する手段と、
前記大陸プレートの運動変化を検出する手段の出力から、大陸プレートの運動変化による膨らみのプレート境界に沿った広がりの長さLを検出する手段と
海洋プレート上にある領域で、時系列の要素パラメータcを(E,N,h)とし、前記変位と加速度の振幅位相比較検出手段を必要としない場合、D(c,m)からD(c,t)とする前記変位変換手段と1次差分値の時刻シフト手段に於いて、2次差分を取る間隔kをゼロとした前記変位と速度の振幅位相比較検出手段とで得た変位D(E,t)とD(N,t)と、それら変位に作用している速度V(E,t)とV(N,t)との関係から、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段か、
もしくは、前記変位と速度の振幅位相比較検出手段と前記変位と加速度の振幅位相比較検出手段で得た変位D(E,t)か変位D(N,t)と、それら変位に作用している速度と加速度のV(E,t)とA(E,t)かV(N,t)とA(N,t)との関係からか、海洋プレートの大陸プレートとの境界方向への速度成分が、通常の速度の約2倍以上に急加速し、その直後、急減速、急停止する異常加速運動を検出する海洋プレートの異常加速運動検出手段と
大陸プレートの運動変化による膨らみと海洋プレートの異常加速運動を検出後、もしくは、膨らみか異常加速運動のいずれかを検出後、今にも発生しそうな巨大地震の発生時刻、マグニチュードを予知する手段と
を備えてなる大地震と巨大地震の予知プログラム。 - 請求項5に記載の大地震と巨大地震の予知プログラムにおいて、前記、震源要素信号化手段と時系列化手段により抽出した地震発生時刻の時間間隔(INT)の信号d(INT,j)を地震の任意な発生順序mまで、地震の任意な発生個数をa(a=1、2、・・)とし、a個ずつ、移動平均するか、a個ずつ、累積加算しCI(m,a)とする時間間隔累積検出手段か、もしくは、
m番目の地震とm−a番目の地震の発生時刻の差をCI(m,a)とする時間間隔累積検出手段と、
前記、CI(m,a)値がNCI(m,a)となるように規格化するNCI(m,a)の規格化手段と
前記、NCI(m,a)から大地震や巨大地震発生前の自然現象の一つである地震発生の変化が加速されるAMR現象(臨界状態)を検出し大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知プログラムにおいて
対象とする領域の広さを、順次変える、観測領域の縮小と拡大手段と
前記、震源要素信号化手段と時系列化手段により抽出した地震の震源の深さ(DEP)の信号d(DEP,j)を地震の任意な発生順序mまで、a個ずつ移動平均するか、
a個ずつ、累積加算しCD(m,a)とする震源の深さDEPの累積検出手段と
前記、累積値CD(m,a)がNCD(m,a)となるように規格化するNCD(m,a)の規格化手段と
観測領域の縮小と拡大手段から出力された各領域のCI(m,a)とCD(m,a)のそれぞれ規格化したNCI(m,a)とNCD(m,a)とを、領域に蓄積と開放される規格化された歪エネルギー密度とする歪エネルギー変換手段と
領域の歪エネルギー変換手段から出力した地殻の蓄積歪エネルギー密度が巨大地震、大地震発生直前の臨界状態に到達したかを検出する地殻の臨界状態検出手段と、
前記、地殻の臨界状態検出手段で臨界状態を検出後大地震や巨大地震発生時期を予知する手段と
を備えてなる大地震と巨大地震の予知プログラム。 - 請求項7の発明は、請求項5ないし請求項6のいずれかに記載の大地震と巨大地震の予知プログラムを記録した記録媒体。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012243002A JP5798545B2 (ja) | 2012-11-02 | 2012-11-02 | 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2012243002A JP5798545B2 (ja) | 2012-11-02 | 2012-11-02 | 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2014092454A JP2014092454A (ja) | 2014-05-19 |
JP5798545B2 true JP5798545B2 (ja) | 2015-10-21 |
Family
ID=50936621
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2012243002A Active JP5798545B2 (ja) | 2012-11-02 | 2012-11-02 | 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5798545B2 (ja) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891884A (zh) * | 2016-06-28 | 2016-08-24 | 中国石油大学(北京) | 一种微地震震源机制反演方法及装置 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP6666023B2 (ja) * | 2015-07-02 | 2020-03-13 | 東京瓦斯株式会社 | 感震センサ及び地震判定方法 |
CN112241561A (zh) * | 2020-06-30 | 2021-01-19 | 同恩(上海)工程技术有限公司 | 一种构件宏观应力指标的监测方法,系统以及存储介质 |
CN113153282A (zh) * | 2021-04-28 | 2021-07-23 | 中国地质大学(北京) | 覆岩水力及岩石力学监测方法及系统 |
CN115343751A (zh) * | 2022-08-04 | 2022-11-15 | 中国地震局地震预测研究所 | 一种活动地块边界带断层震级确定方法及系统 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2849297B2 (ja) * | 1993-01-12 | 1999-01-20 | 鹿島建設株式会社 | 震度予測システム |
JP2858426B2 (ja) * | 1995-09-18 | 1999-02-17 | 耕一 岩田 | 地震前兆の観測方法及び装置 |
JP2948527B2 (ja) * | 1995-11-13 | 1999-09-13 | セラタ ショーセイ | 直下型巨大地震の定量予知方法及びその装置 |
JP4608643B2 (ja) * | 2003-01-17 | 2011-01-12 | 株式会社武田エンジニアリング・コンサルタント | 地震の予知方法、地震の予知システム、地震の予知プログラム及び記録媒体 |
WO2006103929A1 (ja) * | 2005-03-29 | 2006-10-05 | National University Corporation The University Of Electro-Communications | 地震予知システム、地震予知方法、地震予知プログラム及び地震予知プログラムを記録した記録媒体 |
JP4721945B2 (ja) * | 2006-04-14 | 2011-07-13 | 電気化学工業株式会社 | 地殻変動検知装置および地殻変動検知方法 |
-
2012
- 2012-11-02 JP JP2012243002A patent/JP5798545B2/ja active Active
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105891884A (zh) * | 2016-06-28 | 2016-08-24 | 中国石油大学(北京) | 一种微地震震源机制反演方法及装置 |
Also Published As
Publication number | Publication date |
---|---|
JP2014092454A (ja) | 2014-05-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Gischig et al. | Numerical investigation of seismically induced rock mass fatigue as a mechanism contributing to the progressive failure of deep-seated landslides | |
JP5798545B2 (ja) | 大地震と巨大地震の予知方法、予知装置、予知プログラム及び記録媒体 | |
Wesnousky et al. | Integration of geological and seismological data for the analysis of seismic hazard: A case study of Japan | |
Tilmann et al. | The updip seismic/aseismic transition of the Sumatra megathrust illuminated by aftershocks of the 2004 Aceh-Andaman and 2005 Nias events | |
Geist et al. | Implications of the 26 December 2004 Sumatra–Andaman earthquake on tsunami forecast and assessment models for great subduction-zone earthquakes | |
Uchida et al. | A decade of lessons learned from the 2011 Tohoku‐Oki earthquake | |
Crowell et al. | G‐FAST earthquake early warning potential for great earthquakes in Chile | |
de Joussineau et al. | Segmentation along strike-slip faults revisited | |
Ganas et al. | The January 2010 Efpalion earthquakes (Gulf of Corinth, Central Greece): earthquake interactions and blind normal faulting | |
Zhu et al. | Support Vector Machine‐Based Rapid Magnitude Estimation Using Transfer Learning for the Sichuan–Yunnan Region, China | |
JP4608643B2 (ja) | 地震の予知方法、地震の予知システム、地震の予知プログラム及び記録媒体 | |
CN111443379A (zh) | 基于地震波分解提取前兆信息的地震预测时域分析方法 | |
Ritsema et al. | The 2011 tohoku earthquake | |
Berrino | Detection of vertical ground movements by sea-level changes in the Neapolitan volcanoes | |
Bevilacqua et al. | Data analysis of the unsteadily accelerating GPS and seismic records at Campi Flegrei caldera from 2000 to 2020 | |
JP2022087776A (ja) | 地震波を時間領域に分解して地震前兆現象を捉える地震予測時間領域解析システム | |
Tonini et al. | The 29 September 2009 Samoa Islands tsunami: Simulations based on the first focal mechanism solutions and implications on tsunami early warning strategies | |
Chiba | Stress state inferred from b value and focal mechanism distributions in the aftershock area of the 2005 West Off Fukuoka prefecture earthquake | |
Ogiso et al. | Ocean‐wave gradiometry: Visualizing and extracting propagation features of the 15 January 2022 Tsunami wavefield with dense ocean‐bottom pressure gauge arrays | |
SOBIESIAK et al. | Recent seismicity in northeastern Venezuela and tectonic implications | |
Cannavò et al. | Modeling ground deformation associated with the destructive earthquakes occurring on Mt. Etna's southeastern flank in 1984 | |
Groos | Broadband seismic noise: classification and Green's function estimation | |
Hino et al. | Ocean bottom pressure records of the 2011 Tohoku-Oki earthquake | |
Choudhury | The Institute of Seismological Research-A Premier Research Center for Seismology in India | |
Keon et al. | Protecting our shorelines: Modeling the effects of tsunamis and storm waves |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20121228 |
|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20140226 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20141127 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20141216 |
|
A521 | Request for written amendment filed |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20150109 |
|
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: 20150602 |
|
R155 | Notification before disposition of declining of application |
Free format text: JAPANESE INTERMEDIATE CODE: R155 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20150821 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5798545 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |