JP3514564B2 - Mirage occurrence prediction method and prediction device - Google Patents
Mirage occurrence prediction method and prediction deviceInfo
- Publication number
- JP3514564B2 JP3514564B2 JP31231395A JP31231395A JP3514564B2 JP 3514564 B2 JP3514564 B2 JP 3514564B2 JP 31231395 A JP31231395 A JP 31231395A JP 31231395 A JP31231395 A JP 31231395A JP 3514564 B2 JP3514564 B2 JP 3514564B2
- Authority
- JP
- Japan
- Prior art keywords
- mirage
- meteorological
- temperature
- data
- occurrence
- 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.)
- Expired - Fee Related
Links
- 208000030984 MIRAGE syndrome Diseases 0.000 title claims description 77
- TVLSRXXIMLFWEO-UHFFFAOYSA-N prochloraz Chemical compound C1=CN=CN1C(=O)N(CCC)CCOC1=C(Cl)C=C(Cl)C=C1Cl TVLSRXXIMLFWEO-UHFFFAOYSA-N 0.000 title claims description 69
- 238000000034 method Methods 0.000 title claims description 22
- 238000004458 analytical method Methods 0.000 claims description 8
- 238000012545 processing Methods 0.000 claims description 3
- 238000004364 calculation method Methods 0.000 description 5
- 230000003442 weekly effect Effects 0.000 description 4
- 238000011002 quantification Methods 0.000 description 3
- 238000011160 research Methods 0.000 description 3
- 238000007619 statistical method Methods 0.000 description 3
- 238000012795 verification Methods 0.000 description 3
- 238000011161 development Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 239000000155 melt Substances 0.000 description 2
- 239000013535 sea water Substances 0.000 description 2
- 238000012876 topography Methods 0.000 description 2
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000002349 favourable effect Effects 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000000491 multivariate analysis Methods 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
Landscapes
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Description
【0001】[0001]
【発明の属する技術分野】この発明は蜃気楼の発生予測
方法に係り、特に明日以降少なくとも7日先までの蜃気
楼発生の有無を統計的手法をかりて精度よく予測する方
法及び装置に関する。BACKGROUND OF THE INVENTION 1. Field of the Invention The present invention relates to a mirage occurrence prediction method, and more particularly to a method and apparatus for accurately predicting the occurrence of a mirage from tomorrow to at least 7 days ahead by using a statistical method.
【0002】[0002]
【従来の技術】文献によると蜃気楼は、古来より北海道
オホーツク海沿岸、東北から山陰地方にいたる日本海沿
岸、或いは千葉県、和歌山県、鹿児島県をはじめ太平洋
沿岸各地で観察されていた現象であり、その発生原因に
ついてはかなり以前から調査され議論されていたが、蜃
気楼の発生を予め予知するまでには至っていない。中で
も毎年3月上旬から6月下旬にかけて出現する富山県魚
津市の蜃気楼は世界的にも有名であり、魚津市の蜃気楼
研究会では過去20年間にわたって当地の蜃気楼現象を
観察しており、発生時の気象条件等を記載した貴重な記
録が大切に保存されている。しかし、蜃気楼の発生と気
象条件との統計的な関係はまだ解明されておらず、その
出現の有無は過去の経験と勘に頼らざるを得ないのが現
状である。2. Description of the Related Art According to the literature, mirage is a phenomenon that has been observed since ancient times along the coast of the Sea of Okhotsk in Hokkaido, along the coast of the Sea of Japan from Tohoku to the San'in region, or around the Pacific coast including Chiba, Wakayama and Kagoshima prefectures. , The cause of its occurrence has been investigated and discussed for a long time, but it has not been predicted in advance about the occurrence of mirage. The mirage in Uozu City, Toyama Prefecture, which appears every year from the beginning of March to the end of June, is famous worldwide, and the mirage research group in Uozu has observed the mirage phenomenon here for the past 20 years. Valuable records that describe the meteorological conditions, etc. are kept carefully. However, the statistical relationship between the occurrence of mirages and meteorological conditions has not yet been clarified, and it is the current situation that the occurrence of mirages must be relied on past experience and intuition.
【0003】[0003]
【発明が解決しようとする課題】しかしながら上記のよ
うに、蜃気楼出現の予測を人の経験と勘に依存する手法
では、精度の高い適切な予測を期待することは不可能で
あり、蜃気楼の発生を地方開発、街おこしの観光資源と
して利用するとしても、信頼性に欠けるばかりでなく、
発生の日時を予測する人に大きな負担となることは必定
である。However, as described above, it is impossible to expect a highly accurate and appropriate prediction with the method that depends on the experience and intuition of a person to predict the appearance of a mirage. Not only is it lacking in reliability even if it is used as a tourism resource for rural development and town development,
It is inevitable that a person who predicts the date and time of occurrence will be a heavy burden.
【0004】この発明は、このような事情に鑑み、蜃気
楼の発生と気象条件との統計的な関係を明らかにし、蜃
気楼出現の有無を正確かつ容易に予測できるようにする
ことを目的とする。In view of such circumstances, an object of the present invention is to clarify the statistical relationship between the occurrence of a mirage and meteorological conditions, and to accurately and easily predict the presence or absence of a mirage.
【0005】[0005]
【課題を解決するための手段】本発明者は、上記の目的
を達成するため鋭意研究を行った結果、長年に亘る蜃気
楼観察記録データの解明に統計的手法を適用し、かつ、
その解析用アイテムに特定の気象データを選択・使用す
ることにより、所期の目的が達成されることを見出だ
し、この発明を完成するに至った。Means for Solving the Problems As a result of intensive research to achieve the above-mentioned object, the present inventor has applied a statistical method to elucidate mirage observation record data for many years, and
The inventors have found that the intended purpose can be achieved by selecting and using specific meteorological data for the analysis item, and have completed the present invention.
【0006】すなわち、この発明は、過去の蜃気楼発生
に係わる観測気象データを解析して蜃気楼発生の有無を
判別するための予測式を導出し、これを用いて対象地域
における予測気象データを解析し、その解析値を過去の
観測気象データの解析値と照合することにより少なくと
も7日先までの蜃気楼発生の有無を判別することを特徴
とするものであり、また、上記の方法を実施するに当た
り、予め過去の蜃気楼発生に係る観測気象データを入力
して蜃気楼発生の有無を解析・判別できるよう学習制御
した演算処理システムに、対象地域における現在以降の
予測気象データを入力することにより、少なくとも7日
先までの蜃気楼の発生予測情報を出力させることを特徴
とする蜃気楼の発生予測方法を要旨とするものである。That is, the present invention analyzes the observed meteorological data relating to the occurrence of a mirage in the past, derives a prediction formula for determining the presence or absence of a mirage, and uses this to analyze the predicted meteorological data in the target area. By comparing the analysis value with the analysis value of the past observed meteorological data, it is possible to determine whether or not a mirage has occurred at least 7 days ahead. Further, in carrying out the above method, By inputting the forecasted meteorological data after the present in the target area into the arithmetic processing system that has been learning-controlled so that the observed meteorological data related to the past mirage occurrences can be input and analyzed / determined whether or not the mirage occurred The gist is a method of predicting the occurrence of a mirage, which is characterized by outputting the prediction information of the occurrence of a mirage.
【0007】またさらに、蜃気楼発生に係わる過去の観
測気象データを解析して蜃気楼発生の有無を判別するた
めの予測式を導出する演算装置と、この演算装置に過去
の観測気象データ及び対象地域における予測気象データ
を入力する入力装置と、入力データを予測式で解析して
得られた蜃気楼発生の有無の判別結果を出力する出力装
置とからなる蜃気楼の発生予測装置を要旨とするもので
ある。Furthermore, a computing device for analyzing the past observed meteorological data relating to the occurrence of the mirage and deriving a prediction formula for determining the presence or absence of the mirage, and the computing device for observing the past observed meteorological data and the target area The gist of a mirage generation prediction device is an input device that inputs predicted weather data, and an output device that outputs a determination result of the presence or absence of a mirage obtained by analyzing the input data by a prediction formula.
【0008】[0008]
【発明の実施の態様】本発明を実施するに当っては、ま
ず蜃気楼発生の有無と関係が深い気象項目を選択する必
要がある。そこで、過去の文献にみられる蜃気楼発生記
録や気象関係資料をつぶさに検討したところ、各対象地
域に共通する気象条件として、海面または地上に接する
気層の気温変化率が高くて密度の異なる大気層が形成さ
れており、大気が安定していてこの状態が崩れにくく、
晴天で見通しがよい状態にあるとき蜃気楼が発生し易く
また観測しやすいことが確認された。BEST MODE FOR CARRYING OUT THE INVENTION In carrying out the present invention, it is first necessary to select a meteorological item closely related to the presence or absence of a mirage. Therefore, a close examination of the mirage occurrence records and meteorological data found in the past literature shows that the common atmospheric conditions in each target area are atmospheric layers with high temperature change rates and different densities in the air layer contacting the sea surface or the ground. Is formed, the atmosphere is stable and this state is difficult to collapse,
It was confirmed that the mirage is likely to occur and is easy to observe when the weather is fine and the visibility is good.
【0009】上記の気象条件に係る具体的な気象項目と
しては、通常、移動性高気圧中心位置、サーマルリッジ
存在エリアおよびサーマルトラフ存在エリア等の広域気
象データと、海面または地表温度と気温との差、大気低
層の気温減率、ショワルターのスタビリティ・インデッ
クス、最高気温、天気、日照時間、雨量、風速および風
向、冬季降雪量等の狭域気象データがこれに当る。As specific meteorological items related to the above-mentioned meteorological conditions, generally, wide area meteorological data such as a mobile high pressure center position, a thermal ridge existing area and a thermal trough existing area, and a difference between sea surface or surface temperature and air temperature. This includes narrow-range meteorological data such as the temperature decrease rate in the lower atmosphere, the stability index of Showerter, maximum temperature, weather, sunshine hours, rainfall, wind speed and direction, and winter snowfall.
【0010】本発明を実施するに当っては、上記気象項
目の全てを使用することも可能であるが、有効性の少な
い項目や内部相関の高い項目を探して削減し、できるだ
け少ない気象項目データを使用することが望ましい。こ
れらの各気象項目の重みは常に一定ではなく、対象地域
の地形など地理的条件により変動するものなので、この
作業が適切であれば計算所要時間が短縮されるばかりで
なく、計算精度を高めることも可能である。In practicing the present invention, it is possible to use all of the above-mentioned meteorological items, but the items that are less effective or have a high internal correlation are searched for and reduced, and as few meteorological item data as possible are used. Is preferred. The weight of each of these meteorological items is not always constant and varies depending on geographical conditions such as the topography of the target area. Therefore, if this work is appropriate, not only the time required for calculation but also the accuracy of calculation should be improved. Is also possible.
【0011】例えば、気温より海面または地表温度の方
が低くなる季節に観察される蜃気楼(以下、春夏型蜃気
楼と略称する)については、移動性高気圧中心位置とサ
ーマルリッジ存在エリアの2ケの広域気象データと、海
面または地表温度と気温との差、最高気温、ショワルタ
ーのスタビリティ・インデックス、天気、日照時間、最
大風速時の風向および冬季降雪量の中から選択される2
〜7ケの狭域気象データとからなる計4〜9ケの気象項
目データを使用するのが好適である。For example, for mirages (hereinafter abbreviated as spring-summer mirages) that are observed in the season when the sea surface temperature or the surface temperature is lower than the air temperature, there are two locations, a mobile high pressure center position and a thermal ridge existence area. Selected from wide area meteorological data, difference between sea surface or surface temperature and temperature, maximum temperature, Shoalter's stability index, weather, sunshine hours, wind direction at maximum wind speed and winter snowfall 2
It is preferable to use a total of 4 to 9 meteorological item data consisting of 7 to 7 narrow area meteorological data.
【0012】また、気温よりも海面または地表温度の方
が高くなる季節に観察される蜃気楼(以下、冬型蜃気楼
と略称する)については、移動性高気圧中心位置、サー
マルトラフ存在エリアおよびサーマルリッジ存在エリア
からなる広域気象データと、海面または地表温度と気温
の差、ショワルターのスタビリティ・インデックス、天
気、日照時間、および最大風速生起時の風向からなる狭
域気象データの中から計6〜8ケの気象項目データを使
用するのが好適である。Regarding mirages (hereinafter abbreviated as winter mirages) observed in the season when the sea surface temperature or the surface temperature is higher than the air temperature, the movable high pressure central position, the thermal trough existing area, and the thermal ridge existing area. 6-8 in total from regional wide area meteorological data and difference in sea surface or surface temperature and temperature, stability index of shoalter, weather, sunshine hours, and wind direction at maximum wind speed. It is preferable to use the weather item data of.
【0013】前記のごとき気象項目を選定し使用する理
由は下記の通りである。
1.海面または地表に接する気層の気温変化率:海面ま
たは地表の温度が低くて気温が高いか、その逆の場合
は、気層鉛直方向の気温変化率が大となり蜃気楼が出現
し易い。すなわち、前記の春夏型蜃気楼は日中の気温が
上昇して海水温度より高くなる季節に発生し、しかも、
最高気温が高くて海面に接する気層の気温変化率が大き
くなる日ほど発生確率が高い。対象地域の近辺にサーマ
ルリッジ(温度の峰であり、天気図では等温線が北に盛
り上がっている部分)が存在することはかかる気象状態
の出現に極めて好都合である。また逆に、冬型蜃気楼の
発生時には気温が下降して海水温度よりかなり低くなる
ことが条件である。対象地域の近辺にサーマルトラフ
(温度の谷であり天気図では等温線が南に垂れ下がって
いる場所を意味する)が存在することは、気温低下の目
安となり冬型蜃気楼の出現にとって好ましいことであ
る。なお、対象地域によっては海面温度に替えて、雪解
け水となり海面温度を左右する積雪やこれを溶かす雨量
を使用することもできる。The reasons for selecting and using the above-mentioned meteorological items are as follows. 1. Rate of temperature change in the air layer in contact with the sea surface or surface: If the temperature of the sea surface or the surface of the earth is low and the temperature is high, or vice versa, the temperature change rate in the vertical direction of the air layer becomes large and a mirage tends to appear. That is, the above-mentioned spring-summer mirage occurs in the season when the daytime temperature rises and becomes higher than the seawater temperature, and
The probability of occurrence is higher on days when the maximum temperature is higher and the temperature change rate of the air layer in contact with the sea surface is higher. The presence of a thermal ridge (a temperature peak, where the isotherm rises to the north in the weather map) near the target area is extremely convenient for the appearance of such weather conditions. On the other hand, when the winter mirage occurs, the temperature drops and becomes considerably lower than the seawater temperature. The presence of a thermal trough near the target area (a valley of temperature that means that the isotherm hangs to the south in the weather map) is a good indicator of the temperature drop and is preferable for the emergence of a winter mirage. Depending on the target area, instead of the sea surface temperature, it is also possible to use snowfall that melts snow and affects the sea surface temperature and the amount of rainfall that melts the snow.
【0014】2.大気の安定度:一般に高気圧圏内では
大気の状態は安定している。大気上層の安定度を示す尺
度としては、通常ショワルターのスタビリティ・インデ
ックス(SSI)が使われる。該SSI値は、850h
Pa面の空気塊を断熱的に500hPaの高度に持ち上
げたときの周囲との温度差を示すものであり、この値が
大きい(つまり、持ち上げた空気の方が冷たい)ときは
高気圧圏内にあり上層の大気が安定しているといえる。
また、大気下層の安定度を図る尺度としては、下式で求
まる低層気温減率を使用する。
低層気温減率(℃/100m)=(T1000−Ts )/Z1000
但し、T1000:1,000hPa面の気温
Ts :地上の気温
Z1000:1,000hPa面の高度
通常、大気中では上空へいくほど温度が低く、その減率
は100メートルにつき0.6度から1.0度ぐらいで
ある。減率がこれより大きいと大気は不安定である。逆
に上空ほど気温が高い場合を「逆転」と呼び、この場合
は大気は非常に安定である。2. Atmospheric stability: Atmospheric conditions are generally stable in high pressure areas. The Schwarter Stability Index (SSI) is usually used as a measure of the stability of the upper atmosphere. The SSI value is 850h
It shows the temperature difference from the surroundings when the air mass on the Pa side is adiabatically lifted to an altitude of 500 hPa. When this value is large (that is, the lifted air is colder), it is in the high pressure area and the upper layer. It can be said that the atmosphere is stable.
In addition, as a measure for the stability of the lower atmosphere, the lower temperature decrease rate obtained by the following equation is used. Lower air temperature decrease rate (℃ / 100m) = (T 1000 −T s ) / Z 1000 However, T 1000 : 1,000 hPa surface temperature T s : Ground temperature Z 1000 : 1,000 hPa surface altitude Normally, in the atmosphere Then, the temperature becomes lower toward the sky, and the rate of decrease is about 0.6 to 1.0 degrees per 100 meters. If the rate of decrease is higher than this, the atmosphere is unstable. Conversely, when the temperature is higher in the sky, it is called "reversal", and in this case, the atmosphere is very stable.
【0015】3.天候:まず天気は、蜃気楼の発生及び
観測のためには曇天よりも晴天が望ましい。また、地域
によっては当日だけでなく、前日および前々日の天気が
蜃気楼の発生と関係する。風速は、あまり強くなく、あ
まり弱くもない風速が適している。風向は、地域によっ
て異なるが、地形との相関関係において好ましい風向が
存在する。海または湖沼に面した地方では海面または水
面の温度が低く、陸上は晴天で気温が高い状態のとき蜃
気楼の出現確率が高い。このような時には局地的な海風
または湖沼風が吹くことになる。風向・風速はこのよう
な観点からも蜃気楼発生を予測するための指標となる。3. Weather: First of all, it is preferable that the weather is clearer than cloudy for the generation and observation of the mirage. In addition, in some regions, the weather of the day before and the day before the day is related to the occurrence of the mirage. The wind speed is not so strong and not so weak. The wind direction differs depending on the region, but there is a favorable wind direction in correlation with the topography. In regions facing the sea or lakes, the temperature of the sea or water is low, and on land, there is a high probability that a mirage will appear when the weather is sunny and the temperature is high. At such times, a local sea breeze or lake breeze will blow. From such a viewpoint, the wind direction and speed are also indicators for predicting the occurrence of a mirage.
【0016】本発明を実施するに当っては、明らかにし
たい事柄(蜃気楼の出現の有無)とそれに関係する多く
の事柄(前記の各気象項目データ)をともに質的データ
に変換し、数量化II類と称される多変量解析法を使用す
るのが便利である。分析を実行するに当っては、まず過
去の観測気象データについて、蜃気楼出現の有無を目的
変数とし、特定の気象項目データを説明変数としてその
間の関係を調べ、目的変数に強い影響を及ぼしている説
明変数を発見し、蜃気楼発生の有無を判別するための予
測式を導出する。ついで、この予測式を用いて明日以降
の予測気象データを解析し、得られた結果を過去の観測
気象データの解析値と照合して蜃気楼発生の有無を判別
することになる。In carrying out the present invention, both the matters to be clarified (whether the mirage appears or not) and many matters related thereto (the above-mentioned meteorological item data) are converted into qualitative data and quantified. It is convenient to use a multivariate analysis method called class II. In carrying out the analysis, first, regarding the past observed meteorological data, the presence or absence of a mirage is used as the objective variable, and the specific meteorological item data is used as the explanatory variable, and the relationship between them is investigated, and the objective variable is strongly influenced. Discover the explanatory variables and derive the prediction formula for determining the presence or absence of a mirage. Next, using this prediction formula, the forecasted weather data from tomorrow onward will be analyzed, and the obtained results will be collated with the analyzed values of the past observed weather data to determine whether or not a mirage has occurred.
【0017】上記において使用する過去の観測気象デー
タは、主として日本気候表、印刷天気図、高層気象観測
データ、地方気象台気象月報等の気象庁所管の各種気象
関係資料より得られ、現在以降の予測気象データは気象
庁のGPV、週間予報並びに予想天気図などの狭域ない
し広域気象データを基に作成されるものである。この場
合、判別精度を高めるには上記GPVを5キロメッシュ
の予想値に加工して更に狭域の気象データとして使用す
ることが望ましい。The past observed meteorological data used in the above is mainly obtained from various weather-related materials under the jurisdiction of the Meteorological Agency, such as the Japan climate table, printed weather maps, high-level meteorological observation data, and local monthly weather report. The data is created on the basis of narrow-range or wide-area weather data such as the Japan Meteorological Agency's GPV, weekly forecasts, and forecast weather maps. In this case, in order to improve the discrimination accuracy, it is desirable to process the GPV into an expected value of 5 km mesh and use it as meteorological data in a narrow area.
【0018】以上の通りこの発明においては、統計的手
法を用いて蜃気楼の発生を予測するようにしたので、従
来人の勘と経験に依存していた方法に比べて高精度の予
測が可能となり、しかも明日の予報は勿論、少なくとも
7日先までの予報を速やかに提供できるようになった。
これにより蜃気楼が観光サービス資源として十分役立つ
日も近くなることであろう。As described above, in the present invention, since the occurrence of the mirage is predicted by using the statistical method, it is possible to make the prediction with higher accuracy than the method which has conventionally depended on the intuition and experience of the person. Moreover, it is now possible to provide forecasts for at least 7 days ahead, as well as tomorrow's forecasts.
As a result, the mirage will soon be fully utilized as a tourism service resource.
【0019】[0019]
【実施例】つぎに、実施例を挙げて本発明をより詳細に
説明する。
実施例1
表1は明日の蜃気楼発生の有無を予測するに際し、適用
される目的変数と8ケの説明変数を示した変数表であ
る。これを数量化II類による解析モデルとし、富山湾に
おける既存の春夏型蜃気楼の観察記録(1976年4月
から1994年6月までの間の観察データ)282サン
プル(うち10ケは検証用独立データとして使用)につ
き解析を行った。EXAMPLES Next, the present invention will be described in more detail with reference to examples. Example 1 Table 1 is a variable table showing objective variables and eight explanatory variables that are applied in predicting the presence or absence of a mirage tomorrow. Using this as a quantification type II analysis model, 282 samples of observation records of existing spring-summer mirages in Toyama Bay (observation data from April 1976 to June 1994) (of which 10 are independent for verification) (Used as data) was analyzed.
【0020】[0020]
【表1】 [Table 1]
【0021】上記8ケの説明変数、すなわち、最高気温
の平年偏差、日照時間、SSI、降雪の深さなどの数量
データ、天気、風向などの文字データ、および移動性高
気圧中心位置、サーマルリッジ存在エリアなどの図形デ
ータを全てカテゴリーデータに変換した。図1は移動性
高気圧中心位置のカテゴリー化を、図2は850hPa
面サーマルリッジ存在エリアのカテゴリー化を示すもの
である。数量データと文字データのカテゴリー区分は表
2と表3に纏めて示した。The above eight explanatory variables, namely, normal temperature deviation of maximum temperature, sunshine hours, SSI, numerical data such as snowfall depth, character data such as weather and wind direction, mobile high pressure center position, presence of thermal ridge All graphic data such as areas were converted into category data. Fig. 1 shows the categorization of the mobile high pressure center position, and Fig. 2 shows 850 hPa.
It shows the categorization of the surface thermal ridge existing area. The categorical classification of quantitative data and character data is summarized in Tables 2 and 3.
【0022】上記のとおり、カテゴリーデータに変換さ
れた変数をコンピュータ演算処理システムに入力するこ
とにより、各説明変数の影響の強さを表す偏相関係数と
レンジおよび、説明変数についてそれぞれのカテゴリー
・スコアと相関比(説明変数全体での影響の強さを表す
係数)を算出した。偏相関係数は表1に、カテゴリース
コアは表2と表3にそれぞれ示されている。As described above, by inputting the variables converted into categorical data into the computer arithmetic processing system, the partial correlation coefficient and range representing the strength of influence of each explanatory variable and the category Scores and correlation ratios (coefficients representing the strength of influence on all explanatory variables) were calculated. The partial correlation coefficients are shown in Table 1, and the category scores are shown in Tables 2 and 3, respectively.
【0023】[0023]
【表2】 [Table 2]
【0024】[0024]
【表3】 [Table 3]
【0025】説明するまでもなく数量化II類は、目的変
数の群の判別に貢献する度合いに応じて、説明変数の各
カテゴリーに得点を与える手法であり、この得点が上記
のカテゴリースコアである。サンプルごとに、当て嵌ま
るカテゴリーのカテゴリースコアを足し合わせた値をサ
ンプルスコアといい、この値から各々のサンプルが目的
変数のどの群に属するか(ここでは蜃気楼の発生ありと
蜃気楼の発生なし)を判別する。Needless to say, the quantification type II is a method of giving a score to each category of the explanatory variables according to the degree of contribution to the discrimination of the group of objective variables, and this score is the above category score. . The value obtained by adding the category scores of the applicable categories for each sample is called the sample score, and from this value it is determined which group each sample belongs to in the objective variable (here, mirage occurrence and mirage occurrence). Determine.
【0026】図3は、上記270ケのサンプルについて
計算されたサンプルスコアの度数分布を表すグラフであ
る。図において右側の山は目的変数の第1群(蜃気楼発
生あり)を表し、左側の山は第2群(蜃気楼発生なし)
を表している。このように判別グラフが2つの山にはっ
きりと分かれるということは、計算精度の良さを裏付け
ている。FIG. 3 is a graph showing the frequency distribution of sample scores calculated for the above 270 samples. In the figure, the mountain on the right side represents the first group of objective variables (with mirage occurrence), and the mountain on the left side is the second group (without mirage occurrence).
Is represented. The fact that the discrimination graph is clearly divided into two peaks in this way supports the good calculation accuracy.
【0027】なお、表2、表3および図3から明らかな
ように、計算精度については相関比0.7952,判別
的中率は91.5%という良好な結果が得られた。As is clear from Tables 2 and 3, and FIG. 3, good results were obtained in terms of calculation accuracy, with a correlation ratio of 0.7952 and a discriminant predictive ratio of 91.5%.
【0028】ついで、前記のとおり12ケの検証用独立
データにより検証したところ、的外れは12サンプルの
中1例であり、判別精度は0.92(=11/12)で
あった。この値は前記した判別的中率とほぼ同じであ
る。Then, as described above, when 12 pieces of independent data for verification were used for verification, one of the 12 samples had a miss, and the discrimination accuracy was 0.92 (= 11/12). This value is almost the same as the above-mentioned discriminant predictive value.
【0029】実施例2
7日先までの蜃気楼の発生予測を行う手法について検討
した。予測用データとしては、ルーチンに入手可能なデ
ータを使用することとし、説明変数は表4に示す4ケを
使用した。各説明変数のカテゴリーデータへの変換は実
施例1と同様である。Example 2 A method for predicting the occurrence of a mirage up to 7 days ahead was examined. As the predictive data, the data available in the routine were used, and the four explanatory variables shown in Table 4 were used. The conversion of each explanatory variable into categorical data is the same as in the first embodiment.
【0030】[0030]
【表4】 [Table 4]
【0031】ここでも解析手法に数量化II類を使用し、
実施例1と同じ春夏型蜃気楼の観察記録(1976年4
月〜1994年6月までの間の観察データ)282サン
プルにつき解析を行った。この解析モデルについて得ら
れた相関比は0.7678、判別的中率は88.2%で
あった。このように変数が少ないわりには良い精度が得
られたといえる。Here again, the quantification type II is used for the analysis method,
Observation record of the same spring-summer mirage as in Example 1 (April 1976)
282 samples of observation data from the month to June 1994) were analyzed. The correlation ratio obtained for this analytical model was 0.7678, and the discriminant predictive value was 88.2%. It can be said that good accuracy was obtained despite the small number of variables.
【0032】上記の計算結果を予報用に使用する場合の
手順は次の通りである。
1)図4に示す予報用ワークシートに7日先までの毎日
の予報値を記入する。天気と最高気温は地方気象台の富
山週間予報をベースに、独自の修正を加えた値を使用す
る。移動性高気圧中心位置とサーマルリッジ存在エリア
は週間予報支援図をみて、カテゴリーを決める。
2)表5により予報値に対するカテゴリーを決め、その
カテゴリーのカテゴリースコアを読取り、ワークシート
に記入する。
3)ワークシート上で変数4ケのカテゴリースコアを合
計してサンプルスコアを求める。
4)サンプルスコアにより、(判別的中点を参考とし
て)蜃気楼の発生有無を決める。The procedure for using the above calculation result for forecasting is as follows. 1) Enter the daily forecast values up to 7 days ahead on the forecast worksheet shown in FIG. The weather and maximum temperature are based on the Toyama Weekly Forecast of the Local Meteorological Observatory and use their own modified values. The central location of the mobile high pressure and the area where the thermal ridge exists are determined by looking at the weekly forecast support map. 2) Determine the category for the forecast value from Table 5, read the category score of that category, and enter it on the worksheet. 3) Sum up the category scores of four variables on the worksheet to obtain the sample score. 4) Determine the presence or absence of mirage (with reference to the discriminant midpoint) based on the sample score.
【0033】[0033]
【表5】 [Table 5]
【0034】[0034]
【表6】 [Table 6]
【0035】なお、ワークシートを用いる手法に限ら
ず、演算装置を用いて蜃気楼発生の有無を自動出力させ
ることができるのは勿論である。It is needless to say that it is possible to automatically output the presence or absence of a mirage by using an arithmetic device, not limited to the method using the worksheet.
【0036】ついで、1995年4月と5月の2カ月間
(計61サンプル)について、上記の手法を適用し、検
証を行った。結果は表6の通りであリ、外れは61サン
プルの中8例であり、判別精度は0.87(=(6+4
7)/61)であった。Next, the above method was applied and verified for two months (total 61 samples) in April and May 1995. The results are shown in Table 6, and the deviation was 8 out of 61 samples, and the discrimination accuracy was 0.87 (= (6 + 4
It was 7) / 61).
【0037】冬型蜃気楼の発生予測を行う手法について
も、上記と同様な手法を用いて判定できる。As for the method of predicting the occurrence of the winter mirage, the same method as described above can be used for the determination.
【0038】[0038]
【発明の効果】以上説明したようにこの発明は、過去の
気象データを統計的に解析し、これに基づいて蜃気楼の
発生を判別する方法及び装置であるから、信頼性のよい
予測を速やかに行うことができると共に、少なくとも7
日先までの蜃気楼発生の情報を通信サービス網を通じて
何時でもどこからでもアクセスできるようになるという
顕著な効果を有する。As described above, the present invention is a method and apparatus for statistically analyzing past meteorological data and discriminating the occurrence of a mirage based on the statistical data. Therefore, reliable prediction can be made quickly. Can do at least 7
It has a remarkable effect that the information about the occurrence of mirage up to the future can be accessed anytime and anywhere through the communication service network.
【図1】移動性高気圧中心位置のカテゴリー区分を示す
図である。FIG. 1 is a diagram showing category division of a mobile high pressure center position.
【図2】850hPa面サーマルリッジ存在エリアのカ
テゴリー区分を示す図である。FIG. 2 is a diagram showing category division of a thermal ridge existing area of 850 hPa surface.
【図3】サンプルスコアの度数分布を表すグラフの一例
である。FIG. 3 is an example of a graph showing a frequency distribution of sample scores.
【図4】週間予報用ワークシートの一例である。FIG. 4 is an example of a weekly forecast worksheet.
───────────────────────────────────────────────────── フロントページの続き (56)参考文献 特開 平9−49884(JP,A) 特開 平7−225284(JP,A) 特開 平6−300860(JP,A) 特開 平2−221997(JP,A) 特開 平7−260946(JP,A) 特開 平5−134055(JP,A) 特開 平4−253282(JP,A) 吹山直樹、望月東、浅見泰造,“気象 因子を用いた落雷予測の検討”,電気学 会放電高電圧合同研究会資料,日本,社 団法人電気学会,1995年10月26日,ED −95−197〜207,p.23−32 (58)調査した分野(Int.Cl.7,DB名) G01W 1/00 - 1/18 G06F 17/60 - 19/00 JICSTファイル(JOIS)─────────────────────────────────────────────────── ─── Continuation of the front page (56) References JP-A-9-49884 (JP, A) JP-A-7-225284 (JP, A) JP-A-6-300860 (JP, A) JP-A-2- 221997 (JP, A) JP 7-260946 (JP, A) JP 5-134055 (JP, A) JP 4-253282 (JP, A) Fukiyama Naoki, Mochizuki Higashi, Asami Taizo, "Weather Of Lightning Prediction Using Factors ", The Japan Society for Electrical Discharge and High Voltage Joint Research Workshop, Japan, The Institute of Electrical Engineers of Japan, October 26, 1995, ED-95-197-207, p. 23-32 (58) Fields surveyed (Int.Cl. 7 , DB name) G01W 1/00-1/18 G06F 17/60-19/00 JISST file (JOIS)
Claims (9)
タを解析して蜃気楼発生の有無を判別するための予測式
を導出し、これを用いて対象地域における予測気象デー
タを解析し、その解析値を過去の観測気象データの解析
値と照合することにより少なくとも7日先までの蜃気楼
発生の有無を判別することを特徴とする蜃気楼の発生予
測方法。1. A forecast formula for determining the presence or absence of a mirage is derived by analyzing past observed meteorological data related to the occurrence of a mirage, and the forecasted meteorological data in the target area is analyzed using this, and the analysis value is obtained. A mirage generation prediction method characterized by determining whether or not a mirage has occurred at least 7 days ahead by comparing with the analysis value of past observed meteorological data.
データを入力して蜃気楼発生の有無を解析・判別できる
よう学習制御した演算処理システムに、対象地域におけ
る現在以降の予測気象データを入力することにより、少
なくとも7日先までの蜃気楼の発生予測情報を出力させ
ることを特徴とする蜃気楼の発生予測方法。2. The forecasted meteorological data after the present in the target area is input to an arithmetic processing system that is learning-controlled so that past observed meteorological data related to the occurrence of a mirage can be input and analyzed / determined whether or not the mirage occurs. According to the above, the mirage occurrence prediction method is characterized by outputting the mirage occurrence prediction information for at least 7 days ahead.
対象地域の海面または地表に接する気層の気温変化率、
大気の安定度および天候に関与する気象項目データの組
合わせからなる蜃気楼の発生予測方法。3. The meteorological data according to claim 1 or 2,
Rate of temperature change in the air layer in contact with the sea surface or surface of the target area,
A method for predicting the occurrence of mirages, which consists of a combination of meteorological item data related to atmospheric stability and weather.
高気圧中心位置、サーマルリッジ存在エリアおよびサー
マルトラフ存在エリア、海面または地表温度と気温との
差、大気低層の気温減率、ショワルターのスタビリティ
・インデックス、天気、日照時間、雨量、風速および風
向である蜃気楼の発生予測方法。4. The meteorological item data according to claim 3 includes a mobile high pressure center position, a thermal ridge existing area and a thermal trough existing area, a difference between sea surface or surface temperature and air temperature, a temperature decrease rate in the lower atmosphere, and a shoaler's data. Stability index, weather, daylight hours, rainfall, wind speed and wind direction mirage occurrence prediction method.
面または地表温度と気温との差、最高気温の平年偏差、
大気低層の気温減率、ショワルターのスタビリティ・イ
ンデックス、天気、日照時間、最大風速生起時の風向お
よび冬季降雪量の狭域気象データより選ばれた少なくと
も2ケの気象項目データと、移動性高気圧中心位置、サ
ーマルリッジ存在エリアおよびサーマルトラフ存在エリ
アの広域気象データより選ばれた少なくとも2ケの気象
項目データの組合わせからなる蜃気楼の発生予測方法。5. The meteorological item data according to claim 3, wherein the difference between the sea surface or surface temperature and the temperature, the normal deviation of the maximum temperature,
At least two meteorological item data selected from low-temperature deceleration of the lower atmosphere, stability index of shoalter, weather, sunshine duration, wind direction at maximum wind speed and narrow range meteorological data of winter snowfall, and mobility A mirage generation prediction method comprising a combination of at least two meteorological item data selected from wide area meteorological data of a high pressure center position, a thermal ridge existing area, and a thermal trough existing area.
タを解析して蜃気楼発生の有無を判別するための予測式
を導出する演算装置と、この演算装置に過去の観測気象
データ及び対象地域における予測気象データを入力する
入力装置と、入力データを予測式で解析して得られた蜃
気楼発生の有無の判別結果を出力する出力装置とからな
る蜃気楼の発生予測装置。6. An arithmetic device for deriving a prediction formula for analyzing whether or not a mirage has occurred by analyzing past observed meteorological data relating to the occurrence of a mirage, and a past observation meteorological data and a prediction in a target area in this arithmetic device. A mirage generation prediction device comprising an input device for inputting meteorological data and an output device for outputting a determination result of presence or absence of a mirage obtained by analyzing the input data by a prediction formula.
の海面または地表に接する気層の気温変化率、大気の安
定度および天候に関与する気象項目データの組合わせか
らなる蜃気楼の発生予測装置。7. A mirage outbreak prediction in which the meteorological data according to claim 6 is a combination of a temperature change rate of an air layer contacting the sea surface or the surface of the target area, atmospheric stability, and meteorological item data related to weather. apparatus.
高気圧中心位置、サーマルリッジ存在エリアおよびサー
マルトラフ存在エリア、海面または地表温度と気温との
差、大気低層の気温減率、ショワルターのスタビリティ
・インデックス、天気、日照時間、雨量、風速および風
向である蜃気楼の発生予測装置。8. The meteorological item data set forth in claim 6 includes a mobile high pressure center position, a thermal ridge existing area and a thermal trough existing area, a difference between sea surface or surface temperature and air temperature, a temperature decrease rate of the lower atmosphere, and a shoaler value. Stability index, weather, sunshine hours, rainfall, wind speed and wind direction mirage occurrence prediction device.
面または地表温度と気温との差、最高気温の平年偏差、
大気低層の気温減率、ショワルターのスタビリティ・イ
ンデックス、天気、日照時間、最大風速生起時の風向お
よび冬季降雪量の狭域気象データより選ばれた少なくと
も2ケの気象項目データと、移動性高気圧中心位置、サ
ーマルリッジ存在エリアおよびサーマルトラフ存在エリ
アの広域気象データより選ばれた少なくとも2ケの気象
項目データの組合わせからなる蜃気楼の発生予測装置。9. The meteorological item data according to claim 6, wherein the difference between the sea surface temperature or the surface temperature and the temperature, the average deviation of the maximum temperature,
At least two meteorological item data selected from low-temperature deceleration of the lower atmosphere, stability index of shoalter, weather, sunshine duration, wind direction at maximum wind speed and narrow range meteorological data of winter snowfall, and mobility A mirage generation prediction device comprising a combination of at least two meteorological item data selected from wide area meteorological data of a high pressure center position, a thermal ridge existing area, and a thermal trough existing area.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31231395A JP3514564B2 (en) | 1995-11-30 | 1995-11-30 | Mirage occurrence prediction method and prediction device |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| JP31231395A JP3514564B2 (en) | 1995-11-30 | 1995-11-30 | Mirage occurrence prediction method and prediction device |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| JPH09152491A JPH09152491A (en) | 1997-06-10 |
| JP3514564B2 true JP3514564B2 (en) | 2004-03-31 |
Family
ID=18027750
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| JP31231395A Expired - Fee Related JP3514564B2 (en) | 1995-11-30 | 1995-11-30 | Mirage occurrence prediction method and prediction device |
Country Status (1)
| Country | Link |
|---|---|
| JP (1) | JP3514564B2 (en) |
Families Citing this family (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN104280791B (en) * | 2014-09-26 | 2017-03-22 | 宁波绮耘软件股份有限公司 | Meteorological data quality control and processing method and system |
| CN110631566B (en) * | 2019-10-10 | 2025-03-14 | 敦煌研究院 | A special device for measuring the height of desert phantom sea |
| CN111707012B (en) * | 2020-06-19 | 2024-04-09 | 中国地质大学(武汉) | Device for manufacturing mirage based on hot air prism |
-
1995
- 1995-11-30 JP JP31231395A patent/JP3514564B2/en not_active Expired - Fee Related
Non-Patent Citations (1)
| Title |
|---|
| 吹山直樹、望月東、浅見泰造,"気象因子を用いた落雷予測の検討",電気学会放電高電圧合同研究会資料,日本,社団法人電気学会,1995年10月26日,ED−95−197〜207,p.23−32 |
Also Published As
| Publication number | Publication date |
|---|---|
| JPH09152491A (en) | 1997-06-10 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Haiden et al. | Evaluation of ECMWF forecasts, including the 2018 upgrade | |
| Yang et al. | Accuracy of NWS 8" standard nonrecording precipitation gauge: Results and application of WMO intercomparison | |
| Parthasarathy et al. | Prediction of all‐India summer monsoon rainfall with regional and large‐scale parameters | |
| US6535817B1 (en) | Methods, systems and computer program products for generating weather forecasts from a multi-model superensemble | |
| Bedia et al. | Sensitivity of fire weather index to different reanalysis products in the Iberian Peninsula | |
| Weckwerth et al. | Radar refractivity retrieval: Validation and application to short-term forecasting | |
| CN106772697B (en) | Sea of clouds natural landscape forecasting procedure and system | |
| CN115907201A (en) | A short-term rainfall prediction method, device and system based on a bidirectional LSTM network | |
| Hdidou et al. | Impact of the variational assimilation of ground-based GNSS zenith total delay into AROME-Morocco model | |
| CN117172554A (en) | Ice disaster risk prediction methods, devices, equipment and storage media | |
| Yang et al. | High-resolution GEM-LAM application in marine fog prediction: Evaluation and diagnosis | |
| JP3514564B2 (en) | Mirage occurrence prediction method and prediction device | |
| CN116451088A (en) | Preferred station substituting method based on multi-element feature similarity and geographic region clustering | |
| Gringorten | Methods of objective weather forecasting | |
| de Souza et al. | Impact of radar data assimilation on the simulation of a heavy rainfall event over manaus in the Central Amazon | |
| Kumar et al. | Reflectivity-Rain Rate Relationship for Orographic Rainfall at Mahabaleshwar Over the Indian Western Ghats: AK Srivastava et al. | |
| McCLAIN | On the relation of satellite viewed cloud conditions to vertically integrated moisture fields | |
| Müller et al. | Avalanche Problem Solver (APS)–a decision support system for forecasters (part 1) | |
| CN114723120B (en) | Near-ground wind speed forecasting method and device | |
| Tracy et al. | Mesoscale Influences of Land Use, Topography, Antecedent Rainfall, and Atmospheric Conditions on Summertime Convective Storm Initiation under Weak Synoptic-Scale Forcing | |
| Rahmat et al. | Implementation of Bayesian model averaging on the weather data forecasting applications utilizing open weather map | |
| Lewis | Forecasting advective sea fog with the use of classification and regression tree analyses for Kunsan Air Base | |
| Ghanbari et al. | Spatial analysis of fars province climatic drought from the period 1990 to 2014 | |
| Liu et al. | Application Method of Multi-model Fusion for Heavy Rainfall Forecast | |
| Null | A climatology of San Francisco rainfall, 1849-1991 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20040113 |
|
| R250 | Receipt of annual fees |
Free format text: JAPANESE INTERMEDIATE CODE: R250 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090123 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20090123 Year of fee payment: 5 |
|
| FPAY | Renewal fee payment (event date is renewal date of database) |
Free format text: PAYMENT UNTIL: 20100123 Year of fee payment: 6 |
|
| LAPS | Cancellation because of no payment of annual fees |