JP5511196B2 - Back turbulence detector - Google Patents
Back turbulence detector Download PDFInfo
- Publication number
- JP5511196B2 JP5511196B2 JP2009036369A JP2009036369A JP5511196B2 JP 5511196 B2 JP5511196 B2 JP 5511196B2 JP 2009036369 A JP2009036369 A JP 2009036369A JP 2009036369 A JP2009036369 A JP 2009036369A JP 5511196 B2 JP5511196 B2 JP 5511196B2
- Authority
- JP
- Japan
- Prior art keywords
- template
- unit
- wind speed
- turbulence
- predetermined
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Expired - Fee Related
Links
Images
Classifications
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
Landscapes
- Radar Systems Or Details Thereof (AREA)
Description
この発明は、大気中に対する電磁波の送受信によって、飛行する航空機の後方に発生する後方乱気流を検出し、検出した後方乱気流に関する情報を表示する後方乱気流検出装置に関するものである。 The present invention relates to a wake turbulence detector that detects wake turbulence generated behind a flying aircraft by transmitting and receiving electromagnetic waves to and from the atmosphere and displays information about the detected wake turbulence.
遠隔点に存在する物体の位置を計測するものとして、レーダ装置が知られている。レーダ装置は、電磁波や音波などの波動を空間に放射し、対象となる物体で反射された波動を受信し、その信号を解析することにより、レーダ装置から物体までの距離や角度を計測する。 A radar device is known as a device for measuring the position of an object existing at a remote point. The radar apparatus radiates waves such as electromagnetic waves and sound waves to the space, receives the waves reflected by the target object, and analyzes the signal to measure the distance and angle from the radar apparatus to the object.
レーダ装置の中でも、大気中に浮遊する微小な液体または固体の粒子(エアロゾル)を対象とし、反射された波動の位相回転量からエアロゾルの動く速度、すなわち風の速度を知ることができる気象レーダ装置が知られている。 Among radar devices, a meteorological radar device that targets minute liquid or solid particles (aerosol) floating in the atmosphere and can know the velocity of the aerosol, that is, the velocity of the wind, from the amount of phase rotation of the reflected wave It has been known.
気象レーダ装置の中でも特に電磁波として光を用いるレーザーレーダ装置は、放射するビームの広がりが極めて小さく、高い角度分解能で物体を観測することが可能であり、風向風速レーダ装置として使用されている(例えば、非特許文献1参照)。 Laser radar devices that use light as electromagnetic waves, among other weather radar devices, have a very small beam spread and can observe objects with high angular resolution, and are used as wind direction and wind velocity radar devices (for example, Non-Patent Document 1).
国内で過去に生じた大型機による人的被害を伴う航空事故の約60%が乱気流を事故の要因としている。また、世界的にもニューヨークで乱気流が事故の一因とみられている墜落事故が起きるなど、乱気流の危険性は重大視されている。 About 60% of aviation accidents involving human damage caused by large aircraft in the past in Japan are caused by turbulence. In addition, the danger of turbulence is regarded as critical, such as the occurrence of a crash in New York, where turbulence is considered to be a cause of the accident.
これらの乱気流を後方乱気流といい、この後方乱気流は、図8に示すように、航空機の飛行に伴い航空機の両翼端を中心として航空機の後方に発生する二組の渦状の乱気流である。上記のような事故抑止のために、航空管制当局では離発着にはある程度の時間間隔を空けるといった対処を行っているが、その時間間隔の定量的な値は明らかになっていない。このため、航空機の離着陸において、後方乱気流を適切に捉えることは、飛行計画の効率化や、後続機の安全を確保する上で重要である。これを解決するために、近年、上記のようなレーザーレーダ装置による風向風速の観測技術を用いて、航空機の後方乱気流を観測する例がある(例えば、非特許文献2参照)。 These turbulences are referred to as wake turbulences. These wake turbulences are two sets of vortex turbulences generated at the rear of the aircraft around the wing tips of the aircraft as the aircraft flies, as shown in FIG. In order to prevent accidents as described above, the air traffic control authorities take measures such as leaving a certain time interval for takeoff and landing, but the quantitative value of the time interval is not clear. For this reason, it is important to capture the wake turbulence appropriately in taking off and landing of an aircraft in order to improve the efficiency of the flight plan and to ensure the safety of the following aircraft. In order to solve this problem, in recent years, there is an example in which the wake turbulence of an aircraft is observed using the above-described technology for observing wind direction and wind speed by a laser radar device (see, for example, Non-Patent Document 2).
例として、図8の観測環境におけるレーザーレーダ装置10の観測状況を図9に示す。従来のレーザーレーダ装置10は、図9に示すように、レーザーレーダ装置10が航空機の側面に向けられることにより、大気中における航空機の飛行方向(進行方向)に直交する断面(以下、これを空間断面と呼ぶ)を観測し、当該観測結果に基づいて航空機の通過によって生じた後方乱気流の位置及び強度を検出するようになされている。また、一定の距離分解能、及びクロスレンジ分解能(角度分解能)を持つレーザーレーダ装置10のビームをビーム走査範囲内において走査させ、処理単位(所定の距離分解能及び所定の角度分解能からなる単位)毎に、視線方向の風速を得ることができる。
As an example, FIG. 9 shows an observation state of the
レーザーレーダ装置10によって後方乱気流を観測した場合の模擬図を図10に示す。この図10は、観測領域の空間断面における後方乱気流を示している。図10に示すように、レーザーレーダ装置10に対して向かい風を正とした場合、一様に存在する背景風に対して、渦上においては背景風に対して小さい風と大きい風が観測される。
FIG. 10 shows a simulation diagram when wake turbulence is observed by the
この観測模擬図を図11に示す。図11(1)に示すように、地上付近に航空機が飛行し、レーザーレーダ装置10の仰角が低い状況下において、レーザーレーダ装置10による観測結果よりドップラ風速の分布図を作成すると、航空機の通過時に発生する後方乱気流は、背景風を基準とすると正と負のマトリクス構造を持つ後方乱気流が観測される。
This observation simulation diagram is shown in FIG. As shown in FIG. 11 (1), when an aircraft flies near the ground and the elevation angle of the
また、等レンジビン上で風速と高度のプロットを作成すると、図11(2)のように、極大値と極小値を持つ距離分解能毎(等レンジビン)の風速プロットが作成できる。このように、従来では後方乱気流は正と負のマトリクス構造を持つ結果を利用して、レーザーレーダ装置10による乱気流観測において得られるドップラ風速の分布図から、正と負のマトリクスで構築されるテンプレートを用いたテンプレートマッチング処理によって、後方乱気流渦の検出が行われてきた(例えば、特許文献1参照)。
When a plot of wind speed and altitude is created on the equi-range bin, a wind speed plot for each distance resolution (equal range bin) having a maximum value and a minimum value can be created as shown in FIG. Thus, conventionally, the wake turbulence uses a result having a positive and negative matrix structure, and is a template constructed with a positive and negative matrix from a distribution map of Doppler wind speed obtained in turbulence observation by the
しかしながら、従来技術には、以下のような課題がある。すなわち、航空機の飛行高度が高い場合、その後方乱気流を観測しようとすると、レーザーレーダ装置10の仰角を上げる必要がある。その際に得られるドップラ風速の分布図は図12のようになる。仰角が低い図11の場合とは違い、仰角が高い図12の場合は渦がレンジをまたぎやすくなる。またがない場合には、従来の正と負で構成されるマトリクスのテンプレートを用いる。渦がレンジをまたいでレンジ内に両渦が観測された場合は、レーザーレーダ装置10の距離分解能の制約により、渦と渦の間付近を通るレーザービームでは両渦の風速を観測するため、顕著な正の渦風速を得ることができるが、レンジからはみ出した渦領域はレンジに対してその領域は小さいため、領域内の背景風の風速の方が支配的となり、得られるデータは渦の風速ではなく、背景風の値が得られる。このため、従来の正と負で構成されるマトリクスのテンプレートを用いたテンプレートマッチングでは渦の検出は不可能となる。
However, the prior art has the following problems. That is, when the flight altitude of the aircraft is high, it is necessary to increase the elevation angle of the
本発明は、前記のような課題を解決するためになされたものであり、航空機の飛行高度が高いため、仰角が高い状況で観測している場合においても後方乱気流を検出することができる後方乱気流検出装置を得ることを目的とする。 The present invention has been made in order to solve the above-described problems, and since the flight altitude of the aircraft is high, the wake turbulence can detect the wake turbulence even when observing in a situation where the elevation angle is high. The object is to obtain a detection device.
本発明に係る後方乱気流検出装置は、航空機の飛行方向に直交する方向に大気中へ電磁波を所定仰角の範囲内で発射した後、その大気からの反射波を一空間断面の信号として受信する電磁波送受信部と、前記電磁波送受信部からの受信信号に基づき、空間断面のドップラ風速分布を算出し、算出したドップラ風速分布から空間断面における後方乱気流の位置を算出する信号処理部と、前記信号処理部により算出された空間断面における後方乱気流の位置を表示する表示部とを備え、前記信号処理部は、前記電磁波送受信部からの受信信号に基づき、所定の距離分解能及び所定の角度分解能からなる処理単位毎に風速値を算出し、それらの結果より空間断面におけるドップラ風速分布を算出するドップラ風速分布算出部と、前記電磁波送受信部から発射されるビームの仰角に応じたテンプレートを使用して、前記ドップラ風速分布算出部により算出されたドップラ風速分布にテンプレートマッチング処理を行うことにより、後方乱気流を検出する検出処理部と、前記テンプレートマッチング処理の相関値が所定の値である空間断面の所定の位置から、空間断面における後方乱気流の位置を算出する乱気流位置算出部とを有し、前記テンプレートは、所定の角度以下である仰角が低い場合は、正と負のマトリクス構造を持つマトリクス型テンプレートを使用し、前記所定の角度より仰角が高い場合には、正とゼロのマトリクス構造を持ちかつ前記マトリクス型テンプレートを仰角に応じた所定量傾けた傾斜マトリクス型テンプレート、正と負のマトリクス構造を持ちかつ前記マトリクス型テンプレートを仰角に応じた所定量傾けた傾斜マトリクス型テンプレートの少なくとも1つを使用するものである。 The wake turbulence detection device according to the present invention emits an electromagnetic wave into the atmosphere in a direction perpendicular to the flight direction of the aircraft within a predetermined elevation angle, and then receives the reflected wave from the atmosphere as a signal of one spatial section. A signal processing unit that calculates a Doppler wind speed distribution in a spatial section based on a reception signal from the transmission / reception unit, the electromagnetic wave transmission / reception unit, and calculates a position of wake turbulence in the spatial section from the calculated Doppler wind speed distribution, and the signal processing unit A display unit that displays the position of wake turbulence in the spatial cross section calculated by the above, and the signal processing unit is a processing unit having a predetermined distance resolution and a predetermined angle resolution based on a received signal from the electromagnetic wave transmitting / receiving unit A wind speed value is calculated every time, and a Doppler wind speed distribution calculating section for calculating a Doppler wind speed distribution in a spatial section from the results, and the electromagnetic wave transmitting / receiving section Using the template in accordance with the elevation angle of the beam et fired, said by performing template matching in Doppler velocity distribution calculated by the Doppler velocity distribution calculating section, and a detection processing unit for detecting a wake turbulence, the template A turbulence position calculation unit that calculates a position of wake turbulence in the spatial section from a predetermined position of the spatial section where the correlation value of the matching process is a predetermined value, and the template has an elevation angle that is equal to or less than a predetermined angle If it is low, a matrix template having a positive and negative matrix structure is used. If the elevation angle is higher than the predetermined angle, the matrix template has a positive and zero matrix structure and the matrix template is adapted to the elevation angle. Tilt matrix template with a fixed amount of tilt, positive and negative matrix structure and the matrix The template is to use at least one of a predetermined amount tilted inclined matrix template in accordance with the elevation angle.
本発明に係る後方乱気流検出装置によれば、航空機の飛行高度が高いため、仰角が高い状況で観測している場合においても後方乱気流を検出することができる。 According to the backward turbulence detection device according to the present invention, since the flight altitude of the aircraft is high, the backward turbulence can be detected even when the elevation angle is high.
以下、本発明の後方乱気流検出装置の好適な実施の形態につき図面を用いて説明する。 Hereinafter, a preferred embodiment of a backward turbulence detection device of the present invention will be described with reference to the drawings.
実施の形態1.
この発明の実施の形態1に係る後方乱気流検出装置について図1及び図2を参照しながら説明する。この発明の実施の形態1に係る後方乱気流検出装置は、ドップラ風速分布モデルを使用して後方乱気流を検出するものである。図1は、この発明の実施の形態1に係る後方乱気流検出装置の構成を示すブロック図である。なお、各図中、同一符号は同一又は相当部分を示す。
A backward turbulence detector according to
図1において、この発明の実施の形態1に係る後方乱気流検出装置は、大気中に対する電磁波の送受信によって、飛行する航空機の後方に発生する後方乱気流を検出し、検出した後方乱気流に関する情報を表示するもので、航空機の飛行方向に直交する方向に大気中へ電磁波を所定仰角の範囲内で発射した後、その大気からの反射波を一空間断面の信号として受信する電磁波送受信部1と、この電磁波送受信部1からの受信信号に基づき、空間断面のドップラ風速分布を算出し、算出したドップラ風速分布から空間断面における後方乱気流の位置を算出する信号処理部2と、この信号処理部2により算出された空間断面における後方乱気流の位置を表示する表示部3とが設けられている。
In FIG. 1, the wake turbulence detection apparatus according to
また、信号処理部2は、電磁波送受信部1からの受信信号に基づき、所定の距離分解能及び所定の角度分解能からなる処理単位毎に風速値を算出し、それらの結果より空間断面におけるドップラ風速分布を算出するドップラ風速分布算出部21と、このドップラ風速分布算出部21により算出されたドップラ風速分布と、後方乱気流を模擬するためのパラメータにより算出されたドップラ風速分布モデルとの相互相関処理を行うことにより、後方乱気流を検出する検出処理部22と、相互相関処理の相関値が所定の値である空間断面の所定の位置から、空間断面における後方乱気流の位置を算出する乱気流位置算出部23とが設けられている。
Further, the
図2は、この発明の実施の形態1に係る後方乱気流検出装置の検出処理部の構成を示すブロック図である。
FIG. 2 is a block diagram showing a configuration of a detection processing unit of the backward turbulence detection device according to
図2において、検出処理部22は、後方乱気流を模擬するための航空機に関するパラメータを使用して後方乱気流モデルを含んだドップラ風速分布モデルを算出するドップラ風速分布モデル算出部101と、このドップラ風速分布モデル算出部101により算出したドップラ風速分布モデルとドップラ風速分布算出部21により算出されたドップラ風速分布との相互相関処理により相関値を求める相互相関処理部102と、この相互相関処理部102により求められた相関値のうち最も高い相関値と所定の閾値との比較を行って空間断面における後方乱気流を検出する相関値比較部103と、最も高い相関値が閾値を超えていない場合は、パラメータを変更してドップラ風速分布モデル算出部101へ出力するパラメータ出力部104とが設けられている。
In FIG. 2, the
つぎに、この実施の形態1に係る後方乱気流検出装置の動作について図面を参照しながら説明する。
Next, the operation of the wake turbulence detector according to
まず、電磁波送受信部1は、図8及び図9に示したように、航空機の飛行方向(進行方向)に直交する方向に大気中へ電磁波である光パルスを所定仰角の範囲内で発射した後、その大気からの反射光を受信し、当該受信した反射光を一空間断面の受信信号として信号処理部2へ送出する。この電磁波送受信部1が受信した反射光には、図10に示すように、反射光の反射位置における風速に応じてドップラ効果が生じている。
First, as shown in FIGS. 8 and 9, the electromagnetic wave transmission /
次に、信号処理部2は、電磁波送受信部1からの受信信号に基づき、空間断面のドップラ風速分布を算出し、算出したドップラ風速分布から空間断面における後方乱気流の位置を算出する。この信号処理部2の詳細な処理については以下に説明する。
Next, the
信号処理部2内のドップラ風速分布算出部21は、電磁波送受信部1からの受信信号に基づき、所定の距離分解能及び所定の角度分解能からなる処理単位毎に、風速値を算出する。それらの結果より空間断面におけるドップラ風速分布を算出する。
The Doppler wind speed
その後、信号処理部2内の検出処理部22は、観測により得られたドップラ風速分布と、パラメータにより算出されたドップラ風速分布モデルとの相互相関処理を行うことにより、後方乱気流を検出する。この検出処理部22の詳細な処理については以下に説明する。
Thereafter, the
検出処理部22内のドップラ風速分布モデル算出部101は、パラメータ出力部104から、観測を行う航空機の種類(翼の幅)、当該後方乱気流検出装置前を通過する航空機の高度とその速度、大気密度等のパラメータを取得し、このパラメータを使用して処理単位毎にモデル値を求めて、渦モデル(後方乱気流モデル)を含んだドップラ風速分布モデルを算出する。航空機の翼の幅は後方乱気流渦の大きさを与え、航空機の高度は渦が発生すると想定される場所(高度)の参考値となる。航空機の速度や大気密度は渦の強度を算出させることができる。このドップラ風速分布モデルを作成できるものならば上記のパラメータ以外を使用してもよい。
The Doppler wind speed distribution
次に、検出処理部22内の相互相関処理部102は、算出したドップラ風速分布モデルと観測により得られたドップラ風速分布との相互相関処理を行う。その際、処理単位毎に、相互相関処理を行い、分布の処理単位毎に相関係数(相関値)を求め、その結果を相関値比較部103へ送出する。
Next, the
次に、検出処理部22内の相関値比較部103は、得られた最も高い相関係数(相関値)と閾値との比較を行い、最も高い相関係数が閾値を超えていない場合は、後方乱気流が検出できない場合として、パラメータを変更すべくパラメータ出力部104へ移行する。パラメータを種々変えても、最も高い相関係数が閾値を超えていない場合は、その分布上に後方乱気流がないと判断する。最も高い相関係数が閾値を超えている場合(両者が等しい場合を含む)は、後方乱気流を検出できた場合として、乱気流位置算出部23へ移行する。
Next, the correlation value comparison unit 103 in the
検出処理部22内のパラメータ出力部104は、当初は予め認識している、観測を行う航空機の種類(翼の幅)、当該後方乱気流検出装置前を通過する航空機の高度とその速度、大気密度等のパラメータをドップラ風速分布モデル算出部101に出力するが、最も高い相関係数が閾値を超えていない場合は、観測を行う航空機の種類(翼の幅)、当該後方乱気流検出装置前を通過する航空機の高度とその速度、大気密度等のパラメータのいずれか、あるいはいくつかを変更してドップラ風速分布モデル算出部101に出力する。このドップラ風速分布モデル算出部101は、変更されたパラメータに基づき新しいドップラ風速分布モデルを作成する。
The
信号処理部2内の乱気流位置算出部23は、最も高い相関係数が閾値を超えている場合は、その結果から後方乱気流の渦構造の場所を算出し、さらに渦の移動予測を行う。相関位置(処理単位の位置)から、実空間での後方乱気流の位置を算出する。乱気流位置算出部23は、相関係数が高い処理単位の位置から、空間断面における後方乱気流の位置を算出する。後方乱気流の位置は、2つの渦のそれぞれの中心や、2つの渦の間の中心などである。
When the highest correlation coefficient exceeds the threshold value, the turbulence
そして、表示部3は、信号処理部2内の乱気流位置算出部23により算出された、実空間での後方乱気流の位置を表示する。空間断面における後方乱気流の位置をディスプレイに表示する。
And the
以上のように、実施の形態1によれば、航空機の飛行高度が高いため、仰角が高い状況で観測している場合においても後方乱気流を検出することができる。 As described above, according to the first embodiment, since the flight altitude of the aircraft is high, the wake turbulence can be detected even when observation is performed in a situation where the elevation angle is high.
実施の形態2.
この発明の実施の形態2に係る後方乱気流検出装置について図3から図7までを参照しながら説明する。この発明の実施の形態2に係る後方乱気流検出装置は、テンプレートを使用して後方乱気流を検出するものである。図3は、この発明の実施の形態2に係る後方乱気流検出装置の検出処理部の構成を示すブロック図である。なお、この発明の実施の形態2に係る後方乱気流検出装置の全体構成は、上記の実施の形態1と同様である。
A backward turbulence detection apparatus according to
図3において、検出処理部22は、空間断面におけるドップラ風速分布に後方乱気流の検出範囲である所定のフレームを設定するフレーム設定部201と、電磁波送受信部1から発射されるビームの仰角に応じた複数種類のテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を行うテンプレートマッチング部202と、テンプレートマッチング処理により得られた特徴量のうち、所定の特徴量を後方乱気流候補として選出する特徴量算出部203と、背景風の風速値を算出し、前記後方乱気流候補の風速値と前記背景風の風速値の差分から背景風比較値を算出する背景風比較値算出部204と、この背景風比較値算出部204により算出された背景風比較値と所定の閾値との比較を行って空間断面における後方乱気流を検出する背景風比較値比較部205と、背景風比較値が閾値を超えていない場合は、前記テンプレートを変更してテンプレートマッチング部202へ出力するテンプレート出力部206とが設けられている。
In FIG. 3, the
つぎに、この実施の形態2に係る後方乱気流検出装置の動作について図面を参照しながら説明する。
Next, the operation of the wake turbulence detector according to
まず、検出処理部22内のフレーム設定部201は、空間断面におけるドップラ風速分布に後方乱気流の検出範囲である、航空機の翼の幅で決まるフレーム(ウィンドウとも言う)幅を設定する。後述するテンプレートマッチング部202において、その処理を行わせる範囲を決定するためである。つまり、航空機が通過した場所または後方乱気流があると想定される領域に任意の大きさのフレームを設定する。このフレームは、テンプレートの大きさを示す。なお、このフレーム設定部201をなくし、空間断面におけるドップラ風速分布全体にテンプレートマッチング処理を行わせてもよい。
First, the
次に、検出処理部22内のテンプレートマッチング部202は、テンプレート出力部206から取得した仰角に応じたテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を行う。その際、テンプレートマッチング処理では、航空機の飛行高度によってテンプレートを図4のように変更する。図4に示すように、後方乱気流検出装置10Aの電磁波送受信部1から発射される光パルス(ビーム)の仰角や、後方乱気流の渦構造によってテンプレートを変化させる。例えば、仰角が高い場合には、テンプレートAやテンプレートBを使用し、仰角が低い場合には、テンプレートCやテンプレートDを使用する。
Next, the
次に、検出処理部22内の特徴量算出部203は、設定された航空機の翼の幅で決まるフレームをドップラ風速分布上で任意に動かし、テンプレートマッチング(相互相関)処理することにより、空間断面内における後方乱気流が存在する位置が強調され、すなわち、処理単位における後方乱気流が存在する位置に対応する積分値が積み上がり、実際に後方乱気流が存在する場合には、そのフレーム部分には渦があるものとして後方乱気流を検出することができる。この特徴量算出部203は、処理単位毎に得られた積分値(特徴量)の中から最も大きい特徴量、または所定の相関値(特徴量を算出するための閾値)を超えた特徴量を後方乱気流候補(渦候補)として選出する。フレームサイズは任意で、上記以外の手法を用いてもよい。
Next, the feature amount calculation unit 203 in the
その後、検出処理部22内の背景風比較値算出部204は、まず、処理単位毎に背景風の風速値を算出する。背景風の風速値の算出手法については、各レンジビンの風速の平均や、その他過去の観測によって得られた渦構造のないドップラ風速分布、または風速の高度分布モデルといった背景風の風速値を得られるものならばどのような手法でもよい。上述したように、テンプレートマッチング部202によってテンプレートマッチング処理を行って、特徴量算出部203により選出された特徴量を後方乱気流候補(渦候補)とする。
Thereafter, the background wind comparison
背景風比較値算出部204は、次に、その後方乱気流候補(渦候補)と背景風の風速値から、次の式(1)を用いて処理単位毎に背景風比較値(m/s)を算出する。
Next, the background wind comparison
ここで、Vcompは背景風比較値、Vobsは観測された後方乱気流候補(渦候補)の風速値、Vbgは背景風の風速値をそれぞれ表す。 Here, Vcomp represents the background wind comparison value, Vobs represents the wind speed value of the observed wake turbulence candidate (vortex candidate), and Vbg represents the wind speed value of the background wind.
この背景風比較値の結果より誤検出を除去し、高精度に渦構造を抽出する。なお、背景風比較値の評価の詳細については後述する。 The false detection is removed from the result of the background wind comparison value, and the vortex structure is extracted with high accuracy. Details of the background wind comparison value evaluation will be described later.
ここで、テンプレート出力部206から出力されるテンプレートについて説明する。テンプレート例として、フレーム高度が低い、すなわち当該後方乱気流検出装置10Aの仰角が低い場合は、図5(1)に示すような、正と負のマトリクス構造を持つマトリクス型テンプレートを用いる。図5(2)及び(3)は、上記のマトリクス型テンプレートにおいて、後方乱気流の2つの渦が隣り合っている場合で、渦風速変化が高い領域に大きい重みがかかり、渦構造の端、つまり正と負の境界において渦風速変化が小さいため小さい重みがかかるような重み関数をかけている。さらに、図5(4)は、渦構造が離れている場合を想定したもので、正と負のマトリクス構造を持つマトリクス型テンプレートであるが、2つに分かれているマトリクス型テンプレートである。
Here, the template output from the
図5のようなマトリクス型テンプレートに対して、フレーム高度が高い、すなわち当該後方乱気流検出装置10Aの仰角が高い場合には、図6(1)−(4)に示すような、正とゼロ(0)のマトリクス構造を持ち、かつ所定の傾きを持つ傾斜マトリクス型テンプレートを用いる。この所定の傾きは、図6に示す傾きだけでなく、種々の傾きを持つ傾斜マトリクス型テンプレートを用意する。また、正と負のマトリクス構造を持ち、かつ所定の傾きを持つ傾斜マトリクス型テンプレートを用いる。このテンプレートは、図6(1)−(4)に示すゼロ(0)の位置に負を割り当てたものである。本実施の形態2は、当該後方乱気流検出装置10Aの仰角によって、図5及び図6に示すテンプレートを変えて使用することを特徴としている。
When the frame altitude is high with respect to the matrix template as shown in FIG. 5, that is, when the elevation angle of the wake
テンプレート出力部206から出力されるテンプレートで、当該後方乱気流検出装置10Aの仰角が高い場合は、図6(1)のように、正の部分のみを特徴量として抽出を行う。また、観測状況、渦の構造によっては2つの正の部分の位置関係は変化することも考慮すべきであるため、図6(2)のように、2つの正の部分を任意にずらした傾斜マトリクス型テンプレートも用意する。さらに、図6(3)及び(4)のように、抽出精度向上のため、渦構造の中心部分に大きい重みがかかるような重み関数をかけた、傾斜マトリクス型を用意する。これにより、観測されたデータの渦構造の不要部分による影響を抑止できる。
In the template output from the
テンプレート出力部206は、図6のような汎用型である、傾斜マトリクス型テンプレートの他にも、仰角に応じた過去の観測例または観測される渦のモデルをテンプレートとして出力することも可能である。図7(2)のように得られる等レンジビンプロットの過去の渦観測または渦モデルを用意し、観測結果に対して等レンジビンプロットの相互相関をとらせることによって、渦を検出することも可能である。また、渦を観測している図7(2)のような等レンジビンプロットの渦部分のみを切り出し、それを図6(1)−(4)の正の部分にテンプレートとして使用し、相互相関と同様な効果を得ることも可能である。その際、上述した重みをかけてもよい。図7(3)のような相互相関型テンプレートは、図6に示す傾斜マトリクス型テンプレートだけでなく、図5に示すマトリクス型テンプレートにも適用することが可能である。
The
背景風比較値算出部204は、上述したように、式(1)を用いて背景風比較値を算出する。図5に示すマトリクス型テンプレートの正領域/負領域や、図6に示す傾斜マトリクス型テンプレートの正領域においては、背景風比較値が大きければ大きい程有意となる。後述する背景風比較値比較部205は、この背景風比較値が閾値を超えている、または最も大きい場合は、それを後方乱気流の渦として検出する。ただし、当該後方乱気流検出装置10Aの仰角が高い場合には、すなわち図6及び図7の背景風比較領域においては、背景風と同等の値が得られていることが後方乱気流の条件であるため、背景風比較値が0に近いほど後方乱気流であると判断する。
The background wind comparison
背景風比較値比較部205は、算出された背景風比較値(m/s)と所定の閾値(m/s)を比較する。背景風比較値が閾値を超えていない場合は、後方乱気流が検出できない場合として、テンプレートを変更すべくテンプレート出力部206へ移行する。テンプレートを種々変えても、背景風比較値が閾値を超えていない場合は、その空間断面上に後方乱気流がないと判断する。背景風比較値が閾値を超えている場合(両者が等しい場合を含む)は、後方乱気流を検出できた場合として、乱気流位置算出部23へ移行する。
The background wind comparison value comparison unit 205 compares the calculated background wind comparison value (m / s) with a predetermined threshold value (m / s). If the background wind comparison value does not exceed the threshold value, the process proceeds to the
検出処理部22内のテンプレート出力部206は、当初は予め認識している、電磁波送受信部1から発射される光パルス(ビーム)の仰角に応じたテンプレートをテンプレートマッチング部202へ出力するが、背景風比較値が閾値を超えていない場合は、当初と異なる仰角に対応した異なるテンプレートや、観測した後方乱気流の渦構造によって異なるテンプレートに変更してテンプレートマッチング部202へ出力する。このテンプレートマッチング部202は、変更されたテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を再び行う。
The
なお、背景風比較値が閾値を超えていない場合おいては、フレーム設定部201によりフレームを広げたりして、フレームを変更して再処理を行ってもよい。
If the background wind comparison value does not exceed the threshold, the
乱気流位置算出部23は、背景風比較値が閾値を超えている場合は、その結果から後方乱気流の渦構造の場所を算出し、さらに渦の移動予測を行う。当該背景風比較値の処理単位の位置から、実空間での後方乱気流の位置を算出する。乱気流位置算出部23は、当該背景風比較値に対応する処理単位の位置から、空間断面における後方乱気流の位置を算出する。後方乱気流の位置は、2つの渦のそれぞれの中心や、2つの渦の間の中心などである。
When the background wind comparison value exceeds the threshold value, the turbulence
そして、表示部3は、信号処理部2内の乱気流位置算出部23により算出された、実空間での後方乱気流の位置を表示する。空間断面における後方乱気流の位置をディスプレイに表示する。
And the
以上のように、実施の形態2によれば、航空機の飛行高度が高いため、仰角が高い状況で観測している場合においても後方乱気流を検出することができる。 As described above, according to the second embodiment, since the flight altitude of the aircraft is high, the wake turbulence can be detected even when observation is performed in a situation where the elevation angle is high.
1 電磁波送受信部、2 信号処理部、3 表示部、10A 後方乱気流検出装置、21 ドップラ風速分布算出部、22 検出処理部、23 乱気流位置算出部、101 ドップラ風速分布モデル算出部、102 相互相関処理部、103 相関値比較部、104 パラメータ出力部、201 フレーム設定部、202 テンプレートマッチング部、203 特徴量算出部、204 背景風比較値算出部、205 背景風比較値比較部、206 テンプレート出力部。
DESCRIPTION OF
Claims (2)
前記電磁波送受信部からの受信信号に基づき、空間断面のドップラ風速分布を算出し、算出したドップラ風速分布から空間断面における後方乱気流の位置を算出する信号処理部と、
前記信号処理部により算出された空間断面における後方乱気流の位置を表示する表示部とを備え、
前記信号処理部は、
前記電磁波送受信部からの受信信号に基づき、所定の距離分解能及び所定の角度分解能からなる処理単位毎に風速値を算出し、それらの結果より空間断面におけるドップラ風速分布を算出するドップラ風速分布算出部と、
前記電磁波送受信部から発射されるビームの仰角に応じたテンプレートを使用して、前記ドップラ風速分布算出部により算出されたドップラ風速分布にテンプレートマッチング処理を行うことにより、後方乱気流を検出する検出処理部と、
前記テンプレートマッチング処理の相関値が所定の値である空間断面の所定の位置から、空間断面における後方乱気流の位置を算出する乱気流位置算出部とを有し、
前記テンプレートは、所定の角度以下である仰角が低い場合は、正と負のマトリクス構造を持つマトリクス型テンプレートを使用し、前記所定の角度より仰角が高い場合には、正とゼロのマトリクス構造を持ちかつ前記マトリクス型テンプレートを仰角に応じた所定量傾けた傾斜マトリクス型テンプレート、正と負のマトリクス構造を持ちかつ前記マトリクス型テンプレートを仰角に応じた所定量傾けた傾斜マトリクス型テンプレートの少なくとも1つを使用する
ことを特徴とする後方乱気流検出装置。 An electromagnetic wave transmitting / receiving unit that emits an electromagnetic wave into the atmosphere in a direction perpendicular to the flight direction of the aircraft within a range of a predetermined elevation angle, and then receives a reflected wave from the atmosphere as a signal of one spatial section;
Based on the received signal from the electromagnetic wave transmission / reception unit, calculate the Doppler wind speed distribution of the spatial cross section, a signal processing unit that calculates the position of the wake turbulence in the spatial cross section from the calculated Doppler wind speed distribution;
A display unit for displaying the position of wake turbulence in the spatial cross section calculated by the signal processing unit,
The signal processing unit
A Doppler wind speed distribution calculation unit that calculates a wind speed value for each processing unit having a predetermined distance resolution and a predetermined angle resolution based on a reception signal from the electromagnetic wave transmission / reception unit, and calculates a Doppler wind speed distribution in a spatial section from the results. When,
A detection processing unit that detects wake turbulence by performing a template matching process on the Doppler wind speed distribution calculated by the Doppler wind speed distribution calculation unit using a template corresponding to the elevation angle of the beam emitted from the electromagnetic wave transmission / reception unit When,
A turbulence position calculation unit that calculates a position of wake turbulence in the spatial section from a predetermined position of the spatial section where the correlation value of the template matching process is a predetermined value;
The template uses a matrix template having a positive and negative matrix structure when the elevation angle is lower than a predetermined angle, and a positive and zero matrix structure when the elevation angle is higher than the predetermined angle. At least one of has and the matrix template a predetermined amount tilted inclined matrix template in accordance with the elevation angle, the positive and negative has a matrix structure and inclined matrix template the matrix template is inclined a predetermined amount corresponding to the elevation angle A wake turbulence detection device characterized by using.
空間断面におけるドップラ風速分布に後方乱気流の検出範囲である所定のフレームを設定するフレーム設定部と、
前記電磁波送受信部から発射されるビームの仰角に応じた複数種類のテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を行うテンプレートマッチング部と、
前記テンプレートマッチング処理により得られた特徴量のうち、所定の特徴量を後方乱気流候補として選出する特徴量算出部と、
背景風の風速値を算出し、前記後方乱気流候補の風速値と前記背景風の風速値の差分から背景風比較値を算出する背景風比較値算出部と、
前記背景風比較値算出部により算出された背景風比較値と所定の閾値との比較を行って空間断面における後方乱気流を検出する背景風比較値比較部と、
背景風比較値が閾値を超えていない場合は、前記テンプレートを変更して前記テンプレートマッチング部へ出力するテンプレート出力部とを含む
ことを特徴とする請求項1記載の後方乱気流検出装置。 The detection processing unit
A frame setting unit that sets a predetermined frame that is a detection range of wake turbulence in the Doppler wind speed distribution in the spatial section;
A template matching unit that performs a template matching process on a set frame using a plurality of types of templates according to the elevation angle of the beam emitted from the electromagnetic wave transmission / reception unit;
Of the feature quantities obtained by the template matching process, a feature quantity calculation unit that selects a predetermined feature quantity as a posterior turbulence candidate; and
A background wind comparison value calculating unit that calculates a wind speed value of a background wind and calculates a background wind comparison value from a difference between the wind speed value of the wake turbulence candidate and the wind speed value of the background wind;
A background wind comparison value comparing unit that compares the background wind comparison value calculated by the background wind comparison value calculating unit with a predetermined threshold and detects wake turbulence in the spatial section;
The rear turbulence detection device according to claim 1, further comprising: a template output unit that changes the template and outputs the template to the template matching unit when a background wind comparison value does not exceed a threshold value.
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009036369A JP5511196B2 (en) | 2009-02-19 | 2009-02-19 | Back turbulence detector |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
JP2009036369A JP5511196B2 (en) | 2009-02-19 | 2009-02-19 | Back turbulence detector |
Publications (2)
Publication Number | Publication Date |
---|---|
JP2010190772A JP2010190772A (en) | 2010-09-02 |
JP5511196B2 true JP5511196B2 (en) | 2014-06-04 |
Family
ID=42816963
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
JP2009036369A Expired - Fee Related JP5511196B2 (en) | 2009-02-19 | 2009-02-19 | Back turbulence detector |
Country Status (1)
Country | Link |
---|---|
JP (1) | JP5511196B2 (en) |
Families Citing this family (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US9037319B2 (en) | 2013-09-24 | 2015-05-19 | Honeywell International Inc. | System and method for processing and displaying wake turbulence |
JP6257845B2 (en) | 2015-05-12 | 2018-01-10 | 三菱電機株式会社 | Laser radar apparatus and wind speed observation method |
JP7012907B2 (en) * | 2019-07-16 | 2022-01-28 | 三菱電機株式会社 | Eddy detection device and eddy detection method |
JP7470413B2 (en) | 2020-11-05 | 2024-04-18 | 国立研究開発法人宇宙航空研究開発機構 | Information processing device, information processing method, and program |
TWI784625B (en) * | 2021-07-20 | 2022-11-21 | 國立中央大學 | Method for detecting atmospheric turbulence and detection device of atmospheric turbulence |
Family Cites Families (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP3322214B2 (en) * | 1998-06-17 | 2002-09-09 | 三菱電機株式会社 | Airflow detection method and laser radar device |
JP3633363B2 (en) * | 1999-04-28 | 2005-03-30 | 三菱電機株式会社 | Back turbulence detection device and back turbulence detection method |
JP3570353B2 (en) * | 2000-07-25 | 2004-09-29 | 三菱電機株式会社 | Wake turbulence observation device |
JP3638862B2 (en) * | 2000-08-01 | 2005-04-13 | 三菱電機株式会社 | Turbulence detection device and turbulence detection method |
JP3636076B2 (en) * | 2001-01-19 | 2005-04-06 | 三菱電機株式会社 | Turbulence detector |
JP3628621B2 (en) * | 2001-03-06 | 2005-03-16 | 三菱電機株式会社 | Turbulence detection device and turbulence detection method |
JP3633499B2 (en) * | 2001-03-23 | 2005-03-30 | 三菱電機株式会社 | Back turbulence detection device and back turbulence detection method |
JP3958036B2 (en) * | 2001-12-10 | 2007-08-15 | 三菱電機株式会社 | Turbulence detection device and turbulence detection method |
JP3933926B2 (en) * | 2001-12-10 | 2007-06-20 | 三菱電機株式会社 | Turbulence detection device and turbulence detection method |
JP2005114416A (en) * | 2003-10-03 | 2005-04-28 | Mitsubishi Electric Corp | Disturbance detection device |
-
2009
- 2009-02-19 JP JP2009036369A patent/JP5511196B2/en not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
JP2010190772A (en) | 2010-09-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CA2934636C (en) | Method for increasing the situational awareness and the location detection of obstacles in the presence of aerosol clouds | |
US11422255B2 (en) | Slope stability LiDAR | |
EP3070501B1 (en) | Prediction of ice crystal presence in a volume of airspace | |
US5523759A (en) | In flight doppler weather radar wind shear detection system | |
Sabatini et al. | A laser obstacle warning and avoidance system for manned and unmanned aircraft | |
JP5511196B2 (en) | Back turbulence detector | |
JP2010181406A (en) | System and method for displaying radar-measured turbulence intensity on vertical display | |
JP4806949B2 (en) | Laser radar equipment | |
Barbaresco et al. | Radar monitoring of a wake vortex: Electromagnetic reflection of wake turbulence in clear air | |
CN111399535A (en) | Unmanned aerial vehicle obstacle avoidance method and device, unmanned aerial vehicle and storage medium | |
US9075074B2 (en) | Flow determination method | |
CN103983975B (en) | Method and system for detecting atmospheric vertical movement speed based on two types of radars | |
Kaikwong et al. | Aircraft wake fortex observations in Hong Kong | |
JP7386136B2 (en) | Cloud height measurement device, measurement point determination method, and cloud type determination method | |
JP2020038135A (en) | Wind tunnel facility and wind tunnel testing method | |
Schwarz et al. | A full-scale rotor-wake investigation of a free-flying helicopter in ground effect using BOS and PIV | |
EP3339184B1 (en) | Apparatus and method for detecting stall condition | |
JP2007232695A (en) | Method of detecting turbulence | |
EP0695426B1 (en) | In flight doppler weather radar wind shear detection system | |
JP2016090409A (en) | Bird observation method | |
Li et al. | Circulation retrieval of simulated wake vortices under rainy condition with a side-looking scanning radar | |
JP3933926B2 (en) | Turbulence detection device and turbulence detection method | |
Jacob et al. | Assessment of Lockheed Martin's aircraft wake vortex circulation estimation algorithms using simulated lidar data | |
Shen et al. | Aircraft wake recognition and strength classification based on deep learning | |
Wassaf et al. | Wake Vortex Tangential Velocity Adaptive Spectral (TVAS) algorithm for pulsed Lidar systems. |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
A621 | Written request for application examination |
Free format text: JAPANESE INTERMEDIATE CODE: A621 Effective date: 20110914 |
|
A977 | Report on retrieval |
Free format text: JAPANESE INTERMEDIATE CODE: A971007 Effective date: 20121114 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20130108 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20130301 |
|
A131 | Notification of reasons for refusal |
Free format text: JAPANESE INTERMEDIATE CODE: A131 Effective date: 20131210 |
|
A521 | Written amendment |
Free format text: JAPANESE INTERMEDIATE CODE: A523 Effective date: 20140203 |
|
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: 20140225 |
|
A61 | First payment of annual fees (during grant procedure) |
Free format text: JAPANESE INTERMEDIATE CODE: A61 Effective date: 20140325 |
|
R150 | Certificate of patent or registration of utility model |
Ref document number: 5511196 Country of ref document: JP Free format text: JAPANESE INTERMEDIATE CODE: R150 |
|
LAPS | Cancellation because of no payment of annual fees |