JP5090955B2

JP5090955B2 - 発雷予測装置およびその方法、プログラム

Info

Publication number: JP5090955B2
Application number: JP2008031989A
Authority: JP
Inventors: 浩史佐藤; 聡杉山
Original assignee: Nippon Telegraph and Telephone Corp
Current assignee: Nippon Telegraph and Telephone Corp
Priority date: 2008-02-13
Filing date: 2008-02-13
Publication date: 2012-12-05
Anticipated expiration: 2028-02-13
Also published as: JP2009192312A

Description

本発明は、機械学習を用いて気象データを元に発雷を直前に予測する方法、特に発雷予測のための機械学習の訓練データを作成して用いる発雷予測装置およびその方法、プログラムに関する。

従来の発雷予測では、雷センサを用いて直前の発雷位置を知り、気象レーダを用いて雷雲（積乱雲）の動きを知り、それらを組み合わせて次の発雷位置を予測する方法や（非特許文献１）、気象レーダにより解析した雷雲パターンを過去の発雷時のパターンと照らし合わせ、その類似度が高い場合に発雷するとみなす方法がとられている（非特許文献２）。

また、数値予報による気圧面毎の温度や湿度などから大気安定度を計算し、不安定な場合に発雷可能性が高いと判断する方法もある（非特許文献３）。
株式会社フランクリンジャパン：「Lightning Scope」（登録商標）http://www.fjc.co.jp/service/cs/ 木瀬若桜ほか「複数のセンサデータを用いた事例検索による発雷予測支援方式」信学論Ｂ．Ｖｏｌ．Ｊ８３−Ｂ，Ｎｏ．１１，ｐｐ．１５９８−１６０７（２０００−１１）大野久雄著「雷雨とメソ気象」（東京堂出版）第５章第４節

しかしながら、従来の技術における雷センサを使う方法は当然ながら雷センサを必要とするが、一機数千万円程度する上、標定には最低３機が必要となるので、簡単に導入できるものではなかった。

また、雷情報を提供する商用サービスもあるが、リアルタイムでそれらを受けるためには月に百万円単位の費用がかかるので、いずれにせよ安価ではない。

また、発雷が始まらないと次が予測できないため、積乱雲が発達した直後の最初の発雷を予測することが出来ないという問題もある。

また、過去の雷雲パターンと照らし合わせる方法については、上述の「最初の発雷を予測できない」という問題はないが、独自に気象レーダを運用しなければ入手できない特殊なデータを必要としているため、実施が難しく、あくまでも「予測支援」なので、人間による雷発生予測の手助けとしてしか使えなかった。

また、大気安定度を使う方法は、数時間単位での大ざっぱな発雷予測であり、さらにあくまでもポテンシャルを測る手段でしかなく、発雷というより積乱雲の出来やすさを見ているに過ぎないため、発雷を直前に正確に予測することはできなかった。さらに、元となる数値予測自体も観測結果ではなく予報結果であり、位置精度も粗いので、結果もおおまかなものでしかなかった。

また、コンピュータに学習させる「機械学習（Machine learning）」を使った分類問題という観点で雷予測を考えると、雷は全体から見れば頻度の非常に小さい気象現象であり、コンピュータによる機械学習をするために必要な訓練データのバランスが悪くなりがち、すなわち、発雷した事例と比較して発雷しない事例が圧倒的に多くなるので、そのまま学習しても分類器が発雷しない事例に強く引きずられ、結果として高い精度が出ない、という問題があった。

また、精度良く予測できたとしても、予測の元データとなるリアルタイムの気象データの受信が遅れた場合、予測のリードタイムが短くなり、予測情報の価値が著しく下がるという問題もあった。

本発明は上述の課題に鑑みてなされたものであり、その目的は、機械学習後の予測において雷データを必要とせず、高価な雷センサの設置の必要がなく、また、雷現況情報の受信契約を不要にすることにある。

また、機械学習後の予測において雷データを参考にせずに、その雷雲の最初の発雷を予測することにある。

また、一般的に安価に入手可能な気象データのみを入力として、独自のレーダなどの設置を不要にすることにある。

また、気象レーダによるエコー値を予測に用いて、雷も同様の精度で直前に正確に予測することにある。

また、訓練データのバランスをとるための手段を提供して、機械学習時に正例と負例をバランスよく学習させ、予測精度を向上させることにある。

また、訓練データのバランスをとる際に、積乱雲の発達の度合いという発雷に直接的なデータに沿ってエントリ数を調整して、学習に重要な判断の際どいエントリを重点的に学習させ、結果として分類器の精度を向上させることにある。

また、機械学習手法としてＳＶＭを採用して、汎化性の高い分類を可能とすることにある。

また、発雷を予測する前にパターン変化予測を行って、予測のリードタイムをより稼ぐことにある。

課題を解決するために、請求項１に記載の本発明は、メッシュ状に区切られた複数の領域ごとに発雷を直前に予測するための発雷予測装置において、前記領域のそれぞれごとに過去の時系列に沿って気象状況を記録した気象データに基づいて生成されたエントリと、前記エントリの総数において発雷の有無の割合が所定の割合となるように発雷の無かったエントリを間引いて生成した訓練データと、前記訓練データを用いてコンピュータにより機械学習を実行して得た学習結果に基づいて構成され、現在の気象データに対して発雷を予測するための分類器と、を備え、現在の気象データの一部もしくは全部にパターン変化予測を適用し、その予測結果を前記分類器へ入力する。

また、請求項２に記載の本発明は、請求項１において、前記訓練データにおいて、積乱雲の発達の度合いと相関のある値を参照し、発達の度合いが低い時は前記間引きを多くし、発達の度合いが高いときは前記間引きを少なくする。

また、請求項３に記載の本発明は、請求項１または２において、前記機械学習の手法としてＳＶＭ（Support Vector Machine）を用いる。

また、請求項４に記載の本発明は、メッシュ状に区切られた複数の領域ごとに発雷を直前に予測するための発雷予測方法において、前記領域のそれぞれごとに過去の時系列に沿って気象状況を記録した気象データに基づいたエントリを生成するステップと、前記エントリの総数において発雷の有無の割合が所定の割合となるように発雷の無かったエントリを間引いて訓練データを生成するステップと、前記訓練データを用いてコンピュータにより機械学習を実行して得た学習結果に基づいて、現在の気象データに対して発雷を予測するための分類器を構成するステップと、を有し、現在の気象データの一部もしくは全部にパターン変化予測を適用し、その予測結果を前記分類器へ入力する。

また、請求項５に記載の本発明は、請求項４において、前記訓練データにおいて、積乱雲の発達の度合いと相関のある値を参照し、発達の度合いが低い時は前記間引きを多くし、発達の度合いが高いときは前記間引きを少なくする。

また、請求項６に記載の本発明は、請求項４または５において、前記機械学習の手法としてＳＶＭ（Support Vector Machine）を用いる。

また、請求項７に記載の本発明は、メッシュ状に区切られた複数の領域ごとに発雷を直前にコンピュータを用いて予測させるための発雷予測プログラムであって、前記領域のそれぞれごとに過去の時系列に沿って気象状況を記録した気象データに基づいたエントリを生成するステップと、前記エントリの総数において発雷の有無の割合が所定の割合となるように発雷の無かったエントリを間引いて訓練データを生成するステップと、前記訓練データを用いてコンピュータにより機械学習を実行して得た学習結果に基づいて、現在の気象データに対して発雷を予測するための分類器を構成するステップと、をコンピュータに実行させ、現在の気象データの一部もしくは全部にパターン変化予測を適用し、その予測結果を前記分類器へ入力する。

本発明によれば、機械学習後の予測において雷データを必要としないので、高価な雷センサの設置が必要がなく、また、雷現況情報の受信契約の必要がなくなる。

また、機械学習後の予測において雷データを参考にしないので、その雷雲の最初の発雷を予測できる。

また、一般的に安価に入手可能な気象データのみを入力としているので、独自のレーダなどを設置する必要がなくなる。

また、気象レーダによるエコー値を予測に用いた場合は、このデータが現在１０分間隔で配信されており位置精度も１ｋｍ単位なので、雷も同様の精度で直前に正確に予測することができる。

また、訓練データのバランスをとるための手段を提供しているので、機械学習時に正例と負例をバランスよく学習することができ、予測精度を向上させることができる。

また、訓練データのバランスをとる際に、積乱雲の発達の度合いという発雷に直接的なデータに沿ってエントリ数を調整するので、学習に重要な判断の際どいエントリを重点的に学習することができ、結果として分類器の精度を向上させることができる。

また、機械学習手法としてＳＶＭを採用しているので、汎化性の高い分類が可能になる。

また、発雷を予測する前にパターン変化予測を行っているので、予測のリードタイムをより稼ぐことができる。

＜第１の実施の形態＞
本発明の実施の形態では、気象業務支援センタが提供している過去の一般的な気象データ（降水データおよび数値予報データなど）および、雷情報提供会社などが提供する過去の雷情報（発雷の位置、時刻、特性などが収録されている）を用いて、正解ラベル付き訓練データを作成し、「教師あり機械学習（Supervised learning）」の手法を適用して発雷を予測するための分類器を作成する。

この分類器にしたがって、気象業務支援センタなどよりリアルタイムで配信される気象データ、および、その気象データをパターン変化予測にかけた結果を入力として用い、発雷予測を行う。

なお、「機械学習」とはサンプリングした多数のデータについて解析を行い、その解析結果に基づいて対象となるデータ群における規則性を抽出するための技術であり、コンピュータを用いて人間が備える学習能力と同様の機能を実現するためのシステムや方法、プログラムである。

最初に図１に示すのは、本発明の実施の形態を説明するための全体の処理フローである。

図１のように、まず学習を行い、その後、予測を繰り返す、といった順になる。必要に応じて学習をし直したりすることもできる。発雷の予測のための機械学習が開始されると（Ｓ１）、図２を参照して後述する学習フローが実行される（Ｓ２）。この学習結果に基づいて図３を参照して後述する発雷の予測フローを実行する（Ｓ３）。予測フローが終了するまで繰り返され、終了後は機械学習も終了する（Ｓ４、Ｓ５）。

次の図２は、図１中に示したＳ２の学習フローを説明するためのフローチャートを示している。

この図２に示す「教師あり機械学習」の学習フローでは、教師すなわち正解ラベルのついた訓練データが必要となる。本実施の形態においては、ある領域のある時刻の、周辺および前後のメッシュの気象データ群と、その領域でその時刻に発雷があったかどうかのペアが訓練データの１エントリとなる。

まず学習フローが開始されると（Ｓ１０）、データの入力がされる（Ｓ１１）。この際に入力されるデータは気象データ群である。次にＳ１１にて入力された気象データ群に対応する正解ラベルが付与される（Ｓ１２）。ここで付与される正解ラベルは気象データ群におけるメッシュで区切られた領域ごとの発雷の有無と時刻のデータである。

次に、間引き処理が行われる（Ｓ１３）。この間引き処理は正解ラベルが付与された訓練データの中から、後述するように、エコー頂高度が低い場合を多めに間引き、エコー頂高度が高い場合を少なめに間引く処理を行う。次にＳＶＭ（サポートベクターマシン：Support vector machine）学習を行う（Ｓ１４）。そして学習結果が出力され（Ｓ１５）、学習フローが終了する（Ｓ１６）。

また、図３は、図１中に示した予測フローを説明するためのフローチャートを示している。

この図３において予測フローが開始されると（Ｓ２０）、データの入力がされる（Ｓ２１）。ここで入力されるデータは発雷の予測対象であり観測された気象データ群である。次にパターン変化予測が入力された気象データ群に対して実行される（Ｓ２２）。パターン変化予測の方法としては、たとえば単位時間前のデータと現在のデータの変化を認識し、その変化量を現在のデータに適用する、といった技術である。

次に、後述の分類器にて発雷予測を行う（Ｓ２３）。発雷の予測結果が出力され（Ｓ２４）、予測フローが終了する（Ｓ２５）。

次の図４は、本実施の形態の処理フローにおける、各データと各処理との関連を説明するための説明図を示している。この図４には、説明のための構成として、過去気象データＤＢ１と、過去の雷データ２と、過去の気象データ３と、正解ラベル付き訓練データ４と、機械学習のステップＳ３０と、学習結果５と、気象データ６と、分類器７と、予測結果８と、が示されている。

この構成において、雷データ２には、一般に時刻と位置（緯度経度）および雷の電流データなどが羅列してある。従って、この雷データを参照しながら対象のメッシュ（地理的な領域）において発雷があったかどうかを判定すれば良い。

また、図４に示す学習結果５は、コンピュータ上のファイルデータとして得られる。これを分類器（ソフトウェア）７に読み込ませ、分類処理を行うこととなる。

本実施の形態では、気象データ３として図５に示す「エコー強度（気象レーダーで観測される換算降水強度）」と「エコー頂高度（気象レーダーで観測される降水エコーの高さ）」の２つを採用する。３つ以上になっても考え方は同様である。これらは両者とも安価にリアルタイムに入手することができる。

入力する気象データ３としては、対象となるメッシュＸのエコー強度、エコー頂高度、メッシュＸの周辺メッシュのエコー強度の平均値、エコー頂高度の平均値、それぞれの値の１０分前からの差分、を用いることとする。なお、上で「時刻」と書いているが、実際は瞬間ではなくある単位時間（たとえば１０分間）を意味する。以下簡単のため「時刻」と表現するが同様である。

次に、過去の気象データ３より、ある領域、ある時刻の周辺メッシュの各値を抽出し整列させる。その領域でその時刻に発雷があれば正例（＋１）、なければ負例（−１）のラベルを貼り付ける。同様の処理を繰り返すことで、正解ラベル付きの訓練データ４の候補となる。

また、発雷はデータ全体から見ればまれな事象なので、一般には負例（−１）が圧倒的に多くなる。訓練データに含まれる正例（＋１）と負例（−１）の比率は同程度の方が良いとされており、このままでは効果的な学習に障害をきたすので、訓練データの間引き処理を行う。

一般に積乱雲が発達すればするほど発雷が起きる可能性が高い。たとえば、エコー頂高度（降水粒子が存在する最高度すなわち積乱雲の高さ）の値が大きければ発雷が起きる可能性が高い。そのため、負例のエントリのうち学習する立場から価値が高いと言えるのは、エコー頂高度が高いものである。これは、負例と正例を見分ける上で際どいところに位置するエントリと言えるからである。

そこで、訓練データから負例のエントリを間引く際に、エコー頂高度が低い場合を多めに間引き、エコー頂高度が高い場合を少なめに間引く。たとえば、訓練データ候補において、正例エントリが１，２００件、負例エントリが１，０００，８００件あったとする。これについて正例：負例を仮に６：４にしたいと考えた場合、負例エントリを１，０００，０００件間引く必要がある。

ここで負例エントリのエコー頂高度を量子化する。たとえば０以上４未満（ｚ１件）、４以上７未満（ｚ２件）、７以上１０未満（ｚ３件）、１０以上（ｚ４件）、の４段階に量子化する。ｚ１〜ｚ４は十分大きいとする。このｚ１〜ｚ４に対し、単調減少の数列、たとえば、４，３，２，１を順に割り振る。この数列は０以上の整数からなる列とし、狭義の単調減少である必要はない。

次に、ｙ１：ｙ２：ｙ３：ｙ４＝４：３：２：１かつ、ｙ１＋ｙ２＋ｙ３＋ｙ４＝１，０００，０００といった連立方程式を解く。このｙ１〜ｙ４がそれぞれｚ１〜ｚ４から間引くべき数となる。あとはランダムにそれぞれのエントリから指定された件数だけ間引けば良い。この結果が図６に示す訓練データとなる。以上が間引き処理（Ｓ１３）である。

次に、このように間引きを行った訓練データをＳＶＭ学習により処理する。ＳＶＭ学習により、正例と負例を二分する超平面すなわち分類器７が学習され構築される。図７に示すのは、この分類器７を構成する超平面の概念図である（「超」平面のため、正確な図は描けない）。なお、ＳＶＭ学習にはソフトマージンという考え方があり、必ずしも超平面が正例と負例を完全には二分せず、そこに多少の妥協を入れることができる。本実施の形態においては、これについてはどちらを選択してもよい。以上が学習フェーズ（Ｓ２、Ｓ３０）である。

次に、予測対象となる領域に対し、その周辺や前後も含めたメッシュ領域にそれぞれ気象データ６を入力する。たとえば１０分後を予測する場合、これらデータの１０分後の値をパターン変化予測を用いて求める。その際、現在と過去のデータを入力とする。たとえば１０分前と現在のデータから、１０分後のデータを予測する、などができる。図８はこの概念図であって、四角い升目の一つが領域を指している。また、ハッチングされた領域が発雷の予測される領域である。なお、パターン予測のために入力する気象データ６にはエコー強度とエコー頂高度を含んでおり、その両方か、あるいはいずれか一方を用いてもよい。

ここで実行するパターン変化予測は、一般に、単位時間前のデータと現在のデータの変化を認識し、その変化量を現在のデータに適用するものが多い。高度なものだと、流体力学の方程式により変化を予測するものもある。たとえば、Dynamics Texture法（文献：境野英朋ほか、“流体モデルによる時系列レーダパターン変化の予測方法”、電子情報通信学会論文、ｖｏｌ．Ｊ８２−Ａ，ｎｏ．４，ｐｐ．５５０−５６２，１９９９）などが適用できる。なお、ここで実行するパターン予測は、「パターン予測した結果の変化量を得る」場合と「変化量をパターン予測する」場合のどちらをも含んでいる。

次に、この１０分後を予測したデータ組に対し、分類を行う。この分類には事前に学習しておいた分類器７を適用する。分類器７では、そのデータ組が正例と判断されるか負例と判断されるかに応じて、その領域で雷が発生するかどうかを判断する。以上で、雷予測のフローの説明を終える（Ｓ３、Ｓ２３）。

＜第２の実施の形態＞
第１の実施の形態にて説明した分類器７で構成された予測エンジンたる発雷予測部１１が中心となっている。しかし実際のサービス運用においては、図９に示すように気象データ受信部１０、発雷予測部１１、予測結果解釈部１２、予測結果変換部１３、予測結果配信部１４なども合わせて運用されることを想定している。

まず、気象データ受信部１０は、気象業務支援センタなどから気象データの配信を受け、必要に応じてデータ形式の変換などを行う。

次に、予測結果解釈部１２は、分類器７が予測した結果を元に発雷による災害発生の可能性や規模を計算し、警報などを出す必要性を判断する。

予測結果変換部１３は、予測結果や前記の予測結果解釈部１３の解釈結果を必要に応じて画像や文章に変換する。

予測結果配信部１４は、前記の予測結果変換部１３の変換結果をウェブサーバや警報メール送信サーバなどを通じて配信する。

また、予測結果配信部１４からの配信を受けて、ユーザ側での処理として、雷サージに敏感なコンピュータシステムを停止したり、電源系統をバッテリーや自家発電に切り替えるなどの処理およびその自動化も想定している。

また、本実施の形態は発雷の予測としているが、雷のように全体からの頻度が小さい現象、たとえば雹や地震などにも適用可能である。その場合、積乱雲の発達の度合いと相関のある値に対応するものがあればそれを利用すればよい。

また、本実施の形態では、概して「発雷」予測と称したが、訓練データへの正解ラベル付与において、「発雷」を雲放電（いわゆる雲間雷）や、対地放電（いわゆる対地雷／落雷）、または、正極雷、負極雷などに限定することで、それらに限定した予測も可能となる。

訓練データおよび予測時の入力データでは気象データをその要素としたが、実際はそのデータを観測した地域の特徴や、その時刻、日時または季節なども要素となりうる。

同様に、その時刻、その領域、および時間軸での前後、その周辺地域のデータに加えて、それらデータ間でのなんらかの演算結果（たとえば差分など）も要素となりうる。

このように本発明の実施の形態においては、運用を安価に抑えるために、安価に入手可能なデータのみを用いて予測を行っている。また、発雷の始まりを予測するために、雷現況情報に頼らずに気象データのみから予測を行っている。

また、高い精度で予測するために、時間精度、位置精度の高いデータを中心に入力として使用している。

また、分類器の精度を向上させるために、訓練データの前処理において、積乱雲の発達の度合いという発雷に直接的なデータに沿ってエントリ数を調整している。

また、汎化性能を上げるためにＳＶＭを採用し、そしてリードタイムをかせぐために、パターン変化予測を合わせて使用している。

以上説明した本発明の実施の形態によれば、機械学習後の予測において雷データを必要としないので、高価な雷センサの設置が必要がなく、また、雷現況情報の受信契約をする必要がない。

また、機械学習後の予測において雷データを参考にしないので、観測した発雷を頼りに次を予測する手法と違い、その雷雲の最初の発雷を予測できる可能性がある。

また、一般的に安価に入手可能な気象データのみを入力としているので、独自のレーダなどを設置する必要がない。

また、発雷予測に用いるデータの代表的なものとして気象レーダによるエコー値があるが、このデータは現在１０分間隔で配信されており、位置精度も１ｋｍ単位なので、雷も同様の精度で直前に正確に予測することができる。

また、訓練データのバランスをとるための手段を提供しているので、機械学習時に正例と負例をバランスよく学習することができ、結果として精度が向上する。

また、訓練データのバランスをとる際に、積乱雲の発達の度合いという発雷に直接的なデータに沿ってエントリ数を調整するので、学習に重要な判断の際どいエントリを重点的に学習することができ、結果として分類器の精度が向上する。

本発明の第１の実施の形態を説明するための全体の処理フローを示す。本発明の第１の実施の形態を説明するための学習フローを示す。本発明の第１の実施の形態を説明するための予測フローを示す。本実施の第１の形態の処理フローにおける、各データと各処理との関連を説明するための説明図を示す。気象データの一例を示す。訓練データの一例を示す。分類器を構成する超平面の概念図を示す。メッシュ領域を説明するための説明図を示す。本発明の第２の実施の形態を説明するための全体の処理フローを示す。

符号の説明

１…過去気象データＤＢ
２…過去の雷データ
３…過去の気象データ
４…訓練データ
５…学習結果
６…気象データ
７…分類器
８…予測結果

Claims

メッシュ状に区切られた複数の領域ごとに発雷を直前に予測するための発雷予測装置において、
前記領域のそれぞれごとに過去の時系列に沿って気象状況を記録した気象データに基づいて生成されたエントリと、
前記エントリの総数において発雷の有無の割合が所定の割合となるように発雷の無かったエントリを間引いて生成した訓練データと、
前記訓練データを用いてコンピュータにより機械学習を実行して得た学習結果に基づいて構成され、現在の気象データに対して発雷を予測するための分類器と、
を備え、
現在の気象データの一部もしくは全部にパターン変化予測を適用し、その予測結果を前記分類器へ入力する
ことを特徴とする発雷予測装置。
前記訓練データにおいて、
積乱雲の発達の度合いと相関のある値を参照し、発達の度合いが低い時は前記間引きを多くし、発達の度合いが高いときは前記間引きを少なくすること
を特徴とする請求項１に記載の発雷予測装置。
前記機械学習の手法としてＳＶＭ（Support Vector Machine）を用いること
を特徴とする請求項１または２に記載の発雷予測装置。
メッシュ状に区切られた複数の領域ごとに発雷を直前に予測するための発雷予測方法において、
前記領域のそれぞれごとに過去の時系列に沿って気象状況を記録した気象データに基づいたエントリを生成するステップと、
前記エントリの総数において発雷の有無の割合が所定の割合となるように発雷の無かったエントリを間引いて訓練データを生成するステップと、
前記訓練データを用いてコンピュータにより機械学習を実行して得た学習結果に基づいて、現在の気象データに対して発雷を予測するための分類器を構成するステップと、
を有し、
現在の気象データの一部もしくは全部にパターン変化予測を適用し、その予測結果を前記分類器へ入力する
ことを特徴とする発雷予測方法。
前記訓練データにおいて、
積乱雲の発達の度合いと相関のある値を参照し、発達の度合いが低い時は前記間引きを多くし、発達の度合いが高いときは前記間引きを少なくすること
を特徴とする請求項４に記載の発雷予測方法。
前記機械学習の手法としてＳＶＭ（Support Vector Machine）を用いること
を特徴とする請求項４または５に記載の発雷予測方法。
メッシュ状に区切られた複数の領域ごとに発雷を直前にコンピュータを用いて予測させるための発雷予測プログラムであって、
前記領域のそれぞれごとに過去の時系列に沿って気象状況を記録した気象データに基づいたエントリを生成するステップと、
前記エントリの総数において発雷の有無の割合が所定の割合となるように発雷の無かったエントリを間引いて訓練データを生成するステップと、
前記訓練データを用いてコンピュータにより機械学習を実行して得た学習結果に基づいて、現在の気象データに対して発雷を予測するための分類器を構成するステップと、
をコンピュータに実行させ、
現在の気象データの一部もしくは全部にパターン変化予測を適用し、その予測結果を前記分類器へ入力する
ことを特徴とする発雷予測プログラム。