2019-11-11 マップ差替えます！ 地震の予測マップと時系列予測グラフ 12日の地震列島は茨城沖でM4.5, 長野,大阪,秋田沖でM3.X！ 解説：新しいARIMA時系列予測方式を構築中です！ 本日の予測マップ&グラフは？
オペミスでマップが変わっていなかったので、差替えます m(__)m 16:28
常に最新記事なら [こちら最新！] をアクセスし、ブックマーク！
・ 防災科研さんから午前0時に2日前の詳細データが公開され、もって1年分のデータ解析を行なっています、題名先頭にある日付が解析データ1年分の最終日です、記事は毎日、午前1時〜3時頃にアップされます
・ 地震の予測マップは、毎日更新です
・ 地震の時系列予測グラフは、2019-11-01に更新されました、次回の更新は2019-11-14で、約14日間隔で更新されます
[こよみの計算 - 国立天文台暦計算室] 、[月の地心座標 - 国立天文台暦計算室] より：
・ 11月11日は月齢13.5日 東進Day、月距離 39.8556万km
・ 11月12日は月齢14.5日 東進Day、月距離 39.5425万km
・ 11月13日は月齢15.5日 西進Day 満月です、月距離 39.2094万km 右肩下がりライン上
東進Dayとは新月から満月前日までの日々、西進Dayとは満月から新月前日まで
月距離とは地球から月までの距離で、約35〜40万kmの間を振動しています、38.7万kmが大体の滞在時間（日数）中間となる境界ラインで、上部を凸部、下部を凹部と称しています、この振動の周期は月の公転周期27.3日となっており、1年で13周期あり中央38.7万kmで分けた凸部凹部の数は26個/年（これを1年26期と言う）になり、1期の日数は約14日となります
月の地球に対する公転周期は27.3日で、新月から次の新月前日までの月齢一月は29.5日と、月齢の方が2日少々長いのは、地球の太陽に対する公転運動を月の公転運動が追いかけている為です
経験上、千島海溝、日本海溝、でM6.5以上の西進破壊による巨大地震は：
・ 赤マークで覆われている領域で、西進Dayに発生します
・ 青マークで覆われている領域で、東進Dayに発生します
＝ 最新地震情報１１月１２日（M3.0以上かつ震度1以上）です ＝
Yahooさん [4] より掲載（元データは [気象庁] さん）、マップ上★が震源位置
★ 11月12日08時25分、長野南部でM3.5、深さ10km、震度3
★ 11月12日12時53分、大阪南部でM3.1、深さはごく浅く、震度1
★ 11月12日15時36分、秋田沖でM3.7、深さ10km、震度1
★ 11月12日21時55分、茨城沖でM4.5、深さ50km、震度3
＊ 12日は、大阪南部でM3.1でした、場所的に見て和歌山との県境で、中央構造線直下の断層による地震と思われます、ですが、深さがごく浅くでM3.1で震度が1、という事は揺れが小さく、地盤はシッカリしているものと思われます
解説： 新しいARIMA時系列予測方式を構築中です！
＊ 今日は、新しいARIMA予測方式についてご説明致したく、お付き合い頂けますよう、よろしくお願い致します
＊ その前に、現在の西域における過去予測実績にバグがありました、西域におけるM5.0以上地震のカウントにバグがあり、M5.0以上が実態より多く見えていていました
修正した結果、西域についても>=M4.8条件にてARIMAを掛ける事になり、前回まではHIt_P予測で、Hit_P予測実績が72.0%でしたが、今回の修正でHit_N予測となり、Hit_N予測実績が69.2%となっています、尚、本日より予測本文を修正してあります
また>=M4.8とした結果、地震数は31から17と少なくなりました
＊ それにしても、現在の予測方式で、東域の予測実績を見るとHit_P予測実績が16で空振りKaraburi数が10で、合わせて26で、予測段階ではこの1年26期すべてHit_P予測をしていた事になります
同様にバグ修正をした西域の予測実績を見ると、Hit_N予測実績が18で見逃しMinogashi数が8と、足して26で、これは予測段階ではすべてHit_N予測をしていた事になります
私は、現在のARIMA予測方式に限界を感じていまして、新しい但しARIMAに基づく方式を検討しています
＊ 従来のARIMA方式は1年26期におけるM4.8以上の地震数を食わせて予測させていましたが、新しい方式では1年26期の各期を年として1年26期26年のデータを作成しARIMAに食わせます
従来の方式をARIMA26期方式、新しい方式をARIMA26期26年方式と称する事にします、新しい方式では約2年分のデータを食わせる事になるので、まぁいいか、と思っています
＊ 新しいARIMA26期26年方式による2019-11-01段階での27期と28期の予測は出ていて、
東域では：
中域では：
西域では：
となっています
従来と決定論的に異なるのはゼロが無いという事です、ここで横軸は1〜26迄あるのですが、各々が各期を最終とする1年26期の抽象化されたON数の総和となっているからで、最大で26となり（26を越える事はない！）、最小でも西域で直近の8となっています
＊ まだ各種メッセージを組み込んでいないのと、このARIMAグラフ結果からHit_P予測かHit_N予測かを抽出し、過去予測実績のプログラムに組み込んで実績％を得る、等、未だこれからの作業がかなり有りますが、次回の更新2019-11-14までには間に合わせたい、と考えております
これにより適切なHit_P予測とHit_N予測の混在が実現出来るか？という所です
尚、Rでpng出力が出来ました、ので、多少ですが工数が少なくなっています
以上、分かりづらい内容にお付き合い頂きまして、誠にありがとう御座いました、感謝です
そして被災地は今... [happy-ok3の日記] 地震・豪雨・台風と、被災地の現状をレポートするhappy-ok3 さんの考えさせられるブログです、関心を持ち続けて欲しい と
＝ 以下、2019-11-11 迄データ1年分による本日の地震の予測マップ&時系列予測グラフ ＝
赤マークは東進圧力、青マークは西進圧力を示す圧力方向の解析表示です
救急マークがM5.0以上の地震を予測する注意ポイントで、6kmマップにあります
[防災科学技術研究所 Hi-net 高感度地震観測網]、[気象庁｜震源データ] を参照しています
＝ 地震の予測マップ・ピッチ36kmマップです ＝
東進西進圧力表示・ピッチ36km予測マップです
ピンクの小さな★マークは、南海トラフ巨大地震発生ヶ所で、西から、1854安政南海M8.4、1946昭和南海M8.4、1707宝永M8.6、1944昭和東南海M8.2、1854安政東海M8.4
＝ 地震の予測マップ・ピッチ6kmマップとポイント予測マップと月距離地震グラフと時系列地震予測グラフです ＝
東進西進圧力表示・ピッチ6km予測マップ東域です 凡例は36kmマップと同じ
次がポイント予測・東域 救急マーク がM5.0以上の発生予測注意ポイントです
根室沖M7.8〜8.5の確率80%、青森東方沖及び岩手沖北部M7.0〜7.5の確率90% 、宮城沖M7.0〜M7.5の確率90%、福島沖M7.0〜M7.5の確率50% [海溝で起こる地震 | 地震本部] 2019年1月1日算定基準日の30年間確率、以下同様
凸部と凹部比率： East ["47.8%", "52.2%", 46] Days ["49.3%", "50.7%", 365]
時系列地震予測グラフ（Rのforecast機能）東域は、2019-11-01 の予測結果で、
ミニ解説・東域：
38.7万kmを境目とし1年の凸部凹部を26期で示し、続く27期と28期を予測、予測値実数は整数化し判断、topは凸部、btmはボトムで凹部
Forecast East M>=4.8: next period 27th-> 2.0, next of the next 28th-> 1.7
27th(top)-> 2019/11/02 to 2019/11/14, 28th(btm)-> 2019/11/15 to 2019/11/29
author selected ARIMA model order-> (0,0,1) out of 6 trials
東域では11月2日〜11月14日迄の間にM4.8以上が2回、11月15日から11月29日の間で1回発生するだろう、とHit_P予測しています、下記の過去直近1年の予測実績より、Hit_P予測された期でM5.0以上が発生する確率は61.5%です、次の更新は 2019-11-14 です
過去1年26期のM5.0発生実績と予測を比較する抽象化サマリーは、
Hit_P抽象化とは「来ると予測してM5.0以上が来た」期をONとする数え方で個数を問わない、Hit_P Rate = Hit_P数 /（Hit_P数 + 空振りKaraburi数）即ちスウィングして当たった%
過去直近１年の予測実績を示す抽象化サマリー東域は、M5.0以上について：
Hit_P : 16 periods Matched with Hit_P condition
Hit_N : Nothing, not detected
Karaburi : 10 periods Matched with Karaburi condition
Minogashi : Nothing, not detected
Hit_P Rate of the East in the year 2018/11/09 ... 2019/11/01-> 16/(16+10) = 61.5%
東進西進圧力表示・ピッチ6km予測マップ中域です 凡例は36kmマップと同じ
次がポイント予測・中域 救急マーク がM5.0以上の発生予測注意ポイントです
青森東方沖及び岩手沖北部M7.0〜7.5の確率90% 、宮城沖M7.0〜M7.5の確率90%、福島沖M7.0〜M7.5の確率50%、茨城沖M7.0〜M7.5の確率80%、相模トラフ南関東直下地震M6.7〜7.3の確率70% [海溝で起こる地震 | 地震本部]
凸部と凹部比率： Midl ["45.7%", "54.3%", 46] Days ["49.3%", "50.7%", 365]
時系列地震予測グラフ（Rのforecast機能）中域は、2019-11-01 の予測結果で、
ミニ解説・中域：
38.7万kmを境目とし1年の凸部凹部を26期で示し、続く27期と28期を予測、予測値実数は整数化し判断、topは凸部、btmはボトムで凹部
Forecast Midl M>=4.8: next period 27th-> 1.0, next of the next 28th-> 1.0
27th(top)-> 2019/11/02 to 2019/11/14, 28th(btm)-> 2019/11/15 to 2019/11/29
auto selected ARIMA model order-> (0,1,1) out of 18 trials
中域では11月2日〜11月14日迄の間にM4.8以上が1回、11月15日〜11月29日迄の間でも1回発生するだろう、とHIt_P予測しています、下記の過去直近1年の予測実績より、Hit_P予測された期でM5.0以上が発生する確率は68.0%です、次の更新は 2019-11-14 です
過去1年26期のM5.0発生実績と予測を比較する抽象化サマリーは、
Hit_P抽象化とは「来ると予測してM5.0以上が来た」期をONとする数え方で個数を問わない、Hit_P Rate = Hit_P数 /（Hit_P数 + 空振りKaraburi数）即ちスウィングして当たった%
過去直近１年の予測実績を示す抽象化サマリー中域は、M5.0以上について：
Hit_P : 17 periods Matched with Hit_P condition
Hit_N : 1 periods Matched with Hit_N condition
Karaburi : 8 periods Matched with Karaburi condition
Minogashi : Nothing, not detected
Hit_P Rate of the Midl in the year 2018/11/09 ... 2019/11/01-> 17/(17+8) = 68.0%
東進西進圧力表示・ピッチ6km予測マップ西域です 凡例は36kmマップと同じ
次がポイント予測・西域 救急マーク がM5.0以上の発生予測注意ポイントです
南海トラフM8〜M9クラスの確率70%〜80%、日向灘M7.1前後の確率70〜80% [海溝で起こる地震 | 地震本部]
ピンクの小さな★マークは、南海トラフ巨大地震発生ヶ所で、西から、1854安政南海M8.4、1946昭和南海M8.4、1707宝永M8.6、1944昭和東南海M8.2、1854安政東海M8.4
凸部と凹部比率： West ["47.1%", "52.9%", 17] Days ["49.3%", "50.7%", 365]
時系列地震予測グラフ（Rのforecast機能）西域は、2019-11-01 の予測結果で、
ミニ解説・西域：
38.7万kmを境目とし1年の凸部凹部を26期で示し、続く27期と28期を予測、予測値実数は整数化し判断、topは凸部、btmはボトムで凹部
Forecast West M>=4.8: next period 27th-> 0.6, next of the next 28th-> 0.6
27th(btm)-> 2019/10/19 to 2019/11/01, 28th(top)-> 2019/11/02 to 2019/11/14
author selected ARIMA model order-> (0,0,1) out of 6 trials
西域では11月2日から11月14日の間でM4.8以上は０回、11月15日から11月29日の間でも０回、とHit_N予測しています、Hit_N予測された期でM5.0以上が発生しなかった確率は69.2%です、次の更新は 2019-11-14 です
過去1年26期のM5.0発生実績と予測を比較する抽象化サマリーは、
Hit_N抽象化とは「来ないと予測してM5.0以上が来なかった」期をONとする数え方、Hit_N Rate = Hit_N数 /（Hit_N数 + 見逃しMinogashi数）即ち見送り判断して当たった%、見逃し数では抽象化しているので個数は問わない
過去直近１年の予測実績を示す抽象化サマリー西域は、M5.0以上について：
Hit_P : 0 periods Matched with Hit_P condition
Hit_N : 18 periods Matched with Hit_N condition
Karaburi : Nothing, not detected
Minogashi : 8 periods Matched with Minogashi condition
Hit_N Rate of the West in the year 2018/11/09 ... 2019/11/01-> 18/(18+8) = 69.2%
＝ まとめです ＝
＊ 2017年の放出エネルギーは過去最低 、2018年は増加、2019年も現在増加しています [2019-10 地震の予測マップ ここ２６年間の北西太平洋地域・地震放出エネルギー推移を測定する！ - 地震の予測マップ]
＊ ここで赤い救急マークは、M5.0以上の地震ヶ所をピンポイントで予測しています
＊ 時系列予測では、各域において直近将来約14日間の地震発生を予測しています
＊ どちらも、2019山形沖地震M6.8、2018胆振地震M6.7、2018大阪北部地震M6.1、2016熊本地震の前震M6.5、2011年東北太平洋沖地震の前震M7.3、を正しく予測しています
政府の地震調査委員「超巨大地震：北海道沖、M９地震予測 本州にも被害の恐れ」といった報道をするのであれば、注意喚起領域でM５．５クラスが来た時には「注意喚起情報」を出されるべきではないか、と思います。
政府の地震調査委員会
最後まで読んで頂き、ありがとう御座いました
・ 東進西進の原理、東進西進の識別方法、等の説明はこちら [テクニカル事項]
・「地震の予測マップ」のデータ更新タイミングの説明はこちら [データ更新タイミング]
・ 国土地理院さん提供の地殻変動マップはこちら [地殻変動情報] javaがインストールされている必要があります
・「太陽黒点数の推移を追う！」は別ブログへ [太陽黒点数の推移を追う]
＝ 以上です ＝
謝辞： 本予測は「気象庁・文部科学省が協力してデータを処理した結果」の「気象庁一元化処理震源要素」データ一年分（暫定）を「防災科学技術研究所」サイトよりダウンロードして解析しています。 [2] このデータによって初めて一般にリアルタイム解析が可能となったもので、構築にご尽力頂きました各国立大学、各官庁と関連する機関、都道府県と関連する機関、等の関係各位殿に深く謝意を述べさせて頂きます。
免責： 本予測は個人の推論によるもので、プログラムバグやデータ解釈ミス等も含め、ここで表示された結果について何ら責任を負うものではありません。
引用：
[1] スロースリップ - Wikipedia プレートがゆっくりと移動し大きな破壊を伴わずにエネルギー解放する現象ですが、プレート深部には応力歪が伝搬され、これが原因で深部のスラブ内領域で通常の地震が発生します。
[3] 気象庁|震源データ
[6] 国立天文台 太陽観測科学プロジェクト 三鷹太陽地上観測
