JPH0972965A - 気象擾乱の予測方法 - Google Patents

Abstract

(57)【要約】 【目的】 気象擾乱のセルの発達フェーズの違いを考慮
して、レーダー雨量計による計測値の短時間先予測、及
び、予測結果を表示する方法を提供することにある。 【構成】 降雨強度などの第１の物理量であるｗを計測
するとともに、ｗ以外の雲頂高度などの第２の物理量で
あるｚを計測し、計測されたｚから移動ベクトルを算出
し（１１）、算出された移動ベクトルとｚからから気象
擾乱の発達の特性を表す発達特性曲線を同定する（１
２）とともに、ｚから気象擾乱の相関関係を表す相関曲
線を同定し（１３）、さらに移動ベクトルを予測し（１
４）、移動ベクトル、発達特性曲線、相関関係曲線およ
びｗからｗの所定時間先の状態を予測する（２１）。

Description

【発明の詳細な説明】

【０００１】

【産業上の利用分野】本発明は、気象擾乱の所定時間先
の状態を予測する気象擾乱の予測方法に係る。より具体
的に説明すれば、気象擾乱の発達・移動の短時間先予測
方法に関し、特に、空間スケール数ｋｍ〜数百ｋｍ、生
成から消滅までの寿命がおよそ１時間〜１日程度の気象
擾乱の状態の予測方法に関する。

【０００２】また、予測した状態の表示方法にも関す
る。

【０００３】

【従来の技術】都市部では、市街地の整備に伴い、舗装
率が向上し、雨水が地下浸透せずに流出することによる
洪水、浸水の危険性が高くなっている。この危険を回避
するために、管渠網、排水ポンプ、貯溜池、浸透ますな
どを整備し、これら施設の適切な運転、運用による雨水
排水、雨水貯溜、雨水浸透などを行なっている。この適
切な運転、運用のためには、豪雨をもたらす気象擾乱の
短時間先予測をする必要がある。

【０００４】このような気象擾乱の短時間先予測方法と
して、従来、レーダー雨量計によりレーダー反射因子分
布（または反射因子を降雨強度に換算したものの分布）
を測定、算出し、この分布の動きを追跡する方法などが
用いられている。

【０００５】例えば、下水道協会誌論文集、vol.28、N
o.330、P1-P13、「雨水排水システムのための短時間降
雨量予測法の構築」には、（１）雨域の移動は雨域の重
心点の移動として外挿し、（２）雨域の発達は雨域の平
均の変化として外挿する方法がある。

【０００６】また、第２８回水理講演会論文集、P423-P
428、「移流モデルによる短時間降雨予測手法の検討」
には、簡単な移流モデルを導入し、計測値からこのモデ
ルパラメータを同定し、このパラメータを使って、移流
モデルにより雨域の発達・移動を予測する方法がある。

【０００７】

【発明が解決しようとする課題】一般に、雨域、もしく
は、雨域をもたらす雲は、いくつかのセルに分けて考え
ることができるが、これらセルは発達中のものもあれば
減衰中のものもある。また、その発達、減衰速度も急激
なときもあれば緩やかなときもある。このような複雑な
様相を呈する気象擾乱に対して、過去のトレンドから移
動分や発達分を単純に外挿するだけでは、短時間先予測
には不十分である。特に、発達分は移動分と比較して短
時間での影響が大きいため、発達の物理過程を十分考慮
する必要がある。

【０００８】雲の発達の物理過程の特徴をよく表す物理
量として、雲頂高度（または温度）、雲量、電位差、上
昇気流などがあり、これら物理量を測定する計測器（セ
ンサ）として、気象衛星赤外画像、気象衛星可視画像、
雷探知器、ドップラーレーダー、風速計などがある。

【０００９】本発明の目的は、上記のような発達の物理
過程の特徴をよく表す物理量の計測値とレーダー雨量計
による計測値とを統合することにより、気象擾乱のセル
に対する発達過程をとらえ、レーダー雨量計による計測
値の短時間先予測を行なう方法を提供することにある。
また、予測された結果を表示する方法を提供することに
ある。

【００１０】

【課題を解決するための手段】上記の短時間先予測の目
的を達成するために、本発明は、以下の構成を採用し
た。

【００１１】今、雲頂高度（または温度）、雲量、電位
差、上昇気流の大きさなど、発達過程の特徴を表す物理
量をｚと表記する。また、予測が行なわれるレーダー反
射因子、または、降雨強度をｗと表記する。ｚ、ｗは空
間的な分布として計測される。

【００１２】（１）ある方法で算出した移動ベクトル
（ｄｘ、ｄｙ）に基づいて、２つの時刻で得られたｚ、
ｚ’を対応する点で比較照合し、（ｚ，ｚ’）のペアを
求め、横軸をｚ’−ｚ，縦軸をｚとした平面にプロット
し、この平面にプロットされた点の分布を近似する１つ
の発達特性曲線を求める処理 （２）ほぼ同時刻に得られたｚとｗの対応する点を比較
照合し、（ｚ，ｗ）のペアを求め、横軸をｗ、縦軸をｚ
とした平面にプロットし、この平面にプロットされた点
の分布を近似する１つの相関関係曲線を求める処理 （３）ある地点（ｘｉ，ｙｉ）を選択し、この地点にお
けるｗｉを求め、ｗｉに対応するｚｉを相関関係曲線よ
り求め、このｚｉを基点ｚｉ（０）として、発達特性曲
線に従って、ｚｉの短時間先予測値ｚｉ（ｔ）を求め、
ｚｉ（ｔ）に対応するｗｉ（ｔ）を相関関係曲線より求
め、ある方法で算出した予測移動ベクトル（ｄｘ
（ｔ）、ｄｙ（ｔ））とｗｉ（ｔ）とに基づいて、地点
（ｘｉ，ｙｉ）の予測移動先でのｗｉの予測値を求める
処理を有効な範囲の地点全てに対して行なう処理 また、上記の予測結果の表示の目的を達成するために、
以下の３つの処理のいずれかを有することを特徴とす
る。

【００１３】上記（１）（２）の処理によって短時間先
予測を行なったとき、発達特性曲線に基づいて、ｗの予
測値を発達期、成熟期、減衰期などの発達フェーズのい
ずれかに分類した上で、 （４）隣接するｗの予測値の発達フェーズが異なると
き、この隣接するｗの境界に線を入れる処理 （５）発達フェーズごとに、特定の記号を定めておき、
この記号をｗの予測結果の表示に重ねて表示する処理 （６）発達フェーズごとに、色の系統を定めておき、こ
の色の系統に基づいて、ｗの予測結果を表示する処理

【００１４】

【作用】上記処理（１）の新しさは、２つの時刻から得
られた情報だけから雲の発達の一連の過程を同定するこ
とにある。これは、同じ時刻の雲でも、それぞれのセル
の発達フェーズはばらばらであることを利用している。
もし、従来のように雲全体を１つのかたまりとみなし
て、一連の発達の過程を捉えようとすれば、雲が生成を
始めてから雲全体が消滅するまで追跡しなければならな
くなる。

【００１５】上記処理（１）によれば、発達期にあるも
のは、平面の右側に分布し、発達速度の速いものほど、
より右側になる。逆に、減衰期にあるものは、平面の左
側に分布し、減衰速度の速いものほど、より左側にな
る。成熟期にあるものは平面の中央（縦軸）付近に分布
する。これら、分布を１つの曲線で近似すると、それが
発達特性曲線である。この発達特性曲線に対する分布の
ばらつきが小さいほど、対象としている気象擾乱の発達
過程が一様であると言え、発達の予測精度が高くなる。
この発達特性曲線を求めておけば、任意の値から発達、
成熟、減衰の過程をたどりながら予測できる。このこと
は、より長い時間先の予測に効果が現れる。

【００１６】上記処理（２）の新しさも、処理（１）の
ものと同様で、１つの時刻の比較で雲の発達の一連の過
程に対応する降水量を同定することにある。

【００１７】上記処理（２）によれば、気象擾乱ごとに
ｚとｗの相関関係がわかる。求めた相関関係曲線に対す
る分布のばらつきが小さいほど、ｚとｗの相関関係は強
いと言え、変換の精度が高くなり、発達に特徴的な物理
量ｚを介して、ｗの発達を予測できる。

【００１８】上記処理（３）によれば、移動ベクトル、
相関関係曲線、発達特性曲線を介して、任意の地点の短
時間先の移動地点、及び、その移動地点におけるｗの予
測値を算出できる。有効な範囲の地点全体に対してこの
処理を行なえば、ｗの短時間先予測値の空間分布が求ま
る。

【００１９】一方、ｗの予測値をメッシュごとに、色の
階調、または、数値で表示したとして、上記処理（４）
によれば、発達期にあるメッシュ、減衰期にあるメッシ
ュなどを取り囲むようにに線が表示される。この表示に
より気象擾乱において発達フェーズが異なる境界がわか
る。

【００２０】上記処理（５）によれば、例えば、発達期
にあるメッシュには○、減衰期にあるメッシュにはＸと
いうような記号をつければ、そのメッシュがどの発達フ
ェーズにあるのかがわかる。なお、処理（４）と処理
（５）を合わせて、境界線に囲まれている領域に１つだ
け記号を付すこともできる。

【００２１】上記処理（６）によれば、例えば、発達期
にあるメッシュには赤系統、減衰期にあるメッシュには
青系統というように色の系統を定めておけば、そのメッ
シュがどの発達フェーズにあるのかがわかる。この場合
には、処理（４）を利用しなくても発達フェーズが異な
る境界がわかる。

【００２２】

【実施例】本発明の説明においては、発達の物理過程の
特徴を表す物理量ｚとして、気象衛星の赤外画像による
雲頂高度を用いる。ｚとしては、その他に、雲頂温度、
雲量、電位差、気流などを採用してもよい。

【００２３】図１は、本発明の気象擾乱の発達・移動の
把握による短時間先予測方法全体の処理とデータの流れ
を示したものである。

【００２４】処理は、ｚの空間分布を計測する時間間隔
ｔｚごとに行われるブロック１０１と、レーダー雨量計
による計測値ｗの空間分布を計測する時間間隔ｔｗごと
に行われるブロック１０２とに分かれる。ここで、計測
値ｗは、レーダー反射因子Ｄでも、それを次式のＢ−β
変換により降雨強度Ｒに換算したものでもよい。

【００２５】

【数１】

【００２６】ここで、Ｂ，βはパラメータである。以下
ではｗをレーダー反射因子として説明する。また、時間
間隔ｔｚとｔｗとは一致している必要はなく、例えば、
気象衛星赤外画像ならば３０分、レーダー雨量計ならば
５分が標準的である。さらに、それぞれの空間分解能も
一致している必要はない。

【００２７】はじめに、ブロック１０１について説明す
る。ブロック１０１は、移動ベクトルを算出する処理１
１と、ｚの発達過程を同定する処理１２と、ｚとｗとの
相関を同定する処理１３と、次の移動ベクトルを予測す
る処理１４とからなる。

【００２８】移動ベクトル（ｄｘ，ｄｙ）を算出する処
理１１は、従来より用いられている相互相関係数法や重
心法などの方法をそのまま適用してもよく、また、上層
計測による風向・風速を移動ベクトルとしてもよい。こ
こでは、赤外画像のしきい値処理した画像によるパター
ンマッチング処理による移動ベクトルの算出方法につい
て説明する。

【００２９】図２は、移動ベクトルの一算出方法の説明
図である。処理２１１は、連続する２時刻で得られた気
象衛星による赤外画像２０１、２０２の雲頂高度をしき
い値処理し、２値画像２０３、２０４に変換する。図で
は０、１に２値化している。処理２１２は、２０３の画
像から１に符号化されたある領域２０５を完全に含む閉
領域２０６を取り出し、２０４の画像から、閉領域２０
６と全く同じ大きさの閉領域２０７を任意に選択する。
処理２１３は、２つの閉領域を比較照合し、０と０が対
応している点、及び、１と１が対応している点の数をカ
ウントする。処理２１２，２１３をくり返し、処理２１
４は、カウントが最大になるように選んだ閉領域２０７
を、閉領域２０６の移動とみなし、移動ベクトル（ｄ
ｘ，ｄｙ）２０８を算出する。この方法は、相互相関係
数法による比較照合を２値化により簡素化したものとみ
なせる。なお、このままでは閉領域２０７を選択する候
補が多く、処理に時間がかかるので、１時刻前（ここで
はｔｚ前）に算出した移動ベクトルを参照するなどし
て、選択する候補をある範囲にしぼり、全体の処理時間
を短くするのが望ましい。

【００３０】次に、ｚの発達過程を同定する処理１２を
説明する。図３は本発明の物理量の発達過程を同定する
方法の説明図である。

【００３１】移動ベクトルを算出する処理１１によって
算出された移動ベクトル（ｄｘ，ｄｙ）２０８を用い
て、処理３１１は、連続する２時刻で得られた気象衛星
による赤外画像２０１、２０２の対応する点を比較照合
する。ここで、赤外画像２０１の雲頂高度をｚ、赤外画
像２０２の雲頂高度をｚ’と表記する。例えば、赤外画
像２０１の点（ｘ１、ｙ１）の雲頂高度がｚ１ならば、
比較照合されるのは、赤外画像２０２の点（ｘ１＋ｄ
ｘ、ｙ１＋ｄｙ）の雲頂高度ｚ１’である。このように
して、いくつかのｚとｚ’のペアが求まる。処理３１２
は、横軸にｚ’−ｚ、縦軸にｚをとり、この平面に比較
照合したｚとｚ’のペア３０１をプロットし、これらの
ペアを最もよく近似する１つの発達特性曲線３０２を求
める。発達特性曲線３０２を求めるには、最小２乗法な
どを活用すればよい。この方法の新しさは、２つの時刻
から得られた情報だけから雲の発達の一連の過程を同定
することにある。これは、同じ時刻の雲でも、それぞれ
のセルの発達フェーズはばらばらであることを利用して
いる。もし、従来のように雲全体を１つのかたまりとみ
なして、一連の発達の過程を捉えようとすれば、雲が生
成を始めてから雲全体が消滅するまで追跡しなければな
らなくなる。

【００３２】次に、ｚとｗとの相関を同定する処理１３
を説明する。図４は本発明の２つの物理量の相関関係を
同定する方法の説明図である。

【００３３】処理４１１は、赤外画像２０２が得られた
時刻とほぼ同じ時刻のレーダー雨量計によるレーダー反
射因子ｗ４００と、赤外画像２０２の雲頂高度ｚのの対
応する点を比較照合する。ここで、ｚとｗとの空間分解
能が異なる場合もある。例えば、ｚが気象衛星の赤外画
像は５ｋｍメッシュ、レーダー雨量計は５００ｍメッシ
ュのように異なることがある。このような場合には、１
つの点のｚに対して、比較対象となる複数のｗを対応さ
せてもよいし、比較対象となるｗの平均値を対応させて
もよい、例えば、赤外画像２０２の点（ｘ１、ｙ１）の
雲頂高度（または温度）がｚ１ならば、レーダー雨量計
の点（ｘ１、ｙ１）のいくつかの近傍の点のレーダー反
射因子ｗ１１、ｗ１２、ｗ１３．．．、または、これら
の平均値を対応させる。このようにして、いくつかのｚ
とｗのペアが求まる。処理４１２は、横軸にｗ、縦軸に
ｚをとり、この平面に比較照合したｗとｚのペア４０１
をプロットし、これらのペアを最もよく近似する１つの
相関関係曲線４０２を求める。相関関係曲線４０２を求
めるには、最小２乗法などを活用すればよい。

【００３４】次に、次の移動ベクトルを予測する処理１
４を説明する。移動ベクトルを算出する処理１１が算出
した最新の移動ベクトル（ｄｘ、ｄｙ）２０８を、過去
に処理１１が算出した移動ベクトルを用いて補正して、
次の移動ベクトル（ｄｘ’、ｄｙ’）５０１を予測す
る。一般的には、移動ベクトルの大きさの変化分と移動
ベクトルの方向の変化分のトレンドに基づいて、それぞ
れの次の変化分を推定し、これに基づいて最新の移動ベ
クトル２０８を補正する。

【００３５】次に、ブロック１０２について説明する。
ブロック１０２は、ｗの短時間先予測する処理２１から
なる。

【００３６】図５、６は、本発明のレーダー反射因子、
または、降雨強度の短時間先予測方法の第１の説明図お
よび第２の説明図である。

【００３７】処理６１１は、地点（ｘｉ，ｙｉ）を選択
し、この地点におけるレーダー反射因子（または降雨強
度）ｗｉを求める。処理６１２は、ｗｉ７０１に対応す
る相関関係曲線４０２上のｚｉを求める。処理６１３
は、ｚｉ７０２に対応する発達特性曲線３０２上の点ｚ
ｉ（０）を求める。なお、ｚｉ、ｚｉ（０）が相関関係
曲線４０２、発達特性曲線３０２上のどこになるのか
は、地点（ｘｉ，ｙｉ）のｚ，または、１時刻前（ここ
ではｔｗ前）に処理６１７で求めたｔｗ後の地点（ｘ
ｉ，ｙｉ）のｚの予測値を蓄積したデータベース６１８
を参照して求める。

【００３８】処理６１４は、発達特性曲線３０２に基づ
いて、ｚｉ（０）を基点として、ｚｉ（５），ｚｉ（１
０），．．．．のように雲頂高度の短時間先予測を行な
う。処理６１５は、雲頂高度の予測値のそれぞれに対応
する相関関係曲線４０２上のｗｉ（５），ｗｉ（１
０），．．．．のようにレーダー反射因子（または降雨
強度）の短時間先予測値を求める。処理６１６は、処理
１４で予測した次の移動ベクトル（ｄｘ’、ｄｙ’）５
０１と、レーダー反射因子（または降雨強度）の短時間
先予測値とに基づいて、対応する移動先の短時間先予測
値を求める。つまり、地点（ｘｉ，ｙｉ）のレーダー反
射因子ｗｉは、以下のように発達・移動すると予測され
る。

【００３９】５分後、地点（ｘｉ＋ｄｘ’（５），ｙｉ
＋ｄｙ’（５）） レーダー反射因子ｗｉ（５） １０分後、地点（ｘｉ＋ｄｘ’（１０），ｙｉ＋ｄｙ’
（１０）） レーダー反射因子ｗｉ（１０） 処理６１７は、処理１４で予測した次の移動ベクトル
（ｄｘ’、ｄｙ’）５０１と、処理６１４で求めた雲頂
高度の短時間先予測値とに基づいて、対応する移動先の
ｔｗ後の雲頂高度の予測値をデータベース６１８へセー
ブし、処理６１２、６１３のための参照データとする。

【００４０】処理６１１から処理６１７を有効な範囲の
地点に対して全て行なうと、レーダー反射因子の分布が
短時間先予測される。

【００４１】以上、図１に示した本発明の気象擾乱の発
達・移動の把握による短時間先予測方法の処理について
説明した。さらに、若干の補足説明をする。

【００４２】一般に、気象衛星による計測では、レーダ
ー雨量計にくらべて、より広範囲な分布を計測すること
ができる。また、降雨が降り始める前に雲が生成されて
いるのがほとんどであるので、降雨をもたらすような気
象擾乱をレーダー雨量計の計測よりも先に発見できる。

【００４３】したがって、レーダー雨量計の計測範囲外
の地点や、レーダー雨量計のデータが０の地点に対して
も、雲に関する物理量を基点ｚｉ（０）として、処理６
１４以降を行なうことで、その地点に対してのレーダー
反射因子の予測値を得ることができる。ただし、この場
合には、気象衛星の計測の空間分解能よりも細かい空間
分解能で予測することは困難である。

【００４４】また、時間間隔ｔｚとｔｗの違い、及び、
計測時刻のずれに伴うブロック１０１、ブロック１０２
の処理のタイミング、または、利用する情報の時刻につ
いて図１を用いて補足説明する。この説明に特に関係す
るのは、ブロック１０１の処理１３とブロック１０２の
処理２１である。まず、処理１３であるが、これは、相
関関係について算出するものなので、原則的には同時刻
で計測されたものを用いる必要があるが、ｔｗの時間間
隔がある程度小さければ、最新で得られた第１の物理量
と比較照合することにすれば、最大ずれていたとしても
ｔｗ分である。

【００４５】処理２１は、予測の初期情報から追跡する
だけなので、ブロック１０１の最新の結果を反映しさえ
すればよく、時刻のずれに対する問題点は生じない。

【００４６】次に、本発明の気象擾乱の短時間先予測結
果の表示方法について説明する。

【００４７】上記で説明した短時間先予測方法を用いれ
ば、図６に示したとおり、予測したものが、発達期７０
１、成熟期７０２、減衰期７０３などいくつかの発達フ
ェーズに分類することができる。このことを利用して、
通常のレーダー反射因子の分布の予測結果の表示に工夫
を施す。

【００４８】図７は、本発明のレーダー反射因子の予測
結果の表示方法を示したものである。

【００４９】通常、レーダー反射因子の分布の予測結果
の表示は、レーダー反射因子の値を量子化し、量子化し
た値にグレーレベル８０２、または、色を割り当て、表
示装置の画面８０１に、各レーダー反射因子の予測値を
対応する箇所に割り当てたグレーレベルを表示して、予
測分布を表示するのが一般的である。対応する箇所に数
値をそのまま表示することもある。また、予測移動ベク
トル５０１を表示することもある。このような表示への
追加表示として、分類した発達フェーズに基づいて、隣
接するメッシュにおいて、発達フェーズが異なるとき、
メッシュの境界に太線などの線を挿入すると、これらの
線をつなぎあわせたもの８０３は発達フェーズが同じ領
域を囲む。さらに、発達フェーズに記号８０４を割り当
て、領域上にこの記号を表示する。これによって、気象
擾乱のどの領域がどの発達フェーズにあるのかが容易に
わかる。

【００５０】なお、以上の説明においては、線と記号で
発達フェーズを識別表示する方法であるが、発達フェー
ズごとに色の系統、例えば、青系統、緑系統、赤系統な
どを定めておき、それぞれの系統に対して、量子化した
値に色のレベルを割り当てて発達フェーズを識別表示す
る方法も考えられる。

【００５１】

【発明の効果】本発明によれば、レーダー雨量計のデー
タだけでなく、より発達の物理過程に特徴的な物理量デ
ータを活用して、気象擾乱のセルごとの発達の予測が可
能となり、気象擾乱の各セルの発達のフェーズの違いを
取り扱って、レーダー雨量計のデータの短時間先予測を
することができる。

【００５２】また、本発明によれば、レーダー雨量計よ
りも広範囲な計測手段を用いれば、レーダー雨量計の計
測範囲外から移動してくる気象擾乱に対して予測できる
し、降雨が降り出す前に値を持つ物理量の計測手段を用
いれば、レーダー雨量計が計測値をもたない地点に対す
る予測ができ、降雨の降始めを捉えることができる。さ
らに、本発明によれば、分布の予測結果の表示におい
て、発達フェーズの異なる領域を識別することができ
る。

【図面の簡単な説明】

【図１】本発明の気象擾乱の発達・移動の把握による短
時間先予測方法の全体の処理フローを示す図である。

【図２】移動ベクトルの一算出方法の説明図である。

【図３】本発明の物理量の発達過程を同定する方法の説
明図である。

【図４】本発明の２つの物理量の相関関係を同定する方
法の説明図である。

【図５】本発明のレーダー反射因子、降雨強度の短時間
先の予測の仕方を説明する第１の説明図である。

【図６】本発明のレーダー反射因子、降雨強度の短時間
先の予測の仕方を説明する第２の説明図である。

【図７】本発明のレーダー反射因子の予測結果の表示方
法の説明図である。

【符号の説明】

５０１…予測移動ベクトル、３０２…発達特性曲線、４
０２…相関関係曲線

Claims (9)

【特許請求の範囲】
1. 【請求項１】気象擾乱の所定の状態値である第１の物理
   量の分布を計測し、 前記第１の物理量以外の気象擾乱の所定の状態値である
   第２の物理量の分布を計測し、 前記気象擾乱の移動ベクトルを前記第２の物理量を用い
   て算出し、 前記移動ベクトルに基づいて、２時刻で計測された第２
   の物理量の分布を比較照合し、第２の物理量の時間によ
   る変化の特性を示す発達特性曲線を同定する処理と、 所定時刻における前記第１の物理量の分布と前記第２の
   物理量の分布とを比較照合し、前記第１の物理量と前記
   第２の物理量の関係を示す相関関係曲線を同定する処理
   と、 前記移動ベクトル、前記発達特性曲線および前記相関関
   係曲線に基づいて、所定時間先の前記第１の物理量の状
   態を予測することにより、気象擾乱の状態を予測するこ
   とを特徴とする気象擾乱の予測方法。
2. 【請求項２】請求項１に記載の気象擾乱の予測方法にお
   いて、 前記発達特性曲線は、前記移動ベクトルを用いて、２つ
   の時刻で得られた第２の物理量の対応する点を比較照合
   し、第２の物理量のペアを求め、一方の軸を該ペアの差
   分、他方の軸を該ペアのいずれかとした平面にプロット
   し、該平面にプロットされた点の分布を近似する１つの
   発達特性曲線を求めることを特徴とする気象擾乱の予測
   方法。
3. 【請求項３】請求項１に記載の気象擾乱の予測方法にお
   いて、 前記相関関係曲線処理は、前記第１の物理量と前記第２
   の物理量の対応する点を比較照合し、一方の軸を前記第
   １の物理量、他方の軸を前記第２の物理量とした平面に
   プロットし、該平面にプロットされた点の分布を近似す
   る１つの相関関係曲線を求めることを特徴とする気象擾
   乱の予測方法。
4. 【請求項４】前記所定時間先の前記第１の物理量の状態
   は、所定の地点を選択し、該地点における第１の物理量
   を求め、該第１の物理量に対応する第２の物理量を上記
   相関関係曲線より求める処理と、 上記第２の物量量を基点として、上記発達特性曲線に従
   って、第２の物理量の短時間先予測値を求める処理と、 上記第２の物理量の短時間先予測値に対応する第１の物
   理量の短時間先予測値を上記相関関係曲線より求める処
   理と、 上記移動ベクトルに基づく予測移動ベクトルと、上記第
   １の物理量の短時間先予測値とに基づいて、先に選択し
   た地点の予測移動先での第１の物理量の所定時間先予測
   値を求める処理とを有し、 以上の処理を予め定めた複数の地点に対して行なうこと
   を特徴とする気象擾乱の予測方法。
5. 【請求項５】請求項１に記載の気象擾乱の予測方法にお
   いて、 気象衛星赤外画像の計測による雲頂高度、または、雲頂
   温度を第２の物理量の分布として用いることを特徴とす
   る気象擾乱の予測方法。
6. 【請求項６】請求項１に記載の気象擾乱の予測方法にお
   いて、 気象衛星可視画像の計測による雲量を第２の物理量の分
   布として用いることを特徴とする気象擾乱の予測方法。
7. 【請求項７】請求項１に記載の気象擾乱の予測方法にお
   いて、 雷探知器の計測による空間電位、または、空間電位差、
   または、電荷を第２の物理量の分布として用いることを
   特徴とする気象擾乱の発達・移動予測方法。
8. 【請求項８】請求項１に記載の気象擾乱の予測方法にお
   いて、 風速計、または、ドップラーレーダーの計測による上昇
   気流、または、下降気流の大きさを第２の物理量の分布
   として用いることを特徴とする気象擾乱の予測方法。
9. 【請求項９】気象擾乱の所定の状態値である第１の物理
   量の分布を計測し、 前記第１の物理量以外の気象擾乱の所定の状態値である
   第２の物理量の分布を計測し、 前記気象擾乱の移動ベクトルを前記第２の物理量を用い
   て算出し、 前記移動ベクトルに基づいて、２時刻で計測された第２
   の物理量の分布を比較照合し、第２の物理量の時間によ
   る変化の特性を示す発達特性曲線を同定する処理と、 所定時刻における前記第１の物理量の分布と前記第２の
   物理量の分布とを比較照合し、前記第１の物理量と前記
   第２の物理量の関係を示す相関関係曲線を同定する処理
   と、 前記移動ベクトル、前記発達特性曲線および前記相関関
   係曲線に基づいて、所定時間先の前記第１の物理量の状
   態を予測し、 前記発達特性曲線に基づいて、予測した第１の物理量を
   予め定めた発達フェーズの所定のものに分類し、 空間的に隣接する第１の物理量の発達フェーズが異なる
   とき、該隣接する第１の物理量の境界を示す線を表示す
   ることを特徴とする気象擾乱の表示方法。