JP4511979B2 - Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program - Google Patents

Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program Download PDF

Info

Publication number
JP4511979B2
JP4511979B2 JP2005062708A JP2005062708A JP4511979B2 JP 4511979 B2 JP4511979 B2 JP 4511979B2 JP 2005062708 A JP2005062708 A JP 2005062708A JP 2005062708 A JP2005062708 A JP 2005062708A JP 4511979 B2 JP4511979 B2 JP 4511979B2
Authority
JP
Japan
Prior art keywords
rainfall
attenuation
distribution
probability
rain
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
Application number
JP2005062708A
Other languages
Japanese (ja)
Other versions
JP2006246375A (en
Inventor
史明 峯松
陽一 鈴木
和義 正源
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.)
Japan Broadcasting Corp
Original Assignee
Japan Broadcasting 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 Japan Broadcasting Corp filed Critical Japan Broadcasting Corp
Priority to JP2005062708A priority Critical patent/JP4511979B2/en
Publication of JP2006246375A publication Critical patent/JP2006246375A/en
Application granted granted Critical
Publication of JP4511979B2 publication Critical patent/JP4511979B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Monitoring And Testing Of Transmission In General (AREA)
  • Mobile Radio Communication Systems (AREA)
  • Radio Relay Systems (AREA)

Description

本発明は、降雨による電波の減衰を補償する技術に係り、特に降雨による電波の減衰を補償するために、電波の減衰量を算出する技術に関する。   The present invention relates to a technique for compensating for radio wave attenuation due to rain, and more particularly to a technique for calculating radio wave attenuation in order to compensate for radio wave attenuation due to rain.

一般に、衛星放送や、人工衛星(以下、衛星という)から送信される放送波(以下、電波という)を用いた通信では、送受信される電波の周波数帯が通常数GHz以上となるため、降雨による電波の減衰(以下、降雨減衰という)が生じる。この降雨減衰によって、受信機や送受信局(以下、単に送受信局という)で受信される番組や通信が一時的に途絶えることがあった。特に、周波数が大きいほど降雨減衰が大きくなるため、将来利用可能となる20GHz程度の準ミリ波以上の周波数帯域では、放送・通信サービスに更に大きな影響を及ぼす。   In general, in communications using satellite broadcasts and broadcast waves (hereinafter referred to as radio waves) transmitted from artificial satellites (hereinafter referred to as satellites), the frequency band of the transmitted and received radio waves is usually several GHz or more, and therefore due to rain. Radio wave attenuation (hereinafter referred to as rain attenuation) occurs. Due to this rain attenuation, programs and communications received by receivers and transmitting / receiving stations (hereinafter simply referred to as transmitting / receiving stations) may be temporarily interrupted. In particular, since the rain attenuation increases as the frequency increases, the frequency band of quasi-millimeter waves of about 20 GHz or more that can be used in the future has a greater influence on broadcasting / communication services.

そのため、従来、降雨減衰を統計的に推定し、送受信局から電波を送信する送信電力を予め高めに設定して、送信電力に一定の降雨マージンを与えることで、降雨減衰の補償(降雨減衰補償)を行う機能を有する衛星が開示されている(特許文献1参照)。この衛星は、強い降雨の発生している地域が比較的狭い範囲であることに着目し、強い降雨が発生している期間にその地域だけに、より大きな送信電力を設定して電波を送信することで、衛星の規模の増大を抑えつつ、効果的な降雨減衰補償を行うものである。   Therefore, conventionally, the rain attenuation is statistically estimated, the transmission power for transmitting radio waves from the transmitting and receiving stations is set high in advance, and a certain rain margin is given to the transmission power to compensate for the rain attenuation (rain attenuation compensation). ) Has been disclosed (see Patent Document 1). This satellite pays attention to the fact that the area where heavy rainfall is occurring is a relatively narrow area, and during the period when heavy rain occurs, it sets a larger transmission power and transmits radio waves only to that area Thus, effective rain attenuation compensation is performed while suppressing an increase in the scale of the satellite.

また、降雨減衰を統計的に推定する様々な手法が提案され、その大半は、1年間である閾値を超える降雨減衰が発生する確率を推定するものである。更に、複数の観測地点について1分間ごとの降雨強度(1分間降雨強度)を観測し、この1分間降雨強度の確率分布に基づいて、観測された降雨による電波の減衰量を推定する方法が開示されている(非特許文献1参照)。この方法では、所定の観測地点における10年間の平均と、特に降雨減衰の大きい月及び年との降雨強度分布及びこの降雨による減衰量について検討している。   Various methods for statistically estimating rainfall attenuation have been proposed, most of which estimate the probability of occurrence of rainfall attenuation exceeding a threshold of one year. Furthermore, a method of observing the rainfall intensity per minute (one minute rainfall intensity) at a plurality of observation points and estimating the attenuation of radio waves due to the observed rainfall based on the probability distribution of the one minute rainfall intensity is disclosed. (See Non-Patent Document 1). In this method, the average of 10 years at a predetermined observation point, the rainfall intensity distribution especially in the month and year when the rain attenuation is large, and the attenuation due to this rain are examined.

更に、近年、気象予報技術の発達により、この先に発生する降雨エリアを高い地理的分解能で予測することが可能になっている。例えば、日本国内全域を、5km四方の正方形の単位で分割し、分割された区分(以下、メッシュという)ごとの1時間あたりの雨量(1時間雨量)の予測値が、気象庁から降水短時間予報として発表されている。
特開2004−193963号公報(段落0009〜0042) 森田 和夫、「衛星通信回線における伝搬特性の推定法(準ミリ〜ミリ波帯の場合)」、電気通信研究所研究実用化報告、日本電信電話研究開発本部、1979年8月、第28巻、第8号、p.1661〜1676
Furthermore, in recent years, with the development of weather forecasting technology, it is possible to predict the rainfall area that will occur in the future with high geographical resolution. For example, the entire Japan area is divided into units of 5 km squares, and the predicted value of rainfall per hour (rainfall per hour) for each divided section (hereinafter referred to as mesh) is predicted by the Japan Meteorological Agency for a short-term precipitation forecast. It has been announced as.
JP 2004-193963 A (paragraphs 0009 to 0042) Kazuo Morita, “Estimation of propagation characteristics in satellite communication links (in the case of quasi-millimeter to millimeter-wave band)”, Research report on practical application of Research Institute of Electrical Communication, Nippon Telegraph and Telephone Research and Development Headquarters, August 1979, Vol. 28, No. 8, p. 1661 to 1676

しかしながら、特許文献1に記載の衛星を用いて降雨減衰補償を行うためには、強い降雨が発生している地域と、その地域における降雨減衰量を予測する必要があるが、従来の年間を通じた発生確率によって降雨減衰量を推定する手法では、衛星が電波を送信する時点においてどの地域にどれだけの大きさの降雨マージンを与える必要があるかを予測することはできない。   However, in order to perform rainfall attenuation compensation using the satellite described in Patent Document 1, it is necessary to predict a region where strong rainfall occurs and the amount of rainfall attenuation in that region. In the method of estimating the rain attenuation amount based on the probability of occurrence, it is impossible to predict how much rain margin needs to be given to which area when the satellite transmits radio waves.

また、降水短時間予報の1時間雨量の予測値は、降雨マージンを与える必要がある降雨減衰が発生する確率の高い地域を予測することには有益であるが、この1時間の間にどの程度の確率で閾値を超える降雨減衰量となるかは分からなかった。つまり、例えば、1時間の間、降雨強度が一定であると仮定して降雨による減衰量を算出し、降雨マージンを決定することはできても、実際にはこの1時間の間において降雨強度は一定でなく変動するため、実際の降雨減衰を効果的に補償することができない。   In addition, the predicted value of 1-hour rainfall in the short-term precipitation forecast is useful for predicting an area where there is a high probability of rain attenuation that needs to be given a rain margin. It was not known if the rain attenuation amount exceeded the threshold with the probability of. In other words, for example, it is assumed that the rainfall intensity is constant for one hour, and the amount of rain attenuation can be calculated and the rain margin can be determined. Because it is not constant and fluctuates, it cannot effectively compensate for actual rain attenuation.

また、非特許文献1の方法では、複数の観測地点について1分間降雨強度を測定する必要があり、また、例えば1時間のような所定時間内における1分間降雨強度の確率分布は、所定時間ごとの雨量の大きさによって異なる。そのため、非特許文献1に記載されるように、限られた地点について、1分間降雨強度が観測された降雨による電波の減衰量を、観測後に観測値から推定することは可能であるが、例えば、1時間ごと雨量の予測値に基づいて降雨による減衰量を日本全域について解析する場合には、予め複数の観測地点において、様々な1時間ごとの雨量について、1分間降雨強度を観測しておく必要が生じる。そして、降雨強度の確率分布は、観測地点によっても異なるため、誤差を小さくするためには膨大な観測地点について観測データを長期にわたって収集する必要があるが、このような観測データを収集するためには多大なコストと時間がかかるため、入手が困難である。   Further, in the method of Non-Patent Document 1, it is necessary to measure the rainfall intensity for one minute at a plurality of observation points, and the probability distribution of the rainfall intensity for one minute within a predetermined time such as one hour is, for example, every predetermined time. Depending on the amount of rainfall. Therefore, as described in Non-Patent Document 1, it is possible to estimate the amount of radio wave attenuation due to rainfall in which rainfall intensity is observed for one minute from a observed value after observation. When analyzing the amount of rain attenuation for the entire region of Japan based on the predicted value of rainfall every hour, the rainfall intensity for 1 minute is observed in advance for various hourly rainfalls at a plurality of observation points in advance. Need arises. And since the probability distribution of rainfall intensity varies from observation point to observation point, it is necessary to collect observation data for a huge amount of observation points over a long period of time in order to reduce the error. In order to collect such observation data, Is expensive and time consuming and difficult to obtain.

そのため、メッシュごとの1時間雨量の予測値のみでは、降雨による減衰量の予測値を算出することができず、降雨減衰を補償するために、どの地域に対してどの程度の降雨マージンを与えればよいかを予測することができないという問題があった。   Therefore, it is not possible to calculate the predicted value of attenuation due to rainfall only with the predicted value of rainfall for one hour for each mesh, and what amount of rainfall margin should be given to which region in order to compensate for rain attenuation There was a problem that it could not be predicted.

本発明は、前記従来技術の問題を解決するために成されたもので、時間とともに変化する降雨による電波の減衰量を解析することができる降雨減衰解析装置、区分別減衰量解析装置、降雨減衰解析方法及び区分別減衰量解析プログラムを提供することを目的とする。   The present invention has been made to solve the above-described problems of the prior art, and is a rain attenuation analysis apparatus, a divisional attenuation analysis apparatus, a rain attenuation, which can analyze the attenuation of radio waves due to rainfall that varies with time. It is an object to provide an analysis method and a category-specific attenuation analysis program.

前記課題を解決するため、請求項1に記載の降雨減衰解析装置は、ある観測地点について観測された10分間ごとの降雨強度のうち、間における雨量が所定値となる前記降雨強度を示す地点別雨量データに基づいて、前記観測地点における前記時間あたりの雨量の降雨による電波の減衰量の確率分布を解析する降雨減衰解析装置であって、降雨強度分布解析手段と、減衰解析手段とを備える構成とした。 To solve the above problems, rain attenuation analysis apparatus according to claim 1, of the rain rate every 10 minutes observed for an observation point, the rainfall intensity rainfall definitive between 1 hour becomes a predetermined value based on the locational rainfall data indicating, a rain attenuation analyzer for analyzing a probability distribution of a radio wave attenuation due to rainfall rainfall per the hour in the observation point, the rainfall intensity distribution analysis unit, attenuation analyzing means It was set as the structure provided with.

かかる構成によれば、降雨減衰解析装置は、降雨強度分布解析手段によって、地点別雨量データによって示される降雨強度の確率分布を対数正規分布で近似して、この確率分布を示す降雨強度分布データを生成する。また、降雨減衰解析装置は、減衰解析手段によって、降雨強度分布データである10分間降雨強度の確率分布の平均値及び標準偏差を、1分間降雨強度の確率分布の平均値及び標準偏差とみなして、1時間あたりの雨量の降雨における減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する。 According to such a configuration, the rain attenuation analysis apparatus approximates the probability distribution of the rainfall intensity indicated by the point-by-point rainfall data by the log normal distribution by the rain intensity distribution analysis means, and obtains the rain intensity distribution data indicating this probability distribution. Generate. Further, the rain attenuation analysis apparatus regards the average value and standard deviation of the probability distribution of 10-minute rainfall intensity, which is the rainfall intensity distribution data , as the average value and standard deviation of the probability distribution of 1-minute rainfall intensity by the attenuation analysis means. Attenuation amount distribution data indicating the probability distribution is generated by analyzing the probability distribution of the attenuation amount during one hour of rainfall.

ここで、本発明者らは、10分間ごとの降雨強度の実測値と、1分間ごとの降雨強度の実測値との確率分布について検討を行った。その結果、10分間ごとの降雨強度の実測値及び1分間ごとの降雨強度の確率分布は、観測地点ごとの所定時間あたりの雨量に依存し、観測地点及び所定時間あたりの雨量が同じ場合には、どちらの確率分布も対数正規分布に近い分布となるとともに、2つの分布がほぼ一致することを見出した。そのため、降雨減衰解析装置は、10分間ごとの降雨強度を示す地点別雨量データに基づいて、観測地点ごとの所定時間あたりの雨量が所定値となる降雨の10分間ごとの降雨強度の確率分布を対数正規分布で近似し、この確率分布を1分間ごとの降雨強度の確率分布とみなすことで、所定時間内における減衰量の確率分布を解析することが可能になる。   Here, the present inventors examined the probability distribution between the measured value of the rainfall intensity every 10 minutes and the measured value of the rainfall intensity every minute. As a result, the measured value of rainfall intensity every 10 minutes and the probability distribution of rainfall intensity every minute depend on the rainfall per predetermined time for each observation point, and when the rainfall per observation point and the predetermined time are the same Both probability distributions were found to be close to lognormal distributions, and the two distributions almost matched. Therefore, the rain attenuation analysis device calculates the probability distribution of the rainfall intensity every 10 minutes of the rainfall where the rainfall per predetermined time at each observation point becomes a predetermined value based on the rainfall data by point indicating the rainfall intensity every 10 minutes. By approximating with a lognormal distribution and regarding this probability distribution as a probability distribution of rainfall intensity every minute, it becomes possible to analyze the probability distribution of attenuation within a predetermined time.

また、請求項2に記載の区分別減衰量解析装置は、ある観測地点について観測された10分間ごとの降雨強度のうち、所定時間内における雨量が所定値となる前記降雨強度を示す地点別雨量データによって示される前記降雨強度の確率分布を対数正規分布で近似して、前記確率分布を示す降雨強度分布データを生成する降雨強度分布解析手段と、この降雨強度分布解析手段によって生成された降雨強度分布データに基づいて、前記所定時間あたりの雨量の降雨における前記減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する減衰解析手段とを備える降雨減衰解析装置によって生成された、ある観測地点において所定時間あたりの雨量が所定値となる降雨による電波の減衰量の確率分布を示す減衰量分布データに基づいて、前記減衰量の解析対象となる地域内の複数の区分の各々について降雨による前記電波の減衰量を解析する区分別減衰量解析装置であって、減衰量分布データ記憶手段と、雨量データ入力手段と、減衰量分布データ読み出し手段と、電波減衰解析手段とを備える構成とした。 Further, the classified attenuation amount analysis device according to claim 2 is a point-by-point rainfall amount indicating the rain intensity at a predetermined value within a predetermined time period, among the rainfall intensity values every 10 minutes observed at a certain observation point. A rainfall intensity distribution analyzing means for generating a rainfall intensity distribution data indicating the probability distribution by approximating the probability distribution of the rain intensity indicated by the data with a logarithmic normal distribution, and a rainfall intensity generated by the rainfall intensity distribution analyzing means Based on the distribution data, generated by a rain attenuation analysis device comprising attenuation analysis means for analyzing the probability distribution of the attenuation amount in the rainfall of the rainfall per predetermined time and generating attenuation distribution data indicating the probability distribution. In addition, based on attenuation distribution data indicating the probability distribution of the attenuation of radio waves due to rainfall with a predetermined amount of rainfall per predetermined time at a certain observation point. An attenuation analysis device for each category for analyzing the attenuation of the radio wave due to rainfall for each of a plurality of categories in the area to be analyzed for the attenuation, comprising attenuation distribution data storage means, and rainfall data input Means, attenuation amount distribution data reading means, and radio wave attenuation analysis means.

かかる構成によれば、区分別減衰量解析装置は、減衰量分布データ記憶手段によって、減衰量分布データを観測地点及び雨量と対応付けて、複数の観測地点及び複数の雨量について記憶する。また、雨量データ入力手段によって、区分ごとの所定時間あたりの雨量を示す雨量データを入力し、減衰量分布データ読み出し手段によって、減衰量の確率分布を示す減衰量分布データを記憶する減衰量分布データ記憶手段から、区分ごとに、当該区分に予め対応付けられた観測地点と、雨量データによって示される雨量とに対応する減衰量分布データを読み出す。そして、電波減衰解析手段によって、減衰量分布データ読み出し手段で読み出された減衰量分布データによって示される減衰量の確率分布と、この確率分布において減衰量がある値より大きな値をとる確率である超過確率の許容値を示す許容値データとに基づいて、区分ごとに超過確率がこの許容値となる減衰量を解析する。   According to such a configuration, the attenuation analysis apparatus for each category stores the attenuation distribution data by associating the attenuation distribution data with the observation points and the rainfall, and stores the plurality of observation points and the plurality of rainfalls. Moreover, the rain amount data input means inputs rain amount data indicating the rainfall per predetermined time for each section, and the attenuation amount distribution data reading means stores the attenuation amount distribution data indicating the attenuation probability distribution. For each section, attenuation distribution data corresponding to the observation point associated with the section in advance and the rainfall indicated by the rainfall data are read from the storage means. Then, the probability distribution of the attenuation indicated by the attenuation distribution data read by the attenuation distribution data reading means by the radio wave attenuation analysis means, and the probability that the attenuation takes a value larger than a certain value in this probability distribution. Based on the tolerance value data indicating the tolerance value of the excess probability, the attenuation amount at which the excess probability becomes the tolerance value is analyzed for each section.

これによって、区分別減衰量解析装置は、区分ごとに、雨量データによって示される所定時間あたりの雨量に対応する減衰量分布データを減衰量分布データ記憶手段から読み出し、超過確率が許容値データによって示される許容値となる減衰量を解析することができる。   As a result, the category-by-category attenuation amount analysis apparatus reads the attenuation amount distribution data corresponding to the rainfall per predetermined time indicated by the rainfall amount data for each category from the attenuation amount distribution data storage means, and the excess probability is indicated by the allowable value data. It is possible to analyze the amount of attenuation that is an allowable value.

更に、請求項3に記載の区分別減衰量解析装置は、ある観測地点について観測された所定時間内における雨量が所定値となる降雨の10分間ごとの降雨強度の確率分布を対数正規分布で近似して生成された、当該確率分布を示す降雨強度分布データを、前記観測地点及び前記雨量と対応付けて、複数の前記観測地点及び複数の前記雨量について記憶する降雨強度分布データ記憶装置から前記降雨強度分布データを読み出して、減衰量の解析対象となる地域内の複数の区分の各々について降雨による電波の減衰量を解析する区分別減衰量解析装置であって、雨量データ入力手段と、降雨強度分布データ読み出し手段と、減衰解析手段と、電波減衰解析手段とを備える構成とした。   Furthermore, the attenuation analysis device according to claim 3 approximates the probability distribution of the rainfall intensity every 10 minutes of the rain when the rainfall within a predetermined time observed at a certain observation point becomes a predetermined value by a lognormal distribution. The rainfall intensity distribution data indicating the probability distribution generated in this manner is associated with the observation point and the rainfall amount, and stored for the plurality of observation points and the plurality of rainfall amounts. A category-specific attenuation analysis device that reads out intensity distribution data and analyzes the attenuation of radio waves due to rainfall for each of a plurality of categories in the area subject to attenuation analysis. A distribution data reading unit, an attenuation analysis unit, and a radio wave attenuation analysis unit are provided.

かかる構成によれば、区分別減衰量解析装置は、雨量データ入力手段によって、区分ごとの所定時間あたりの雨量を示す雨量データを入力し、降雨強度分布データ読み出し手段によって、10分間ごとの降雨強度の確率分布を示す降雨強度分布データを記憶する降雨強度分布データ記憶装置から、区分ごとに、各々の区分に予め対応付けられた観測地点と、雨量データによって示される雨量とに対応する降雨強度分布データを読み出す。また、減衰解析手段によって、降雨強度分布データ読み出し手段で読み出された降雨強度分布データに基づいて、所定時間あたりの雨量の降雨による電波の減衰量の確率分布を解析し、電波減衰解析手段によって、この確率分布と、確率分布において減衰量がある値より大きな値をとる確率である超過確率の許容値を示す許容値データとに基づいて区分ごとに超過確率が許容値となる減衰量を解析する。   According to this configuration, the category-by-category attenuation amount analysis apparatus inputs the rainfall data indicating the rainfall per predetermined time for each category by the rainfall data input unit, and the rainfall intensity every 10 minutes by the rainfall intensity distribution data reading unit. From the rainfall intensity distribution data storage device that stores the rainfall intensity distribution data indicating the probability distribution of each, the rainfall intensity distribution corresponding to the observation point previously associated with each category and the rainfall indicated by the rainfall data for each category Read data. Further, the attenuation analysis means analyzes the probability distribution of the attenuation of the radio wave due to the rainfall per predetermined time based on the rainfall intensity distribution data read by the rainfall intensity distribution data reading means, and the radio wave attenuation analysis means Based on this probability distribution and the tolerance value data indicating the tolerance value of the excess probability, which is the probability that the attenuation amount takes a value larger than a certain value in the probability distribution, the attenuation amount in which the excess probability becomes the tolerance value for each category is analyzed. To do.

これによって、区分別減衰量解析装置は、区分ごとに、雨量データによって示される所定時間あたりの雨量に対応する、10分間ごとの降雨強度の確率分布を示す降雨強度分布データを降雨強度分布データ記憶装置から読み出し、10分間ごとの降雨強度の確率分布を1分間ごとの降雨強度の確率分布とみなして減衰量分布データを生成して、超過確率が許容値データによって示される許容値となる減衰量を解析することができる。   Accordingly, the category-specific attenuation amount analysis apparatus stores the rainfall intensity distribution data indicating the probability distribution of the rainfall intensity every 10 minutes corresponding to the rainfall per predetermined time indicated by the rainfall data for each category. Attenuation amount in which the probability distribution of the excess intensity becomes an allowable value indicated by the allowable value data by generating the attenuation amount distribution data by considering the probability distribution of the rainfall intensity every 10 minutes as the probability distribution of the rainfall intensity every minute. Can be analyzed.

また、請求項4に記載の降雨減衰解析方法は、ある観測地点について観測された10分間ごとの降雨強度のうち、間における雨量が所定値となる前記降雨強度を示す地点別雨量データに基づいて、前記観測地点における前記時間あたりの雨量の降雨による電波の減衰量の確率分布を解析する降雨減衰解析方法であって、降雨強度分布解析ステップと、減衰解析ステップとを含むことを特徴とする。 Also, rain attenuation analysis method according to claim 4, of the rain rate every 10 minutes observed for an observation point, locational rainfall data indicating the rain intensity rainfall definitive between 1 hour becomes a predetermined value based on the a rain attenuation analyzing method for analyzing a probability distribution of a radio wave attenuation due to rain rainfall per hour in the observation point, the rainfall intensity distribution analysis step, that comprises a damping analysis step Features.

この方法によれば、降雨強度分布解析ステップによって、地点別雨量データによって示される降雨強度の確率分布を対数正規分布で近似して、この確率分布を示す降雨強度分布データを生成する。続いて、減衰解析ステップによって、降雨強度分布データである10分間降雨強度の確率分布の平均値及び標準偏差を、1分間降雨強度の確率分布の平均値及び標準偏差とみなして、1時間あたりの雨量の降雨における減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する。 According to this method, the rainfall intensity distribution analysis step generates the rainfall intensity distribution data indicating the probability distribution by approximating the probability distribution of the rainfall intensity indicated by the point-by-point rainfall data with the lognormal distribution. Subsequently, in the attenuation analysis step, the average value and standard deviation of the probability distribution of 10-minute rainfall intensity, which is the rainfall intensity distribution data , are regarded as the average value and standard deviation of the probability distribution of 1-minute rainfall intensity per hour. Attenuation distribution data representing the probability distribution is generated by analyzing the probability distribution of the attenuation in the rainfall.

これによって、10分間ごとの降雨強度を示す地点別雨量データに基づいて、観測地点ごとの所定時間あたりの所定の雨量の降雨の10分間ごとの降雨強度の確率分布を対数正規分布で近似することで、所定時間内における減衰量の確率分布を解析することができる。   Thus, based on the point-by-point rainfall data indicating the rainfall intensity every 10 minutes, the probability distribution of the rainfall intensity every 10 minutes of the rainfall of the predetermined rainfall per predetermined time for each observation point is approximated by a lognormal distribution. Thus, the probability distribution of the attenuation amount within a predetermined time can be analyzed.

また、請求項5に記載の区分別減衰量解析プログラムは、ある観測地点について観測された10分間ごとの降雨強度のうち、所定時間内における雨量が所定値となる前記降雨強度を示す地点別雨量データによって示される前記降雨強度の確率分布を対数正規分布で近似して、前記確率分布を示す降雨強度分布データを生成する降雨強度分布解析手段と、この降雨強度分布解析手段によって生成された降雨強度分布データに基づいて、前記所定時間あたりの雨量の降雨における前記減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する減衰解析手段とを備える降雨減衰解析装置によって生成された、ある観測地点において所定時間あたりの雨量が所定値となる降雨による電波の減衰量の確率分布を示す減衰量分布データを、前記観測地点及び前記雨量と対応付けて、複数の前記観測地点及び複数の前記雨量について記憶する減衰量分布データ記憶装置から前記減衰量分布データを読み出して、前記減衰量の解析対象となる地域内の複数の区分の各々について降雨による前記電波の減衰量を解析するためにコンピュータを、雨量データ入力手段、減衰量分布データ読み出し手段、電波減衰解析手段として機能させることとした。 Further, according to the category-specific attenuation amount analysis program, the point-by-point rain amount indicating the rain intensity at a predetermined value within a predetermined time out of the 10-minute rainfall intensity observed at a certain observation point. A rainfall intensity distribution analyzing means for generating a rainfall intensity distribution data indicating the probability distribution by approximating the probability distribution of the rain intensity indicated by the data with a logarithmic normal distribution, and a rainfall intensity generated by the rainfall intensity distribution analyzing means Based on the distribution data, generated by a rain attenuation analysis device comprising attenuation analysis means for analyzing the probability distribution of the attenuation amount in the rainfall of the rainfall per predetermined time and generating attenuation distribution data indicating the probability distribution. Attenuation distribution data showing the probability distribution of the attenuation of radio waves due to rainfall with a certain amount of rainfall per given time at a certain observation point The attenuation distribution data is read out from the attenuation distribution data storage device that stores the plurality of observation points and the plurality of rainfalls in association with the observation point and the rainfall, and the area to be analyzed for the attenuation amount In order to analyze the attenuation amount of the radio wave due to rain for each of the plurality of sections, the computer is caused to function as a rain data input means, an attenuation distribution data reading means, and a radio attenuation analysis means.

かかる構成によれば、区分別減衰量解析プログラムは、雨量データ入力手段によって、区分ごとの所定時間あたりの雨量を示す雨量データを入力し、減衰量分布データ読み出し手段によって、減衰量の確率分布を示す減衰量分布データを記憶する減衰量分布データ記憶装置から、区分ごとに、当該区分に予め対応付けられた観測地点と、雨量データによって示される雨量とに対応する減衰量分布データを読み出す。また、電波減衰解析手段によって、減衰量分布データ読み出し手段で読み出された減衰量分布データによって示される減衰量の確率分布と、この確率分布において減衰量がある値より大きな値をとる確率である超過確率の許容値を示す許容値データとに基づいて、区分ごとに超過確率が許容値となる減衰量を解析する。   According to this configuration, the category-specific attenuation amount analysis program inputs the rainfall amount data indicating the rainfall amount per predetermined time for each category by the rainfall amount data input unit, and the attenuation amount distribution data reading unit calculates the probability distribution of the attenuation amount. From the attenuation distribution data storage device that stores the attenuation distribution data shown, the attenuation distribution data corresponding to the observation point previously associated with the section and the rain indicated by the rain data is read for each section. In addition, the probability distribution of the attenuation indicated by the attenuation distribution data read by the attenuation distribution data reading means by the radio wave attenuation analysis means, and the probability that the attenuation is greater than a certain value in this probability distribution. Based on the tolerance data indicating the tolerance value of the excess probability, the attenuation amount at which the excess probability becomes the tolerance value is analyzed for each category.

これによって、区分別減衰量解析プログラムは、区分ごとに、雨量データによって示される所定時間あたりの雨量に対応する減衰量分布データを減衰量分布データ記憶装置から読み出し、超過確率が許容値データによって示される許容値となる減衰量を解析することができる。   As a result, the attenuation analysis program for each category reads the attenuation distribution data corresponding to the rainfall per predetermined time indicated by the rainfall data for each category from the attenuation distribution data storage device, and the excess probability is indicated by the allowable value data. It is possible to analyze the amount of attenuation that is an allowable value.

本発明に係る降雨減衰解析装置、区分別減衰量解析装置、降雨減衰解析方法及び区分別減衰量解析プログラムでは、以下のような優れた効果を奏する。   The rain attenuation analysis device, the category attenuation analysis device, the rainfall attenuation analysis method, and the category attenuation analysis program according to the present invention have the following excellent effects.

請求項1又は請求項4に記載の発明によれば、10分間ごとの降雨強度に基づいて、所定時間あたりの雨量が所定値となる降雨による電波の減衰量の確率分布を解析できる。そして、例えば、気象庁の地域気象観測システム(AMeDAS;Automated Meteorological Data Acquisition System)の各々の雨量計によって観測された全国約1300箇所の観測地点の観測データのように、10分間ごとの雨量を観測した観測データは容易に入手できるため、この観測データを用いて電波の減衰量の確率分布を解析することが可能になる。   According to the invention described in claim 1 or claim 4, it is possible to analyze the probability distribution of the attenuation amount of the radio wave due to rainfall in which the rainfall per predetermined time becomes a predetermined value based on the rainfall intensity every 10 minutes. And, for example, the rainfall for every 10 minutes was observed like the observation data of about 1300 observation points nationwide observed by each rain gauge of the Meteorological Agency's regional meteorological observation system (AMeDAS; Automated Meteorological Data Acquisition System) Since the observation data can be easily obtained, it is possible to analyze the probability distribution of the attenuation of the radio wave using this observation data.

請求項2、請求項3又は請求項5に記載の発明によれば、10分間ごとの降雨強度から解析された減衰量の確率分布に基づいて、区分ごとに、雨量データによって示される所定時間あたりの雨量の降雨による、超過確率が許容値データによって示される許容値となる減衰量を解析することができる。そして、例えば、この減衰量を補償する降雨マージンを送信電力に与えて電波を衛星から送信すると、放送や通信が降雨によって遮断される確率は、この許容値によって示される確率となる。この許容値を十分小さい値に設定して解析された減衰量に基づいて、降雨マージンを与えて電波を送信することで、時間とともに変化する降雨強度に応じた降雨減衰補償を行うことができる。   According to the invention of claim 2, claim 3 or claim 5, based on the probability distribution of attenuation analyzed from the rainfall intensity every 10 minutes, for each predetermined time indicated by the rainfall data for each section. It is possible to analyze an attenuation amount at which the excess probability becomes an allowable value indicated by the allowable value data due to the rainfall of the current amount of rain. For example, if a rain margin that compensates for this attenuation is applied to the transmission power and radio waves are transmitted from the satellite, the probability that the broadcast or communication is interrupted by the rain is the probability indicated by this allowable value. By setting the allowable value to a sufficiently small value and transmitting a radio wave with a rain margin based on the analyzed attenuation amount, it is possible to perform rain attenuation compensation according to the rainfall intensity that changes with time.

また、区分ごとに減衰量を解析するので、降雨の集中によって減衰量が特に大きくなる範囲を特定でき、例えば、衛星放送の電波をフェーズドアレーアンテナのような、電波を送受信する地域の特定の範囲のみに高い電力で電波を送信できるアンテナによって送信することで、減衰量が特に大きくなる範囲に対して局所的に電波の強度を高くすることができるため、効率よく降雨減衰の補償ができる。   In addition, since the attenuation is analyzed for each category, it is possible to specify the range where the attenuation is particularly large due to the concentration of rainfall.For example, the specific range of the area where radio waves are transmitted and received, such as satellite broadcast waves, such as phased array antennas By transmitting with an antenna that can transmit radio waves with high power only, the intensity of radio waves can be locally increased with respect to a particularly large range of attenuation, so that rain attenuation can be compensated efficiently.

以下、本発明の実施の形態について図面を参照して説明する。
[降雨マージン算出装置による降雨減衰補償の概要]
まず、図1及び図2を参照して、本発明の実施の形態である降雨マージン算出装置1(1A)による降雨減衰補償の概要について説明する。図1は、本発明における降雨マージン算出装置による降雨減衰補償の概要を模式的に示す説明図である。図2は、本発明における降雨マージン算出装置によって算出された降雨マージンを与えて、降雨減衰補償を行う例を示す説明図、(a)は、フェーズドアレーアンテナによって、電波を送受信する地域に降雨マージンを与えた電波を送信し、降雨の集中する一部の範囲に特に高い強度の電波を送信する例を示す模式図、(b)は、(a)において送信される電波の強さを表すグラフである。
Hereinafter, embodiments of the present invention will be described with reference to the drawings.
[Outline of Rain Attenuation Compensation by Rain Margin Calculation Device]
First, with reference to FIG.1 and FIG.2, the outline | summary of the rain attenuation compensation by the rain margin calculation apparatus 1 (1A) which is embodiment of this invention is demonstrated. FIG. 1 is an explanatory view schematically showing an outline of rain attenuation compensation by a rain margin calculating apparatus according to the present invention. FIG. 2 is an explanatory diagram showing an example in which a rain margin is calculated by giving a rain margin calculated by the rain margin calculating apparatus according to the present invention, and (a) shows a rain margin in an area where radio waves are transmitted and received by a phased array antenna. The schematic diagram which shows the example which transmits the radio wave which gave the wave, and transmits the especially high intensity | strength radio wave to the one part range where rainfall concentrates, (b) is a graph showing the strength of the radio wave transmitted in (a) It is.

図1に示すように、降雨マージン算出装置1(1A)は、所定の観測地点において観測された10分間ごとの雨量(10分間雨量)を示す地点別雨量データに基づいて地点ごとの降雨減衰量の確率分布を示す減衰量を解析し、降雨の発生時における衛星放送の電波の減衰量を解析して降雨マージンを算出するものである。   As shown in FIG. 1, the rain margin calculating apparatus 1 (1A) is configured to calculate the rain attenuation amount for each point based on the point-by-point rainfall data indicating the rainfall for every 10 minutes (rainfall for 10 minutes) observed at a predetermined observation point. The rain margin is calculated by analyzing the attenuation amount indicating the probability distribution of the satellite broadcast and analyzing the attenuation amount of the radio wave of the satellite broadcast when the rain occurs.

なお、ここでは、降雨マージン算出装置1(1A)は、降雨減衰の解析対象となる地域である日本国内全域を、5km四方の正方形の単位で分割した区分ごとの1時間雨量の予測値である雨量データを入力し、メッシュごとに減衰量を解析して、降雨マージンを算出することとした。この雨量データには、例えば、気象庁によって発表される、メッシュごとの1時間雨量の予測値を示す降水短時間予報のデータを用いることができる。なお、この区分は、前記のものに限定されることなく、例えば、10km四方の正方形の単位で分割したものであってもよいし、また、人口の多い所定の地域等としてもよい。   Here, the rainfall margin calculation device 1 (1A) is a predicted value of the rainfall for one hour for each section obtained by dividing the entire area of Japan, which is an analysis target of rain attenuation, into square units of 5 km square. The rainfall margin was calculated by inputting the rainfall data and analyzing the attenuation for each mesh. For this rainfall data, for example, data of a short-term precipitation forecast that is announced by the Japan Meteorological Agency and indicates a predicted value of one-hour rainfall for each mesh can be used. In addition, this division is not limited to the above-mentioned ones, and may be divided, for example, in units of a 10 km square, or may be a predetermined area with a large population.

また、図2(a)に示すように、この電波を送信するアンテナを、フェーズドアレーアンテナAとすることで、電波を送受信する地域の一部の範囲に特に高い強度の電波(増力ビーム)を送信できる。このフェーズドアレーアンテナAは、電波を送信する複数のホーンh、h、…と、このホーンh、h、…によって送信された電波を反射して、電波を地上の送受信局に送信する反射鏡Rとを備え、各々のホーンhから送信される電波の位相や振幅を調整することで、送信される電波の強度を所望の放射パターンにすることができる。   Further, as shown in FIG. 2 (a), the antenna for transmitting the radio wave is a phased array antenna A, so that a particularly high intensity radio wave (intensifying beam) is transmitted to a part of the area where the radio wave is transmitted and received. Can be sent. This phased array antenna A reflects a plurality of horns h, h,... That transmit radio waves, and a reflector R that reflects the radio waves transmitted by the horns h, h,. By adjusting the phase and amplitude of the radio wave transmitted from each horn h, the intensity of the radio wave transmitted can be set to a desired radiation pattern.

そのため、フェーズドアレーアンテナAが、本発明における降雨マージン算出装置1(1A)によって算出された降雨マージンを送信電力に与えて、降雨による減衰量の大きい範囲D1、D2に増力ビームB1、B2を送信することで、効率よく降雨減衰補償を行うことができる。図2の例では、降雨が集中する範囲D1、D2の減衰量が特に大きいため、図2(b)に示すように、フェーズドアレーアンテナAによって、この範囲D1、D2には、晴天時に必要となる送信電力値Icに、更に降雨マージンIb1、Ib2を加えた増力ビームB1、B2を送信する。 Therefore, the phased array antenna A gives the rain margin calculated by the rain margin calculating device 1 (1A) according to the present invention to the transmission power, and the boost beam B 1 , in the ranges D 1 and D 2 where the attenuation due to rain is large, By transmitting B 2 , rain attenuation compensation can be performed efficiently. In the example of FIG. 2, because the attenuation amount in the range D 1, D 2 which rainfall is concentrated particularly large, as shown in FIG. 2 (b), the phased array antenna A, in this range D 1, D 2 are, Booster beams B 1 and B 2 are transmitted by adding rain margins Ib 1 and Ib 2 to the transmission power value Ic required in fine weather.

[降雨マージン算出装置の構成(第一の実施の形態)]
次に、図3を参照して、本発明における第一の実施の形態である降雨マージン算出装置1の構成について説明する。図3は、本発明における第一の実施の形態である降雨マージン算出装置の構成を示したブロック図である。降雨マージン算出装置1は、外部から地点別雨量データを入力し、この観測地点ごとの降雨減衰量の確率分布を解析するとともに、この降雨減衰量の確率分布を示す減衰量分布データと外部から入力された雨量データと、1時間あたりで放送が遮断されない時間を示す要求放送時間データとに基づいて、降雨マージンを算出するものである。降雨マージン算出装置1は、降雨減衰解析装置3と、メッシュ減衰量解析装置5とを備える。
[Configuration of Rainfall Margin Calculation Device (First Embodiment)]
Next, the configuration of the rain margin calculating apparatus 1 according to the first embodiment of the present invention will be described with reference to FIG. FIG. 3 is a block diagram showing the configuration of the rain margin calculating apparatus according to the first embodiment of the present invention. The rainfall margin calculation device 1 inputs rainfall data for each point from the outside, analyzes the probability distribution of rainfall attenuation at each observation point, and inputs the attenuation distribution data indicating the probability distribution of this rainfall attenuation from the outside. The rain margin is calculated on the basis of the rain data and the requested broadcast time data indicating the time during which broadcasting is not interrupted per hour. The rainfall margin calculation device 1 includes a rain attenuation analysis device 3 and a mesh attenuation amount analysis device 5.

(降雨減衰解析装置の構成)
降雨減衰解析装置3は、地点別雨量データに基づいて、各々の観測地点における衛星放送の電波の減衰量の確率分布を解析するものである。ここでは、降雨減衰解析装置3は、地点別雨量データ入力手段31と、降雨強度分布解析手段32と、解析データ記憶手段33と、減衰解析手段34と、減衰量分布データ出力手段35とを備える。そして、ここでは、降雨減衰解析装置3は、複数の観測地点において継続的に観測された10分間雨量を示す地点別雨量データを入力し、観測地点ごとに衛星放送の電波の減衰量の確率分布を解析して、この減衰量の確率分布を示す減衰量分布データを、観測地点と、この地点別雨量データによって示される10分間雨量を1時間分加算した1時間雨量とに対応させてメッシュ減衰量解析装置5の減衰量分布データ記憶手段52に記憶することとした。
(Configuration of rain attenuation analyzer)
The rain attenuation analysis device 3 analyzes the probability distribution of the attenuation amount of the radio wave of the satellite broadcast at each observation point based on the rainfall data for each point. Here, the rain attenuation analysis device 3 includes point-by-point rainfall data input means 31, rainfall intensity distribution analysis means 32, analysis data storage means 33, attenuation analysis means 34, and attenuation distribution data output means 35. . Here, the rain attenuation analysis device 3 inputs point-by-point rainfall data indicating 10-minute rainfall continuously observed at a plurality of observation points, and the probability distribution of the attenuation of radio waves of satellite broadcasting at each observation point. The attenuation distribution data indicating the probability distribution of this attenuation amount is correlated with the observation point and the 1-hour rainfall amount obtained by adding the 10-minute rainfall amount indicated by this point-specific rainfall data for one hour. The attenuation amount distribution data storage means 52 of the amount analysis device 5 is stored.

地点別雨量データ入力手段31は、外部から地点別雨量データを入力するものである。ここでは、地点別雨量データ入力手段31は、全国約1300箇所の観測地点に設置されたAMeDASの雨量計によって過去に観測された複数の10分間雨量のデータを、地点別雨量データとして入力することとした。ここで入力された地点別雨量データは、降雨強度分布解析手段32に出力される。   The point-by-point rainfall data input means 31 inputs point-by-point rainfall data from the outside. Here, the point-by-point rainfall data input means 31 inputs a plurality of 10-minute rainfall data observed in the past by AMeDAS rain gauges installed at about 1,300 observation points nationwide as point-by-point rain data. It was. The point-by-point rainfall data input here is output to the rainfall intensity distribution analysis means 32.

降雨強度分布解析手段32は、地点別雨量データ入力手段31から入力された地点別雨量データに基づいて、観測地点ごとに継続して観測された6つの10分間雨量の合計である1時間(所定時間)あたりの雨量を算出し、この1時間雨量が所定値となる時の、当該10分間雨量によって示されるこの10分間における降雨強度(10分間降雨強度)の確率分布を対数正規分布で近似して、当該確率分布を示す降雨強度分布データを生成するものである。ここで生成された降雨強度分布データは、減衰解析手段34の減衰分布算出手段342に出力される。   The rainfall intensity distribution analysis means 32 is based on the point-by-point rainfall data input from the point-by-point rainfall data input means 31 for one hour (predetermined total of six 10-minute rainfalls observed at each observation point. The amount of rainfall per hour) is calculated, and the probability distribution of the 10-minute rainfall intensity (10-minute rainfall intensity) indicated by the 10-minute rainfall when the 1-hour rainfall reaches a predetermined value is approximated by a lognormal distribution. Thus, rainfall intensity distribution data indicating the probability distribution is generated. The rainfall intensity distribution data generated here is output to the attenuation distribution calculating means 342 of the attenuation analyzing means 34.

この所定値となる1時間雨量は、例えば、5mm、10mm、15mm、…のように5mmおきの雨量を設定することとしてもよいし、2mm、5mm、10mm、20mm、…のように雨量が小さいほど間隔を狭くして設定することとしてもよいし、また、降雨減衰の生じやすい大きい雨量の間隔を狭くして設定することとしてもよい。そして、降雨強度分布解析手段32は、入力された10分間雨量を1時間分加算して1時間雨量を算出し、この1時間雨量が所定値となる当該6つの10分間雨量の組を複数取り出して、観測地点及び1時間雨量ごとに、10分間雨量によって示される10分間降雨強度の確率分布を対数正規分布で近似する。そして、この確率分布を示す降雨強度分布データを生成する。   The predetermined one-hour rainfall may be set at 5 mm intervals such as 5 mm, 10 mm, 15 mm,..., And the rainfall is small as 2 mm, 5 mm, 10 mm, 20 mm,. The interval may be set to be narrower, or the interval between large rainfalls where rain attenuation is likely to occur may be set to be narrower. Then, the rainfall intensity distribution analysis means 32 calculates the 1-hour rainfall by adding the input 10-minute rainfall for one hour, and takes out a plurality of sets of the six 10-minute rainfalls for which the one-hour rainfall is a predetermined value. Then, the probability distribution of the 10-minute rainfall intensity indicated by the 10-minute rainfall is approximated by a lognormal distribution for each observation point and one-hour rainfall. Then, rainfall intensity distribution data indicating this probability distribution is generated.

ここで、図4を参照(適宜図3参照)して、降雨強度分布解析手段32が、10分間降雨強度の確率分布を対数正規分布で近似する理由について説明する。図4は、東京及び鹿児島の2つの観測地点において観測された、1時間雨量が5mmとなる降雨の1分間降雨強度及び10分間降雨強度の確率分布を対数正規確率紙にプロットしたグラフである。なお、ここでは、東京における降雨の1分間降雨強度を三角で示し、10分間降雨強度の近似直線を点線で示した。また、鹿児島における降雨の1分間降雨強度を丸で示し、10分間降雨強度の近似直線を一点鎖線で示した。   Here, with reference to FIG. 4 (refer to FIG. 3 as appropriate), the reason why the rainfall intensity distribution analysis means 32 approximates the probability distribution of the 10-minute rainfall intensity with a lognormal distribution will be described. FIG. 4 is a graph obtained by plotting a probability distribution of a 1-minute rainfall intensity and a 10-minute rainfall intensity of 1-minute rainfall observed at two observation points in Tokyo and Kagoshima on a log-normal probability paper. Here, the rainfall intensity for 1 minute of rainfall in Tokyo is indicated by a triangle, and the approximate straight line of the rainfall intensity for 10 minutes is indicated by a dotted line. Moreover, the 1-minute rainfall intensity in Kagoshima is indicated by a circle, and the approximate straight line of the 10-minute rainfall intensity is indicated by a one-dot chain line.

図4に示すように、2箇所の観測地点(東京及び鹿児島)において1分間雨量及び10分間雨量を計測し、1時間雨量が5mmとなる1分間雨量及び10分間雨量に基づいて1分間降雨強度及び10分間降雨強度を算出して、当該1分間降雨強度及び10分間降雨強度の各々の値を超える確率である超過確率を算出して、対数正規確率紙に打点したところ、それぞれの観測地点の1分間降雨強度及び10分間降雨強度の点がほぼ直線上に配列された。そのため、1分間降雨強度及び10分間降雨強度の確率分布は対数正規分布で近似できることが分かった。また、各々の観測地点の1分間降雨強度及び10分間降雨強度の点がほぼ同一直線上に配列されたため、1分間降雨強度及び10分間降雨強度がほぼ同一の確率分布を有することが分かった。そこで、降雨強度分布解析手段32は、地点別雨量データ入力手段31から入力される地点別雨量データに基づいて観測地点及び1時間降雨強度ごとに、10分間雨量から10分間降雨強度を算出し、この10分間降雨強度の確率分布を対数正規分布で近似することとした。   As shown in FIG. 4, 1-minute rainfall and 10-minute rainfall are measured at two observation points (Tokyo and Kagoshima), and 1-minute rainfall intensity based on 1-minute rainfall and 10-minute rainfall with 1-hour rainfall of 5 mm. And 10-minute rainfall intensity, and the excess probability, which is the probability of exceeding the 1-minute rainfall intensity and 10-minute rainfall intensity, is calculated on the log-normal probability paper. The points of 1 minute rainfall intensity and 10 minute rainfall intensity were arranged almost on a straight line. Therefore, it was found that the probability distribution of 1-minute rainfall intensity and 10-minute rainfall intensity can be approximated by a lognormal distribution. Moreover, since the points of 1-minute rainfall intensity and 10-minute rainfall intensity at each observation point were arranged on substantially the same straight line, it was found that the 1-minute rainfall intensity and the 10-minute rainfall intensity had almost the same probability distribution. Therefore, the rainfall intensity distribution analyzing means 32 calculates the rainfall intensity for 10 minutes from the rainfall for 10 minutes for each observation point and 1 hour rainfall intensity based on the rainfall data for each point input from the rainfall data input means 31 for each point. The probability distribution of this 10-minute rainfall intensity was approximated by a lognormal distribution.

図3に戻って説明を続ける。以下、降雨強度分布解析手段32が、地点別雨量データに基づいて、1時間雨量が所定値となる10分間降雨強度の確率分布を対数正規分布で近似して、当該確率分布を示す降雨強度分布データを生成する方法の例を説明する。降雨強度分布解析手段32は、観測地点及び所定値の1時間雨量ごとに、10分間雨量を1時間分加算した1時間雨量が所定値となる6個の10分間雨量の組を複数取り出す。そして、降雨強度分布解析手段32は、10分間雨量の各々について、値を6倍して10分間降雨強度を算出する。そして、降雨強度分布解析手段32は、各々の組の6個の10分間降雨強度の値を昇順に並び替え、値が小さいものから順に5つの10分間降雨強度の超過確率を83.33%、66.67%、50%、33.33%、16.67%とする。   Returning to FIG. 3, the description will be continued. Hereinafter, the rainfall intensity distribution analysis means 32 approximates the probability distribution of the 10-minute rainfall intensity with a predetermined value of the one-hour rainfall based on the point-by-point rainfall data with a lognormal distribution, and shows the rainfall intensity distribution indicating the probability distribution. An example of a method for generating data will be described. The rainfall intensity distribution analyzing means 32 takes out a plurality of sets of six 10-minute rainfalls, each of which has a predetermined value of the one-hour rainfall obtained by adding the 10-minute rainfall for one hour for each observation point and a predetermined value of one-hour rainfall. Then, the rainfall intensity distribution analysis means 32 calculates the 10-minute rainfall intensity by multiplying the value by 6 for each 10-minute rainfall. Then, the rainfall intensity distribution analysis means 32 sorts the six 10-minute rainfall intensity values of each set in ascending order, and sets the five 10-minute rainfall intensity excess probabilities in the order of increasing values from 83.33%, 66.67%, 50%, 33.33%, and 16.67%.

更に、降雨強度分布解析手段32は、複数の組の10分間降雨強度とその超過確率とに基づいて、超過確率p(pは83.33%、66.67%、50%、33.33%、16.67%のいずれか)における10分間降雨強度Rpを、以下の式(1)によって算出する。ここで、Median{}は括弧内の数値の中央値を示す。また、Rp[k]は、n個の組のうちのk番目の組において超過確率がpとなる10分間降雨強度を示す。
Rp=Median{Rp[1]、Rp[2]、…、Rp[n]} …(1)
Further, the rainfall intensity distribution analyzing means 32 determines the excess probability p (p is 83.33%, 66.67%, 50%, 33.33%) based on a plurality of sets of 10-minute rainfall intensity and the excess probability. , 16.67%) is calculated by the following formula (1). Here, Median {} indicates the median value in the parentheses. Further, Rp [k] represents the 10-minute rainfall intensity in which the excess probability is p in the kth group among the n groups.
Rp = Median {Rp [1], Rp [2],..., Rp [n]} (1)

そして、降雨強度分布解析手段32は前記の式(1)によって算出された、超過確率が83.33%、66.67%、50%、33.33%、16.67%の10分間降雨強度の平均値と標準偏差とを算出し、降雨強度分布データとする。なお、前記したように、10分間降雨強度と1分間降雨強度はほぼ一致するため、ここで算出された10分間降雨強度の確率分布の平均値と標準偏差を1分間降雨強度の確率分布の平均値と標準偏差とし、後記する減衰分布算出手段342によって減衰量を算出する際に用いることとした。   Then, the rainfall intensity distribution analysis means 32 calculates the intensity of precipitation for 10 minutes with the probability of excess being 83.33%, 66.67%, 50%, 33.33%, 16.67% calculated by the above equation (1). The average value and the standard deviation are calculated and used as rainfall intensity distribution data. As described above, since the 10-minute rainfall intensity and the 1-minute rainfall intensity substantially coincide, the average value and standard deviation of the 10-minute rainfall intensity calculated here are the average of the 1-minute rainfall intensity probability distribution. The value and the standard deviation are used when calculating the attenuation amount by the attenuation distribution calculating means 342 described later.

解析データ記憶手段33は、後記する減衰解析手段34によって減衰量の確率分布を算出するために必要となるデータを予め記憶するもので、半導体メモリ、ハードディスク等の一般的な記憶手段である。ここでは、解析データ記憶手段33に、後記する降雨区間長算出手段341によって、衛星放送の電波の伝送路上において降雨が発生している区間の長さである降雨区間長を算出する際に必要となる、観測地点位置データ、衛星軌道位置データ、及び、雨滴層高度データと、後記する減衰分布算出手段342によって減衰量の確率分布を算出する際に必要となる周波数パラメータとを記憶することとした。   The analysis data storage means 33 stores in advance data necessary for calculating the probability distribution of the attenuation amount by the attenuation analysis means 34 described later, and is a general storage means such as a semiconductor memory or a hard disk. Here, it is necessary for the analysis data storage means 33 to calculate the rain section length, which is the length of the section where the rain is occurring on the satellite broadcast radio wave transmission path, by the rain section length calculating means 341 described later. The observation point position data, the satellite orbit position data, the raindrop height data, and the frequency parameters necessary for calculating the probability distribution of the attenuation amount by the attenuation distribution calculating means 342 described later are stored. .

なお、観測地点位置データは、観測地点の位置(緯度と経度)を示すものである。また、衛星軌道位置データは、衛星の軌道位置を示すものである。更に、雨滴層高度データは、観測地点ごとの、降雨雲が雨滴から氷晶へと変化する気温0℃高度である雨滴層高度[km]を示すものである。   The observation point position data indicates the position (latitude and longitude) of the observation point. The satellite orbit position data indicates the orbit position of the satellite. Further, the raindrop layer altitude data indicates the raindrop layer altitude [km] at a temperature of 0 ° C. at which the rain clouds change from raindrops to ice crystals for each observation point.

また、周波数パラメータは、降雨区間長の単位長さあたりの減衰量を示す降雨減衰係数A0が、以下の式(2)によって近似される場合において、電波の周波数と偏波とに依存するパラメータγ、nと、周波数とを対応させたものである。ここで、Rは1分間降雨強度[mm/h]である。
0=γ・Rn[dB/km] …(2)
The frequency parameter is a parameter that depends on the frequency and polarization of the radio wave when the rain attenuation coefficient A 0 indicating the attenuation per unit length of the rain section length is approximated by the following equation (2). γ, n and frequency are associated with each other. Here, R is the rainfall intensity [mm / h] for 1 minute.
A 0 = γ · R n [dB / km] (2)

減衰解析手段34は、降雨強度分布解析手段32によって生成された降雨強度分布データに基づいて、降雨強度が解析された観測地点及び1時間雨量に対応する減衰量の確率分布を示す減衰量分布データを生成するものである。ここでは、減衰解析手段34は、降雨区間長算出手段341と、減衰分布算出手段342とを備える。   The attenuation analysis means 34 is attenuation distribution data indicating the probability distribution of the attenuation corresponding to the observation point where the rainfall intensity is analyzed and the one hour rainfall based on the rainfall intensity distribution data generated by the rainfall intensity distribution analysis means 32. Is generated. Here, the attenuation analysis unit 34 includes a rainfall section length calculation unit 341 and an attenuation distribution calculation unit 342.

降雨区間長算出手段341は、観測地点ごとに降雨区間長を算出するものである。ここで算出された降雨区間長は、減衰分布算出手段342に出力される。なお、この降雨区間長は、例えば、非特許文献1に記載されている方法等の一般的な方法で算出することができる。ここで、図5を参照(適宜図3参照)して、降雨区間長算出手段341が、降雨区間長を算出する例について説明する。図5は、雨滴層高度と、衛星放送の電波の降雨区間長とを模式的に示した説明図である。   The rainfall section length calculation means 341 calculates the rain section length for each observation point. The rainfall section length calculated here is output to the attenuation distribution calculating means 342. In addition, this rainfall section length is computable by general methods, such as the method described in the nonpatent literature 1, for example. Here, with reference to FIG. 5 (refer to FIG. 3 as appropriate), an example in which the rainfall section length calculation means 341 calculates the rain section length will be described. FIG. 5 is an explanatory diagram schematically showing the height of the raindrop layer and the length of the rain section of satellite broadcast radio waves.

図5に示すように、仰角θ[°]の衛星放送の電波の伝送路Tにおいて、観測地点Oから気温0℃高度までの長さD[km]が、降雨区間長となる。そして、降雨区間長算出手段341は、この降雨区間長Dを、以下の式(3)に基づいて算出する。なお、H[km]は、各々の観測地点Oにおける雨滴層高度であり、また、伝送路Tの仰角θは、解析データ記憶手段33に記憶された、観測地点位置データと衛星軌道位置データとに基づいて算出することができる。
D=H・cosecθ[km] …(3)
As shown in FIG. 5, in a satellite broadcast radio wave transmission path T at an elevation angle θ [°], the length D [km] from the observation point O to the temperature of 0 ° C. is the rainfall section length. Then, the rainfall section length calculation means 341 calculates the rain section length D based on the following equation (3). H [km] is the raindrop height at each observation point O, and the elevation angle θ of the transmission line T is the observation point position data and satellite orbit position data stored in the analysis data storage means 33. Can be calculated based on
D = H · cosecθ [km] (3)

図3に戻って説明を続ける。減衰分布算出手段342は、降雨区間長算出手段341から入力される降雨区間長と、降雨強度分布解析手段32から入力される降雨強度分布データとに基づいて、観測地点ごとの降雨による減衰量の確率分布を解析し、この確率分布を示す減衰量分布データを生成するものである。ここで生成された減衰量分布データは、観測地点と、1時間雨量とに対応付けられて減衰量分布データ出力手段35に出力される。   Returning to FIG. 3, the description will be continued. The attenuation distribution calculating means 342 is configured to calculate the amount of attenuation due to rain for each observation point based on the rainfall section length input from the rainfall section length calculating means 341 and the rain intensity distribution data input from the rainfall intensity distribution analyzing means 32. The probability distribution is analyzed, and attenuation distribution data indicating the probability distribution is generated. The attenuation distribution data generated here is output to the attenuation distribution data output means 35 in association with the observation point and the one-hour rainfall.

なお、ここでは、減衰量の確率分布が対数正規分布で近似できることとし、減衰分布算出手段342は、降雨強度分布解析手段32によって解析された10分間降雨強度の確率分布の平均値と標準偏差を、1分間降雨強度の確率分布の平均値と標準偏差とみなして、この平均値と標準偏差とを、この1分間降雨強度の確率分布を有する降雨による減衰量の確率分布を示す平均値と標準偏差とに変換して、減衰量分布データとすることとした。ここで、図6を参照(適宜図3参照)して、減衰量の確率分布を対数正規分布で近似する理由について説明する。図6は、東京において観測された、1時間雨量が5mm及び10mmとなる降雨の1分間降雨強度及び減衰量の確率分布を対数正規確率紙にプロットしたグラフである。なお、ここでは、1時間雨量が5mmの時の1分間降雨強度を三角、減衰量を白抜きの三角、1時間雨量が10mmの時の1分間降雨強度を丸、減衰量を白抜きの丸で示した。また、各々の分布を示す点の近似直線も示した。   Here, it is assumed that the probability distribution of attenuation can be approximated by a lognormal distribution, and the attenuation distribution calculating means 342 calculates the average value and standard deviation of the probability distribution of 10-minute rainfall intensity analyzed by the rainfall intensity distribution analyzing means 32. Considering the average value and standard deviation of the probability distribution of 1-minute rainfall intensity, the average value and standard deviation are used as the average value and standard indicating the probability distribution of attenuation due to rainfall having the probability distribution of 1-minute rainfall intensity. It was decided that it would be converted into a deviation and used as attenuation amount distribution data. Here, the reason why the probability distribution of the attenuation amount is approximated by a lognormal distribution will be described with reference to FIG. 6 (refer to FIG. 3 as appropriate). FIG. 6 is a graph obtained by plotting probability distributions of 1-minute rainfall intensity and attenuation of rainfall with an hourly rainfall of 5 mm and 10 mm observed in Tokyo on a lognormal probability paper. In addition, here, the rainfall intensity for 1 minute when the hourly rainfall is 5 mm is triangular, the attenuation is white triangle, and the rainfall intensity for 1 minute when the hourly rainfall is 10 mm is round, and the attenuation is white. It showed in. Moreover, the approximate straight line of the point which shows each distribution was also shown.

図6に示すように、1時間雨量が5mmと、10mmの場合について、1分間雨量と12GHz帯の電波の降雨による減衰量とを計測し、1分間雨量から算出された1分間降雨強度及び減衰量の各々の値を超える確率である超過確率を算出して、対数正規確率紙に打点したところ、それぞれの点がほぼ直線上に配列された。そのため、1分間降雨強度と同様に、減衰量の確率分布も対数正規分布で近似できることが分かった。   As shown in FIG. 6, when the hourly rainfall is 5 mm and 10 mm, the one-minute rainfall and the attenuation due to the rain of the radio wave in the 12 GHz band are measured, and the one-minute rainfall intensity and attenuation calculated from the one-minute rainfall. When the excess probability, which is the probability of exceeding each value of the quantity, was calculated and scored on a lognormal probability paper, each point was arranged on a substantially straight line. Therefore, it was found that the probability distribution of attenuation can be approximated by a lognormal distribution as well as the 1-minute rainfall intensity.

図3に戻って説明を続ける。ここで、減衰分布算出手段342は、降雨強度分布解析手段32によって解析された降雨強度分布データによって示される平均値と標準偏差とを、例えば、非特許文献1に記載されている方法等によって、減衰量の確率分布を示す平均値と標準偏差とに変換することができる。   Returning to FIG. 3, the description will be continued. Here, the attenuation distribution calculating unit 342 calculates the average value and the standard deviation indicated by the rainfall intensity distribution data analyzed by the rainfall intensity distribution analyzing unit 32 by, for example, a method described in Non-Patent Document 1. It can convert into the average value and standard deviation which show the probability distribution of attenuation amount.

以下、減衰分布算出手段342が降雨強度分布データによって示される平均値と標準偏差とを、減衰量の確率分布を示す平均値と標準偏差とに変換する方法について説明する。減衰量Aは、以下の式(4)で表される。ここで、γ、nは、解析データ記憶手段33に記憶された周波数パラメータ、Rは1分間降雨強度[mm/h]、積分は伝送路にわたっての線積分である。
A=γ・∫Rndx[dB] …(4)
Hereinafter, a method in which the attenuation distribution calculating unit 342 converts the average value and the standard deviation indicated by the rainfall intensity distribution data into the average value and the standard deviation indicating the probability distribution of the attenuation amount will be described. The attenuation amount A is expressed by the following formula (4). Here, γ and n are frequency parameters stored in the analysis data storage means 33, R is the rainfall intensity [mm / h] for 1 minute, and the integral is a line integral over the transmission line.
A = γ · ∫R n dx [dB] (4)

そして、1分間降雨強度Rが対数正規分布に従い、log10Rの平均値をm、標準偏差をsとすると、減衰量Aの分布は、対数正規分布で近似でき、log10Aの平均値me、標準偏差seは、式(5)、式(6)に示すようになる。なお、Kは、log10e(=0.4343)である。 Then, if the rainfall intensity R for one minute follows a lognormal distribution, the average value of log 10 R is m, and the standard deviation is s, the distribution of the attenuation A can be approximated by a log normal distribution, and the average value m of log 10 A e and standard deviation s e are as shown in equations (5) and (6). Note that K is log 10 e (= 0.4343).

Figure 0004511979
Figure 0004511979

ここで、μA、σA 2は、減衰量Aの平均値及び分散で、それぞれ、以下の式(7)、式(8)によって示される。なお、Dは、前記の式(3)によって求められた降雨区間長、αは、1分間降雨強度Rの空間相関を示すパラメータ(α≒0.35)である。 Here, μ A and σ A 2 are the average value and the variance of the attenuation A, and are represented by the following equations (7) and (8), respectively. Here, D is the length of the rainfall section obtained by the above equation (3), and α is a parameter (α≈0.35) indicating the spatial correlation of the rainfall intensity R for 1 minute.

Figure 0004511979
Figure 0004511979

なお、減衰解析手段34が降雨強度分布データに基づいて減衰量の確率分布を解析する方法は前記の方法に限定されることなく、1分間降雨強度から減衰量を算出する他の方法を適用してもよい。また、ここでは、減衰解析手段34が、降雨強度分布データによって示される平均値と標準偏差を減衰量の確率分布を示す平均値と標準偏差とに変換することとしたが、例えば、減衰解析手段34が、所定の超過確率となる降雨強度を複数算出し、各々の降雨強度を1分間降雨強度として減衰量に変換することとしてもよい。   The method for the attenuation analysis means 34 to analyze the probability distribution of the attenuation based on the rainfall intensity distribution data is not limited to the above method, and other methods for calculating the attenuation from the rainfall intensity for 1 minute are applied. May be. Here, the attenuation analysis unit 34 converts the average value and the standard deviation indicated by the rainfall intensity distribution data into the average value and the standard deviation indicating the probability distribution of the attenuation amount. For example, the attenuation analysis unit 34 34 may calculate a plurality of rainfall intensities that have a predetermined excess probability, and convert each of the rainfall intensities as a one-minute rainfall intensity into an attenuation amount.

減衰量分布データ出力手段35は、減衰解析手段34から入力された減衰量分布データを外部へ出力するものである。なお、ここでは、減衰量分布データ出力手段35は、減衰量分布データを観測地点と1時間雨量とに対応付けてメッシュ減衰量解析装置5に出力することとした。   The attenuation amount distribution data output means 35 outputs the attenuation amount distribution data input from the attenuation analysis means 34 to the outside. Here, the attenuation amount distribution data output means 35 outputs the attenuation amount distribution data to the mesh attenuation amount analysis device 5 in association with the observation point and the one hour rainfall.

(メッシュ減衰量解析装置の構成)
次に、メッシュ減衰量解析装置5の構成について説明する。メッシュ減衰量解析装置5は、外部から入力される、メッシュごとの1時間雨量の予測値を示す雨量データと、降雨減衰解析装置3によって生成された減衰量分布データとに基づいて、衛星放送の電波の減衰量を解析し、降雨による減衰を補償する降雨マージンを算出するものである。ここでは、メッシュ減衰量解析装置5は、減衰量分布データ入力手段51と、減衰量分布データ記憶手段52と、雨量データ入力手段53と、減衰量分布データ読み出し手段54と、要求放送時間データ入力手段55と、電波減衰解析手段56と、算出データ記憶手段57と、降雨マージン算出手段58と、降雨マージン出力手段59とを備える。
(Configuration of mesh attenuation analysis device)
Next, the configuration of the mesh attenuation amount analysis apparatus 5 will be described. The mesh attenuation amount analysis device 5 is based on the rainfall data indicating the predicted value of the hourly rainfall for each mesh and the attenuation amount distribution data generated by the rain attenuation analysis device 3 input from the outside. It analyzes the amount of attenuation of radio waves and calculates a rain margin that compensates for attenuation due to rain. Here, the mesh attenuation amount analysis apparatus 5 includes an attenuation amount distribution data input unit 51, an attenuation amount distribution data storage unit 52, a rain amount data input unit 53, an attenuation amount distribution data reading unit 54, and a requested broadcast time data input. Means 55, radio wave attenuation analysis means 56, calculation data storage means 57, rainfall margin calculation means 58, and rainfall margin output means 59 are provided.

減衰量分布データ入力手段51は、降雨減衰解析装置3の減衰解析手段34によって生成された減衰量分布データを減衰量分布データ出力手段35から入力し、減衰量分布データ記憶手段52に記憶するものである。   The attenuation distribution data input means 51 receives the attenuation distribution data generated by the attenuation analysis means 34 of the rainfall attenuation analysis device 3 from the attenuation distribution data output means 35 and stores it in the attenuation distribution data storage means 52. It is.

減衰量分布データ記憶手段(減衰量分布データ記憶装置)52は、降雨減衰解析装置3によって解析された減衰量分布データを観測地点と1時間雨量とに対応付けて記憶するもので、半導体メモリ、ハードディスク等の一般的な記憶手段である。この減衰量分布データは、後記する減衰量分布データ読み出し手段54によって読み出され、電波減衰解析手段56によって、メッシュごとの衛星放送の電波の減衰量を解析する際に参照されて用いられる。   The attenuation distribution data storage means (attenuation distribution data storage device) 52 stores the attenuation distribution data analyzed by the rain attenuation analysis device 3 in association with the observation point and the one-hour rainfall. This is a general storage means such as a hard disk. This attenuation amount distribution data is read by an attenuation amount distribution data reading unit 54 described later, and is used by being referred to when the attenuation amount of satellite broadcast radio waves for each mesh is analyzed by the radio wave attenuation analysis unit 56.

雨量データ入力手段53は、メッシュごとの1時間雨量を示す雨量データを外部から入力するものである。なお、ここでは、この雨量データをメッシュごとの1時間雨量の予測値とするが、この雨量データは、予め設定された所定時間あたりの雨量や、AMeDAS等によって所定の単位のメッシュごとに観測された所定時間あたりの雨量の実測値であってもよい。ここで入力された雨量データは、減衰量分布データ読み出し手段54に出力される。   The rainfall data input means 53 inputs rain data indicating the rainfall for one hour for each mesh from the outside. Here, this rainfall data is used as a predicted value for one hour of rainfall for each mesh, but this rainfall data is observed for each predetermined unit of mesh by a predetermined amount of rainfall per predetermined time, AMeDAS, or the like. Alternatively, it may be an actual measurement value of rainfall per predetermined time. The rain data input here is output to the attenuation distribution data reading means 54.

減衰量分布データ読み出し手段54は、雨量データ入力手段53から入力された雨量データに基づいて、減衰量分布データ記憶手段52から減衰量分布データを読み出すものである。ここで、すべてのメッシュはその近傍の観測地点と予め対応付けられており、減衰量分布データ読み出し手段54は、メッシュごとに、各々のメッシュに予め対応付けられた観測地点における、雨量データによって示される1時間雨量に対応する減衰量分布データを読み出す。ここで読み出された減衰量分布データは、電波減衰解析手段56に出力される。   The attenuation distribution data reading unit 54 reads the attenuation distribution data from the attenuation distribution data storage unit 52 based on the rainfall data input from the rainfall data input unit 53. Here, all the meshes are associated in advance with the observation points in the vicinity thereof, and the attenuation amount distribution data reading means 54 is indicated by the rainfall data at the observation points associated in advance with each mesh for each mesh. The attenuation distribution data corresponding to the one hour rainfall is read out. The attenuation amount distribution data read here is output to the radio wave attenuation analysis means 56.

要求放送時間データ入力手段55は、1時間の放送時間のうち送受信局の受信機における映像の再生が可能となる時間として設定される値を示す要求放送時間データを外部から入力するものである。この要求放送時間データ(許容値データ)は、減衰量分布データにおける超過確率の許容値を示し、例えば、この要求放送時間データが1時間の放送時間のうちの59分を示す場合には、超過確率の許容値として98.33%(59/60)を示している。ここで入力された要求放送時間データは、電波減衰解析手段56に出力される。   The requested broadcast time data input means 55 inputs request broadcast time data indicating a value set as a time during which video can be reproduced in the receiver of the transmission / reception station within one hour of broadcast time. This required broadcast time data (allowable value data) indicates an allowable value of the excess probability in the attenuation amount distribution data. For example, when the required broadcast time data indicates 59 minutes of one hour of broadcast time, it exceeds 98.33% (59/60) is shown as the allowable value of the probability. The requested broadcast time data input here is output to the radio wave attenuation analysis means 56.

電波減衰解析手段56は、メッシュごとに、減衰量分布データ読み出し手段54によって読み出された減衰量分布データによって示される減衰量の確率分布において、要求放送時間データ入力手段55から入力される要求放送時間データによって示される超過確率の許容値に対応する減衰量を解析するものである。ここで、例えば、要求放送時間データが、1時間の放送時間のうちの59分を示す場合には、電波減衰解析手段56は、減衰量の確率分布において超過確率が98.33%となる減衰量を求める。ここで解析された減衰量は、放送時間中に変動する降雨による減衰量が98.33%の確率で超過する値を示している。この減衰量は、降雨マージン算出手段58に出力される。   The radio wave attenuation analysis means 56 is a request broadcast input from the request broadcast time data input means 55 in the probability distribution of the attenuation indicated by the attenuation distribution data read by the attenuation distribution data reading means 54 for each mesh. The attenuation corresponding to the allowable value of the excess probability indicated by the time data is analyzed. Here, for example, when the requested broadcast time data indicates 59 minutes of one hour of broadcast time, the radio wave attenuation analysis means 56 attenuates such that the excess probability is 98.33% in the probability distribution of attenuation. Find the amount. The amount of attenuation analyzed here indicates a value at which the amount of attenuation due to rainfall that fluctuates during the broadcast time exceeds with a probability of 98.33%. This attenuation amount is output to the rain margin calculating means 58.

算出データ記憶手段57は、後記する降雨マージン算出手段58による降雨マージンの算出に必要となるデータを記憶するもので、半導体メモリ、ハードディスク等の一般的な記憶手段である。ここでは、算出データ記憶手段57に、送受信局データと、衛星データと、伝送条件データとを記憶することとした。   The calculation data storage means 57 stores data necessary for calculating the rain margin by the rain margin calculation means 58 described later, and is a general storage means such as a semiconductor memory or a hard disk. Here, the calculation data storage means 57 stores transmission / reception station data, satellite data, and transmission condition data.

送受信局データとは、衛星放送の電波を受信あるいは送受信する地上の送受信局のパラメータであり、例えば、送受信局の受信アンテナの利得値、送信アンテナの実効等方放射電力値(EIRP;Effective Isotropic Radiated Power)、送受信局の受信装置のノイズフィギュア値、アンテナのポインティングロス値、カップリングロス値等である。また、衛星データとは、例えば、衛星の受信アンテナの利得値、送信アンテナの実効等方放射電力値、衛星の受信装置のノイズフィギュア値等である。   The transmission / reception station data is a parameter of a ground transmission / reception station that receives or transmits satellite broadcast radio waves. For example, the reception antenna gain value of the transmission / reception station, the effective isotropic radiated power value (EIRP) of the transmission antenna, and the like. Power), the noise figure value of the receiving device of the transmitting / receiving station, the pointing loss value of the antenna, the coupling loss value, and the like. The satellite data includes, for example, a gain value of a satellite reception antenna, an effective isotropic radiated power value of a transmission antenna, a noise figure value of a satellite reception device, and the like.

伝送条件データは、変調方式や誤り訂正方式等の、CN比(career to noise ratio;搬送波電力対雑音電力比)に影響を与える因子に、衛星放送の再生に必要となるCN比を対応させたデータである。この送受信局データ、衛星データ及び伝送条件データは、降雨マージン算出手段58によって降雨マージンを算出する際に参照されて用いられる。   In the transmission condition data, the CN ratio required for satellite broadcast reproduction is made to correspond to factors affecting the CN ratio (carrier to noise ratio) such as modulation scheme and error correction scheme. It is data. The transmission / reception station data, satellite data, and transmission condition data are referred to and used when the rain margin is calculated by the rain margin calculating means 58.

降雨マージン算出手段58は、電波減衰解析手段56から入力された減衰量に基づいて、降雨マージンを算出するものである。ここで算出された降雨マージンは降雨マージン出力手段59に出力される。   The rain margin calculating unit 58 calculates a rain margin based on the attenuation amount input from the radio wave attenuation analyzing unit 56. The rain margin calculated here is output to the rain margin output means 59.

ここで、降雨マージン算出手段58は、算出データ記憶手段57に記憶された送受信局データによって示される送受信局側の利得と、衛星データによって示される衛星側の利得と、酸素吸収による減衰と、電波減衰解析手段56から入力された減衰量となる降雨による減衰の影響をと加えた際に、送受信局の受信機において、伝送条件データによって示される再生に必要なCN比となるために必要な送信電力を、メッシュごとに算出する。送信電力は以下の式(9)によって算出される。ここで、衛星の送信電力をPt[dBm]、送受信局における受信電力をPr[dBm]、送受信局側の利得をGt[dB]、衛星側の利得をGr[dB]、電波の波長をλ[m]、伝送路長をr[m]、酸素吸収による減衰量をa[dB/m]、降雨による減衰量をA[dB]とする。
Pt=Pr−Gt−Gr−20log(λ/4πr)+ar+A …(9)
Here, the rain margin calculating means 58 is a transmission / reception station side gain indicated by the transmission / reception station data stored in the calculation data storage means 57, a satellite side gain indicated by the satellite data, an attenuation due to oxygen absorption, a radio wave. When the influence of attenuation due to rain, which is the amount of attenuation input from the attenuation analysis means 56, is added, the transmission necessary for achieving the CN ratio necessary for reproduction indicated by the transmission condition data in the receiver of the transmission / reception station. The power is calculated for each mesh. The transmission power is calculated by the following equation (9). Here, the transmission power of the satellite is Pt [dBm], the reception power at the transmission / reception station is Pr [dBm], the gain at the transmission / reception station side is Gt [dB], the gain at the satellite side is Gr [dB], and the wavelength of the radio wave is λ [M], the transmission line length is r [m], the attenuation due to oxygen absorption is a [dB / m], and the attenuation due to rainfall is A [dB].
Pt = Pr−Gt−Gr−20 log (λ / 4πr) + ar + A (9)

そして、降雨マージン算出手段58は、図2(a)によって示されたフェーズドアレーアンテナAのホーンh、h、…によって電波が送信される各々の範囲に対応する複数のメッシュにおいて、前記の式(9)によって算出された複数の送信電力の中の最大値を、当該ホーンhから送出する電波の送信電力とし、降雨マージンを算出する。これによって、降雨マージン算出手段58は、各々のホーンh、h、…によって電波が送信される範囲においてすべての送受信局の受信機によって、1時間の放送時間のうち要求放送時間データによって示される時間の再生が可能となる降雨マージンを算出することができる。   Then, the rain margin calculating means 58 uses the above expression (in the plurality of meshes corresponding to the respective ranges in which radio waves are transmitted by the horns h, h,... Of the phased array antenna A shown in FIG. The maximum value among the plurality of transmission powers calculated in 9) is set as the transmission power of the radio wave transmitted from the horn h, and the rain margin is calculated. As a result, the rain margin calculating means 58 performs the time indicated by the requested broadcast time data within the broadcast time of one hour by the receivers of all the transmitting and receiving stations within the range in which radio waves are transmitted by the respective horns h, h,. It is possible to calculate the rain margin that makes it possible to regenerate.

降雨マージン出力手段59は、降雨マージン算出手段58によって算出された降雨マージンを外部に出力するものである。   The rain margin output means 59 outputs the rain margin calculated by the rain margin calculation means 58 to the outside.

以上のように降雨マージン算出装置1を構成することで、降雨マージン算出装置1は、降雨減衰解析装置3によって、観測地点ごとに、衛星放送の電波の降雨による減衰量の確率分布を解析し、減衰量分布データを生成することができる。そして、メッシュ減衰量解析装置5によって、この減衰量分布データと、外部から入力された雨量データと、要求放送時間データに基づいて、送受信局の受信機において要求放送時間データによって示される時間の再生が可能となるような送信電力を解析し、フェーズドアレーアンテナAのホーンh、h、…(図2参照)の各々の降雨マージンを算出することができる。これによって、フェーズドアレーアンテナAによって、降雨減衰の補償を効果的に行うことができる。   By configuring the rain margin calculating device 1 as described above, the rain margin calculating device 1 analyzes the probability distribution of the attenuation amount due to the rain of the satellite broadcast radio wave for each observation point by the rain attenuation analyzing device 3. Attenuation distribution data can be generated. Based on the attenuation distribution data, the rainfall data input from the outside, and the requested broadcast time data, the mesh attenuation amount analysis device 5 reproduces the time indicated by the requested broadcast time data at the receiver of the transmission / reception station. Can be calculated, and the rain margin of each of the horns h, h,... (See FIG. 2) of the phased array antenna A can be calculated. Accordingly, the rain attenuation can be effectively compensated by the phased array antenna A.

なお、特許請求の範囲に記載の区分別減衰量解析装置は、少なくとも減衰量分布データ記憶手段52と、雨量データ入力手段53と、減衰量分布データ読み出し手段54と、電波減衰解析手段56とを備えていればよい。   It should be noted that the category-specific attenuation amount analyzing apparatus described in the claims includes at least an attenuation amount distribution data storage unit 52, a rainfall amount data input unit 53, an attenuation amount distribution data reading unit 54, and a radio wave attenuation analyzing unit 56. It only has to have.

更に、ここでは、雨量データを、メッシュごとの1時間雨量の予測値としたが、降雨マージン算出装置1は、予め設定された所定時間あたりの雨量や、気象レーダ等によって観測された過去の雨量を雨量データとして入力し、この降雨によって生じた減衰量を解析して降雨マージンを算出することとしてもよい。また、所定時間は、例えば、1時間や2時間等であり、10分間より長い時間であればよい。   Furthermore, here, the rainfall data is the predicted value of the hourly rainfall for each mesh. However, the rainfall margin calculation device 1 can calculate the rainfall per predetermined time set in advance or the past rainfall observed by a weather radar or the like. May be input as rain data, and the rain margin may be calculated by analyzing the attenuation caused by the rain. The predetermined time is, for example, 1 hour or 2 hours, and may be a time longer than 10 minutes.

また、メッシュ減衰量解析装置5は、要求放送時間データ入力手段55によって、外部から要求放送時間データを入力することとしたが、超過確率の許容値を示すデータを入力すればよい。また、メッシュ減衰量解析装置5は、図示しない記憶手段に要求放送時間データを記憶し、電波減衰解析手段56によって記憶手段から要求放送時間データを読み出すこととしてもよい。   Further, although the mesh attenuation amount analyzing apparatus 5 inputs the required broadcast time data from the outside by the required broadcast time data input means 55, it is only necessary to input data indicating the allowable value of the excess probability. Further, the mesh attenuation amount analyzing apparatus 5 may store the required broadcast time data in a storage unit (not shown), and the radio wave attenuation analysis unit 56 may read out the required broadcast time data from the storage unit.

更に、メッシュ減衰量解析装置5は、ここでは降雨減衰解析装置3から減衰量分布データを入力し、当該メッシュ減衰量解析装置5内の減衰量分布データ記憶手段52に記憶することとしたが、メッシュ減衰量解析装置5は、外部に減衰量分布データを記憶する減衰量分布データ記憶装置(図示せず)を備えることとしてもよいし、また、減衰量分布データ記憶手段52に減衰量分布データを予め記憶し、降雨減衰解析装置3とは接続されない構成としてもよい。また、降雨減衰解析装置3及びメッシュ減衰量解析装置5は、衛星からの電波の降雨による減衰量を解析することとしたが、本発明の降雨減衰解析装置及び区分別減衰量解析装置は、衛星からの電波に限定されず、様々な電波の降雨による減衰量を解析することができる。   Further, here, the mesh attenuation amount analysis device 5 inputs the attenuation amount distribution data from the rain attenuation amount analysis device 3 and stores it in the attenuation amount distribution data storage means 52 in the mesh attenuation amount analysis device 5. The mesh attenuation amount analysis device 5 may include an attenuation amount distribution data storage device (not shown) that stores attenuation amount distribution data outside, and the attenuation amount distribution data storage unit 52 stores the attenuation amount distribution data. May be stored in advance and not connected to the rain attenuation analysis device 3. Further, the rain attenuation analysis device 3 and the mesh attenuation amount analysis device 5 analyze the attenuation amount of the radio wave from the satellite, but the rain attenuation analysis device and the divisional attenuation amount analysis device according to the present invention are the satellites. It is not limited to radio waves from, but can analyze the attenuation of various radio waves due to rainfall.

また、降雨減衰解析装置3及びメッシュ減衰量解析装置5は、コンピュータにおいて各手段を各機能プログラムとして実現することも可能であり、各機能プログラムを結合して、降雨減衰解析プログラム及びメッシュ別減衰量解析プログラムとして動作させることも可能である。   Further, the rain attenuation analysis device 3 and the mesh attenuation amount analysis device 5 can also realize each means as a function program in a computer, and by combining the function programs, the rain attenuation analysis program and the mesh-specific attenuation amount. It is also possible to operate as an analysis program.

[降雨マージン算出装置の動作]
次に、図7及び図8(適宜図3参照)を参照して、本発明の第一の実施の形態である降雨マージン算出装置1の降雨減衰解析装置3とメッシュ減衰量解析装置5との動作について説明する。
[Operation of rainfall margin calculation device]
Next, referring to FIG. 7 and FIG. 8 (refer to FIG. 3 as appropriate), the rain attenuation analysis device 3 and the mesh attenuation amount analysis device 5 of the rain margin calculation device 1 according to the first embodiment of the present invention. The operation will be described.

(降雨減衰解析装置の動作)
まず、図7を参照(適宜図3参照)して、降雨減衰解析装置3が、観測地点ごとに所定の1時間雨量に対して衛星放送の電波の減衰量の確率分布を示す減衰量分布データを生成する動作について説明する。図7は、本発明の第一の実施の形態における降雨減衰解析装置の動作を示したフローチャートである。
(Operation of rain attenuation analyzer)
First, referring to FIG. 7 (refer to FIG. 3 as appropriate), the attenuation attenuation data 3 indicates the attenuation distribution data indicating the probability distribution of the attenuation of the radio wave of the satellite broadcast with respect to a predetermined one-hour rainfall for each observation point. The operation of generating FIG. 7 is a flowchart showing the operation of the rain attenuation analysis apparatus according to the first embodiment of the present invention.

(降雨強度分布解析ステップ)
降雨減衰解析装置3は、地点別雨量データ入力手段31によって、外部から地点別雨量データを入力する(ステップS31)。続いて、降雨減衰解析装置3は、降雨強度分布解析手段32によって、ステップS31において入力された地点別雨量データに基づいて、観測地点ごとに1時間雨量を算出し、1時間雨量が所定値となる10分間降雨強度の確率分布を対数正規分布で近似して降雨強度分布データを生成する(ステップS32)。
(Rainfall intensity distribution analysis step)
The rain attenuation analysis device 3 inputs the rainfall data for each point from the outside by the rainfall data input means 31 for each point (step S31). Subsequently, the rainfall attenuation analysis device 3 calculates the rainfall for one hour for each observation point based on the point-by-point rainfall data input in step S31 by the rainfall intensity distribution analysis unit 32, and the one-hour rainfall is a predetermined value. The 10-minute rainfall intensity probability distribution is approximated by a lognormal distribution to generate rainfall intensity distribution data (step S32).

(減衰解析ステップ)
そして、降雨減衰解析装置3は、減衰解析手段34の降雨区間長算出手段341によって、ステップS32において10分間降雨強度が解析された観測地点の降雨区間長を算出する(ステップS33)。続いて、降雨減衰解析装置3は、減衰分布算出手段342によって、ステップS32において生成された降雨強度分布データと、ステップS33において算出された降雨区間長とに基づいて、降雨による減衰量の確率分布を解析し、平均値及び標準偏差を観測地点と1時間あたりの雨量とに対応させた減衰量分布データを生成する(ステップS34)。
(Attenuation analysis step)
Then, the rain attenuation analysis apparatus 3 calculates the rain section length of the observation point where the rain intensity is analyzed in step S32 by the rain section length calculation means 341 of the attenuation analysis means 34 (step S33). Subsequently, the rain attenuation analysis device 3 uses the attenuation distribution calculation means 342 to calculate the probability distribution of the attenuation due to the rain based on the rain intensity distribution data generated in step S32 and the rain section length calculated in step S33. Is generated, and attenuation distribution data in which the average value and the standard deviation are associated with the observation point and the rainfall per hour is generated (step S34).

更に、降雨減衰解析装置3は、減衰量分布データ出力手段35によって、ステップS34において生成された減衰量分布データをメッシュ減衰量解析装置5に出力する(ステップS35)。そして、降雨減衰解析装置3は、降雨強度分布解析手段32によって、すべての観測地点及び所定の1時間雨量について、減衰量分布データの生成を行ったかを判断する(ステップS36)。   Further, the rain attenuation analysis device 3 outputs the attenuation amount distribution data generated in step S34 to the mesh attenuation amount analysis device 5 by the attenuation amount distribution data output means 35 (step S35). Then, the rain attenuation analysis device 3 determines whether or not the attenuation distribution data has been generated for all the observation points and the predetermined one-hour rainfall by the rainfall intensity distribution analysis means 32 (step S36).

そして、すべての観測地点及び所定の1時間雨量について終了していない場合(ステップS36でNo)には、ステップS32に戻って、降雨減衰解析装置3が、降雨強度分布解析手段32によって、次の1時間雨量あるいは観測地点について、降雨強度分布データを生成する動作以降の動作を行う。また、すべての観測地点及び所定の1時間雨量について終了した場合(ステップS36でYes)には、動作を終了する。   And when it is not complete | finished about all the observation points and predetermined 1 hour rainfall (it is No at step S36), it returns to step S32 and the rain attenuation analysis apparatus 3 performs the following by rain intensity distribution analysis means 32. For the hourly rainfall or the observation point, the operation after the operation for generating the rainfall intensity distribution data is performed. In addition, when all the observation points and the predetermined one-hour rainfall are finished (Yes in step S36), the operation is finished.

(メッシュ減衰量解析装置の動作)
次に、図8を参照(適宜図3参照)して、メッシュ減衰量解析装置5が、外部から入力された雨量データに基づいて、衛星放送の電波の送信電力の降雨マージンを算出する動作について説明する。図8は、本発明の第一の実施の形態におけるメッシュ減衰量解析装置の動作を示したフローチャートである。
(Operation of mesh attenuation analyzer)
Next, referring to FIG. 8 (refer to FIG. 3 as appropriate), the mesh attenuation amount analysis apparatus 5 calculates the rain margin of the transmission power of the radio wave of the satellite broadcast based on the rainfall data input from the outside. explain. FIG. 8 is a flowchart showing the operation of the mesh attenuation amount analyzing apparatus according to the first embodiment of the present invention.

なお、ここでは、メッシュ減衰量解析装置5が、減衰量分布データ入力手段51によって、複数の観測地点及び1時間雨量についての減衰量分布データを降雨減衰解析装置3から入力して、既に減衰量分布データ記憶手段52に記憶していることとし、これ以降の動作について説明する。   Here, the mesh attenuation amount analysis device 5 inputs attenuation amount distribution data for a plurality of observation points and one-hour rainfall from the rain attenuation analysis device 3 by the attenuation amount distribution data input means 51, and the attenuation amount has already been obtained. It is assumed that the data is stored in the distribution data storage means 52, and the subsequent operation will be described.

(雨量データ入力手段ステップ)
メッシュ減衰量解析装置5は、雨量データ入力手段53によって、外部から雨量データを入力する。また、メッシュ減衰量解析装置5は、要求放送時間データ入力手段55によって、外部から要求放送時間データを入力する(ステップS51)。
(Rainfall data input means step)
The mesh attenuation amount analysis device 5 inputs rain data from the outside by the rain data input means 53. Further, the mesh attenuation amount analyzing apparatus 5 inputs the requested broadcast time data from the outside by the requested broadcast time data input means 55 (step S51).

(減衰量分布データ読み出しステップ)
そして、メッシュ減衰量解析装置5は、減衰量分布データ読み出し手段54によって、メッシュごとに、減衰量分布データ記憶手段52から、メッシュに予め対応付けられた観測地点と、ステップS51において入力された雨量データによって示される当該メッシュの1時間雨量とに対応する減衰量分布データを読み出す(ステップS52)。
(Attenuation distribution data reading step)
Then, the mesh attenuation amount analyzing apparatus 5 uses the attenuation amount distribution data reading unit 54 for each mesh, from the attenuation amount distribution data storage unit 52 to the observation point previously associated with the mesh, and the rainfall input in step S51. Attenuation distribution data corresponding to the hourly rainfall of the mesh indicated by the data is read (step S52).

(電波減衰解析ステップ)
更に、メッシュ減衰量解析装置5は、電波減衰解析手段56によって、ステップS51において入力された要求放送時間データによって示される超過確率の許容値に基づいて、メッシュごとに、ステップS52において読み出された減衰量分布データにおいて、超過確率がこの許容値となる減衰量を解析する(ステップS53)。
(Radio attenuation analysis step)
Further, the mesh attenuation amount analyzing device 5 is read out in step S52 for each mesh by the radio wave attenuation analyzing unit 56 based on the allowable value of the excess probability indicated by the requested broadcast time data input in step S51. In the attenuation amount distribution data, an attenuation amount whose excess probability becomes the allowable value is analyzed (step S53).

そして、メッシュ減衰量解析装置5は、電波減衰解析手段56によって、ステップS51において入力された雨量データに示されたすべてのメッシュについて、超過確率の許容値に対応する減衰量の解析を行ったかを判断する(ステップS54)。そして、すべてのメッシュについて終了していない場合(ステップS54でNo)には、ステップS52に戻って、メッシュ減衰量解析装置5が、減衰量分布データ読み出し手段54によって、減衰量分布データ記憶手段52から、次のメッシュについての減衰量分布データを読み出す動作以降の動作を行う。   Then, the mesh attenuation amount analysis device 5 has analyzed whether the attenuation amount corresponding to the allowable value of the excess probability has been analyzed for all the meshes indicated in the rainfall amount data input in step S51 by the radio wave attenuation analysis unit 56. Judgment is made (step S54). If the processing has not been completed for all meshes (No in step S54), the process returns to step S52, and the mesh attenuation amount analysis device 5 causes the attenuation amount distribution data reading unit 54 to use the attenuation amount distribution data storage unit 52. Then, the operation after the operation of reading the attenuation distribution data for the next mesh is performed.

また、すべてのメッシュについて終了した場合(ステップS54でYes)には、メッシュ減衰量解析装置5は、降雨マージン算出手段58によって、ステップS53において解析された減衰量に基づいてメッシュごとに必要な送信電力を算出し、フェーズドアレーアンテナAのホーンh、h、…(図2参照)によって電波が送信される各々の範囲に対応する複数のメッシュに対して算出された複数の送信電力値の中の最大値を、該当する送信範囲に割り当て、電波の送信電力として降雨マージンを算出する(ステップS55)。   When all the meshes have been completed (Yes in step S54), the mesh attenuation amount analysis apparatus 5 uses the rain margin calculating means 58 to transmit necessary transmissions for each mesh based on the attenuation amount analyzed in step S53. Among the plurality of transmission power values calculated for the plurality of meshes corresponding to the respective ranges in which the radio waves are transmitted by the horns h, h,... (See FIG. 2) of the phased array antenna A. The maximum value is assigned to the corresponding transmission range, and the rain margin is calculated as the radio wave transmission power (step S55).

そして、メッシュ減衰量解析装置5は、降雨マージン出力手段59によって、ステップS55において生成された降雨マージンを出力し(ステップS56)、動作を終了する。   Then, the mesh attenuation amount analysis apparatus 5 outputs the rain margin generated in step S55 by the rain margin output means 59 (step S56), and ends the operation.

[降雨マージン算出装置の構成(第二の実施の形態)]
次に、図9(適宜図3参照)を参照して、本発明における第二の実施の形態である降雨マージン算出装置1Aの構成について説明する。図9は、本発明における第二の実施の形態である降雨マージン算出装置の構成を示したブロック図である。図9に示すように、降雨マージン算出装置1Aは、外部から地点別雨量データを入力し、この観測地点ごとの降雨強度の確率分布を解析するとともに、この降雨強度の確率分布を示す降雨強度分布データと、外部から入力された雨量データと、1時間あたりで放送が遮断されない時間を示す要求放送時間データとに基づいて、降雨マージンを算出するものである。降雨マージン算出装置1Aは、降雨減衰解析装置3Aと、メッシュ減衰量解析装置5Aとを備える。
[Configuration of Rain Margin Calculation Device (Second Embodiment)]
Next, with reference to FIG. 9 (refer to FIG. 3 as appropriate), the configuration of the rain margin calculating apparatus 1A according to the second embodiment of the present invention will be described. FIG. 9 is a block diagram showing a configuration of a rain margin calculating apparatus according to the second embodiment of the present invention. As shown in FIG. 9, the rain margin calculating apparatus 1A inputs the rainfall data by point from the outside, analyzes the probability distribution of the rainfall intensity at each observation point, and also displays the rainfall intensity distribution indicating the probability distribution of the rainfall intensity. The rainfall margin is calculated based on the data, the rainfall data input from the outside, and the requested broadcast time data indicating the time during which broadcasting is not interrupted per hour. The rainfall margin calculation device 1A includes a rain attenuation analysis device 3A and a mesh attenuation amount analysis device 5A.

(降雨減衰解析装置の構成)
降雨減衰解析装置3Aは、地点別雨量データに基づいて、各々の観測地点における10分間降雨強度の確率分布を解析するものである。この降雨減衰解析装置3Aは、降雨減衰解析装置3の解析データ記憶手段33と、減衰解析手段34と、減衰量分布データ出力手段35とを備えず、降雨強度分布解析手段32の代わりに降雨強度分布解析手段32Aを備え、更に、降雨強度分布データ出力手段36Aを備える。なお、降雨減衰解析装置3A内の地点別雨量データ入力手段31は、図3に示したものと同一であるので、同一の符号を付し、説明を省略する。
(Configuration of rain attenuation analyzer)
The rain attenuation analysis device 3A analyzes the probability distribution of the 10-minute rainfall intensity at each observation point based on the point-by-point rainfall data. This rain attenuation analysis apparatus 3A does not include the analysis data storage means 33, the attenuation analysis means 34, and the attenuation amount distribution data output means 35 of the rain attenuation analysis apparatus 3, and instead of the rainfall intensity distribution analysis means 32, the rainfall intensity A distribution analysis means 32A is provided, and a rainfall intensity distribution data output means 36A is further provided. Since the point-by-point rainfall data input means 31 in the rain attenuation analysis apparatus 3A is the same as that shown in FIG. 3, the same reference numerals are given and description thereof is omitted.

降雨強度分布解析手段32Aは、地点別雨量データ入力手段31から入力された地点別雨量データに基づいて、1時間雨量が所定値となる、観測地点ごとの10分間降雨強度の確率分布を対数正規分布で近似して、当該確率分布を示す降雨強度分布データを生成するものである。ここでは、10分間降雨強度の確率分布の平均値と標準偏差を降雨強度布データとし、観測地点及び1時間雨量と対応させることとした。このようにして生成された降雨強度分布データは、降雨強度分布データ出力手段36Aに出力される。   The rainfall intensity distribution analysis unit 32A logarithmically normalizes the probability distribution of the 10-minute rainfall intensity for each observation point where the hourly rainfall is a predetermined value based on the point-by-point rainfall data input from the point-by-point rainfall data input unit 31. The rainfall intensity distribution data representing the probability distribution is generated by approximating the distribution. Here, the average value and standard deviation of the probability distribution of the 10-minute rainfall intensity are used as the rainfall intensity distribution data, and are associated with the observation point and the 1-hour rainfall. The rainfall intensity distribution data generated in this way is output to the rainfall intensity distribution data output means 36A.

降雨強度分布データ出力手段36Aは、降雨強度分布解析手段32Aから入力された降雨強度分布データを外部へ出力するものである。なお、ここでは、降雨強度分布データ出力手段36Aは、降雨強度分布データをメッシュ減衰量解析装置5Aに出力することとした。なお、この降雨強度分布データ出力手段36Aは、減衰量分布データ出力手段35(図3)に比べて出力するデータを減衰量分布データから降雨強度分布データとしただけで機能は同じものである。   The rainfall intensity distribution data output unit 36A outputs the rainfall intensity distribution data input from the rainfall intensity distribution analysis unit 32A to the outside. Here, the rain intensity distribution data output means 36A outputs the rain intensity distribution data to the mesh attenuation amount analyzer 5A. The rainfall intensity distribution data output means 36A has the same function as the output data compared with the attenuation distribution data output means 35 (FIG. 3) except that the rainfall intensity distribution data is changed from the attenuation distribution data.

(メッシュ減衰量解析装置の構成)
メッシュ減衰量解析装置5Aは、外部から入力される、メッシュごとの1時間雨量の予測値を示す雨量データと、降雨減衰解析装置3Aによって算出された降雨強度分布データとに基づいて、衛星放送の電波の減衰量を解析し、降雨による減衰を補償する降雨マージンを算出するものである。このメッシュ減衰量解析装置5Aは、減衰量分布データ入力手段51と、減衰量分布データ記憶手段52と、減衰量分布データ読み出し手段54とを備えず、雨量データ入力手段53と、電波減衰解析手段56との代わりに雨量データ入力手段53Aと、電波減衰解析手段56Aとを備え、更に、降雨強度分布データ入力手段61Aと、降雨強度分布データ記憶手段62Aと、降雨強度分布データ読み出し手段64Aと、解析データ記憶手段33Aと、減衰解析手段34Aとを備える。なお、メッシュ減衰量解析装置5A内の降雨強度分布データ入力手段61A、降雨強度分布データ記憶手段62A、雨量データ入力手段53A、降雨強度分布データ読み出し手段64A、解析データ記憶手段33A及び減衰解析手段34A、電波減衰解析手段56A以外の構成は、図3に示したものと同一であるので、同一の符号を付し、説明を省略する。
(Configuration of mesh attenuation analysis device)
The mesh attenuation amount analysis apparatus 5A is based on the rainfall data indicating the predicted value of rainfall for one hour for each mesh and the rainfall intensity distribution data calculated by the rain attenuation analysis apparatus 3A. It analyzes the amount of attenuation of radio waves and calculates a rain margin that compensates for attenuation due to rain. This mesh attenuation amount analyzing apparatus 5A does not include the attenuation amount distribution data input means 51, the attenuation amount distribution data storage means 52, and the attenuation amount distribution data reading means 54, but the rain amount data input means 53 and the radio wave attenuation analysis means. In addition to 56, a rainfall data input means 53A and a radio wave attenuation analysis means 56A are provided. Further, a rainfall intensity distribution data input means 61A, a rainfall intensity distribution data storage means 62A, a rainfall intensity distribution data reading means 64A, Analysis data storage means 33A and attenuation analysis means 34A are provided. It should be noted that the rainfall intensity distribution data input means 61A, the rainfall intensity distribution data storage means 62A, the rainfall data input means 53A, the rainfall intensity distribution data read means 64A, the analysis data storage means 33A, and the attenuation analysis means 34A in the mesh attenuation amount analysis apparatus 5A. Since the configuration other than the radio wave attenuation analysis means 56A is the same as that shown in FIG. 3, the same reference numerals are given and the description thereof is omitted.

降雨強度分布データ入力手段61Aは、降雨減衰解析装置3Aの降雨強度分布解析手段32Aによって生成された降雨強度分布データを降雨強度分布データ出力手段36Aから入力し、降雨強度分布データ記憶手段62Aに記憶するものである。なお、この降雨強度分布データ入力手段61Aは、減衰量分布データ入力手段51(図3)に比べて入力するデータを減衰量分布データから降雨強度分布データとしただけで機能は同じものである。   The rainfall intensity distribution data input means 61A inputs the rain intensity distribution data generated by the rainfall intensity distribution analysis means 32A of the rain attenuation analyzer 3A from the rain intensity distribution data output means 36A and stores it in the rain intensity distribution data storage means 62A. To do. The function of the rainfall intensity distribution data input means 61A is the same as that of the attenuation distribution data input means 51 (FIG. 3) except that the input data is changed from the attenuation distribution data to the rainfall intensity distribution data.

降雨強度分布データ記憶手段(降雨強度分布データ記憶装置)62Aは、降雨減衰解析装置3Aによって解析された降雨強度分布データを観測地点及び1時間雨量と対応付けて記憶するもので、半導体メモリ、ハードディスク等の一般的な記憶手段である。この降雨強度分布データは、後記する降雨強度分布データ読み出し手段64Aによって読み出され、減衰解析手段34Aによって、観測地点ごとの衛星放送の電波の減衰量の確率分布を解析する際に参照されて用いられる。   The rainfall intensity distribution data storage means (rain intensity distribution data storage device) 62A stores the rainfall intensity distribution data analyzed by the rain attenuation analysis device 3A in association with the observation point and the hourly rainfall amount. It is a general storage means such as. This rainfall intensity distribution data is read by the rain intensity distribution data reading means 64A described later, and is used by referring to the attenuation analysis means 34A when analyzing the probability distribution of the attenuation amount of the satellite broadcast radio wave for each observation point. It is done.

雨量データ入力手段53Aは、メッシュごとの1時間雨量を示す雨量データを外部から入力するものである。ここで入力された雨量データは、降雨強度分布データ読み出し手段64Aに出力される。なお、この降雨強度分布データ出力手段36Aは、減衰量分布データ出力手段35(図3)に比べて、雨量データの出力先を減衰量分布データ読み出し手段54から降雨強度分布データ読み出し手段64Aとしただけで機能は同じものである。   The rainfall data input means 53A is for inputting rainfall data indicating an hourly rainfall for each mesh from the outside. The rainfall data input here is output to the rainfall intensity distribution data reading means 64A. In this rain intensity distribution data output means 36A, as compared with the attenuation distribution data output means 35 (FIG. 3), the output destination of the rain data is the attenuation intensity distribution data read means 54 to the rain intensity distribution data read means 64A. Just the function is the same.

降雨強度分布データ読み出し手段64Aは、雨量データ入力手段53Aから入力された雨量データに基づいて、降雨強度分布データ記憶手段62Aから、各々のメッシュに予め対応付けられた観測地点における、雨量データによって示される1時間雨量に対応する降雨強度分布データを読み出すものである。ここで読み出された減衰量分布データは、減衰解析手段34Aに出力される。   The rainfall intensity distribution data reading means 64A is indicated by the rainfall data at the observation points previously associated with the respective meshes from the rainfall intensity distribution data storage means 62A based on the rainfall data input from the rainfall data input means 53A. The rainfall intensity distribution data corresponding to the one hour rainfall is read out. The attenuation distribution data read here is output to the attenuation analysis means 34A.

解析データ記憶手段33Aは、後記する減衰解析手段34Aによって減衰量の確率分布を算出するために必要となるデータを予め記憶するもので、半導体メモリ、ハードディスク等の一般的な記憶手段である。ここでは、解析データ記憶手段33Aに、観測地点位置データ、衛星軌道位置データ、雨滴層高度データ及び周波数パラメータを記憶することとした。なお、この解析データ記憶手段33Aは、解析データ記憶手段33(図3)に比べてメッシュ減衰量解析装置5Aに含まれることとしただけで記憶されるデータは同じものである。   The analysis data storage unit 33A stores in advance data necessary for calculating the probability distribution of the attenuation amount by the attenuation analysis unit 34A described later, and is a general storage unit such as a semiconductor memory or a hard disk. Here, the observation point position data, satellite orbit position data, raindrop height data, and frequency parameters are stored in the analysis data storage means 33A. It should be noted that the analysis data storage means 33A is the same as the data stored only by being included in the mesh attenuation amount analysis device 5A as compared with the analysis data storage means 33 (FIG. 3).

減衰解析手段34Aは、降雨強度分布データ読み出し手段64Aによって読み出された降雨強度分布データに基づいて、減衰量の確率分布を示す減衰量分布データを生成するものである。ここでは、減衰解析手段34Aは、降雨区間長算出手段341と、減衰分布算出手段342Aとを備える。なお、降雨区間長算出手段341は、図3に示したものと同一であるので、同一の符号を付し、説明を省略する。   The attenuation analysis unit 34A generates attenuation amount distribution data indicating the probability distribution of the attenuation amount based on the rainfall intensity distribution data read by the rainfall intensity distribution data reading unit 64A. Here, the attenuation analysis unit 34A includes a rainfall section length calculation unit 341 and an attenuation distribution calculation unit 342A. In addition, since the rainfall section length calculation means 341 is the same as that shown in FIG. 3, the same reference numerals are given and the description thereof is omitted.

減衰分布算出手段342Aは、降雨区間長算出手段341から入力される降雨区間長と、降雨強度分布データ読み出し手段64Aから入力される降雨強度分布データとに基づいて、観測地点ごとの降雨による減衰量の確率分布を解析するものである。ここで解析された確率分布を示す減衰量分布データは、減衰量分布データ出力手段35に出力される。   The attenuation distribution calculating means 342A is based on the rainfall section length input from the rainfall section length calculating means 341 and the rain intensity distribution data input from the rainfall intensity distribution data reading means 64A, and the amount of attenuation due to rain for each observation point. The probability distribution is analyzed. The attenuation distribution data indicating the probability distribution analyzed here is output to the attenuation distribution data output means 35.

電波減衰解析手段56Aは、メッシュごとに、減衰解析手段34Aによって生成された減衰量分布データによって示される減衰量の確率分布において、要求放送時間データ入力手段55から入力される要求放送時間データによって示される超過確率の許容値に対応する減衰量を解析するものである。この減衰量は、降雨マージン算出手段58に出力される。なお、この電波減衰解析手段56Aは、電波減衰解析手段56(図3)に比べて、減衰量分布データの入力元を減衰量分布データ読み出し手段54から減衰解析手段34Aとしただけで機能は同じものである。   The radio wave attenuation analysis means 56A is indicated by the requested broadcast time data input from the requested broadcast time data input means 55 in the probability distribution of attenuation indicated by the attenuation distribution data generated by the attenuation analysis means 34A for each mesh. The amount of attenuation corresponding to the permissible excess probability is analyzed. This attenuation amount is output to the rain margin calculating means 58. The radio wave attenuation analysis unit 56A has the same function as the radio wave attenuation analysis unit 56 (FIG. 3) except that the attenuation distribution data is input from the attenuation distribution data reading unit 54 to the attenuation analysis unit 34A. Is.

これによって、降雨マージン算出装置1Aは、メッシュごとの1時間雨量の予測値に基づいて、1時間の放送時間のうち要求放送時間データによって示される時間の再生を可能とするための降雨マージンを算出することができる。   Thereby, the rain margin calculating apparatus 1A calculates a rain margin for enabling reproduction of the time indicated by the requested broadcast time data out of the broadcast time of 1 hour based on the predicted value of the 1 hour rainfall for each mesh. can do.

なお、特許請求の範囲に記載の区分別減衰量解析装置は、少なくとも雨量データ入力手段53Aと、降雨強度分布データ読み出し手段64Aと、減衰解析手段34Aと、電波減衰解析手段56Aとを備えていればよい。   The classification-specific attenuation amount analyzing apparatus described in the claims is provided with at least rainfall data input means 53A, rainfall intensity distribution data reading means 64A, attenuation analysis means 34A, and radio wave attenuation analysis means 56A. That's fine.

また、降雨減衰解析装置3A及びメッシュ減衰量解析装置5Aは、コンピュータにおいて各手段を各機能プログラムとして実現することも可能であり、各機能プログラムを結合して、降雨減衰解析プログラム及びメッシュ別減衰量解析プログラムとして動作させることも可能である。   Further, the rain attenuation analysis device 3A and the mesh attenuation amount analysis device 5A can also realize each means as each function program in a computer, and combine the function programs to obtain the rain attenuation analysis program and the mesh-specific attenuation amount. It is also possible to operate as an analysis program.

[降雨マージン算出装置の動作(第二の実施の形態)]
次に、図10(適宜図9参照)を参照して、本発明の第二の実施の形態である降雨マージン算出装置1Aの動作について説明する。図10は、本発明の第二の実施の形態におけるメッシュ減衰量解析装置の動作を示したフローチャートである。なお、ここでは、降雨減衰解析装置3Aが、複数の観測地点及び1時間雨量についての降雨強度分布データを生成し、メッシュ減衰量解析装置5Aが、降雨強度分布データ入力手段61Aによってこの降雨強度分布データを入力して、降雨強度分布データ記憶手段62Aに記憶していることとし、これ以降のメッシュ減衰量解析装置5Aの動作について説明する。
[Operation of Rain Margin Calculation Device (Second Embodiment)]
Next, with reference to FIG. 10 (see FIG. 9 as appropriate), the operation of the rain margin calculating apparatus 1A according to the second embodiment of the present invention will be described. FIG. 10 is a flowchart showing the operation of the mesh attenuation amount analyzing apparatus according to the second embodiment of the present invention. Here, the rain attenuation analysis device 3A generates rain intensity distribution data for a plurality of observation points and one-hour rainfall, and the mesh attenuation analysis device 5A uses the rain intensity distribution data input means 61A to generate this rain intensity distribution data. Assume that data is input and stored in the rainfall intensity distribution data storage means 62A, and the subsequent operation of the mesh attenuation amount analysis apparatus 5A will be described.

(雨量データ入力手段ステップ)
メッシュ減衰量解析装置5Aは、雨量データ入力手段53Aによって、外部から雨量データを入力する。また、メッシュ減衰量解析装置5Aは、要求放送時間データ入力手段55によって、外部から要求放送時間データを入力する(ステップS61)。
(Rainfall data input means step)
The mesh attenuation amount analysis apparatus 5A inputs rain amount data from the outside by the rain amount data input means 53A. Further, the mesh attenuation amount analyzing apparatus 5A inputs the requested broadcast time data from the outside by the requested broadcast time data input means 55 (step S61).

(降雨強度分布データ読み出しステップ)
そして、メッシュ減衰量解析装置5Aは、降雨強度分布データ読み出し手段64Aによって、メッシュごとに、降雨強度分布データ記憶手段62Aから、メッシュに予め対応付けられた観測地点と、ステップS61において入力された雨量データによって示される当該メッシュの1時間雨量とに対応する降雨強度分布データを読み出す(ステップS62)。
(Step of reading rainfall intensity distribution data)
Then, the mesh attenuation amount analyzing apparatus 5A uses the rainfall intensity distribution data reading unit 64A to, for each mesh, the observation point previously associated with the mesh from the rain intensity distribution data storage unit 62A and the rainfall input in step S61. Rain intensity distribution data corresponding to the hourly rainfall of the mesh indicated by the data is read (step S62).

(減衰解析ステップ)
そして、メッシュ減衰量解析装置5Aは、減衰解析手段34Aの降雨区間長算出手段341によって、ステップS62において読み出された降雨強度分布データの観測地点における降雨区間長を算出する(ステップS63)。続いて、メッシュ減衰量解析装置5Aは、減衰分布算出手段342Aによって、ステップS62において読み出された降雨強度分布データと、ステップS63において算出された降雨区間長とに基づいて、降雨による減衰量の確率分布を解析し、平均値及び標準偏差を観測地点と1時間あたりの雨量とに対応させた減衰量分布データを生成する(ステップS64)。
(Attenuation analysis step)
Then, the mesh attenuation amount analysis apparatus 5A calculates the rain section length at the observation point of the rain intensity distribution data read in step S62 by the rain section length calculation means 341 of the attenuation analysis means 34A (step S63). Subsequently, the mesh attenuation amount analyzing apparatus 5A uses the attenuation distribution calculating unit 342A to calculate the attenuation amount due to rainfall based on the rainfall intensity distribution data read in step S62 and the rainfall section length calculated in step S63. The probability distribution is analyzed, and attenuation distribution data in which the average value and the standard deviation are associated with the observation point and the rainfall per hour is generated (step S64).

(電波減衰解析ステップ)
更に、メッシュ減衰量解析装置5Aは、電波減衰解析手段56Aによって、ステップS61において入力された要求放送時間データによって示される超過確率の許容値に基づいて、メッシュごとに、ステップS64において生成された減衰量分布データにおいて、超過確率がこの許容値となる減衰量を解析する(ステップS65)。
(Radio attenuation analysis step)
Further, the mesh attenuation amount analyzing apparatus 5A uses the attenuation generated in step S64 for each mesh by the radio wave attenuation analyzing unit 56A based on the allowable excess probability indicated by the requested broadcast time data input in step S61. In the quantity distribution data, an attenuation amount whose excess probability becomes the allowable value is analyzed (step S65).

そして、メッシュ減衰量解析装置5Aは、電波減衰解析手段56Aによって、ステップS61において入力された雨量データに示されたすべてのメッシュについて、超過確率の許容値に対応する減衰量の解析を行ったかを判断する(ステップS66)。そして、すべてのメッシュについて終了していない場合(ステップS66でNo)には、ステップS62に戻って、メッシュ減衰量解析装置5Aが、降雨強度分布データ読み出し手段64Aによって、降雨強度分布データ記憶手段62Aから、次のメッシュについての降雨強度分布データを読み出す動作以降の動作を行う。   Then, the mesh attenuation amount analysis device 5A determines whether the attenuation amount corresponding to the allowable value of the excess probability has been analyzed for all the meshes indicated in the rainfall amount data input in step S61 by the radio wave attenuation analysis unit 56A. Judgment is made (step S66). If all the meshes have not been completed (No in step S66), the process returns to step S62, and the mesh attenuation amount analysis device 5A causes the rainfall intensity distribution data storage unit 62A to perform the rainfall intensity distribution data reading unit 64A. Then, the operation after the operation of reading the rainfall intensity distribution data for the next mesh is performed.

また、すべてのメッシュについて終了した場合(ステップS66でYes)には、メッシュ減衰量解析装置5Aは、降雨マージン算出手段58によって、ステップS65において解析された減衰量に基づいてメッシュごとに必要な送信電力を算出し、フェーズドアレーアンテナAのホーンh、h、…(図2参照)によって電波が送信される各々の範囲に対応する複数のメッシュに対して算出された複数の送信電力の中の最大値を、該当する送信範囲に割り当てる電波の送信電力として、降雨マージンを算出する。(ステップS67)。   When the processing is completed for all the meshes (Yes in step S66), the mesh attenuation amount analyzing apparatus 5A transmits the necessary transmission for each mesh based on the attenuation amount analyzed in step S65 by the rain margin calculating unit 58. The power is calculated, and the maximum of the plurality of transmission powers calculated for the plurality of meshes corresponding to the respective ranges in which the radio waves are transmitted by the horns h, h,... (See FIG. 2) of the phased array antenna A The rain margin is calculated using the value as the transmission power of the radio wave assigned to the corresponding transmission range. (Step S67).

そして、メッシュ減衰量解析装置5Aは、降雨マージン出力手段59によって、ステップS67において生成された降雨マージンを出力し(ステップS68)、動作を終了する。   Then, the mesh attenuation amount analysis apparatus 5A outputs the rain margin generated in step S67 by the rain margin output means 59 (step S68), and ends the operation.

本発明における第一の実施の形態である降雨マージン算出装置による降雨減衰補償の概要を模式的に示す説明図である。It is explanatory drawing which shows typically the outline | summary of the rain attenuation compensation by the rain margin calculation apparatus which is 1st embodiment in this invention. 本発明における第一の実施の形態である降雨マージン算出装置によって算出された降雨マージンを与えて、降雨減衰補償を行う例を示す説明図、(a)は、フェーズドアレーアンテナによって、電波を送受信する地域に降雨マージンを与えた電波を送信し、降雨の集中する一部の範囲に特に高い強度の電波を送信する例を示す模式図、(b)は、(a)において送信される電波の強さを表すグラフである。BRIEF DESCRIPTION OF THE DRAWINGS FIG. 1 is an explanatory diagram showing an example in which a rain margin calculated by a rain margin calculating apparatus according to a first embodiment of the present invention is given to compensate for rain attenuation, (a) transmitting and receiving radio waves by a phased array antenna. A schematic diagram showing an example of transmitting a radio wave with a rain margin to a region and transmitting a particularly high intensity radio wave in a part of a region where rainfall is concentrated, (b) is a diagram showing the strength of the radio wave transmitted in (a) It is a graph showing thickness. 本発明における第一の実施の形態である降雨マージン算出装置の構成を示したブロック図である。It is the block diagram which showed the structure of the rain margin calculation apparatus which is 1st embodiment in this invention. 東京及び鹿児島の2つの観測地点において観測された、1時間雨量が5mmとなる降雨の1分間降雨強度及び10分間降雨強度の確率分布を対数正規確率紙にプロットしたグラフである。It is the graph which plotted the probability distribution of the 1-minute rainfall intensity and the 10-minute rainfall intensity of 1-minute rainfall observed at two observation points in Tokyo and Kagoshima on a lognormal probability paper. 本発明の第一の実施の形態における降雨減衰解析装置の降雨区間長算出手段によって算出される衛星放送の電波の降雨区間長と、雨滴層高度とを模式的に示した説明図である。It is explanatory drawing which showed typically the rainfall area length of the radio wave of a satellite broadcast calculated by the rainfall area length calculation means of the rain attenuation analysis apparatus in 1st embodiment of this invention, and raindrop layer height. 東京において観測された、1時間雨量が5mm及び10mmとなる降雨の1分間降雨強度及び減衰量の確率分布を対数正規確率紙にプロットしたグラフである。It is the graph which plotted the probability distribution of the 1-minute rainfall intensity and the attenuation amount of the rainfall observed in Tokyo in which the hourly rainfall is 5 mm and 10 mm on a logarithmic normal probability paper. 本発明の第一の実施の形態における降雨減衰解析装置の動作を示したフローチャートである。It is the flowchart which showed operation | movement of the rain attenuation analysis apparatus in 1st embodiment of this invention. 本発明の第一の実施の形態におけるメッシュ減衰量解析装置の動作を示したフローチャートである。It is the flowchart which showed operation | movement of the mesh attenuation amount analysis apparatus in 1st embodiment of this invention. 本発明における第二の実施の形態である降雨マージン算出装置の構成を示したブロック図である。It is the block diagram which showed the structure of the rain margin calculation apparatus which is 2nd embodiment in this invention. 本発明の第二の実施の形態におけるメッシュ減衰量解析装置の動作を示したフローチャートである。It is the flowchart which showed operation | movement of the mesh attenuation amount analysis apparatus in 2nd embodiment of this invention.

符号の説明Explanation of symbols

3、3A 降雨減衰解析装置
31 地点別雨量データ入力手段
32、32A 降雨強度分布解析手段
34、34A 減衰解析手段
52 減衰量分布データ記憶手段(減衰量分布データ記憶装置)
53、53A 雨量データ入力手段
54 減衰量分布データ読み出し手段
56、56A 電波減衰解析手段
62A 降雨強度分布データ記憶手段(降雨強度分布データ記憶装置)
64A 降雨強度分布データ読み出し手段
3, 3A Rain attenuation analysis device 31 Point-by-point rainfall data input means 32, 32A Rain intensity distribution analysis means 34, 34A Attenuation analysis means 52 Attenuation distribution data storage means (attenuation distribution data storage device)
53, 53A Rainfall data input means 54 Attenuation distribution data reading means 56, 56A Radio wave attenuation analysis means 62A Rainfall intensity distribution data storage means (rainfall intensity distribution data storage device)
64A Rain intensity distribution data reading means

Claims (5)

ある観測地点について観測された10分間ごとの降雨強度のうち、間における雨量が所定値となる前記降雨強度を示す地点別雨量データに基づいて、前記観測地点における前記時間あたりの雨量の降雨による電波の減衰量の確率分布を解析する降雨減衰解析装置であって、
前記地点別雨量データによって示される前記降雨強度の確率分布を対数正規分布で近似して、前記確率分布を示す降雨強度分布データを生成する降雨強度分布解析手段と、
この降雨強度分布解析手段によって生成された降雨強度分布データである10分間降雨強度の確率分布の平均値及び標準偏差を、1分間降雨強度の確率分布の平均値及び標準偏差とみなして、前記時間あたりの雨量の降雨における前記減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する減衰解析手段と
を備え
前記降雨強度分布解析手段は、
前記観測地点ごとに継続して観測された6つの10分間雨量を加算して1時間雨量を算出し、当該1時間雨量が前記所定値となる前記6つの10分間雨量の組を複数取り出し、この取り出した10分間雨量の各値を6倍して10分間降雨強度を算出し、前記取り出した複数の組の各々に対して10分間降雨強度の値を昇順に並び替え、値が小さいものから順にその値をこえる確率を、第1超過確率、第2超過確率、第3超過確率、第4超過確率、及び、第5超過確率とし、この複数の組にわたる各超過確率における10分間降雨強度の中央値を、当該超過確率における10分間降雨強度とし、この5つの超過確率から導かれる10分間降雨強度分布から当該分布に従う10分間降雨強度の平均値と標準偏差とを算出して、前記降雨強度分布データとする
ことを特徴とする降雨減衰解析装置。
Of rainfall every 10 minutes observed for an observation point, based on the locational rainfall data indicating the rain intensity rainfall definitive between 1 hour becomes a predetermined value, rainfall of the per hour in the observation point A rain attenuation analysis device that analyzes the probability distribution of radio attenuation due to rain
A rainfall intensity distribution analyzing means for approximating the probability distribution of the rainfall intensity indicated by the point-by-point rainfall data with a lognormal distribution to generate rainfall intensity distribution data indicating the probability distribution;
The mean and standard deviation of the probability distribution of the 10 minute rainfall intensity is rainfall intensity distribution data generated by the rainfall intensity distribution analysis unit, is regarded as the mean and standard deviation of the probability distribution of 1 minute rain rate, the 1 Attenuation analysis means for analyzing the probability distribution of the attenuation amount in rainfall of rainfall per hour and generating attenuation distribution data indicating the probability distribution , and
The rainfall intensity distribution analysis means includes
By adding six 10-minute rainfalls observed continuously at each observation point to calculate one-hour rainfall, a plurality of pairs of the six 10-minute rainfalls in which the one-hour rainfall becomes the predetermined value are extracted. The 10-minute rainfall intensity is calculated by multiplying the taken 10-minute rainfall values by 6 times, and the 10-minute rainfall intensity values are rearranged in ascending order with respect to each of the plurality of the extracted sets. The probability exceeding the value is defined as the first excess probability, the second excess probability, the third excess probability, the fourth excess probability, and the fifth excess probability, and the center of the 10-minute rainfall intensity at each excess probability across the plurality of sets. The value is defined as the 10-minute rainfall intensity at the excess probability, and the average value and the standard deviation of the 10-minute rainfall intensity according to the distribution are calculated from the 10-minute rainfall intensity distribution derived from the five excess probabilities, and the rainfall intensity distribution is calculated. De Rain attenuation analyzing apparatus characterized by the data.
ある観測地点について観測された10分間ごとの降雨強度のうち、所定時間内における雨量が所定値となる前記降雨強度を示す地点別雨量データによって示される前記降雨強度の確率分布を対数正規分布で近似して、前記確率分布を示す降雨強度分布データを生成する降雨強度分布解析手段と、
この降雨強度分布解析手段によって生成された降雨強度分布データに基づいて、前記所定時間あたりの雨量の降雨における前記減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する減衰解析手段と
を備える降雨減衰解析装置によって生成された、ある観測地点において所定時間あたりの雨量が所定値となる降雨による電波の減衰量の確率分布を示す減衰量分布データに基づいて、前記減衰量の解析対象となる地域内の複数の区分の各々について降雨による前記電波の減衰量を解析する区分別減衰量解析装置であって、
前記減衰量分布データを前記観測地点及び前記雨量と対応付けて、複数の前記観測地点及び複数の前記雨量について記憶する減衰量分布データ記憶手段と、
前記区分ごとの所定時間あたりの雨量を示す雨量データを入力する雨量データ入力手段と、
前記減衰量分布データ記憶手段から、前記区分ごとに、当該区分に予め対応付けられた前記観測地点と、前記雨量データによって示される前記雨量とに対応する減衰量分布データを読み出す減衰量分布データ読み出し手段と、
この減衰量分布データ読み出し手段によって読み出された前記減衰量分布データによって示される前記減衰量の確率分布と、前記確率分布において前記減衰量がある値より大きな値をとる確率である超過確率の許容値を示す許容値データとに基づいて、前記区分ごとに前記超過確率が前記許容値となる減衰量を解析する電波減衰解析手段と
を備えることを特徴とする区分別減衰量解析装置。
Approximate the probability distribution of the rainfall intensity indicated by the rainfall data according to the point indicating the rainfall intensity within a predetermined time out of the rainfall intensity every 10 minutes observed for a certain observation point with a lognormal distribution. A rainfall intensity distribution analyzing means for generating rainfall intensity distribution data indicating the probability distribution;
Based on the rainfall intensity distribution data generated by the rainfall intensity distribution analyzing means, the attenuation for generating the attenuation distribution data indicating the probability distribution by analyzing the probability distribution of the attenuation in the rainfall of the predetermined amount of rainfall Analysis means and
Produced by rain attenuation analyzer comprising, there rainfall per predetermined time at the observation point is based on the attenuation amount distribution data indicating a probability distribution of a radio wave attenuation due to rainfall which becomes a predetermined value, the analysis target of the attenuation A category-specific attenuation analysis device for analyzing the attenuation of the radio wave due to rainfall for each of a plurality of categories in an area,
Attenuation distribution data storage means for storing the plurality of observation points and the plurality of rain amounts in association with the observation point and the rainfall amount,
Rainfall data input means for inputting rainfall data indicating rainfall per predetermined time for each of the categories;
Attenuation distribution data reading for reading out the attenuation distribution data corresponding to the observation points previously associated with the divisions and the rainfall indicated by the rainfall data from the attenuation distribution data storage means. Means,
Allowance of the probability distribution of the attenuation indicated by the attenuation distribution data read by the attenuation distribution data reading means and an excess probability that is a probability that the attenuation takes a value larger than a certain value in the probability distribution. A category-specific attenuation analysis device, comprising: radio wave attenuation analysis means for analyzing an attenuation amount at which the excess probability becomes the tolerance value for each category based on tolerance value data indicating a value.
ある観測地点について観測された所定時間内における雨量が所定値となる降雨の10分間ごとの降雨強度の確率分布を対数正規分布で近似して生成された、当該確率分布を示す降雨強度分布データを、前記観測地点及び前記雨量と対応付けて、複数の前記観測地点及び複数の前記雨量について記憶する降雨強度分布データ記憶装置から前記降雨強度分布データを読み出して、減衰量の解析対象となる地域内の複数の区分の各々について降雨による電波の減衰量を解析する区分別減衰量解析装置であって、
前記区分ごとの所定時間あたりの雨量を示す雨量データを入力する雨量データ入力手段と、
前記降雨強度分布データ記憶装置から、前記区分ごとに、当該区分に予め対応付けられた前記観測地点と、前記雨量データによって示される前記雨量とに対応する降雨強度分布データを読み出す降雨強度分布データ読み出し手段と、
この降雨強度分布データ読み出し手段によって読み出された降雨強度分布データに基づいて、前記所定時間あたりの雨量の降雨による前記電波の減衰量の確率分布を解析する減衰解析手段と、
この減衰解析手段によって解析された前記減衰量の確率分布と、前記確率分布において前記減衰量がある値より大きな値をとる確率である超過確率の許容値を示す許容値データとに基づいて、前記区分ごとに前記超過確率が前記許容値となる減衰量を解析する電波減衰解析手段と
を備えることを特徴とする区分別減衰量解析装置。
Rain intensity distribution data indicating the probability distribution generated by approximating the probability distribution of the rainfall intensity every 10 minutes of rainfall with a predetermined amount of rainfall observed for a certain observation point with a lognormal distribution. The rainfall intensity distribution data is read out from the rainfall intensity distribution data storage device that stores the plurality of observation points and the plurality of rainfalls in association with the observation point and the rainfall, and is stored in the region to be analyzed for attenuation. A category-specific attenuation analysis device for analyzing the attenuation of radio waves due to rainfall for each of a plurality of categories,
Rainfall data input means for inputting rainfall data indicating rainfall per predetermined time for each of the categories;
Read out rainfall intensity distribution data from the rain intensity distribution data storage device, for each of the sections, to read out rainfall intensity distribution data corresponding to the observation points previously associated with the sections and the rainfall indicated by the rainfall data Means,
Based on the rainfall intensity distribution data read by the rainfall intensity distribution data reading means, an attenuation analysis means for analyzing the probability distribution of the attenuation amount of the radio wave due to the rainfall of the rainfall per predetermined time;
Based on the probability distribution of the attenuation amount analyzed by the attenuation analysis means, and tolerance value data indicating an allowable value of an excess probability that is a probability that the attenuation amount takes a value larger than a certain value in the probability distribution, An attenuation analysis device for each category, comprising: radio wave attenuation analysis means for analyzing an attenuation at which the excess probability becomes the allowable value for each category.
ある観測地点について観測された10分間ごとの降雨強度のうち、間における雨量が所定値となる前記降雨強度を示す地点別雨量データに基づいて、前記観測地点における前記時間あたりの雨量の降雨による電波の減衰量の確率分布を解析する降雨減衰解析方法であって、
前記地点別雨量データによって示される前記降雨強度の確率分布を対数正規分布で近似して、前記確率分布を示す降雨強度分布データを生成する降雨強度分布解析ステップと、
この降雨強度分布解析手段によって生成された降雨強度分布データである10分間降雨強度の確率分布の平均値及び標準偏差を、1分間降雨強度の確率分布の平均値及び標準偏差とみなして、前記時間あたりの雨量の降雨における前記減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する減衰解析ステップと
を含み、
前記降雨強度分布解析ステップは、
前記観測地点ごとに継続して観測された6つの10分間雨量を加算して1時間雨量を算出し、当該1時間雨量が前記所定値となる前記6つの10分間雨量の組を複数取り出しこの取り出した10分間雨量の各値を6倍して10分間降雨強度を算出し、前記取り出した複数の組の各々に対して10分間降雨強度の値を昇順に並び替え、値が小さいものから順にその値を超える確率を、第1超過確率、第2超過確率、第3超過確率、第4超過確率、及び、第5超過確率とし、この複数の組にわたる各超過確率における10分間降雨強度の中央値を、当該超過確率における10分間降雨強度とし、この5つの超過確率から導かれる10分間降雨強度分布から当該分布に従う10分間降雨強度の平均値と標準偏差とを算出して、前記降雨強度分布データとする
ことを特徴とする降雨減衰解析方法。
Of rainfall every 10 minutes observed for an observation point, based on the locational rainfall data indicating the rain intensity rainfall definitive between 1 hour becomes a predetermined value, rainfall of the per hour in the observation point A rain attenuation analysis method for analyzing the probability distribution of radio attenuation due to rain
A rainfall intensity distribution analysis step for generating a rainfall intensity distribution data indicating the probability distribution by approximating the probability distribution of the rainfall intensity indicated by the point-by-point rainfall data with a lognormal distribution;
The mean and standard deviation of the probability distribution of the 10 minute rainfall intensity is rainfall intensity distribution data generated by the rainfall intensity distribution analysis unit, is regarded as the mean and standard deviation of the probability distribution of 1 minute rain rate, the 1 by analyzing the probability distribution of the attenuation in rainfall rainfall per hour saw including a damping analysis step of generating the attenuation distribution data indicating the probability distribution,
The rainfall intensity distribution analysis step includes
Wherein by adding six 10 minute rainfall observed continuously for each observation point is calculated for 1 hour rainfall, taken out a plurality of sets of the six 10 minute rainfall which the 1 hour rainfall becomes the predetermined value, the The 10-minute rainfall intensity is calculated by multiplying the taken 10-minute rainfall values by 6 times, and the 10-minute rainfall intensity values are rearranged in ascending order with respect to each of the plurality of the extracted sets. The probability exceeding the value is defined as the first excess probability, the second excess probability, the third excess probability, the fourth excess probability, and the fifth excess probability, and the center of the 10-minute rainfall intensity at each excess probability across the plurality of sets. The value is defined as the 10-minute rainfall intensity at the excess probability, and the average value and the standard deviation of the 10-minute rainfall intensity according to the distribution are calculated from the 10-minute rainfall intensity distribution derived from the five excess probabilities, and the rainfall intensity distribution is calculated. De Rain attenuation analysis method characterized by the data.
ある観測地点について観測された10分間ごとの降雨強度のうち、所定時間内における雨量が所定値となる前記降雨強度を示す地点別雨量データによって示される前記降雨強度の確率分布を対数正規分布で近似して、前記確率分布を示す降雨強度分布データを生成する降雨強度分布解析手段と、
この降雨強度分布解析手段によって生成された降雨強度分布データに基づいて、前記所定時間あたりの雨量の降雨における前記減衰量の確率分布を解析して当該確率分布を示す減衰量分布データを生成する減衰解析手段と
を備える降雨減衰解析装置によって生成された、ある観測地点において所定時間あたりの雨量が所定値となる降雨による電波の減衰量の確率分布を示す減衰量分布データを、前記観測地点及び前記雨量と対応付けて、複数の前記観測地点及び複数の前記雨量について記憶する減衰量分布データ記憶装置から前記減衰量分布データを読み出して、前記減衰量の解析対象となる地域内の複数の区分の各々について降雨による前記電波の減衰量を解析するためにコンピュータを、
前記区分ごとの所定時間あたりの雨量を示す雨量データを入力する雨量データ入力手段、
前記減衰量分布データ記憶装置から、前記区分ごとに、当該区分に予め対応付けられた前記観測地点と、前記雨量データによって示される前記雨量とに対応する減衰量分布データを読み出す減衰量分布データ読み出し手段、
この減衰量分布データ読み出し手段によって読み出された前記減衰量分布データによって示される前記減衰量の確率分布と、前記確率分布において前記減衰量がある値より大きな値をとる確率である超過確率の許容値を示す許容値データとに基づいて、前記区分ごとに前記超過確率が前記許容値となる減衰量を解析する電波減衰解析手段
として機能させることを特徴とする区分別減衰量解析プログラム。
Approximate the probability distribution of the rainfall intensity indicated by the rainfall data according to the point indicating the rainfall intensity within a predetermined time out of the rainfall intensity every 10 minutes observed for a certain observation point with a lognormal distribution. A rainfall intensity distribution analyzing means for generating rainfall intensity distribution data indicating the probability distribution;
Based on the rainfall intensity distribution data generated by the rainfall intensity distribution analyzing means, the attenuation for generating the attenuation distribution data indicating the probability distribution by analyzing the probability distribution of the attenuation in the rainfall of the predetermined amount of rainfall Analysis means and
Produced by rain attenuation analyzer comprising a rainfall per predetermined time attenuation distribution data indicating a probability distribution of a radio wave attenuation due to rainfall which becomes a predetermined value in a certain observation point, corresponding to the observation point and the rainfall In addition, the attenuation distribution data is read out from the attenuation distribution data storage device that stores the plurality of observation points and the plurality of rain amounts, and the rainfall is calculated for each of the plurality of sections in the area to be analyzed for the attenuation amount. To analyze the attenuation of the radio wave by
Rainfall data input means for inputting rainfall data indicating the rainfall per predetermined time for each section,
Attenuation distribution data reading for reading out the attenuation distribution data corresponding to the observation point previously associated with the section and the rainfall indicated by the rainfall data from the attenuation amount distribution data storage device. means,
Allowance of the probability distribution of the attenuation indicated by the attenuation distribution data read by the attenuation distribution data reading means and an excess probability that is a probability that the attenuation takes a value larger than a certain value in the probability distribution. A category-specific attenuation amount analysis program that functions as a radio wave attenuation analysis unit that analyzes an attenuation amount for which the excess probability becomes the tolerance value for each category based on tolerance value data indicating a value.
JP2005062708A 2005-03-07 2005-03-07 Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program Expired - Fee Related JP4511979B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2005062708A JP4511979B2 (en) 2005-03-07 2005-03-07 Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2005062708A JP4511979B2 (en) 2005-03-07 2005-03-07 Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program

Publications (2)

Publication Number Publication Date
JP2006246375A JP2006246375A (en) 2006-09-14
JP4511979B2 true JP4511979B2 (en) 2010-07-28

Family

ID=37052267

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2005062708A Expired - Fee Related JP4511979B2 (en) 2005-03-07 2005-03-07 Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program

Country Status (1)

Country Link
JP (1) JP4511979B2 (en)

Families Citing this family (9)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP5025359B2 (en) * 2007-05-22 2012-09-12 三菱電機株式会社 Radar equipment
JP6040428B2 (en) * 2013-02-20 2016-12-07 株式会社Kddi総合研究所 Network monitoring device, wireless device, weather prediction system and program
US20140372038A1 (en) * 2013-04-04 2014-12-18 Sky Motion Research, Ulc Method for generating and displaying a nowcast in selectable time increments
JP6166636B2 (en) * 2013-10-08 2017-07-19 日本放送協会 Apparatus, satellite system, and method for determining rainfall attenuation compensation
JP6504562B2 (en) * 2015-07-03 2019-04-24 日本無線株式会社 Communication diagnostic system and communication diagnostic method
JP6582976B2 (en) * 2015-12-28 2019-10-02 富士通株式会社 Estimation apparatus, estimation method, and estimation program
JP7097056B2 (en) * 2018-06-14 2022-07-07 国立研究開発法人宇宙航空研究開発機構 Communication equipment, communication methods and programs
CN112712195B (en) * 2020-12-17 2024-01-05 杭州鲁尔物联科技有限公司 Rainfall attenuation coefficient prediction method, rainfall attenuation coefficient prediction device, rainfall attenuation coefficient prediction equipment and storage medium
JP7170763B2 (en) * 2021-02-18 2022-11-14 ソフトバンク株式会社 Apparatus, method and program for predicting attenuation of electromagnetic waves

Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000091975A (en) * 1998-09-14 2000-03-31 Mitsubishi Electric Corp Transmission power control system
JP2003318795A (en) * 2002-04-18 2003-11-07 Matsushita Electric Ind Co Ltd Radio apparatus

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JPS62199130A (en) * 1986-02-26 1987-09-02 Nec Corp Satellite communication system
JP2879836B2 (en) * 1991-08-07 1999-04-05 ケイディディ株式会社 Transmission power control method for satellite communication and broadcasting
JPH11355220A (en) * 1998-06-12 1999-12-24 Mitsubishi Electric Corp Satellite broadcast receiver

Patent Citations (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2000091975A (en) * 1998-09-14 2000-03-31 Mitsubishi Electric Corp Transmission power control system
JP2003318795A (en) * 2002-04-18 2003-11-07 Matsushita Electric Ind Co Ltd Radio apparatus

Also Published As

Publication number Publication date
JP2006246375A (en) 2006-09-14

Similar Documents

Publication Publication Date Title
JP4511979B2 (en) Rain attenuation analysis device, category attenuation analysis device, rainfall attenuation analysis method and category attenuation analysis program
Shrestha et al. Characterization of rain specific attenuation and frequency scaling method for satellite communication in South Korea
CN109936855B (en) Method, device, equipment and medium for inverting precipitation distribution by using base station signals
Crane et al. ACTS propagation experiment: experiment design, calibration, and data preparation and archival
Abdulrahman et al. A new rain attenuation conversion technique for tropical regions
Ewona et al. Tropospheric influence on low-band very high frequency (VHF) radio waves
Adegbindin et al. Intelligent weather awareness technique for mitigating propagation impairment at SHF and EHF satellite network system in a tropical climate
CN113466969B (en) Rainfall monitoring method, receiving equipment, rainfall monitoring system and storage medium
WO2007064953A2 (en) Hf radio network frequency management
JP2005321235A (en) Radio wave attenuation calculation system, maximum attenuation analysis system, rate of influence of rainfall analysis system and radio wave attenuation calculation program
Ananya et al. Atmospheric propagation impairment effects for wireless communications
Kamoru et al. Development of rain attenuation prediction in south west Nigeria on terrestrial link using artificial neural Network
JP3479484B2 (en) Communication quality calculation method and communication quality calculation device for mobile communication system
Luini et al. A time diversity model for EHF satellite communication systems
CN114007239A (en) Service satellite operation state diagnosis method and related device
Pimienta del Valle Contribution to the millimeter-wave propagation characterization for satellite and 5G wireless links
Khan et al. Extra-high frequency line-of-sight propagation for future urban communications
Abdulrahman et al. Terrestrial Rain Attenuation Transformation Model for KU-band in Tropical Regions Using Path Length Reduction Techniques
Khairolanuar et al. Assessment of ITU-R predictions for Ku-Band rain attenuation in Malaysia
Kaya et al. Rain rate and rain attenuation prediction for satellite communication in Ku band beacon over TURKSAT Golbası
Ojo et al. Influence of wind speed on rain-based attenuation at Ku-band for earth-satellite links in Nigeria
Islam Atmospheric Propagation Impairment Effects for Wireless Communications
Sanyaolu et al. Rainfade analysis at C, Ku and Ka Bands in Nigeria
Ukommi et al. Evaluation of rainfall rates and rain-induced signal attenuation for satellite communication in the South-South Region of Nigeria
Odedina et al. Clear-Air Signal Level Measurement for Microwave Line-of-Sight Link Application in South Africa

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20070409

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20090702

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20090804

A521 Written amendment

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20091001

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

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

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

Free format text: PAYMENT UNTIL: 20130514

Year of fee payment: 3

R150 Certificate of patent or registration of utility model

Free format text: JAPANESE INTERMEDIATE CODE: R150

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

Free format text: PAYMENT UNTIL: 20140514

Year of fee payment: 4

R250 Receipt of annual fees

Free format text: JAPANESE INTERMEDIATE CODE: R250

LAPS Cancellation because of no payment of annual fees