JP2009052976A - Weather prediction data analyzer and weather prediction data analysis method - Google Patents

Weather prediction data analyzer and weather prediction data analysis method Download PDF

Info

Publication number
JP2009052976A
JP2009052976A JP2007218762A JP2007218762A JP2009052976A JP 2009052976 A JP2009052976 A JP 2009052976A JP 2007218762 A JP2007218762 A JP 2007218762A JP 2007218762 A JP2007218762 A JP 2007218762A JP 2009052976 A JP2009052976 A JP 2009052976A
Authority
JP
Japan
Prior art keywords
weather
data
determination
conditional probability
probability
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.)
Granted
Application number
JP2007218762A
Other languages
Japanese (ja)
Other versions
JP4908346B2 (en
Inventor
Fumihiko Mizutani
文彦 水谷
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Toshiba Corp
Original Assignee
Toshiba Corp
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Toshiba Corp filed Critical Toshiba Corp
Priority to JP2007218762A priority Critical patent/JP4908346B2/en
Publication of JP2009052976A publication Critical patent/JP2009052976A/en
Application granted granted Critical
Publication of JP4908346B2 publication Critical patent/JP4908346B2/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Abstract

<P>PROBLEM TO BE SOLVED: To provide a precise index for determining various influences caused by a weather phenomenon. <P>SOLUTION: A conditional probability table formation part 18 determines propriety of adaptability of a plurality of weather determination elements relative to a specific weather phenomenon respectively, and forms a conditional probability table showing a conditional probability in a dependence relation between each element based on a totalization result of the determination, concerning weather prediction data and weather observation data acquired during a fixed period. When new weather prediction data and weather observation data are acquired, a probability calculation part 17 determines each element on the acquired data, and calculates a generation probability of the specific weather phenomenon on reference to the conditional probability table corresponding to a weather determination element based on the determination result. When the calculated probability exceeds a threshold, a data determination/delivery part 19 informs a user of information showing the generation probability of the specific weather phenomenon. <P>COPYRIGHT: (C)2009,JPO&INPIT

Description

この発明は、気象予測システムから出力される予測データを解析する気象予測データ解析装置及び気象予測データ解析方法に関する。   The present invention relates to a weather prediction data analysis apparatus and a weather prediction data analysis method for analyzing prediction data output from a weather prediction system.

従来の気象予測システムでは、気象レーダ等で得られる観測データや気象庁から提供されるGPV(Grid Point Value)データ等を用いて大気の流れを計算することで気象予測を行っている。気象予測情報は、人々にとって身近であると同時に、台風や集中豪雨等のように生命や財産に関わる重要な情報であるため、予測データの信頼性の向上が図られている。例えば、観測データを取得する時間間隔が長い場合でも、予測データの精度を継続的に維持できるようにする手法が提案されている(例えば、特許文献1を参照。)。
特開2003−090888号公報
In a conventional weather prediction system, weather prediction is performed by calculating the atmospheric flow using observation data obtained by a weather radar or the like, or GPV (Grid Point Value) data provided by the Japan Meteorological Agency. Weather forecast information is not only familiar to people but also important information related to life and property, such as typhoons and torrential rains, so that the reliability of forecast data is improved. For example, a method has been proposed that allows the accuracy of prediction data to be continuously maintained even when the time interval for obtaining observation data is long (see, for example, Patent Document 1).
JP 2003-090888 A

ところが、気象予測情報は、あくまでも予測に過ぎず、時間的・空間的ずれを伴う不確実性を有する。しかし、予測情報の価値は存在する。例えば、集中豪雨の発生の可能性が事前に予測されていれば、たとえその可能性が低くとも十分に価値のある情報となる。つまり、気象予測情報の提供の仕方が価値を左右するのであり、提供手法にこそ問題点が存在する。   However, the weather prediction information is only a prediction, and has uncertainties accompanied by temporal and spatial deviations. However, the value of prediction information exists. For example, if the possibility of the occurrence of torrential rain is predicted in advance, the information is sufficiently valuable even if the possibility is low. In other words, the way of providing weather forecast information determines the value, and there is a problem in the providing method.

この発明は上記事情に着目してなされたもので、その目的とするところは、気象現象による様々な影響を判断する上で的確な指標を提供することが可能な気象予測データ解析装置及び気象予測データ解析方法を提供することにある。   The present invention has been made paying attention to the above circumstances, and its object is to provide a weather prediction data analysis apparatus and a weather prediction capable of providing an accurate index for judging various effects due to weather phenomena. It is to provide a data analysis method.

上記目的を達成するためにこの発明に係る気象予測データ解析装置は、気象観測データと、前記気象観測データをもとに気象予測モデルに基づいて演算された気象予測データとを取得する取得手段と、前記取得された気象予測データ及び気象観測データをもとに、特定の気象現象について複数の気象判定要素の適合の可否をそれぞれ判定する判定手段と、前記取得手段により一定期間に取得された気象予測データ及び気象観測データについて、前記判定手段による各気象判定要素の判定結果をそれぞれ集計する集計手段と、前記集計された結果をもとに前記複数の気象判定要素間それぞれについて要素間の依存関係における条件付確率を表す条件付確率表(CPT:Conditional Probability Table)を作成する作成手段と、前記取得手段により前記気象予測データ及び気象観測データが取得されたとき、前記判定手段による判定結果について前記気象判定要素に対応する条件付確率表を参照して、前記特定の気象現象の生起確率を算出する確率算出手段と、前記算出された生起確率が閾値を超えた場合に前記特定の気象現象の確率情報を通知する通知手段とを具備することを特徴とする。   In order to achieve the above object, a weather prediction data analysis apparatus according to the present invention includes weather observation data and acquisition means for acquiring weather prediction data calculated based on a weather prediction model based on the weather observation data; Determination means for determining whether or not a plurality of weather determination elements can be adapted to a specific weather phenomenon based on the acquired weather prediction data and weather observation data, and the weather acquired by the acquisition means for a certain period of time About the prediction data and the weather observation data, a totaling unit that totals the determination results of each weather determination element by the determination unit, and a dependency relationship between the plurality of weather determination elements based on the totaled results Creating means for creating a conditional probability table (CPT) that represents the conditional probability in the method; Probability calculation means for calculating the occurrence probability of the specific weather phenomenon with reference to a conditional probability table corresponding to the weather determination element for the determination result by the determination means when the prediction data and the weather observation data are acquired; And a notification means for notifying the probability information of the specific weather phenomenon when the calculated occurrence probability exceeds a threshold value.

上記構成による気象予測データ解析装置では、一定期間に取得された気象予測データ及び気象観測データについて、特定の気象現象について複数の気象判定要素の適合の可否をそれぞれ判定し、その判定の集計結果をもとに各要素間の依存関係における条件付確率を表す条件付確率表を作成する。そして、新たな気象予測データ及び気象観測データが取得されたとき、取得されたデータについて各要素の判定を行い、この判定結果をもとに気象判定要素に対応する条件付確率表を参照して、特定の気象現象の生起確率を算出している。このようにすることで、例えば、交通手段に影響を与える降雪等の特定の気象現象についての発生確率を通知することができるため、気象現象による様々な影響を判断をする上で、的確な指標を提供することが可能となる。   In the weather prediction data analysis apparatus having the above-described configuration, the weather prediction data and the weather observation data acquired for a certain period are respectively determined as to whether or not a plurality of weather determination elements can be applied to a specific weather phenomenon, and the result of the determination is calculated. A conditional probability table representing the conditional probability in the dependency relationship between each element is created. Then, when new weather forecast data and weather observation data are acquired, each element is determined for the acquired data, and the conditional probability table corresponding to the weather determination element is referred to based on the determination result. The probability of occurrence of a specific meteorological phenomenon is calculated. In this way, for example, it is possible to notify the probability of occurrence of a specific weather phenomenon such as snowfall that affects the means of transportation, so an accurate index for determining various effects due to the weather phenomenon Can be provided.

したがってこの発明によれば、気象現象による様々な影響を判断する上で的確な指標を提供することが可能な気象予測データ解析装置及び気象予測データ解析方法を提供することができる。   Therefore, according to the present invention, it is possible to provide a weather prediction data analysis apparatus and a weather prediction data analysis method capable of providing an accurate index for judging various effects due to weather phenomena.

以下、図面を参照しながら本発明の実施の形態を詳細に説明する。
図1は、本発明に係る気象予測データ解析装置の一実施形態を示す機能ブロック図である。この気象予測データ解析装置は、ネットワークNTを介して気象庁データサーバDS0,レーダサイトサーバDS1,DS1などに接続されている。
Hereinafter, embodiments of the present invention will be described in detail with reference to the drawings.
FIG. 1 is a functional block diagram showing an embodiment of a weather forecast data analyzing apparatus according to the present invention. The meteorological prediction data analysis apparatus is connected to the Japan Meteorological Agency data server DS0, radar site servers DS1, DS1, and the like via a network NT.

図1において気象予測データ解析装置は、ネットワークNTと接続される通信インターフェース11と、通信処理部12と、広域予測データ格納部13と、観測データ格納部14と、局地気象予測モデル演算部15と、局地予測データ格納部16とを備え、局地気象予測モデルによる気象予測が実施される。さらに、局地予測データ格納部16に格納された予測データから、顕著気象現象の生起確率を算出・提供するために、確率算出部17と条件付確率表作成部18とデータ判定/配信部19とを備える。   In FIG. 1, the weather prediction data analysis apparatus includes a communication interface 11 connected to the network NT, a communication processing unit 12, a wide area prediction data storage unit 13, an observation data storage unit 14, and a local weather prediction model calculation unit 15. And a local prediction data storage unit 16 for performing weather prediction using a local weather prediction model. Further, in order to calculate / provide the occurrence probability of the remarkable weather phenomenon from the prediction data stored in the local prediction data storage unit 16, the probability calculation unit 17, the conditional probability table creation unit 18, and the data determination / distribution unit 19 With.

気象予測のもとになる観測データ・予測データは、気象庁データサーバDS0やレーダサイトサーバDS1,DS2からネットワークNTを介して入力される。入力された観測データは観測データ格納部14に格納される。広域予測データ格納部13には、例えば、広域の気象予測データとして気象庁RSM(Regional Spectral Model)が格納される。観測データ格納部14及び広域予測データ格納部13に格納されたデータは、局地気象予測モデル演算部15からの要求に応じて選択的に局地気象予測モデル演算部15に送られる。   Observation data / prediction data as a basis for weather prediction is input from the Japan Meteorological Agency data server DS0 and radar site servers DS1, DS2 via the network NT. The input observation data is stored in the observation data storage unit 14. The wide-area prediction data storage unit 13 stores, for example, a Japan Meteorological Agency RSM (Regional Spectral Model) as wide-area weather prediction data. The data stored in the observation data storage unit 14 and the wide area prediction data storage unit 13 are selectively sent to the local weather prediction model calculation unit 15 in response to a request from the local weather prediction model calculation unit 15.

局地気象予測モデルとは、流体力学の方程式系を時間積分により解くことで、気象場である気圧・風・気温・雨量等を予測する方法である。現在、水平格子間隔数100m程度〜10km程度の範囲で予測可能な「非静力学気象モデル」が開発されている。非静力学モデルには、例えば、CReSS(名大地球水循環研究センターCloud Resolving Storm Simulator)のほかに、気象庁NHM(JMA NonHydrostatic Model)、MM5(PSU/NCAR mesoscale model)等の計算コードが存在する。それぞれ、発達した積乱雲(〜10km)程度の水平スケールでの物理過程を考慮したモデルである。   The local weather prediction model is a method for predicting atmospheric pressure, wind, temperature, rainfall, etc., which is a weather field, by solving a fluid dynamics equation system by time integration. Currently, “non-hydrostatic meteorological models” that can be predicted in the range of about 100 m to 10 km of horizontal grid spacing have been developed. For example, in addition to CReSS (Cloud Resolving Storm Simulator), there are calculation codes such as Japan Meteorological Agency NHM (JMA NonHydrostatic Model), MM5 (PSU / NCAR mesoscale model). Each model takes into account a physical process on a horizontal scale of developed cumulonimbus clouds (about 10 km).

局地気象モデルで計算される物理量は、現在一般に予測データとして配信される気圧・風・気温・雨量などのほかに、降雪量も算出可能である。   In addition to atmospheric pressure, wind, temperature, rainfall, etc., which are generally distributed as forecast data, physical quantities calculated using the local weather model can also calculate snowfall.

局地気象予測モデルは、解像度は高いが計算領域が広くなるに従い、大規模な計算機システムが必要となることから、通常より解像度は低いが、広い計算領域をもつ広域の気象予測データを初期値及び境界値として、限られた地域での気象予測を行う。   The local weather forecast model has a high resolution but requires a large-scale computer system as the calculation area widens.Therefore, although the resolution is lower than usual, a wide range of weather forecast data with a wide calculation area is the initial value. As a boundary value, weather forecast in a limited area is performed.

局地気象予測モデル演算部15は、広域の気象予測データの配信に伴って起動し、広域予測データ格納部13及び観測データ格納部14に格納されたデータをもとに気象予測演算を実行する。例えば、気象庁RSMの場合は12時間毎に配信される。局地気象予測モデル演算部15から出力される予測演算結果は局地予測データ格納部16に記憶される。   The local weather prediction model calculation unit 15 is activated along with the distribution of the wide area weather prediction data, and executes the weather prediction calculation based on the data stored in the wide area prediction data storage unit 13 and the observation data storage unit 14. . For example, in the case of the Japan Meteorological Agency RSM, it is delivered every 12 hours. The prediction calculation result output from the local weather prediction model calculation unit 15 is stored in the local prediction data storage unit 16.

また、観測データ格納部14には、配信間隔の短い例えば10分毎に配信される気象レーダデータなどが格納される。新たな観測データが入力されると、局地気象予測モデル演算部15は再び起動し、12時間毎に局地気象予測モデルにより予測された計算値をたとえば1時間毎に取り出し、その前1時間分の最新の気象レーダデータにより得られる雨量・風速データを「ナッジング法」と呼ばれるデータ同化手法を用いて局地気象予測モデルの再計算により補正し、数時間先の予測精度を向上させる。上記処理を繰り返し行うことによって予測演算を行い、結果である予測値を局地予測データ格納部16に格納しておき、例えば画面表示などの形態で局地予測データを提供する。   In addition, the observation data storage unit 14 stores weather radar data distributed every 10 minutes, for example, with a short distribution interval. When new observation data is input, the local weather prediction model calculation unit 15 starts again, takes out the calculated value predicted by the local weather prediction model every 12 hours, for example, every hour, and 1 hour before that The rain and wind speed data obtained from the latest meteorological radar data of the minute are corrected by recalculating the local weather forecast model using a data assimilation method called “nudge method” to improve the prediction accuracy several hours ahead. A prediction calculation is performed by repeating the above processing, and a prediction value as a result is stored in the local prediction data storage unit 16, and the local prediction data is provided in a form such as a screen display.

さらに、条件付確率表作成部18は、特定の気象現象について複数の気象判定要素の依存関係における条件付確率を表す条件付確率表(CPT:Conditional Probability Table)を作成する。確率算出部17は、広域予測データ格納部13及び観測データ格納部14に記憶された予測データ及び観測データをもとに条件付確率表を参照して、特定の気象現象の生起確率値を算出する。確率算出処理の詳細は後述する。算出された確率はデータ判定/配信部19によって、ある閾値以上の確率が算出された場合のみユーザに通知する。通知方法は、例えば携帯電話端末へのメール配信等が考えられる。   Further, the conditional probability table creation unit 18 creates a conditional probability table (CPT) that represents the conditional probability in the dependency relationship of a plurality of weather determination elements for a specific weather phenomenon. The probability calculation unit 17 refers to the conditional probability table based on the prediction data and the observation data stored in the wide area prediction data storage unit 13 and the observation data storage unit 14, and calculates the occurrence probability value of a specific weather phenomenon. To do. Details of the probability calculation process will be described later. The calculated probability is notified to the user only when the data determination / distribution unit 19 calculates a probability equal to or higher than a certain threshold. As a notification method, for example, mail delivery to a mobile phone terminal can be considered.

次に、この発明の特徴となる確率算出処理について説明する。
確率算出処理には、ベイジアンネットワーク(Bayesian Network)の手法を用いる。ベイジアンネットワークとは、不確かな出来事の連鎖について、確率の相互作用を集計する手法であり、「不確実性を扱うための計算モデル」として、認知科学やデータ・マイニング等に応用が広がっている。図2に、本発明で用いるベイジアンネットワークの例を示す。本実施形態では、特定の気象現象の一例として降雪の発生確率を求めるものとする。複数の気象判定要素をノード1〜5で表し、ノード間に張られたリンクが依存関係の方向を表している。条件付確率表はリンク毎に作成される。なお、図2で示したベイジアンネットワークはあくまで一例で、他のネットワーク構成も考えることができる。
Next, the probability calculation process that is a feature of the present invention will be described.
For the probability calculation process, a Bayesian Network method is used. A Bayesian network is a method of counting the interaction of probabilities for a chain of uncertain events, and its application is spreading to cognitive science, data mining, etc. as a “computation model for handling uncertainty”. FIG. 2 shows an example of a Bayesian network used in the present invention. In the present embodiment, it is assumed that the probability of snowfall is obtained as an example of a specific weather phenomenon. A plurality of weather determination elements are represented by nodes 1 to 5, and a link extending between the nodes represents a direction of dependency. A conditional probability table is created for each link. The Bayesian network shown in FIG. 2 is merely an example, and other network configurations can be considered.

(条件付確率表作成・更新処理)
まず、条件付確率表作成部18の処理について説明する。図3に、条件付確率表の作成処理の手順とその内容を示すフローをDFD(Data Flow Diagram)形式で表す。図5は、リンク毎判定集計結果の一例を示す図である。図6は、図5に示す集計結果をもとに作成された条件付確率表を示す図である。
(Conditional probability table creation / update process)
First, the processing of the conditional probability table creation unit 18 will be described. FIG. 3 shows a procedure for creating a conditional probability table and a flow showing the contents in a DFD (Data Flow Diagram) format. FIG. 5 is a diagram illustrating an example of a determination result for each link. FIG. 6 is a diagram showing a conditional probability table created based on the aggregation results shown in FIG.

図3において、条件付確率表作成部18は、気象観測データと、広域予測データと、局地気象予測モデル演算部15により演算された局地気象予測データとをもとに、ノード単体の気象判定要素について適合の可否を判定する(ステップS3a)。この判定により、例えば図4に示すようなノード判定結果テーブルが作成される。図4は、1時間毎に読み出した上記データをもとに、ノード1〜5それぞれの気象判定要素について適合の可否(True/False)を判定した結果である。   In FIG. 3, the conditional probability table creation unit 18 is based on the weather observation data, the wide-area prediction data, and the local weather prediction data calculated by the local weather prediction model calculation unit 15. It is determined whether or not the determination element is compatible (step S3a). By this determination, for example, a node determination result table as shown in FIG. 4 is created. FIG. 4 shows the result of determining the suitability (True / False) of the weather determination elements of the nodes 1 to 5 based on the data read out every hour.

次に、条件付確率表作成部18は、上記ノード判定結果テーブルをもとにリンク毎の判定結果を集計する(ステップS3b)。図5にリンク毎判定集計結果の一例を示す。図5は、ノード1とノード2との間のリンクの集計結果であり、a〜dは、ノード1とノード2の判定結果の組が一定期間に発生した頻度を表している。   Next, the conditional probability table creation unit 18 aggregates the determination results for each link based on the node determination result table (step S3b). FIG. 5 shows an example of determination results for each link. FIG. 5 shows a result of aggregation of the link between the node 1 and the node 2, and a to d represent the frequency at which a set of determination results of the node 1 and the node 2 occurs in a certain period.

条件付確率表作成部18は、上記集計された判定集計結果をもとにリンク毎に条件付確率表を作成する(ステップS3c)。図6は、図5に示した集計結果をもとに作成された条件付確率表を示す図である。ノード1とノード2の間の依存関係における条件付確率の値がそれぞれ算出される。   The conditional probability table creation unit 18 creates a conditional probability table for each link on the basis of the totaled determination totalization results (step S3c). FIG. 6 is a diagram showing a conditional probability table created based on the aggregation results shown in FIG. A conditional probability value in the dependency relation between the node 1 and the node 2 is calculated.

なお、条件付確率表作成部18は、広域予測データ格納部13、観測データ格納部14及び局地予測データ格納部16に日々蓄積されたデータを用いて、上記作成された条件付確率表の条件付確率の値を更新していく。最新データをもとに条件付確率表を更新することで、確率算出部17で算出される確率値の確からしさを向上させることが可能となる。   The conditional probability table creation unit 18 uses the data accumulated daily in the wide area prediction data storage unit 13, the observation data storage unit 14, and the local prediction data storage unit 16. The conditional probability value is updated. By updating the conditional probability table based on the latest data, the probability of the probability value calculated by the probability calculation unit 17 can be improved.

(確率算出処理)
確率算出部17は、このように作成された条件付確率表を用いて、特定の気象現象の生起確率を算出する。図7に、確率算出処理の作成処理の手順とその内容を示すフローをDFD(Data Flow Diagram)形式で表す。
(Probability calculation process)
The probability calculation unit 17 calculates the occurrence probability of a specific weather phenomenon using the conditional probability table created in this way. FIG. 7 shows a flow of the creation process of the probability calculation process and the contents thereof in a DFD (Data Flow Diagram) format.

確率算出部17は、広域予測データ格納部13、観測データ格納部14、又は局地予測データ格納部16に新たなデータが格納されると、このデータをもとにノード毎の判定を行う。この判定によるノード判定結果と条件付確率表作成部18により作成された条件付確率表とをもとに確率算出処理を行う(ステップS7a)。   When new data is stored in the wide area prediction data storage unit 13, the observation data storage unit 14, or the local prediction data storage unit 16, the probability calculation unit 17 performs determination for each node based on this data. Probability calculation processing is performed based on the node determination result by this determination and the conditional probability table created by the conditional probability table creating unit 18 (step S7a).

図8を参照しながら、ベイジアンネットワークを用いた確率算出手法について説明する。ここでは、ノード2の降雪が発生する確率を求めるものとする。確率算出部17は、広域予測データ格納部13、観測データ格納部14、又は局地予測データ格納部16に格納されたデータをもとにノード1,3,4,5の気象判定要素についての適合の可否をそれぞれ判定する。例えば、ノード1=True,ノード3=True,ノード4=False,ノード5=Falseと判定されたとする。このとき、各ノードの判定結果とそのノードに対応する条件付確率表とをもとに、ノード2がTrueである確率P2Tは、ベイズの定理によって図8に示す式のように算出できる。確率算出部17は、このようにして求められた確率値を出力する。 A probability calculation method using a Bayesian network will be described with reference to FIG. Here, it is assumed that the probability of occurrence of snowfall at node 2 is obtained. The probability calculation unit 17 uses the data stored in the wide area prediction data storage unit 13, the observation data storage unit 14, or the local prediction data storage unit 16 for the weather determination elements of the nodes 1, 3, 4, and 5. Judgment of suitability is made. For example, assume that node 1 = True, node 3 = True, node 4 = False, and node 5 = False. At this time, on the basis of the determination result of each node and the conditional probability table corresponding to that node, the probability P 2T that the node 2 is True can be calculated as shown in FIG. 8 by Bayes' theorem. The probability calculation unit 17 outputs the probability value obtained in this way.

以上述べたように上記実施形態では、条件付確率表作成部18は、一定期間に取得された気象予測データ及び気象観測データについて、特定の気象現象について複数の気象判定要素の適合の可否をそれぞれ判定し、その判定の集計結果をもとに各要素間の依存関係における条件付確率を表す条件付確率表を作成する。確率算出部17は、新たな気象予測データ及び気象観測データが取得されたとき、取得されたデータについて各要素の判定を行い、この判定結果をもとに気象判定要素に対応する条件付確率表を参照して、特定の気象現象の生起確率を算出している。このようにすることで、例えば、交通手段に影響を与える降雪等の特定の気象現象についての発生確率を通知することができるため、気象現象による様々な影響を判断をする上で、的確な指標を提供することが可能となる。さらに、作成された条件付確率表を最新の観測データ及び予測データで更新するようにすることで、より的中率の高い確率値を算出することが可能となる。   As described above, in the above-described embodiment, the conditional probability table creation unit 18 determines whether or not a plurality of weather determination elements can be applied to a specific weather phenomenon with respect to the weather prediction data and the weather observation data acquired in a certain period. A conditional probability table representing the conditional probabilities in the dependency relationship between each element is created based on the result of the determination. When new weather prediction data and weather observation data are acquired, the probability calculation unit 17 determines each element for the acquired data, and a conditional probability table corresponding to the weather determination element based on the determination result. The occurrence probability of a specific meteorological phenomenon is calculated with reference to FIG. In this way, for example, it is possible to notify the probability of occurrence of a specific weather phenomenon such as snowfall that affects the means of transportation, so an accurate index for determining various effects due to the weather phenomenon Can be provided. Furthermore, by updating the created conditional probability table with the latest observation data and prediction data, it is possible to calculate a probability value with a higher probability.

したがって、上記実施形態によれば、予測データの不確実性をできる限り排して、気象予測データから得られる情報を的確に提供することができる。この情報の提供により、住民の避難誘導や、ダム放水の運用・都市下水道のポンプ場運転など、気象現象の危険度判断を必要としているユーザに対し、判断を支援する価値の高い情報を提供することができる。   Therefore, according to the above-described embodiment, it is possible to accurately provide information obtained from weather prediction data by eliminating uncertainty of the prediction data as much as possible. By providing this information, high-value information that supports judgment is provided to users who need to judge the risk of weather phenomena, such as evacuation guidance for residents, dam drainage operation, and operation of urban sewer pump stations. be able to.

なお、この発明は上記実施形態に限定されるものではない。上記実施形態では、一例として豪雨となる確率値の算出法を記述したが、この他にも、例えば、土砂災害等をもたらす積算降雨量等の確率値の算出にも適用することができる。上述した(1)〜(5)の判定は、数値データを用いてそれぞれの閾値と比較して行っている。この閾値の値を変更することで、「大雪」や「みぞれ」等の確率値を簡単に算出することが可能である。   The present invention is not limited to the above embodiment. In the above-described embodiment, the calculation method of the probability value of heavy rain is described as an example, but it can also be applied to the calculation of the probability value of the integrated rainfall amount that causes a sediment disaster, for example. The above-described determinations (1) to (5) are made by using numerical data and comparing with each threshold value. By changing the threshold value, it is possible to easily calculate probability values such as “snowfall” and “sleet”.

また、上記示したベイジアンネットワークはあくまで一例であり、他のベイジアンネットワークの構築により「強風」「大雨」等の気象防災に関わる確率値を算出することももちろん可能である。但し、予測したい事象に対して、その裏づけとなる事象をノードに含める点は同様とする。   The Bayesian network shown above is only an example, and it is of course possible to calculate probability values related to weather disaster prevention such as “strong wind” and “heavy rain” by constructing other Bayesian networks. However, it is the same in that an event that supports the event to be predicted is included in the node.

すなわち、この発明は上記実施形態そのままに限定されるものではなく、実施段階ではその要旨を逸脱しない範囲で構成要素を変形して具体化できる。また、上記実施形態に開示されている複数の構成要素の適宜な組み合せにより種々の発明を形成できる。例えば、実施形態に示される全構成要素から幾つかの構成要素を削除してもよい。さらに、異なる実施形態に亘る構成要素を適宜組み合せてもよい。   That is, the present invention is not limited to the above-described embodiment as it is, and can be embodied by modifying the constituent elements without departing from the scope of the invention in the implementation stage. Further, various inventions can be formed by appropriately combining a plurality of constituent elements disclosed in the embodiment. For example, some components may be deleted from all the components shown in the embodiment. Furthermore, you may combine suitably the component covering different embodiment.

この発明に係る気象予測データ解析装置の一実施形態を示す機能ブロック図。The functional block diagram which shows one Embodiment of the weather forecast data analysis apparatus which concerns on this invention. 図1に示す気象予測データ解析装置において降雪の発生確率を算出するためのベイジアンネットワークの一例を示す図。The figure which shows an example of the Bayesian network for calculating the occurrence probability of snowfall in the weather prediction data analysis apparatus shown in FIG. 条件付確率表の作成処理の手順とその内容を示す図。The figure which shows the procedure of the creation process of a conditional probability table | surface, and its content. ノード判定結果テーブルの一例を示す図。The figure which shows an example of a node determination result table. リンク毎判定集計結果の一例を示す図。The figure which shows an example of the determination total result for every link. 条件付確率表の一例を示す図。The figure which shows an example of a conditional probability table | surface. 確率算出処理の手順とその内容を示す図。The figure which shows the procedure of the probability calculation process, and its content. ベイジアンネットワークを用いた確率算出手法を説明する図。The figure explaining the probability calculation method using a Bayesian network.

符号の説明Explanation of symbols

11…通信インターフェース、12…通信処理部、13…広域予測データ格納部、14…観測データ格納部、15…局地気象予測モデル演算部、16…局地予測データ格納部、17…確率算出部、18…条件付確率表作成部、19…データ判定/配信部、NT…ネットワーク、DS0…気象庁データサーバ、DS1,DS2…レーダサイトデータサーバ。 DESCRIPTION OF SYMBOLS 11 ... Communication interface, 12 ... Communication processing part, 13 ... Wide area prediction data storage part, 14 ... Observation data storage part, 15 ... Local weather prediction model calculating part, 16 ... Local prediction data storage part, 17 ... Probability calculation part , 18 ... Conditional probability table creation unit, 19 ... Data determination / distribution unit, NT ... Network, DS0 ... Japan Meteorological Agency data server, DS1, DS2 ... Radar site data server.

Claims (6)

気象観測データと、前記気象観測データをもとに気象予測モデルに基づいて演算された気象予測データとを取得する取得手段と、
前記取得された気象予測データ及び気象観測データをもとに、特定の気象現象について複数の気象判定要素の適合の可否をそれぞれ判定する判定手段と、
前記取得手段により一定期間に取得された気象予測データ及び気象観測データについて、前記判定手段による各気象判定要素の判定結果をそれぞれ集計する集計手段と、
前記集計された結果をもとに前記複数の気象判定要素間それぞれについて要素間の依存関係における条件付確率を表す条件付確率表(CPT:Conditional Probability Table)を作成する作成手段と、
前記取得手段により前記気象予測データ及び気象観測データが取得されたとき、前記判定手段による判定結果について前記気象判定要素に対応する条件付確率表を参照して、前記特定の気象現象の生起確率を算出する確率算出手段と、
前記算出された生起確率が閾値を超えた場合に前記特定の気象現象の確率情報を通知する通知手段と
を具備することを特徴とする気象予測データ解析装置。
Acquisition means for acquiring weather observation data and weather prediction data calculated based on a weather prediction model based on the weather observation data;
Based on the acquired weather forecast data and weather observation data, determination means for respectively determining whether or not a plurality of weather determination elements can be adapted for a specific weather phenomenon;
For weather forecast data and weather observation data acquired in a certain period by the acquisition means, a totaling means for totaling the determination results of each weather determination element by the determination means;
Creating means for creating a conditional probability table (CPT) representing a conditional probability in a dependency relationship between elements for each of the plurality of weather determination elements based on the aggregated results;
When the weather prediction data and the weather observation data are acquired by the acquisition unit, the occurrence probability of the specific weather phenomenon is determined by referring to a conditional probability table corresponding to the weather determination element for the determination result by the determination unit. A probability calculating means for calculating;
A weather prediction data analysis apparatus comprising: a notification unit that notifies probability information of the specific weather phenomenon when the calculated occurrence probability exceeds a threshold value.
前記条件付確率表は、前記複数の気象判定要素間の依存関係により構成されるベイジアンネットワークをもとに作成されることを特徴とする請求項1記載の気象予測データ解析装置。   The weather prediction data analysis apparatus according to claim 1, wherein the conditional probability table is created based on a Bayesian network configured by dependency relationships between the plurality of weather determination elements. 前記作成手段は、前記取得手段により前記気象予測データ及び気象観測データが取得されたとき、この取得された気象予測データ及び気象観測データをもとに前記作成された条件付確率表それぞれの条件付確率の値を更新することを特徴とする請求項1記載の気象予測データ解析装置。   When the weather forecast data and weather observation data are obtained by the obtaining means, the creating means includes a conditional probability table for each of the created conditional probability tables based on the obtained weather forecast data and weather observation data. The weather prediction data analysis apparatus according to claim 1, wherein the probability value is updated. 気象観測データと、前記気象観測データをもとに気象予測モデルに基づいて演算された気象予測データとを取得し、
前記取得された気象予測データ及び気象観測データをもとに、特定の気象現象について複数の気象判定要素の適合の可否をそれぞれ判定し、
前記取得された一定期間の気象予測データ及び気象観測データについて、前記判定された各気象判定要素の判定結果をそれぞれ集計し、
前記集計された結果をもとに前記複数の気象判定要素間それぞれについて要素間の依存関係における条件付確率を表す条件付確率表(CPT:Conditional Probability Table)を作成し、
前記気象予測データ及び気象観測データが取得されたとき、前記判定された結果について前記気象判定要素に対応する条件付確率表を参照して、前記特定の気象現象の生起確率を算出し、
前記算出された生起確率が閾値を超えた場合に前記特定の気象現象の確率情報を通知することを特徴とする気象予測データ解析方法。
Acquiring weather observation data and weather prediction data calculated based on a weather prediction model based on the weather observation data;
Based on the acquired weather forecast data and weather observation data, determine whether or not a plurality of meteorological determination elements can be applied to a specific weather phenomenon,
For the acquired weather forecast data and weather observation data for a certain period, the determination results of each of the determined weather determination elements are totaled,
A conditional probability table (CPT) representing a conditional probability in a dependency relation between elements for each of the plurality of weather determination elements based on the aggregated results is created,
When the weather forecast data and weather observation data are acquired, referring to a conditional probability table corresponding to the weather determination element for the determined result, the occurrence probability of the specific weather phenomenon is calculated,
Providing probability information of the specific weather phenomenon when the calculated occurrence probability exceeds a threshold value, a weather prediction data analysis method,
前記条件付確率表は、前記複数の気象判定要素間の依存関係により構成されるベイジアンネットワークをもとに作成されることを特徴とする請求項4記載の気象予測データ解析方法。   5. The weather forecast data analysis method according to claim 4, wherein the conditional probability table is created based on a Bayesian network configured by dependency relationships between the plurality of weather determination elements. 前記取得手段により前記気象予測データ及び気象観測データが取得されたとき、この取得された気象予測データ及び気象観測データをもとに前記作成された条件付確率表それぞれの条件付確率の値を更新することを特徴とする請求項4記載の気象予測データ解析方法。   When the weather forecast data and weather observation data are acquired by the acquisition means, the conditional probability values of the prepared conditional probability tables are updated based on the acquired weather prediction data and weather observation data. The meteorological prediction data analysis method according to claim 4, wherein:
JP2007218762A 2007-08-24 2007-08-24 Weather forecast data analysis apparatus and weather forecast data analysis method Active JP4908346B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2007218762A JP4908346B2 (en) 2007-08-24 2007-08-24 Weather forecast data analysis apparatus and weather forecast data analysis method

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2007218762A JP4908346B2 (en) 2007-08-24 2007-08-24 Weather forecast data analysis apparatus and weather forecast data analysis method

Publications (2)

Publication Number Publication Date
JP2009052976A true JP2009052976A (en) 2009-03-12
JP4908346B2 JP4908346B2 (en) 2012-04-04

Family

ID=40504188

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2007218762A Active JP4908346B2 (en) 2007-08-24 2007-08-24 Weather forecast data analysis apparatus and weather forecast data analysis method

Country Status (1)

Country Link
JP (1) JP4908346B2 (en)

Cited By (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011028339A (en) * 2009-07-21 2011-02-10 Nippon Telegr & Teleph Corp <Ntt> Information processor, information processing system and information processing method
JP2015190866A (en) * 2014-03-28 2015-11-02 株式会社東芝 weather prediction error analysis system and weather prediction error method
WO2015182142A1 (en) * 2014-05-28 2015-12-03 国立大学法人京都大学 Weather forecasting system
JP2016057152A (en) * 2014-09-09 2016-04-21 株式会社東芝 Weather prediction error analysis device and weather prediction error analysis method
KR101695797B1 (en) * 2015-10-20 2017-01-13 한국외국어대학교 연구산학협력단 Method for meteorological information analysis by linking a plurality of different classification algorithms
JP2018010015A (en) * 2017-10-18 2018-01-18 株式会社東芝 Weather prediction error analysis system and weather prediction error analysis method
CN107765348A (en) * 2017-10-23 2018-03-06 成都信息工程大学 Precipitation meteorological telegraphic messages data error modification method, glaze landscape forecast method
JP2020533897A (en) * 2017-09-11 2020-11-19 ワールドビュー・サテライツ・リミテッド Satellite systems and methods to deal with rainfall fades
CN112213727A (en) * 2020-10-15 2021-01-12 国家卫星气象中心(国家空间天气监测预警中心) Precipitation correction method of satellite-borne radar based on active and passive microwave combined detection
CN114547951A (en) * 2022-04-24 2022-05-27 浙江远算科技有限公司 Dam state prediction method and system based on data assimilation

Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004355075A (en) * 2003-05-27 2004-12-16 Sony Corp Information exhibition device, information exhibition method and computer program
JP2005055190A (en) * 2003-08-01 2005-03-03 Japan Science & Technology Agency Electromagnetic-wave radiation source detecting method and device by bayesian network
JP2006018688A (en) * 2004-07-02 2006-01-19 Toyota Motor Corp Road environment recognition method and road environment recognition apparatus
JP2006181037A (en) * 2004-12-27 2006-07-13 Fuji Photo Film Co Ltd Diagnosis assisting device, diagnosis assisting method and its program
JP2006259862A (en) * 2005-03-15 2006-09-28 Toshiba Corp Bayesian network creating device and method, and program
JP2007102647A (en) * 2005-10-06 2007-04-19 Central Res Inst Of Electric Power Ind Questionnaire preparation method and system

Patent Citations (6)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2004355075A (en) * 2003-05-27 2004-12-16 Sony Corp Information exhibition device, information exhibition method and computer program
JP2005055190A (en) * 2003-08-01 2005-03-03 Japan Science & Technology Agency Electromagnetic-wave radiation source detecting method and device by bayesian network
JP2006018688A (en) * 2004-07-02 2006-01-19 Toyota Motor Corp Road environment recognition method and road environment recognition apparatus
JP2006181037A (en) * 2004-12-27 2006-07-13 Fuji Photo Film Co Ltd Diagnosis assisting device, diagnosis assisting method and its program
JP2006259862A (en) * 2005-03-15 2006-09-28 Toshiba Corp Bayesian network creating device and method, and program
JP2007102647A (en) * 2005-10-06 2007-04-19 Central Res Inst Of Electric Power Ind Questionnaire preparation method and system

Cited By (14)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2011028339A (en) * 2009-07-21 2011-02-10 Nippon Telegr & Teleph Corp <Ntt> Information processor, information processing system and information processing method
JP2015190866A (en) * 2014-03-28 2015-11-02 株式会社東芝 weather prediction error analysis system and weather prediction error method
WO2015182142A1 (en) * 2014-05-28 2015-12-03 国立大学法人京都大学 Weather forecasting system
US10288769B2 (en) 2014-05-28 2019-05-14 Metroweather Co., Ltd. Weather forecasting system
JP2016057152A (en) * 2014-09-09 2016-04-21 株式会社東芝 Weather prediction error analysis device and weather prediction error analysis method
KR101695797B1 (en) * 2015-10-20 2017-01-13 한국외국어대학교 연구산학협력단 Method for meteorological information analysis by linking a plurality of different classification algorithms
JP7183258B2 (en) 2017-09-11 2022-12-05 ワールドビュー・サテライツ・リミテッド Satellite system and method for dealing with rain fade
JP2020533897A (en) * 2017-09-11 2020-11-19 ワールドビュー・サテライツ・リミテッド Satellite systems and methods to deal with rainfall fades
JP2018010015A (en) * 2017-10-18 2018-01-18 株式会社東芝 Weather prediction error analysis system and weather prediction error analysis method
CN107765348A (en) * 2017-10-23 2018-03-06 成都信息工程大学 Precipitation meteorological telegraphic messages data error modification method, glaze landscape forecast method
CN112213727A (en) * 2020-10-15 2021-01-12 国家卫星气象中心(国家空间天气监测预警中心) Precipitation correction method of satellite-borne radar based on active and passive microwave combined detection
CN112213727B (en) * 2020-10-15 2024-01-02 国家卫星气象中心(国家空间天气监测预警中心) Precipitation correction method of satellite-borne radar based on active and passive microwave combined detection
CN114547951A (en) * 2022-04-24 2022-05-27 浙江远算科技有限公司 Dam state prediction method and system based on data assimilation
CN114547951B (en) * 2022-04-24 2022-07-22 浙江远算科技有限公司 Dam state prediction method and system based on data assimilation

Also Published As

Publication number Publication date
JP4908346B2 (en) 2012-04-04

Similar Documents

Publication Publication Date Title
JP4908346B2 (en) Weather forecast data analysis apparatus and weather forecast data analysis method
Emerton et al. Continental and global scale flood forecasting systems
Arnbjerg-Nielsen et al. Impacts of climate change on rainfall extremes and urban drainage systems: a review
Zhang et al. Predicting Typhoon Morakot’s catastrophic rainfall with a convection-permitting mesoscale ensemble system
JP2008050903A (en) Flood prediction method and flood prediction system
Sam et al. Heterogeneity in topographic control on velocities of Western Himalayan glaciers
TW201807434A (en) Method and system for predicting the impact of forecasted weather, environmental and/or geologic conditions
Tramblay et al. Extreme value modelling of daily areal rainfall over Mediterranean catchments in a changing climate
JP4818079B2 (en) Weather forecast data analysis apparatus and weather forecast data analysis method
Woodworth et al. Towards reliable global allowances for sea level rise
Vieux et al. Statistical evaluation of a radar rainfall system for sewer system management
Verstraten et al. Sensitivity of Australian roof drainage structures to design rainfall variability and climatic change
JP2008008772A (en) Weather forecasting system
Juan et al. Developing a radar-based flood alert system for Sugar Land, Texas
JP2006220445A (en) Sandstorm prediction system and its sandstorm alarm issuance method
US20170102482A1 (en) Forensic weather system
US11493666B2 (en) System and method for forecasting snowfall probability distributions
Rathore et al. Collaborative mechanism for tropical cyclone monitoring and prediction over north Indian ocean
Chiu et al. A watershed-based debris flow early warning system using sensor web enabling techniques in heterogeneous environments
JP2007187478A (en) Weather-forecasting system
Hoblit et al. Integrating radar rainfall estimates with digital elevation models and land use data to create an accurate hydrologic model
Looper et al. Distributed hydrologic forecast reliability using next-generation radar
Kumar et al. Statistical evaluation of snow accumulation and depletion from remotely sensed MODIS snow time series data using the SARIMA model
Rosbjerg et al. Nordic contributions to stochastic methods in hydrology
Biantoro et al. Time of Concentration Estimated Using Some Methods and Application in The Ciliwung Watershed, Jakarta

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20090914

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20110826

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20110913

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20111111

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: 20111213

A01 Written decision to grant a patent or to grant a registration (utility model)

Free format text: JAPANESE INTERMEDIATE CODE: A01

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20120112

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150120

Year of fee payment: 3

R151 Written notification of patent or utility model registration

Ref document number: 4908346

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R151

FPAY Renewal fee payment (event date is renewal date of database)

Free format text: PAYMENT UNTIL: 20150120

Year of fee payment: 3

S111 Request for change of ownership or part of ownership

Free format text: JAPANESE INTERMEDIATE CODE: R313111

R350 Written notification of registration of transfer

Free format text: JAPANESE INTERMEDIATE CODE: R350