JP5511196B2 - Back turbulence detector - Google Patents

Back turbulence detector Download PDF

Info

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
Application number
JP2009036369A
Other languages
Japanese (ja)
Other versions
JP2010190772A (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.)
Mitsubishi Electric Corp
Original Assignee
Mitsubishi Electric 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 Mitsubishi Electric Corp filed Critical Mitsubishi Electric Corp
Priority to JP2009036369A priority Critical patent/JP5511196B2/en
Publication of JP2010190772A publication Critical patent/JP2010190772A/en
Application granted granted Critical
Publication of JP5511196B2 publication Critical patent/JP5511196B2/en
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • YGENERAL 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
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information 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 laser radar device 10 in the observation environment of FIG. As shown in FIG. 9, the conventional laser radar device 10 has a cross section (hereinafter referred to as a space) perpendicular to the flight direction (traveling direction) of the aircraft in the atmosphere when the laser radar device 10 is directed to the side of the aircraft. (Referred to as a cross section), and the position and intensity of wake turbulence generated by the passage of the aircraft is detected based on the observation result. Further, the beam of the laser radar apparatus 10 having a constant distance resolution and a cross range resolution (angular resolution) is scanned within the beam scanning range, and is processed for each processing unit (unit consisting of a predetermined distance resolution and a predetermined angular resolution). The wind speed in the sight line direction can be obtained.

レーザーレーダ装置10によって後方乱気流を観測した場合の模擬図を図10に示す。この図10は、観測領域の空間断面における後方乱気流を示している。図10に示すように、レーザーレーダ装置10に対して向かい風を正とした場合、一様に存在する背景風に対して、渦上においては背景風に対して小さい風と大きい風が観測される。   FIG. 10 shows a simulation diagram when wake turbulence is observed by the laser radar device 10. FIG. 10 shows the backward turbulence in the spatial section of the observation region. As shown in FIG. 10, when the head wind is positive with respect to the laser radar device 10, a small wind and a large wind are observed on the vortex with respect to the uniformly existing background wind. .

この観測模擬図を図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 laser radar device 10 is low, a distribution map of Doppler wind speed is created from the observation results of the laser radar device 10. As for the wake turbulence generated sometimes, the wake turbulence with a positive and negative matrix structure is observed with reference to the background wind.

また、等レンジビン上で風速と高度のプロットを作成すると、図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 laser radar device 10. Detection of wake turbulence vortices has been performed by template matching processing using the method (for example, see Patent Document 1).

特開2003−172779号公報JP 2003-1727779 A

深尾昌一郎、浜津享助、「気象と大気のレーダーリモートセンシング」、京都大学学術出版会、2005年3月30日発行Shoichiro Fukao, Kyosuke Hamatsu, “Radar Remote Sensing of Weather and Atmosphere”, Kyoto University Academic Press, March 30, 2005 小松原健史、加来信之、「出発航空機から発する後方乱気流の観測」電子航法研究所研究発表会、第5回、平成17年(2005年)6月Kenfumi Komatsubara, Nobuyuki Kurai, “Observation of wake turbulence from departure aircraft”, Electronic Navigation Research Institute Research Presentation, 5th, June 2005

しかしながら、従来技術には、以下のような課題がある。すなわち、航空機の飛行高度が高い場合、その後方乱気流を観測しようとすると、レーザーレーダ装置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 laser radar device 10 in order to observe the backward turbulence. The distribution map of the Doppler wind speed obtained at this time is as shown in FIG. Unlike the case of FIG. 11 where the elevation angle is low, the vortex tends to cross the range in the case of FIG. 12 where the elevation angle is high. In the case where there is nothing, a conventional matrix template composed of positive and negative is used. When both vortices are observed in the range across the range, the wind speed of both vortices is observed with the laser beam passing between the vortices due to the limitation of the distance resolution of the laser radar device 10, which is remarkable. However, the vortex area that protrudes from the range is smaller than the range, so the background wind speed in the area is more dominant, and the data obtained is the vortex wind speed. Instead, a background value is obtained. For this reason, the vortex cannot be detected by the conventional template matching using a matrix template composed of positive and negative.

本発明は、前記のような課題を解決するためになされたものであり、航空機の飛行高度が高いため、仰角が高い状況で観測している場合においても後方乱気流を検出することができる後方乱気流検出装置を得ることを目的とする。   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.

この発明の実施の形態1に係る後方乱気流検出装置の構成を示すブロック図である。It is a block diagram which shows the structure of the back turbulence detection apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態1に係る後方乱気流検出装置の検出処理部の構成を示すブロック図である。It is a block diagram which shows the structure of the detection process part of the back turbulence detection apparatus which concerns on Embodiment 1 of this invention. この発明の実施の形態2に係る後方乱気流検出装置の検出処理部の構成を示すブロック図である。It is a block diagram which shows the structure of the detection process part of the back turbulence detection apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る後方乱気流検出装置の電磁波送受信部から発射される光パルス(ビーム)の仰角に応じてテンプレートを変化させる様子を示す図である。It is a figure which shows a mode that a template is changed according to the elevation angle of the light pulse (beam) emitted from the electromagnetic wave transmission / reception part of the back turbulence detection apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る後方乱気流検出装置により使用される正と負のマトリクス構造を持つマトリクス型テンプレート例を示す図である。It is a figure which shows the matrix type template example with the positive and negative matrix structure used by the back turbulence detection apparatus which concerns on Embodiment 2 of this invention. この発明の実施の形態2に係る後方乱気流検出装置により使用される正とゼロ(0)のマトリクス構造を持ち、かつ所定の傾きを持つ傾斜マトリクス型テンプレート例を示す図である。It is a figure which shows the example of an inclination matrix type | mold template which has the matrix structure of positive and zero (0) used by the back turbulence detection apparatus which concerns on Embodiment 2 of this invention, and has a predetermined | prescribed inclination. この発明の実施の形態2に係る後方乱気流検出装置により使用される相互相関型テンプレート例を示す図である。It is a figure which shows the example of a cross correlation type template used by the back turbulence detection apparatus which concerns on Embodiment 2 of this invention. 航空機が発生する後方乱気流とそれを観測するレーザーレーダ装置を示す図である。It is a figure which shows the wake turbulence which an aircraft generate | occur | produces, and the laser radar apparatus which observes it. レーザーレーダ装置を用いた後方乱気流の観測例を示す図である。It is a figure which shows the example of observation of a wake turbulence using a laser radar apparatus. 観測領域の空間断面における後方乱気流の概念を示す図である。It is a figure which shows the concept of the wake turbulence in the space cross section of an observation area | region. レーザーレーダ装置を用いた、仰角が低い場合の後方乱気流の観測例を示す図である。It is a figure which shows the example of an observation of back turbulence in case an elevation angle is low using a laser radar apparatus. レーザーレーダ装置を用いた、仰角が高い場合の後方乱気流の観測例を示す図である。It is a figure which shows the example of an observation of the back turbulence in case an elevation angle is high using a laser radar apparatus.

以下、本発明の後方乱気流検出装置の好適な実施の形態につき図面を用いて説明する。   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に係る後方乱気流検出装置の構成を示すブロック図である。なお、各図中、同一符号は同一又は相当部分を示す。
Embodiment 1 FIG.
A backward turbulence detector according to Embodiment 1 of the present invention will be described with reference to FIGS. 1 and 2. The wake turbulence detection apparatus according to Embodiment 1 of the present invention detects wake turbulence using a Doppler wind speed distribution model. 1 is a block diagram showing a configuration of a wake turbulence detection apparatus according to Embodiment 1 of the present invention. In addition, in each figure, the same code | symbol shows the same or equivalent part.

図1において、この発明の実施の形態1に係る後方乱気流検出装置は、大気中に対する電磁波の送受信によって、飛行する航空機の後方に発生する後方乱気流を検出し、検出した後方乱気流に関する情報を表示するもので、航空機の飛行方向に直交する方向に大気中へ電磁波を所定仰角の範囲内で発射した後、その大気からの反射波を一空間断面の信号として受信する電磁波送受信部1と、この電磁波送受信部1からの受信信号に基づき、空間断面のドップラ風速分布を算出し、算出したドップラ風速分布から空間断面における後方乱気流の位置を算出する信号処理部2と、この信号処理部2により算出された空間断面における後方乱気流の位置を表示する表示部3とが設けられている。   In FIG. 1, the wake turbulence detection apparatus according to Embodiment 1 of the present invention detects wake turbulence generated behind a flying aircraft by transmitting and receiving electromagnetic waves to the atmosphere, and displays information on the detected wake turbulence. An electromagnetic wave transmission / reception unit 1 that 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 a reflected wave from the atmosphere as a signal of one spatial section, and the electromagnetic wave Based on the received signal from the transmission / reception unit 1, the Doppler wind speed distribution in the spatial section is calculated, and the signal processing unit 2 that calculates the position of the wake turbulence in the spatial section from the calculated Doppler wind speed distribution, is calculated by the signal processing unit 2. And a display unit 3 for displaying the position of the wake turbulence in the spatial section.

また、信号処理部2は、電磁波送受信部1からの受信信号に基づき、所定の距離分解能及び所定の角度分解能からなる処理単位毎に風速値を算出し、それらの結果より空間断面におけるドップラ風速分布を算出するドップラ風速分布算出部21と、このドップラ風速分布算出部21により算出されたドップラ風速分布と、後方乱気流を模擬するためのパラメータにより算出されたドップラ風速分布モデルとの相互相関処理を行うことにより、後方乱気流を検出する検出処理部22と、相互相関処理の相関値が所定の値である空間断面の所定の位置から、空間断面における後方乱気流の位置を算出する乱気流位置算出部23とが設けられている。   Further, the signal processing unit 2 calculates a wind speed value for each processing unit having a predetermined distance resolution and a predetermined angle resolution based on the received signal from the electromagnetic wave transmission / reception unit 1, and the Doppler wind speed distribution in the spatial section is obtained from these results. A cross-correlation process is performed between the Doppler wind speed distribution calculating unit 21 for calculating the Doppler wind speed distribution calculated by the Doppler wind speed distribution calculated by the Doppler wind speed distribution calculating unit 21 and a parameter for simulating the backward turbulence. Thus, a detection processing unit 22 that detects wake turbulence, a turbulence position calculation unit 23 that calculates the position of wake turbulence in the spatial section from a predetermined position in the spatial section where the correlation value of the cross-correlation processing is a predetermined value, Is provided.

図2は、この発明の実施の形態1に係る後方乱気流検出装置の検出処理部の構成を示すブロック図である。   FIG. 2 is a block diagram showing a configuration of a detection processing unit of the backward turbulence detection device according to Embodiment 1 of the present invention.

図2において、検出処理部22は、後方乱気流を模擬するための航空機に関するパラメータを使用して後方乱気流モデルを含んだドップラ風速分布モデルを算出するドップラ風速分布モデル算出部101と、このドップラ風速分布モデル算出部101により算出したドップラ風速分布モデルとドップラ風速分布算出部21により算出されたドップラ風速分布との相互相関処理により相関値を求める相互相関処理部102と、この相互相関処理部102により求められた相関値のうち最も高い相関値と所定の閾値との比較を行って空間断面における後方乱気流を検出する相関値比較部103と、最も高い相関値が閾値を超えていない場合は、パラメータを変更してドップラ風速分布モデル算出部101へ出力するパラメータ出力部104とが設けられている。   In FIG. 2, the detection processing unit 22 calculates a Doppler wind speed distribution model calculation unit 101 including a wake turbulence model using parameters relating to an aircraft for simulating wake turbulence, and the Doppler wind speed distribution. A cross-correlation processing unit 102 that obtains a correlation value by a cross-correlation process between the Doppler wind speed distribution model calculated by the model calculation unit 101 and the Doppler wind speed distribution calculation unit 21 and a cross-correlation processing unit 102 The correlation value comparison unit 103 that detects the wake turbulence in the spatial section by comparing the highest correlation value with the predetermined threshold value among the obtained correlation values, and if the highest correlation value does not exceed the threshold value, the parameter is A parameter output unit 104 is provided for changing and outputting to the Doppler wind speed distribution model calculation unit 101. It is.

つぎに、この実施の形態1に係る後方乱気流検出装置の動作について図面を参照しながら説明する。   Next, the operation of the wake turbulence detector according to Embodiment 1 will be described with reference to the drawings.

まず、電磁波送受信部1は、図8及び図9に示したように、航空機の飛行方向(進行方向)に直交する方向に大気中へ電磁波である光パルスを所定仰角の範囲内で発射した後、その大気からの反射光を受信し、当該受信した反射光を一空間断面の受信信号として信号処理部2へ送出する。この電磁波送受信部1が受信した反射光には、図10に示すように、反射光の反射位置における風速に応じてドップラ効果が生じている。   First, as shown in FIGS. 8 and 9, the electromagnetic wave transmission / reception unit 1 emits a light pulse, which is an electromagnetic wave, into the atmosphere in a direction perpendicular to the flight direction (traveling direction) of the aircraft within a predetermined elevation angle range. Then, the reflected light from the atmosphere is received, and the received reflected light is sent to the signal processing unit 2 as a reception signal of one spatial section. The reflected light received by the electromagnetic wave transmitting / receiving unit 1 has a Doppler effect according to the wind speed at the reflection position of the reflected light, as shown in FIG.

次に、信号処理部2は、電磁波送受信部1からの受信信号に基づき、空間断面のドップラ風速分布を算出し、算出したドップラ風速分布から空間断面における後方乱気流の位置を算出する。この信号処理部2の詳細な処理については以下に説明する。   Next, the signal processing unit 2 calculates the Doppler wind speed distribution in the spatial section based on the received signal from the electromagnetic wave transmission / reception unit 1, and calculates the position of the wake turbulence in the spatial section from the calculated Doppler wind speed distribution. Detailed processing of the signal processing unit 2 will be described below.

信号処理部2内のドップラ風速分布算出部21は、電磁波送受信部1からの受信信号に基づき、所定の距離分解能及び所定の角度分解能からなる処理単位毎に、風速値を算出する。それらの結果より空間断面におけるドップラ風速分布を算出する。   The Doppler wind speed distribution calculation unit 21 in the signal processing unit 2 calculates a wind speed value for each processing unit having a predetermined distance resolution and a predetermined angle resolution based on the received signal from the electromagnetic wave transmission / reception unit 1. From these results, the Doppler wind speed distribution in the spatial section is calculated.

その後、信号処理部2内の検出処理部22は、観測により得られたドップラ風速分布と、パラメータにより算出されたドップラ風速分布モデルとの相互相関処理を行うことにより、後方乱気流を検出する。この検出処理部22の詳細な処理については以下に説明する。   Thereafter, the detection processing unit 22 in the signal processing unit 2 detects the backward turbulence by performing a cross-correlation process between the Doppler wind speed distribution obtained by the observation and the Doppler wind speed distribution model calculated by the parameters. Detailed processing of the detection processing unit 22 will be described below.

検出処理部22内のドップラ風速分布モデル算出部101は、パラメータ出力部104から、観測を行う航空機の種類(翼の幅)、当該後方乱気流検出装置前を通過する航空機の高度とその速度、大気密度等のパラメータを取得し、このパラメータを使用して処理単位毎にモデル値を求めて、渦モデル(後方乱気流モデル)を含んだドップラ風速分布モデルを算出する。航空機の翼の幅は後方乱気流渦の大きさを与え、航空機の高度は渦が発生すると想定される場所(高度)の参考値となる。航空機の速度や大気密度は渦の強度を算出させることができる。このドップラ風速分布モデルを作成できるものならば上記のパラメータ以外を使用してもよい。   The Doppler wind speed distribution model calculation unit 101 in the detection processing unit 22 receives, from the parameter output unit 104, the type of aircraft to be observed (wing width), the altitude and speed of the aircraft passing in front of the wake turbulence detector, the atmosphere A parameter such as density is acquired, a model value is obtained for each processing unit using this parameter, and a Doppler wind speed distribution model including a vortex model (a wake turbulence model) is calculated. The width of the aircraft wing gives the size of the wake turbulence vortex, and the altitude of the aircraft is a reference value for the location (altitude) where the vortex is expected to occur. Aircraft speed and atmospheric density can be used to calculate vortex strength. Other parameters than the above may be used as long as the Doppler wind speed distribution model can be created.

次に、検出処理部22内の相互相関処理部102は、算出したドップラ風速分布モデルと観測により得られたドップラ風速分布との相互相関処理を行う。その際、処理単位毎に、相互相関処理を行い、分布の処理単位毎に相関係数(相関値)を求め、その結果を相関値比較部103へ送出する。   Next, the cross-correlation processing unit 102 in the detection processing unit 22 performs cross-correlation processing between the calculated Doppler wind speed distribution model and the Doppler wind speed distribution obtained by observation. At this time, cross-correlation processing is performed for each processing unit, a correlation coefficient (correlation value) is obtained for each processing unit of distribution, and the result is sent to the correlation value comparison unit 103.

次に、検出処理部22内の相関値比較部103は、得られた最も高い相関係数(相関値)と閾値との比較を行い、最も高い相関係数が閾値を超えていない場合は、後方乱気流が検出できない場合として、パラメータを変更すべくパラメータ出力部104へ移行する。パラメータを種々変えても、最も高い相関係数が閾値を超えていない場合は、その分布上に後方乱気流がないと判断する。最も高い相関係数が閾値を超えている場合(両者が等しい場合を含む)は、後方乱気流を検出できた場合として、乱気流位置算出部23へ移行する。   Next, the correlation value comparison unit 103 in the detection processing unit 22 compares the obtained highest correlation coefficient (correlation value) with a threshold value, and when the highest correlation coefficient does not exceed the threshold value, When the wake turbulence cannot be detected, the process proceeds to the parameter output unit 104 to change the parameter. If the highest correlation coefficient does not exceed the threshold value even if the parameters are changed, it is determined that there is no wake turbulence on the distribution. When the highest correlation coefficient exceeds the threshold (including the case where both are equal), the process moves to the turbulence position calculation unit 23 as a case where the backward turbulence can be detected.

検出処理部22内のパラメータ出力部104は、当初は予め認識している、観測を行う航空機の種類(翼の幅)、当該後方乱気流検出装置前を通過する航空機の高度とその速度、大気密度等のパラメータをドップラ風速分布モデル算出部101に出力するが、最も高い相関係数が閾値を超えていない場合は、観測を行う航空機の種類(翼の幅)、当該後方乱気流検出装置前を通過する航空機の高度とその速度、大気密度等のパラメータのいずれか、あるいはいくつかを変更してドップラ風速分布モデル算出部101に出力する。このドップラ風速分布モデル算出部101は、変更されたパラメータに基づき新しいドップラ風速分布モデルを作成する。   The parameter output unit 104 in the detection processing unit 22 initially recognizes the type of aircraft to be observed (wing width), the altitude and speed of the aircraft passing in front of the wake turbulence detector, and the atmospheric density. Are output to the Doppler wind speed distribution model calculation unit 101. If the highest correlation coefficient does not exceed the threshold value, the type of aircraft to be observed (the wing width), the turbulence detection device passes through Any one or some of parameters such as altitude and speed of the aircraft to be operated, atmospheric density, etc. are changed and output to the Doppler wind speed distribution model calculation unit 101. The Doppler wind speed distribution model calculation unit 101 creates a new Doppler wind speed distribution model based on the changed parameters.

信号処理部2内の乱気流位置算出部23は、最も高い相関係数が閾値を超えている場合は、その結果から後方乱気流の渦構造の場所を算出し、さらに渦の移動予測を行う。相関位置(処理単位の位置)から、実空間での後方乱気流の位置を算出する。乱気流位置算出部23は、相関係数が高い処理単位の位置から、空間断面における後方乱気流の位置を算出する。後方乱気流の位置は、2つの渦のそれぞれの中心や、2つの渦の間の中心などである。   When the highest correlation coefficient exceeds the threshold value, the turbulence position calculation unit 23 in the signal processing unit 2 calculates the location of the vortex structure of the backward turbulence from the result, and further performs vortex movement prediction. From the correlation position (processing unit position), the position of the wake turbulence in the real space is calculated. The turbulence position calculation unit 23 calculates the position of the backward turbulence in the spatial section from the position of the processing unit having a high correlation coefficient. The position of the wake turbulence is the center of each of the two vortices or the center between the two vortices.

そして、表示部3は、信号処理部2内の乱気流位置算出部23により算出された、実空間での後方乱気流の位置を表示する。空間断面における後方乱気流の位置をディスプレイに表示する。   And the display part 3 displays the position of the back turbulence in real space calculated by the turbulence position calculation part 23 in the signal processing part 2. FIG. The position of the wake turbulence in the spatial section is displayed on the display.

以上のように、実施の形態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と同様である。
Embodiment 2. FIG.
A backward turbulence detection apparatus according to Embodiment 2 of the present invention will be described with reference to FIGS. The backward turbulence detection device according to Embodiment 2 of the present invention detects backward turbulence using a template. FIG. 3 is a block diagram showing a configuration of a detection processing unit of the backward turbulence detection device according to Embodiment 2 of the present invention. The overall configuration of the wake turbulence detection apparatus according to Embodiment 2 of the present invention is the same as that of Embodiment 1 above.

図3において、検出処理部22は、空間断面におけるドップラ風速分布に後方乱気流の検出範囲である所定のフレームを設定するフレーム設定部201と、電磁波送受信部1から発射されるビームの仰角に応じた複数種類のテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を行うテンプレートマッチング部202と、テンプレートマッチング処理により得られた特徴量のうち、所定の特徴量を後方乱気流候補として選出する特徴量算出部203と、背景風の風速値を算出し、前記後方乱気流候補の風速値と前記背景風の風速値の差分から背景風比較値を算出する背景風比較値算出部204と、この背景風比較値算出部204により算出された背景風比較値と所定の閾値との比較を行って空間断面における後方乱気流を検出する背景風比較値比較部205と、背景風比較値が閾値を超えていない場合は、前記テンプレートを変更してテンプレートマッチング部202へ出力するテンプレート出力部206とが設けられている。   In FIG. 3, the detection processing unit 22 sets a predetermined frame that is a detection range of wake turbulence in the Doppler wind speed distribution in the spatial section, and the angle of the beam emitted from the electromagnetic wave transmission / reception unit 1. A template matching unit 202 that performs a template matching process on a set frame using a plurality of types of templates, and a feature quantity that selects a predetermined feature quantity as a posterior turbulence candidate among the feature quantities obtained by the template matching process A background wind comparison value calculation unit 204 that calculates a wind speed value of the background wind, calculates a background wind comparison value from the difference between the wind speed value of the wake turbulence candidate and the wind speed value of the background wind, and the background wind Back turbulence in a spatial section by comparing the background wind comparison value calculated by the comparison value calculation unit 204 with a predetermined threshold value Background Wind comparison value comparison section 205 for detecting, if the background wind comparison value does not exceed the threshold, the template output unit 206 for outputting to the template matching unit 202 to change the template are provided.

つぎに、この実施の形態2に係る後方乱気流検出装置の動作について図面を参照しながら説明する。   Next, the operation of the wake turbulence detector according to Embodiment 2 will be described with reference to the drawings.

まず、検出処理部22内のフレーム設定部201は、空間断面におけるドップラ風速分布に後方乱気流の検出範囲である、航空機の翼の幅で決まるフレーム(ウィンドウとも言う)幅を設定する。後述するテンプレートマッチング部202において、その処理を行わせる範囲を決定するためである。つまり、航空機が通過した場所または後方乱気流があると想定される領域に任意の大きさのフレームを設定する。このフレームは、テンプレートの大きさを示す。なお、このフレーム設定部201をなくし、空間断面におけるドップラ風速分布全体にテンプレートマッチング処理を行わせてもよい。   First, the frame setting unit 201 in the detection processing unit 22 sets a frame (also referred to as a window) width determined by the width of an aircraft wing, which is a detection range of wake turbulence, in the Doppler wind speed distribution in the spatial section. This is because a template matching unit 202 to be described later determines a range in which the processing is performed. That is, a frame of an arbitrary size is set in a place where an aircraft has passed or an area where there is a wake turbulence. This frame indicates the size of the template. The frame setting unit 201 may be omitted, and the template matching process may be performed on the entire Doppler wind speed distribution in the spatial section.

次に、検出処理部22内のテンプレートマッチング部202は、テンプレート出力部206から取得した仰角に応じたテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を行う。その際、テンプレートマッチング処理では、航空機の飛行高度によってテンプレートを図4のように変更する。図4に示すように、後方乱気流検出装置10Aの電磁波送受信部1から発射される光パルス(ビーム)の仰角や、後方乱気流の渦構造によってテンプレートを変化させる。例えば、仰角が高い場合には、テンプレートAやテンプレートBを使用し、仰角が低い場合には、テンプレートCやテンプレートDを使用する。   Next, the template matching unit 202 in the detection processing unit 22 performs a template matching process on the set frame using a template corresponding to the elevation angle acquired from the template output unit 206. At that time, in the template matching process, the template is changed as shown in FIG. 4 according to the flight altitude of the aircraft. As shown in FIG. 4, the template is changed according to the elevation angle of the light pulse (beam) emitted from the electromagnetic wave transmitting / receiving unit 1 of the backward turbulence detection device 10 </ b> A and the vortex structure of the backward turbulence. For example, when the elevation angle is high, template A or template B is used, and when the elevation angle is low, template C or template D is used.

次に、検出処理部22内の特徴量算出部203は、設定された航空機の翼の幅で決まるフレームをドップラ風速分布上で任意に動かし、テンプレートマッチング(相互相関)処理することにより、空間断面内における後方乱気流が存在する位置が強調され、すなわち、処理単位における後方乱気流が存在する位置に対応する積分値が積み上がり、実際に後方乱気流が存在する場合には、そのフレーム部分には渦があるものとして後方乱気流を検出することができる。この特徴量算出部203は、処理単位毎に得られた積分値(特徴量)の中から最も大きい特徴量、または所定の相関値(特徴量を算出するための閾値)を超えた特徴量を後方乱気流候補(渦候補)として選出する。フレームサイズは任意で、上記以外の手法を用いてもよい。   Next, the feature amount calculation unit 203 in the detection processing unit 22 arbitrarily moves a frame determined by the set aircraft wing width on the Doppler wind speed distribution and performs template matching (cross-correlation) processing, thereby obtaining a spatial cross section. The position where the wake turbulence exists is emphasized, that is, the integral value corresponding to the position where the wake turbulence exists in the processing unit is accumulated, and when the wake turbulence actually exists, vortex is generated in the frame portion. It is possible to detect wake turbulence as being. The feature amount calculation unit 203 calculates a feature amount that exceeds the largest feature amount or a predetermined correlation value (threshold for calculating the feature amount) from the integral values (feature amounts) obtained for each processing unit. Selected as a wake turbulence candidate (vortex candidate). The frame size is arbitrary, and methods other than those described above may be used.

その後、検出処理部22内の背景風比較値算出部204は、まず、処理単位毎に背景風の風速値を算出する。背景風の風速値の算出手法については、各レンジビンの風速の平均や、その他過去の観測によって得られた渦構造のないドップラ風速分布、または風速の高度分布モデルといった背景風の風速値を得られるものならばどのような手法でもよい。上述したように、テンプレートマッチング部202によってテンプレートマッチング処理を行って、特徴量算出部203により選出された特徴量を後方乱気流候補(渦候補)とする。   Thereafter, the background wind comparison value calculation unit 204 in the detection processing unit 22 first calculates the wind speed value of the background wind for each processing unit. Regarding the method of calculating the wind speed value of the background wind, it is possible to obtain the wind speed value of the background wind such as the average wind speed of each range bin, the Doppler wind speed distribution without vortex structure obtained by past observations, or the altitude distribution model of the wind speed. Any method can be used. As described above, the template matching process is performed by the template matching unit 202, and the feature quantity selected by the feature quantity calculation unit 203 is set as a posterior turbulence candidate (vortex candidate).

背景風比較値算出部204は、次に、その後方乱気流候補(渦候補)と背景風の風速値から、次の式(1)を用いて処理単位毎に背景風比較値(m/s)を算出する。   Next, the background wind comparison value calculation unit 204 uses the following equation (1) to calculate the background wind comparison value (m / s) from the wake turbulence candidate (vortex candidate) and the wind speed value of the background wind. Is calculated.

Figure 0005511196
Figure 0005511196

ここで、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 template output unit 206 will be described. As a template example, when the frame altitude is low, that is, when the elevation angle of the wake turbulence detection device 10A is low, a matrix template having a positive and negative matrix structure as shown in FIG. FIGS. 5 (2) and (3) show a case where two vortices of wake turbulence are adjacent to each other in the matrix template, and a large weight is applied to a region where the vortex wind speed change is high, that is, the end of the vortex structure, Since the change in vortex wind speed is small at the positive and negative boundaries, a weighting function that applies a small weight is applied. Further, FIG. 5 (4) assumes a case where the vortex structures are separated, and is a matrix template having a positive and negative matrix structure, but is a matrix template divided into two.

図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 turbulence detection device 10A is high, positive and zero (as shown in FIGS. 6 (1)-(4)). An inclined matrix template having a matrix structure of 0) and having a predetermined inclination is used. As the predetermined inclination, not only the inclination shown in FIG. 6 but also an inclination matrix template having various inclinations are prepared. An inclined matrix template having a positive and negative matrix structure and a predetermined inclination is used. In this template, a negative value is assigned to the position of zero (0) shown in FIGS. The second embodiment is characterized in that the template shown in FIGS. 5 and 6 is used by changing the elevation angle of the backward turbulence detection device 10A.

テンプレート出力部206から出力されるテンプレートで、当該後方乱気流検出装置10Aの仰角が高い場合は、図6(1)のように、正の部分のみを特徴量として抽出を行う。また、観測状況、渦の構造によっては2つの正の部分の位置関係は変化することも考慮すべきであるため、図6(2)のように、2つの正の部分を任意にずらした傾斜マトリクス型テンプレートも用意する。さらに、図6(3)及び(4)のように、抽出精度向上のため、渦構造の中心部分に大きい重みがかかるような重み関数をかけた、傾斜マトリクス型を用意する。これにより、観測されたデータの渦構造の不要部分による影響を抑止できる。   In the template output from the template output unit 206, when the elevation angle of the backward turbulence detection device 10A is high, only the positive part is extracted as the feature amount as shown in FIG. In addition, it should be taken into account that the positional relationship between the two positive parts changes depending on the observation situation and the structure of the vortex. Therefore, as shown in FIG. A matrix template is also prepared. Furthermore, as shown in FIGS. 6 (3) and 6 (4), in order to improve the extraction accuracy, an inclined matrix type is prepared by applying a weight function that applies a large weight to the central portion of the vortex structure. Thereby, the influence by the unnecessary part of the vortex structure of the observed data can be suppressed.

テンプレート出力部206は、図6のような汎用型である、傾斜マトリクス型テンプレートの他にも、仰角に応じた過去の観測例または観測される渦のモデルをテンプレートとして出力することも可能である。図7(2)のように得られる等レンジビンプロットの過去の渦観測または渦モデルを用意し、観測結果に対して等レンジビンプロットの相互相関をとらせることによって、渦を検出することも可能である。また、渦を観測している図7(2)のような等レンジビンプロットの渦部分のみを切り出し、それを図6(1)−(4)の正の部分にテンプレートとして使用し、相互相関と同様な効果を得ることも可能である。その際、上述した重みをかけてもよい。図7(3)のような相互相関型テンプレートは、図6に示す傾斜マトリクス型テンプレートだけでなく、図5に示すマトリクス型テンプレートにも適用することが可能である。   The template output unit 206 can output a past observation example or an observed vortex model corresponding to the elevation angle as a template in addition to the general-purpose tilt matrix template as shown in FIG. . It is also possible to detect vortices by preparing past eddy observations or eddy models of equi-range bin plots as shown in Fig. 7 (2), and taking the cross-correlation of equi-range bin plots with the observation results. is there. Further, only the vortex portion of the equi-range bin plot as shown in FIG. 7 (2) observing the vortex is cut out and used as a template in the positive portion of FIGS. Similar effects can be obtained. At that time, the above-described weight may be applied. The cross-correlation template as shown in FIG. 7 (3) can be applied not only to the inclined matrix template shown in FIG. 6, but also to the matrix template shown in FIG.

背景風比較値算出部204は、上述したように、式(1)を用いて背景風比較値を算出する。図5に示すマトリクス型テンプレートの正領域/負領域や、図6に示す傾斜マトリクス型テンプレートの正領域においては、背景風比較値が大きければ大きい程有意となる。後述する背景風比較値比較部205は、この背景風比較値が閾値を超えている、または最も大きい場合は、それを後方乱気流の渦として検出する。ただし、当該後方乱気流検出装置10Aの仰角が高い場合には、すなわち図6及び図7の背景風比較領域においては、背景風と同等の値が得られていることが後方乱気流の条件であるため、背景風比較値が0に近いほど後方乱気流であると判断する。   The background wind comparison value calculation unit 204 calculates the background wind comparison value using the equation (1) as described above. In the positive region / negative region of the matrix template shown in FIG. 5 and the positive region of the inclined matrix template shown in FIG. 6, the larger the background wind comparison value, the more significant. The background wind comparison value comparison unit 205 described later detects the background wind comparison value as a vortex of the backward turbulence when the background wind comparison value exceeds or is the maximum. However, when the backward turbulence detection device 10A has a high elevation angle, that is, in the background wind comparison region of FIGS. 6 and 7, it is the condition of the backward turbulence that a value equivalent to the background wind is obtained. When the background wind comparison value is closer to 0, it is determined that the turbulence is backward.

背景風比較値比較部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 template output unit 206 to change the template as a case where the wake turbulence cannot be detected. If the background wind comparison value does not exceed the threshold value even if the template is changed variously, it is determined that there is no wake turbulence on the spatial section. When the background wind comparison value exceeds the threshold (including the case where both are equal), the process proceeds to the turbulence position calculation unit 23 as a case where the backward turbulence can be detected.

検出処理部22内のテンプレート出力部206は、当初は予め認識している、電磁波送受信部1から発射される光パルス(ビーム)の仰角に応じたテンプレートをテンプレートマッチング部202へ出力するが、背景風比較値が閾値を超えていない場合は、当初と異なる仰角に対応した異なるテンプレートや、観測した後方乱気流の渦構造によって異なるテンプレートに変更してテンプレートマッチング部202へ出力する。このテンプレートマッチング部202は、変更されたテンプレートを使用して、設定されたフレームにテンプレートマッチング処理を再び行う。   The template output unit 206 in the detection processing unit 22 outputs a template according to the elevation angle of the light pulse (beam) emitted from the electromagnetic wave transmission / reception unit 1 that is initially recognized to the template matching unit 202. When the wind comparison value does not exceed the threshold value, the template is changed to a different template corresponding to an elevation angle different from the initial angle, or a different template depending on the observed vortex structure of the wake turbulence and output to the template matching unit 202. The template matching unit 202 performs the template matching process again on the set frame using the changed template.

なお、背景風比較値が閾値を超えていない場合おいては、フレーム設定部201によりフレームを広げたりして、フレームを変更して再処理を行ってもよい。   If the background wind comparison value does not exceed the threshold, the frame setting unit 201 may widen the frame, change the frame, and perform reprocessing.

乱気流位置算出部23は、背景風比較値が閾値を超えている場合は、その結果から後方乱気流の渦構造の場所を算出し、さらに渦の移動予測を行う。当該背景風比較値の処理単位の位置から、実空間での後方乱気流の位置を算出する。乱気流位置算出部23は、当該背景風比較値に対応する処理単位の位置から、空間断面における後方乱気流の位置を算出する。後方乱気流の位置は、2つの渦のそれぞれの中心や、2つの渦の間の中心などである。   When the background wind comparison value exceeds the threshold value, the turbulence position calculation unit 23 calculates the location of the vortex structure of the backward turbulence from the result, and further performs vortex movement prediction. The position of the wake turbulence in the real space is calculated from the position of the processing unit of the background wind comparison value. The turbulence position calculation unit 23 calculates the position of the backward turbulence in the spatial section from the position of the processing unit corresponding to the background wind comparison value. The position of the wake turbulence is the center of each of the two vortices or the center between the two vortices.

そして、表示部3は、信号処理部2内の乱気流位置算出部23により算出された、実空間での後方乱気流の位置を表示する。空間断面における後方乱気流の位置をディスプレイに表示する。   And the display part 3 displays the position of the back turbulence in real space calculated by the turbulence position calculation part 23 in the signal processing part 2. FIG. The position of the wake turbulence in the spatial section is displayed on the display.

以上のように、実施の形態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 SYMBOLS 1 Electromagnetic wave transmission / reception part, 2 Signal processing part, 3 Display part, 10A Backward turbulence detection apparatus, 21 Doppler wind speed distribution calculation part, 22 Detection processing part, 23 Turbulence position calculation part, 101 Doppler wind speed distribution model calculation part, 102 Cross correlation process , 103 correlation value comparison unit, 104 parameter output unit, 201 frame setting unit, 202 template matching unit, 203 feature amount calculation unit, 204 background wind comparison value calculation unit, 205 background wind comparison value comparison unit, 206 template output unit.

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.
JP2009036369A 2009-02-19 2009-02-19 Back turbulence detector Expired - Fee Related JP5511196B2 (en)

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)

* Cited by examiner, † Cited by third party
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)

* Cited by examiner, † Cited by third party
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

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