JP7597367B2 - Airflow measurement device, airflow measurement method, and airflow measurement program - Google Patents

Airflow measurement device, airflow measurement method, and airflow measurement program Download PDF

Info

Publication number
JP7597367B2
JP7597367B2 JP2021043041A JP2021043041A JP7597367B2 JP 7597367 B2 JP7597367 B2 JP 7597367B2 JP 2021043041 A JP2021043041 A JP 2021043041A JP 2021043041 A JP2021043041 A JP 2021043041A JP 7597367 B2 JP7597367 B2 JP 7597367B2
Authority
JP
Japan
Prior art keywords
reflected light
light
airflow
virtual surface
reflected
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.)
Active
Application number
JP2021043041A
Other languages
Japanese (ja)
Other versions
JP2022142821A (en
Inventor
正志 原田
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Japan Aerospace Exploration Agency JAXA
Original Assignee
Japan Aerospace Exploration Agency JAXA
Priority date (The priority date is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the date listed.)
Filing date
Publication date
Application filed by Japan Aerospace Exploration Agency JAXA filed Critical Japan Aerospace Exploration Agency JAXA
Priority to JP2021043041A priority Critical patent/JP7597367B2/en
Publication of JP2022142821A publication Critical patent/JP2022142821A/en
Application granted granted Critical
Publication of JP7597367B2 publication Critical patent/JP7597367B2/en
Active 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

  • Optical Radar Systems And Details Thereof (AREA)

Description

本開示は、例えば航空機が遭遇する乱気流の検知に用いられる気流計測装置、気流計測方法及び気流計測プログラムに関する。 The present disclosure relates to an airflow measurement device, an airflow measurement method, and an airflow measurement program that are used, for example, to detect turbulence encountered by an aircraft.

晴天時に航空機が遭遇する強い乱気流である晴天乱気流(clear air turbulence)は、航空機の乗り心地に影響を与えるだけでなく、しばしば機体破壊の原因となる。晴天乱気流はジェット気流、山岳波などによる速度勾配がケルビン・ヘルムホルツ(Kelvin-Helmhollz)不安定を主な原因として乱流へと移行する際に発生する。時間とともに乱流のスケールが小さくなるほど、乱流のエネルギーが小さくなる。よって発生初期の規模の大きい乱れが機体への影響が大きい。このスケール及びエネルギーの大きな乱流を検知してこれを回避して航行することが望まれる。 Clear air turbulence, the strong turbulence that aircraft encounter on clear days, not only affects the comfort of the aircraft, but is often the cause of aircraft damage. Clear air turbulence occurs when the velocity gradient caused by the jet stream, mountain waves, etc., transitions into turbulence, primarily due to Kelvin-Helmholtz instability. As the scale of the turbulence becomes smaller over time, the energy of the turbulence decreases. Therefore, large-scale turbulence at the beginning of its development has a large impact on aircraft. It is desirable to detect turbulence of this scale and large energy and navigate by avoiding it.

航空機が乱気流の中に入ると大きく揺れ、急激な上昇及び下降を繰り返すため機体に大きな荷重が加わり、最悪の場合破壊に至ることがある。また、シートベルトのサインが遅れたために機内の乗客や乗員が重軽傷を負う事故も毎年に何例か発生している。そして、機体への負担の軽減のほか、機内サービスや乗客の乗り心地といった快適性、乗客や乗員の身体や生命に影響を及ぼすため、パイロットが前方に乱気流があることを発見したときは、可能な限り回避しなければならない。 When an aircraft enters turbulence, it shakes violently and repeatedly ascends and descends rapidly, placing a large load on the aircraft, which in the worst case scenario can lead to destruction. There are also several accidents each year in which passengers and crew suffer serious or minor injuries due to delayed seatbelt signs. In addition to reducing the burden on the aircraft, turbulence can also affect in-flight service, passenger comfort, and the physical and physical health and lives of passengers and crew, so when pilots notice turbulence ahead, they must avoid it as much as possible.

特許文献1では、レーザ光を前方に照射しこの反射光を受信することで、時間遅れ及びドップラーシフトから光軸方向の距離及び速度を知る。さらに詳細な前方の情報を得るために10度程度の角度で上下2方向に、レーザ光を照射し反射光を受信する。反射光の時系列などからもっともらしい前方の速度場を推定する。 In Patent Document 1, laser light is emitted forward and the reflected light is received, and the distance and speed in the optical axis direction are determined from the time delay and Doppler shift. In order to obtain more detailed information about the road ahead, laser light is emitted in two directions, up and down, at an angle of about 10 degrees, and the reflected light is received. A plausible forward velocity field is estimated from the time series of the reflected light, etc.

特開2017-67680号公報JP 2017-67680 A

特許文献1の方法で得られる上下方向の気流速度は、光軸方向の速度に上下に10度程度振る角度の正弦(sin)を掛けたものであるため、精度が低いおそれがある。また、これから通過する予定の飛行経路上の気流速度ではなく、これから通過する予定の空域の上方の領域及び下方の領域の気流速度が算出される。このため、航空機が乱気流を回避するために必要な空域の気流速度を適切に算出できないおそれがある。 The vertical airflow speed obtained by the method of Patent Document 1 is calculated by multiplying the speed in the optical axis direction by the sine of an angle of about 10 degrees up and down, so there is a risk of low accuracy. In addition, the airflow speed calculated is not the airflow speed on the flight path that the aircraft is scheduled to pass through, but the airflow speed in the areas above and below the airspace that the aircraft is scheduled to pass through. For this reason, there is a risk that the airflow speed in the airspace necessary for the aircraft to avoid turbulence cannot be properly calculated.

以上のような事情に鑑み、本開示の目的は、航空機が乱気流を回避するために必要な空域の気流速度を精度よく適切に算出することが可能な気流計測装置、気流計測方法及び気流計測プログラムを提供することにある。 In view of the above circumstances, the objective of the present disclosure is to provide an airflow measurement device, an airflow measurement method, and an airflow measurement program that can accurately and appropriately calculate the airflow speed in an airspace required for an aircraft to avoid turbulence.

本開示の一形態に係る気流計測装置は、
照射装置が光を前記光の光軸に交差する大気中の仮想面に照射してから、受光装置が前記光の反射光を受光するまでの時間遅れを算出し、
前記時間遅れに基づき、前記仮想面で反射した反射光を抽出する
算出装置と、
前記仮想面に対応し複数のグリッドエリアを含むマップの各グリッドエリアに、前記反射光の強度をプロットすることで、反射光の強度のパターンを示す反射光マップを作成し、
前記反射光マップと、所定時間経過後に作成された反射光マップとの前記パターンの変化を演算し、
前記パターンの変化に基づき、前記仮想面上の前記光軸に交差する方向の気流の速度を計測する
演算装置と
を具備する。
An airflow measuring device according to an embodiment of the present disclosure includes:
Calculating a time delay from when an irradiating device irradiates light onto a virtual surface in the atmosphere that intersects with an optical axis of the light until when a light receiving device receives reflected light of the light;
A calculation device that extracts the reflected light reflected by the virtual surface based on the time delay;
creating a reflected light map showing a pattern of the intensity of the reflected light by plotting the intensity of the reflected light in each grid area of a map corresponding to the virtual surface and including a plurality of grid areas;
Calculating a change in the pattern between the reflected light map and a reflected light map created after a predetermined time has elapsed;
and a calculation device that measures the speed of the airflow on the virtual surface in a direction intersecting the optical axis based on the change in the pattern.

本実施形態によれば、気流計測装置は、照射から受光までの反射光の時間遅れを算出し、時間遅れに基づき、光の光軸に交差する仮想面で反射した反射光を抽出し、仮想面で反射した反射光の強度のパターンを示す反射光マップを作成する。これにより、速度を検出すべき流れの方向と一致する面方向を有する仮想面で反射した反射光の反射光マップを作成することができる。 According to this embodiment, the airflow measurement device calculates the time delay of the reflected light from irradiation to reception, extracts the reflected light reflected on a virtual surface that intersects with the optical axis of the light based on the time delay, and creates a reflected light map that shows the intensity pattern of the reflected light reflected on the virtual surface. This makes it possible to create a reflected light map of the reflected light reflected on a virtual surface that has a surface direction that matches the direction of the flow whose speed is to be detected.

前記複数の反射光マップに対応する前記仮想面は、大気中の同じ位置にあり、
前記演算装置は、前記光軸に交差する方向として、前記仮想面の面方向の気流の速度を計測してもよい。
the virtual surfaces corresponding to the plurality of reflected light maps are at the same position in the atmosphere;
The calculation device may measure the speed of the airflow in a planar direction of the imaginary surface as a direction intersecting the optical axis.

本実施形態によれば、時間的に異なる複数の、大気中の同じ位置にある仮想面で反射した反射光の反射光マップを作成することで、仮想面の面方向の気流の速度を計測することができる。 According to this embodiment, by creating a reflected light map of light reflected from multiple virtual surfaces at the same position in the atmosphere at different times, it is possible to measure the airflow speed in the surface direction of the virtual surfaces.

前記仮想面は、前記仮想面上の1点に照射される光の光軸に直交してもよい。 The virtual plane may be perpendicular to the optical axis of the light irradiated to a point on the virtual plane.

本実施形態によれば、光の光軸に直交する面方向の気流の速度を計測することができる。 According to this embodiment, it is possible to measure the speed of airflow in a plane direction perpendicular to the optical axis of the light.

前記気流計測装置は、航空機に搭載されてもよい。 The airflow measuring device may be mounted on an aircraft.

本実施形態によれば、航空機から照射される光の光軸に直交する面方向の気流の速度を計測することができる。 According to this embodiment, it is possible to measure the speed of airflow in a plane perpendicular to the optical axis of the light emitted from the aircraft.

前記仮想面は、飛行中である前記航空機が通過する予定である飛行経路を含んでもよい。 The virtual surface may include a flight path along which the aircraft is scheduled to fly.

本実施形態によれば、航空機が通過する予定である飛行経路に直交する面方向の気流の速度を計測することができる。特に、晴天乱気流の上下方向の速度を計測することができる。また、仮想面のサイズを適切に設定することで、航空機が乱気流を回避するために必要な空域の気流速度を精度よく適切に算出することを図れる。 According to this embodiment, it is possible to measure the air current speed in a plane perpendicular to the flight path that the aircraft is scheduled to take. In particular, it is possible to measure the vertical speed of clear air turbulence. In addition, by appropriately setting the size of the virtual surface, it is possible to accurately and appropriately calculate the air current speed in the airspace required for the aircraft to avoid turbulence.

前記仮想面は、前記仮想面上に乱気流が存在する場合に前記航空機が前記乱気流を回避するために飛行する飛行経路を含んでもよい。 The virtual surface may include a flight path that the aircraft will fly to avoid turbulence if turbulence is present on the virtual surface.

本実施形態によれば、仮想面上の予定される飛行経路に乱気流が存在する場合に、航空機は乱気流を回避して仮想面内の安全な(即ち、乱気流が存在しない)空域を飛行することが可能である。 According to this embodiment, when turbulence is present on a planned flight path on a virtual plane, the aircraft can avoid the turbulence and fly in a safe (i.e., turbulence-free) airspace on the virtual plane.

前記算出装置は、前記航空機が飛行中に飛行経路上の異なる位置から光を照射する毎に、抽出する反射光の時間遅れを異ならせてもよい。 The calculation device may vary the time delay of the extracted reflected light each time the aircraft irradiates light from a different position on the flight path while in flight.

遅れ時間の時間長さは、航空機から仮想面までの距離に依存して変更される。言い換えれば、距離は、航空機が仮想面に近づくに従って小さくなるため、遅れ時間の時間長さも短くなる。算出装置は、より短時間で戻って来た反射光を抽出し続けることで、航空機の位置に依存せず、大気中の同じ位置にある仮想面で反射した反射光を抽出し続けることができる。 The length of the delay time varies depending on the distance from the aircraft to the virtual surface. In other words, the distance decreases as the aircraft approaches the virtual surface, and therefore the length of the delay time also decreases. By continuing to extract reflected light that returns in a shorter time, the calculation device can continue to extract reflected light that is reflected by the virtual surface at the same position in the atmosphere, regardless of the position of the aircraft.

前記演算装置は、
前記複数の反射光マップから、前記仮想面上の同じエリアに対応する検査区画をそれぞれ抽出し、
前記複数の反射光マップからそれぞれ抽出した複数の前記検査区画の前記パターンの変化を演算することにより、
前記複数の反射光マップの前記パターンの変化を演算してもよい。
The computing device includes:
Extracting inspection sections from the plurality of reflected light maps, each of which corresponds to a same area on the virtual surface;
Calculating changes in the patterns of the plurality of inspection sections extracted from the plurality of reflected light maps, respectively;
A change in the pattern of the plurality of reflected light maps may be calculated.

反射光の強度がプロットされたグリッドエリアの数を適切に設定することで、複数の検査区画のパターンの変化をより正確に演算することができる。 By appropriately setting the number of grid areas in which the reflected light intensity is plotted, the pattern changes in multiple inspection sections can be calculated more accurately.

前記検査区画に含まれる前記グリッドエリアの数は可変であり、
前記検査区画に含まれる、反射光の強度がプロットされたグリッドエリアの数は、数個以上数十個以下でもよい。
The number of the grid areas included in the inspection section is variable;
The number of grid areas, in which the intensity of reflected light is plotted, included in the inspection section may be from several to several tens.

グリッドエリアの数は、演算負荷やパイロットの人的負荷を軽くすることを重視するならば固定とすればよいし、演算の正確性を重視するならば可変とすればよい。 The number of grid areas can be fixed if the emphasis is on reducing the computational load and the pilot's physical burden, or it can be variable if the emphasis is on calculation accuracy.

前記照射装置は、
前記仮想面を分割したエリアと同形の断面形状を有するレーザ光を発振するレーザ発振装置と、
前記レーザ光の照射方向を制御して、前記レーザ光で前記仮想面の全てのエリアをスキャンするスキャン装置と
を有してもよい。
The irradiation device includes:
a laser oscillation device that oscillates a laser beam having a cross-sectional shape identical to that of an area obtained by dividing the imaginary surface;
and a scanning device that controls an irradiation direction of the laser light and scans all areas of the virtual surface with the laser light.

これにより、レーザ光が仮想面の全域を短時間で万遍無く照射することができる。 This allows the laser light to irradiate the entire virtual surface evenly in a short period of time.

前記算出装置又は前記演算装置は、前記レーザ光で前記仮想面をスキャンすることにより生じる照射タイミングのずれに応じて時間遅れのずれを補正してもよい。 The calculation device or the arithmetic device may correct the time delay deviation according to the deviation in irradiation timing caused by scanning the virtual surface with the laser light.

照射装置は、仮想面の全域を一度に照射するのではなく、仮想面の全域を順にスキャンすることにより仮想面の全域にレーザ光を照射する。このため、最初にスキャンする矩形エリアから、最後にスキャンする矩形エリアまで、厳密には、照射タイミングに時間的なずれがある。従って、最初にスキャンするタイミングから最後にスキャンするタイミングまで、航空機の位置が厳密には異なる。このため、照射タイミングの時間的なずれに依存して、仮想面が、進行方向(即ち、仮想面の中心の光軸方向)に直交ではなく、進行方向の成分を含む方向に歪んだ(即ち、進行方向に対して傾斜した)形になり得る。算出装置又は演算装置は、この様な照射タイミングのずれに応じて時間遅れのずれを補正することで、仮想面を進行方向に直交する面となるように補正してもよい。 The irradiation device does not irradiate the entire virtual surface at once, but irradiates the entire virtual surface with laser light by scanning the entire virtual surface in sequence. Therefore, strictly speaking, there is a time lag in the irradiation timing from the rectangular area scanned first to the rectangular area scanned last. Therefore, the position of the aircraft is strictly different from the timing of the first scan to the timing of the last scan. Therefore, depending on the time lag in the irradiation timing, the virtual surface may not be perpendicular to the traveling direction (i.e., the optical axis direction at the center of the virtual surface), but may be distorted in a direction that includes a component of the traveling direction (i.e., tilted with respect to the traveling direction). The calculation device or arithmetic device may correct the time delay deviation according to such a deviation in the irradiation timing, so that the virtual surface becomes a surface perpendicular to the traveling direction.

前記算出装置は、さらに、照射された光と、前記仮想面で反射した反射光との間の周波数のドップラーシフト量を算出し、
前記演算装置は、さらに、前記ドップラーシフト量に基づき、前記仮想面上の前記光軸の方向の気流の速度を計測してもよい。
The calculation device further calculates a Doppler shift amount of a frequency between the irradiated light and the reflected light reflected by the virtual surface,
The calculation device may further measure a speed of airflow in a direction of the optical axis on the virtual surface based on the amount of Doppler shift.

演算装置は、記憶装置から、仮想面で反射した反射光のドップラーシフト量と、位置情報とを読み出す。演算装置は、ドップラーシフト量に基づき、仮想面上の光軸の方向(即ち、航空機の進行方向)の気流の速度を計測する。本実施形態によれば、仮想面の面方向の気流の速度を計測するのに用いられる反射光を利用して、航空機の進行方向の気流の速度も計測することができる。 The arithmetic device reads out the Doppler shift amount of the light reflected from the virtual surface and the position information from the storage device. The arithmetic device measures the airflow speed in the direction of the optical axis on the virtual surface (i.e., the direction of travel of the aircraft) based on the Doppler shift amount. According to this embodiment, the reflected light used to measure the airflow speed in the planar direction of the virtual surface can also be used to measure the airflow speed in the direction of travel of the aircraft.

気流計測装置は、前記照射装置と、前記受光装置とをさらに具備してもよい。 The airflow measuring device may further include the irradiation device and the light receiving device.

本実施形態では、最小限のハードウェアであるレーザ発振装置及び受光装置を航空機に搭載するだけでよく、処理装置は受光装置が受光した反射光に基づき、ソフトウェア処理により仮想面の面方向の気流の速度を計測する。また、同時にドップラーシフト量に基づく進行方向の気流速度を検出することも可能である。 In this embodiment, the aircraft only needs to be equipped with the minimum amount of hardware, a laser oscillator and a light receiving device, and the processing device uses software processing to measure the airflow speed in the surface direction of the virtual surface based on the reflected light received by the light receiving device. It is also possible to simultaneously detect the airflow speed in the direction of travel based on the amount of Doppler shift.

本開示の一形態に係る気流計測方法は、
照射装置が光を前記光の光軸に交差する大気中の仮想面に照射してから、受光装置が前記光の反射光を受光するまでの時間遅れに基づき、前記仮想面で反射した反射光を抽出し、
前記仮想面に対応し複数のグリッドエリアを含むマップの各グリッドエリアに、前記反射光の強度をプロットすることで、反射光の強度のパターンを示す反射光マップを作成し、
複数の前記反射光マップの前記パターンの変化を演算し、
前記パターンの変化に基づき、前記仮想面上の前記光軸に交差する方向の気流の速度を計測する。
An airflow measuring method according to an embodiment of the present disclosure includes:
extracting light reflected by the virtual surface based on a time delay between when an irradiating device irradiates light onto a virtual surface in the atmosphere that intersects with an optical axis of the light and when a light receiving device receives reflected light of the light;
creating a reflected light map showing a pattern of the intensity of the reflected light by plotting the intensity of the reflected light in each grid area of a map corresponding to the virtual surface and including a plurality of grid areas;
Calculating changes in the pattern of a plurality of the reflected light maps;
Based on the change in the pattern, the speed of the airflow on the virtual surface in a direction intersecting the optical axis is measured.

本開示の一形態に係る気流計測プログラムは、
気流計測装置の演算装置を、
照射装置が光を前記光の光軸に交差する大気中の仮想面に照射してから、受光装置が前記光の反射光を受光するまでの時間遅れに基づき、前記仮想面で反射した反射光を抽出し、
前記仮想面に対応し複数のグリッドエリアを含むマップの各グリッドエリアに、前記反射光の強度をプロットすることで、反射光の強度のパターンを示す反射光マップを作成し、
複数の前記反射光マップの前記パターンの変化を演算し、
前記パターンの変化に基づき、前記仮想面上の前記光軸に交差する方向の気流の速度を計測する
よう動作させる。
An airflow measurement program according to an embodiment of the present disclosure,
A calculation unit of the airflow measuring device,
extracting light reflected by the virtual surface based on a time delay between when an irradiating device irradiates light onto a virtual surface in the atmosphere that intersects with an optical axis of the light and when a light receiving device receives reflected light of the light;
creating a reflected light map showing a pattern of the intensity of the reflected light by plotting the intensity of the reflected light in each grid area of a map corresponding to the virtual surface and including a plurality of grid areas;
Calculating changes in the pattern of a plurality of the reflected light maps;
The optical system is operated to measure the speed of airflow on the imaginary surface in a direction intersecting the optical axis based on the change in the pattern.

本開示によれば、航空機が乱気流を回避するために必要な空域の気流速度を精度よく適切に算出することを図れる。 The present disclosure makes it possible to accurately and appropriately calculate the airflow speed in an airspace required for an aircraft to avoid turbulence.

なお、ここに記載された効果は必ずしも限定されるものではなく、本開示中に記載されたいずれかの効果であってもよい。 Note that the effects described here are not necessarily limited to those described herein and may be any of the effects described in this disclosure.

本開示の一実施形態に係る気流計測装置を搭載した航空機を示す。1 shows an aircraft equipped with an airflow measuring device according to an embodiment of the present disclosure. 気流計測装置の構成を示す。The configuration of the airflow measuring device is shown. 気流計測装置の動作フローを示す。4 shows the operation flow of the airflow measuring device. マップの一例を示す。An example of the map is shown below. 複数の反射光マップからそれぞれ抽出した複数の検査区画のパターンの変化を模式的に示す。13A and 13B show schematic diagrams of changes in patterns of a plurality of inspection sections extracted from a plurality of reflected light maps, respectively.

以下、図面を参照しながら、本開示の実施形態を説明する。 Embodiments of the present disclosure will be described below with reference to the drawings.

1.気流計測装置を搭載した航空機の概要 1. Overview of aircraft equipped with airflow measurement equipment

図1は、本開示の一実施形態に係る気流計測装置を搭載した航空機を示す。 Figure 1 shows an aircraft equipped with an airflow measurement device according to one embodiment of the present disclosure.

気流計測装置100は、航空機1に搭載される。航空機1は、例えば、無人機、小型航空機、大型旅客機、電動航空機等である。気流計測装置100は、照射及び受光装置2を有する。照射及び受光装置2は、航空機1の機首下部に設けられる。照射及び受光装置2は、レーザ光を照射し、反射光を受光する。気流計測装置100は、光を利用して気流の速度を遠隔的に計測するLiDAR(Light Detection And Ranging)装置の一種である。 The airflow measurement device 100 is mounted on an aircraft 1. The aircraft 1 may be, for example, an unmanned aircraft, a small aircraft, a large passenger aircraft, an electric aircraft, etc. The airflow measurement device 100 has an illumination and light receiving device 2. The illumination and light receiving device 2 is provided under the nose of the aircraft 1. The illumination and light receiving device 2 illuminates laser light and receives reflected light. The airflow measurement device 100 is a type of LiDAR (Light Detection and Ranging) device that uses light to remotely measure the speed of airflow.

2.気流計測装置の構成 2. Configuration of the airflow measurement device

図2は、気流計測装置の構成を示す。 Figure 2 shows the configuration of the airflow measurement device.

気流計測装置100は、照射及び受光装置2と、処理装置130と、記憶装置140と、表示装置150とを有する。照射及び受光装置2は、照射装置110と、受光装置120とを有する。照射装置110は、レーザ発振装置111と、スキャン装置112とを有する光学系である。処理装置130は、プロセッサがROMに記憶された情報処理プログラム(気流計測プログラムを含む)をRAMにロードして実行することにより、算出装置131及び演算装置132として動作する。記憶装置140は、大容量の不揮発性の記憶媒体であり、例えば、HDD、SSD等である。表示装置150はコックピット内のディスプレイであり、速度表示画面151を表示する。 The airflow measurement device 100 has an irradiation and light receiving device 2, a processing device 130, a storage device 140, and a display device 150. The irradiation and light receiving device 2 has an irradiation device 110 and a light receiving device 120. The irradiation device 110 is an optical system having a laser oscillation device 111 and a scanning device 112. The processing device 130 operates as a calculation device 131 and an arithmetic device 132 by the processor loading an information processing program (including an airflow measurement program) stored in the ROM into the RAM and executing it. The storage device 140 is a large-capacity non-volatile storage medium, such as an HDD or SSD. The display device 150 is a display in the cockpit, and displays a speed display screen 151.

3.気流計測装置の動作フロー 3. Airflow measurement device operation flow

図3は、気流計測装置の動作フローを示す。 Figure 3 shows the operational flow of the airflow measuring device.

ステップS201:N=1及び遅れ時間Tを設定 Step S201: Set N=1 and delay time T

演算装置132は、N=1と、基準となる遅れ時間Tとを設定する。Nは、照射装置110が光を大気中の仮想面に照射する回数である。N=1は、1回目の照射を意味し、また、1回目の動作フローを意味する。基準となる遅れ時間Tは、N=1回目に照射装置110が光を照射してから受光装置120が反射光を受光するまでの遅れ時間の時間長さである。言い換えれば、N=1のとき(即ち、1回目の照射時)の遅れ時間がTである。 The calculation device 132 sets N=1 and a reference delay time T. N is the number of times that the irradiation device 110 irradiates the virtual surface in the atmosphere with light. N=1 means the first irradiation, and also means the first operation flow. The reference delay time T is the length of time delay from when the irradiation device 110 irradiates light for the N=1st time until the light receiving device 120 receives the reflected light. In other words, the delay time when N=1 (i.e., when the first irradiation is performed) is T.

遅れ時間の時間長さは、航空機1から速度計測対象である大気までの距離Lに依存して変化する。照射から受光までの遅れ時間の時間長さは、2L/Cである。Cは光速であり、Lは航空機1から速度計測対象である大気までの距離である。航空機1は前進し続けるため、航空機1から速度計測対象である大気までの距離Lは徐々に小さくなり、これに伴い、遅れ時間の時間長さも徐々に短くなる。 The length of the delay time varies depending on the distance L from the aircraft 1 to the atmosphere, the object of speed measurement. The length of the delay time from illumination to reception is 2L/C, where C is the speed of light and L is the distance from the aircraft 1 to the atmosphere, the object of speed measurement. As the aircraft 1 continues to move forward, the distance L from the aircraft 1 to the atmosphere, the object of speed measurement, gradually becomes smaller, and accordingly, the length of the delay time gradually becomes shorter.

ステップS202:レーザ光を照射 Step S202: Irradiate laser light

照射装置110は、光を大気中の仮想面6に照射する(図1参照)。大気中の仮想面6には、速度計測対象である大気が存在する。具体的には、照射装置110のレーザ発振装置111は、レーザ光を発振する光源を含む。照射装置110のスキャン装置112は、レーザ発振装置111が発振したレーザ光の照射方向を制御して、レーザ光で仮想面6をスキャンする。照射装置110は、dt時間のスキャン間隔で、定期的に仮想面6をスキャンする。 The irradiation device 110 irradiates light onto a virtual surface 6 in the atmosphere (see FIG. 1). The virtual surface 6 in the atmosphere contains the air, the object of which is to be measured for speed. Specifically, the laser oscillation device 111 of the irradiation device 110 includes a light source that oscillates laser light. The scanning device 112 of the irradiation device 110 controls the irradiation direction of the laser light oscillated by the laser oscillation device 111, and scans the virtual surface 6 with the laser light. The irradiation device 110 periodically scans the virtual surface 6 at a scan interval of dt time.

具体的には、レーザ発振装置111は、仮想面6を隙間なく照射できるように、レーザ光を光学的に矩形に整形し、レーザ光の断面形状を見込み角の矩形エリアと一致させる。スキャン装置112は、航空機1の進行方向前方の四角錐を、横10以上数百以下×縦10以上数百以下の見込み角に格子状に分割し、それぞれの見込み角の矩形エリアをレーザ光でスキャンする。四角錐の底面が仮想面6に相当する。スキャン装置112は、ポリゴンミラー等でレーザ発振装置111を左右に振り、レーザ発振装置111が見込み角の矩形エリアに正確に向いた瞬間に、レーザ光の強度をパルス状に変化させる。これを繰り返すことで、スキャン装置112は、レーザ光で仮想面6の全ての矩形エリアをスキャンする。これにより、レーザ光が仮想面6の全域を短時間で万遍無く照射することができる。 Specifically, the laser oscillator 111 optically shapes the laser light into a rectangle so that the virtual surface 6 can be irradiated without gaps, and the cross-sectional shape of the laser light matches the rectangular area of the angle of view. The scanning device 112 divides the pyramid in front of the aircraft 1 in the direction of travel into a grid of angles of view of 10 to several hundred horizontally and 10 to several hundred vertically, and scans the rectangular area of each angle of view with the laser light. The bottom of the pyramid corresponds to the virtual surface 6. The scanning device 112 swings the laser oscillator 111 left and right with a polygon mirror or the like, and changes the intensity of the laser light in a pulsed manner at the moment when the laser oscillator 111 is accurately directed toward the rectangular area of the angle of view. By repeating this, the scanning device 112 scans the entire rectangular area of the virtual surface 6 with the laser light. This allows the laser light to irradiate the entire area of the virtual surface 6 evenly in a short time.

四角錐の頂角θ、θは、例えば、0.1度以上5度以下、より具体的には、いずれも1度である。スキャン装置112は、左右角θ及び上下角θの範囲でレーザ光をスキャンする。例えば、スキャン装置112は、仮想面6の左上の頂点から右上の頂点に左右角θの範囲でレーザ光をスキャンする。続いて、スキャン装置112は、上下角θを格子の1段分下げて、左から右へと左右角θの範囲でレーザ光をスキャンする。スキャン装置112は、これを繰り返すことで仮想面6の全域に光を照射する。航空機1から仮想面6までの距離はLである。 The apex angles θ 1 and θ 2 of the quadrangular pyramid are, for example, 0.1 degrees or more and 5 degrees or less, more specifically, each is 1 degree. The scanning device 112 scans the laser light in the range of the left-right angle θ 1 and the up-down angle θ 2. For example, the scanning device 112 scans the laser light from the top left vertex of the virtual surface 6 to the top right vertex in the range of the left-right angle θ 1. Next, the scanning device 112 lowers the up-down angle θ 2 by one lattice step and scans the laser light from left to right in the range of the left-right angle θ 1. The scanning device 112 repeats this process to irradiate the light over the entire virtual surface 6. The distance from the aircraft 1 to the virtual surface 6 is L.

仮想面6は、照射される全ての光の光軸に交差する。中でも、仮想面6は、仮想面6上の1点(例えば、中心点)に照射される光の光軸に直交する。仮想面6は、飛行中である航空機1が通過する予定である飛行経路を含む。例えば、航空機1が通過する予定である飛行経路が、仮想面6上の中心点又は中心点付近に直交する。仮想面6は、仮想面6上に乱気流が存在する場合に航空機1が乱気流を回避するための飛行経路を含む。例えば、仮想面6は、1辺が10kmの正方形である。仮想面6が1辺10kmの正方形であれば、仮想面6の中心(即ち、予定される飛行経路)に乱気流が存在する場合に、航空機1は乱気流を回避して仮想面6内の安全な(即ち、乱気流が存在しない)空域を飛行することが可能である。 The virtual plane 6 intersects with the optical axes of all the irradiated light. In particular, the virtual plane 6 is perpendicular to the optical axis of the light irradiated to one point (e.g., the center point) on the virtual plane 6. The virtual plane 6 includes a flight path that the aircraft 1 is scheduled to pass through during flight. For example, the flight path that the aircraft 1 is scheduled to pass is perpendicular to the center point or near the center point on the virtual plane 6. The virtual plane 6 includes a flight path for the aircraft 1 to avoid turbulence when turbulence exists on the virtual plane 6. For example, the virtual plane 6 is a square with one side of 10 km. If the virtual plane 6 is a square with one side of 10 km, when turbulence exists at the center of the virtual plane 6 (i.e., the planned flight path), the aircraft 1 can avoid the turbulence and fly in a safe airspace (i.e., no turbulence exists) within the virtual plane 6.

ステップS203:レーザ光を受光 Step S203: Receive laser light

受光装置120は、照射装置110が照射した光の、大気中のエアロゾル(塵、微粒子)からの反射光を受光する。 The light receiving device 120 receives the light emitted by the irradiation device 110 that is reflected from aerosols (dust, fine particles) in the atmosphere.

ステップS204:時間遅れを算出 Step S204: Calculate the time delay

算出装置131は、照射装置110が光を大気中に照射してから、受光装置120が反射光を受光するまでの時間遅れを算出する。照射から受光までの時間遅れは、航空機1から反射源(塵等)までの距離に依存して異なる。 The calculation device 131 calculates the time delay between when the irradiation device 110 irradiates light into the atmosphere and when the light receiving device 120 receives the reflected light. The time delay between when the light is irradiated and when the light is received varies depending on the distance from the aircraft 1 to the reflection source (dust, etc.).

算出装置131は、さらに、照射装置110から照射された光と、仮想面6上で反射し受光装置120が受光した反射光との間の周波数のドップラーシフト量を算出する。 The calculation device 131 further calculates the amount of Doppler shift in frequency between the light irradiated from the irradiation device 110 and the reflected light reflected on the virtual surface 6 and received by the light receiving device 120.

算出装置131は、算出した時間遅れに基づき、仮想面6で反射した反射光(即ち、時間遅れが2L/Cである反射光)だけを抽出する。N=1回目の照射時であれば、算出装置131は、時間遅れがTである反射光だけを抽出する。算出装置131は、抽出した反射光の強度と、ドップラーシフト量と、仮想面6上の位置情報とを、記憶装置140に記録する。位置情報は、この反射光が反射した仮想面6上の見込み角の格子を示す。強度は、例えば輝度である。なお、記憶装置140の記憶容量に余裕があれば、算出装置131は、仮想面6で反射した反射光である反射光だけに限らず、他の時間遅れの反射光の強度と、ドップラーシフト量と、位置情報とをさらに記憶してもよい。 The calculation device 131 extracts only the reflected light reflected on the virtual surface 6 based on the calculated time delay (i.e., the reflected light with a time delay of 2L/C). If N=1st irradiation, the calculation device 131 extracts only the reflected light with a time delay of T. The calculation device 131 records the intensity of the extracted reflected light, the amount of Doppler shift, and the position information on the virtual surface 6 in the storage device 140. The position information indicates the grid of the angle of view on the virtual surface 6 where the reflected light is reflected. The intensity is, for example, luminance. Note that if the storage device 140 has a sufficient storage capacity, the calculation device 131 may further store the intensity, the amount of Doppler shift, and the position information of reflected light with other time delays, in addition to the reflected light reflected on the virtual surface 6.

なお、N回目のフローで、算出装置131は、T-(N-1)dtの時間遅れの反射光を抽出する。動作フローがループする度に、Nはインクリメント(+1)される。従って、算出装置131は、動作フローがループする度に、より短時間で戻って来た反射光を抽出し続ける。遅れ時間の時間長さは、航空機1から仮想面6までの距離Lに依存して変更される。言い換えれば、距離Lは、航空機1が仮想面6に近づくに従って小さくなるため、遅れ時間の時間長さも短くなる。算出装置131は、より短時間で戻って来た反射光を抽出し続けることで、航空機1の位置の変化に依存せず、大気中の同じ位置にある仮想面6で反射した反射光を抽出し続けることができる。この様に、算出装置131は、航空機1が飛行中に飛行経路上の異なる位置から光を照射する毎に、抽出する反射光の時間遅れを異ならせる。 In addition, in the Nth flow, the calculation device 131 extracts reflected light with a time delay of T-(N-1)dt. N is incremented (+1) each time the operation flow loops. Therefore, the calculation device 131 continues to extract reflected light that returns in a shorter time each time the operation flow loops. The length of the delay time is changed depending on the distance L from the aircraft 1 to the virtual surface 6. In other words, the distance L becomes smaller as the aircraft 1 approaches the virtual surface 6, and therefore the length of the delay time also becomes shorter. By continuing to extract reflected light that returns in a shorter time, the calculation device 131 can continue to extract reflected light reflected by the virtual surface 6 at the same position in the atmosphere, regardless of changes in the position of the aircraft 1. In this way, the calculation device 131 varies the time delay of the extracted reflected light each time the aircraft 1 irradiates light from a different position on the flight path during flight.

ステップS205:反射光マップを作成 Step S205: Create a reflected light map

図4は、マップの一例を示す。 Figure 4 shows an example of a map.

マップは、仮想面6に対応し、複数のグリッドエリアを含む。グリッドエリアは、マップを縦横のグリッド線で格子状に区画し、この区画された各エリアを意味する。マップは、例えば、18×18(左右の分割数N×上下の分割数N)個のグリッドエリアを含む。マップの全体が仮想面6に対応する。1個のグリッドエリアのサイズ(横×縦)は、(Lθ/N)×(Lθ/N)である。マップに含まれるグリッドエリアは、仮想面6の見込み角の矩形と1対1で一致してもよいし、しなくてもよい。仮想面6の見込み角の矩形を、マップに含まれるグリッドエリアに変換可能であれば、マップと仮想面6とは如何なる対応関係でもよい。 The map corresponds to the virtual surface 6 and includes a plurality of grid areas. The grid areas refer to each area obtained by dividing the map into a lattice shape using vertical and horizontal grid lines. The map includes, for example, 18×18 (left and right division number N 1 × top and bottom division number N 2 ) grid areas. The entire map corresponds to the virtual surface 6. The size (horizontal × vertical) of one grid area is (Lθ 1 /N 1 ) × (Lθ 2 /N 2 ). The grid areas included in the map may or may not match one-to-one with the rectangle of the view angle of the virtual surface 6. As long as the rectangle of the view angle of the virtual surface 6 can be converted into a grid area included in the map, any correspondence relationship between the map and the virtual surface 6 is acceptable.

演算装置132は、記憶装置140から、仮想面6で反射した反射光の強度と、位置情報とを読み出す。演算装置132は、仮想面6の位置情報に対応するグリッドエリアに、反射光の強度をプロットする。具体的には、演算装置132は、反射光の強度を量子化し、量子化された強度をグリッドエリアにプロットする。同図の例では、濃色のグリッドエリアは反射光の強度が高いことを示し、淡色のグリッドエリアは反射光の強度が低いことを示し、白色のグリッドエリアは閾値以上の強度の反射光が無いことを示す。これにより、演算装置132は、仮想面6で反射した反射光の強度のパターンを示す反射光マップ3を作成する。 The arithmetic device 132 reads out the intensity of the reflected light reflected by the virtual surface 6 and the position information from the storage device 140. The arithmetic device 132 plots the intensity of the reflected light in a grid area corresponding to the position information of the virtual surface 6. Specifically, the arithmetic device 132 quantizes the intensity of the reflected light and plots the quantized intensity in the grid area. In the example of the same figure, a dark grid area indicates a high intensity of the reflected light, a light grid area indicates a low intensity of the reflected light, and a white grid area indicates that there is no reflected light with an intensity above the threshold. In this way, the arithmetic device 132 creates a reflected light map 3 that indicates the pattern of the intensity of the reflected light reflected by the virtual surface 6.

なお、N回目のフローで、演算装置132は、T-(N-1)dtの時間遅れの反射光の強度の反射光マップ3を作成する。動作フローがループする度に、Nはインクリメント(+1)される。従って、演算装置132は、動作フローがループする度に、より短時間で戻って来た反射光の反射光マップ3を作成し続ける。遅れ時間の時間長さは、航空機1から仮想面6までの距離Lに依存して変更される。言い換えれば、距離Lは、航空機1が仮想面6に近づくに従って小さくなるため、遅れ時間の時間長さも短くなる。演算装置132は、より短時間で戻って来た反射光の反射光マップ3を作成し続けることで、航空機1の位置の変化に依存せず、大気中の同じ位置にある仮想面6に対応する反射光マップ3を作成し続けることができる。この様に、演算装置132は、航空機1が飛行中に飛行経路上の異なる位置から照射された光に対応する複数の反射光マップ3毎に、反射光の時間遅れを異ならせることにより、複数の反射光マップ3の光軸方向の位置を補正する。 In the Nth flow, the calculation device 132 creates a reflected light map 3 of the intensity of reflected light with a time delay of T-(N-1)dt. Each time the operation flow loops, N is incremented (+1). Therefore, the calculation device 132 continues to create a reflected light map 3 of reflected light that returns in a shorter time each time the operation flow loops. The length of the delay time is changed depending on the distance L from the aircraft 1 to the virtual surface 6. In other words, the distance L becomes smaller as the aircraft 1 approaches the virtual surface 6, so the length of the delay time also becomes shorter. By continuing to create a reflected light map 3 of reflected light that returns in a shorter time, the calculation device 132 can continue to create a reflected light map 3 corresponding to the virtual surface 6 at the same position in the atmosphere, regardless of changes in the position of the aircraft 1. In this way, the calculation device 132 corrects the positions in the optical axis direction of the multiple reflected light maps 3 by varying the time delay of the reflected light for each of the multiple reflected light maps 3 corresponding to light irradiated from different positions on the flight path while the aircraft 1 is flying.

なお、ステップS204又はステップS205で、算出装置131又は演算装置132は、レーザ光で仮想面6をスキャンすることにより生じる照射タイミングのずれに応じて時間遅れのずれを補正してもよい。上述のように、照射装置110は、仮想面6の全域を一度に照射するのではなく、仮想面6の全域を順にスキャンすることにより仮想面6の全域にレーザ光を照射する。このため、最初にスキャンする(例えば、左上の頂点を含む)見込み角の矩形エリアから、最後にスキャンする(例えば右下の頂点を含む)見込み角の矩形エリアまで、厳密には、照射タイミングに時間的なずれがある。従って、最初にスキャンするタイミングから最後にスキャンするタイミングまで、航空機1の位置が厳密には異なる。このため、照射タイミングの時間的なずれに依存して、仮想面6が、進行方向(即ち、仮想面6の中心の光軸方向)に直交ではなく、進行方向の成分を含む方向に歪んだ(即ち、進行方向に対して傾斜した)形になり得る。算出装置131又は演算装置132は、この様な照射タイミングのずれに応じて時間遅れのずれを補正することで、仮想面6を進行方向に直交する面となるように補正してもよい。ただし、照射装置110のスキャン速度性能、仮想面6までの距離、算出装置131又は演算装置132の演算性能等に拠り、照射タイミングのずれに応じて時間遅れのずれを補正しなくてもよい。 In step S204 or S205, the calculation device 131 or the arithmetic device 132 may correct the time delay deviation according to the deviation in the irradiation timing caused by scanning the virtual surface 6 with the laser light. As described above, the irradiation device 110 does not irradiate the entire virtual surface 6 at once, but irradiates the entire virtual surface 6 with the laser light by sequentially scanning the entire virtual surface 6. Therefore, strictly speaking, there is a time deviation in the irradiation timing from the rectangular area of the angle of view that is scanned first (e.g., including the upper left vertex) to the rectangular area of the angle of view that is scanned last (e.g., including the lower right vertex). Therefore, the position of the aircraft 1 is strictly different from the timing of the first scan to the timing of the last scan. Therefore, depending on the time deviation in the irradiation timing, the virtual surface 6 may become distorted in a direction including a component of the traveling direction (i.e., inclined with respect to the traveling direction) rather than perpendicular to the traveling direction (i.e., the optical axis direction of the center of the virtual surface 6). The calculation device 131 or the arithmetic device 132 may correct the time delay deviation according to such deviation in irradiation timing, so that the virtual surface 6 becomes a surface perpendicular to the direction of travel. However, depending on the scanning speed performance of the irradiation device 110, the distance to the virtual surface 6, the computation performance of the calculation device 131 or the arithmetic device 132, etc., it is not necessary to correct the time delay deviation according to the deviation in irradiation timing.

ステップS206:反射光マップ同士を比較 Step S206: Compare reflected light maps

演算装置132は、前回作成した反射光マップ3と、dt時間経過後、今回作成した反射光マップ3とを時系列的に比較する。即ち、前回(即ち、N-1回目に)作成した反射光マップ3(時間遅れ:T-Ndt)と、今回(即ち、N回目に)作成した反射光マップ3(時間遅れ:T-(N-1)dt)とを比較する。具体的には、演算装置132は、前回作成した反射光マップ3の反射光の強度のパターンと、dt時間経過後、今回作成した反射光マップ3の反射光の強度のパターンとを比較することにより、パターンの変化を演算する。 The calculation device 132 chronologically compares the reflected light map 3 created last time with the reflected light map 3 created this time after dt time has passed. That is, it compares the reflected light map 3 created last time (i.e., the N-1th time) (time delay: T-Ndt) with the reflected light map 3 created this time (i.e., the Nth time) (time delay: T-(N-1)dt). Specifically, the calculation device 132 calculates the change in the pattern by comparing the pattern of the reflected light intensity of the reflected light map 3 created last time with the pattern of the reflected light intensity of the reflected light map 3 created this time after dt time has passed.

具体的には、演算装置132は、前回の反射光マップ3及び今回の反射光マップ3を、それぞれ、複数の検査区画4に分割する。前回の反射光マップ3及び今回の反射光マップ3を複数の検査区画4に分割するレイアウトは、同じである。 Specifically, the calculation device 132 divides each of the previous reflected light map 3 and the current reflected light map 3 into multiple inspection sections 4. The layout for dividing the previous reflected light map 3 and the current reflected light map 3 into multiple inspection sections 4 is the same.

検査区画4に含まれるグリッドエリアの数(言い換えれば、検査区画4のサイズ)は、固定でも可変でもよい。演算負荷やパイロットの人的負荷を軽くすることを重視するならば固定とすればよいし、演算の正確性を重視するならば可変とすればよい。固定の場合、例えば、検査区画4は、8×8個のグリッドエリアを含む矩形でよい。可変の場合、例えば、検査区画4に含まれる反射光の強度がプロットされたグリッドエリア(即ち、白色で無いグリッドエリア)の数が、例えば、数個以上数十個以下であるように、検査区画4を抽出すればよい。検査区画4に含まれる、反射光の強度がプロットされたグリッドエリアの数は、例えば、5個以上数十個以下とすればよい。反射光の強度がプロットされたグリッドエリアの数を適切に設定することで、複数の検査区画4のパターンの変化をより正確に演算することができる。 The number of grid areas included in the inspection section 4 (in other words, the size of the inspection section 4) may be fixed or variable. If emphasis is placed on reducing the computational load and the pilot's human workload, the number may be fixed. If emphasis is placed on the accuracy of the computation, the number may be variable. In the case of a fixed number, the inspection section 4 may be, for example, a rectangle including 8 x 8 grid areas. In the case of a variable number, the inspection section 4 may be extracted so that the number of grid areas (i.e., non-white grid areas) in which the intensity of reflected light included in the inspection section 4 is plotted is, for example, several or more and several tens or less. The number of grid areas in which the intensity of reflected light included in the inspection section 4 is plotted may be, for example, five or more and several tens or less. By appropriately setting the number of grid areas in which the intensity of reflected light is plotted, the change in the pattern of the multiple inspection sections 4 can be calculated more accurately.

演算装置132は、複数の反射光マップ3(即ち、前回の反射光マップ3及び今回の反射光マップ3)から、仮想面6上の同じエリアに対応する検査区画4をそれぞれ抽出する。演算装置132は、複数の反射光マップ3からそれぞれ抽出した複数の検査区画4のパターンの変化を演算することにより、複数の反射光マップ3のパターンの変化を演算する。 The calculation device 132 extracts inspection sections 4 corresponding to the same area on the virtual surface 6 from each of the multiple reflected light maps 3 (i.e., the previous reflected light map 3 and the current reflected light map 3). The calculation device 132 calculates the change in the pattern of the multiple reflected light maps 3 by calculating the change in the pattern of the multiple inspection sections 4 extracted from each of the multiple reflected light maps 3.

図5は、複数の反射光マップからそれぞれ抽出した複数の検査区画のパターンの変化を模式的に示す。 Figure 5 shows a schematic diagram of the changes in the patterns of multiple inspection sections extracted from multiple reflected light maps.

左図は、前回の反射光マップ3から分割された1個の検査区画4を示す(図4の検査区画4と同じ)。右図は、dt時間経過後の反射光マップ3から分割された1個の検査区画4'を示す。検査区画4及び検査区画4'の反射光マップ3内の位置(座標)は、同じである。 The left image shows one inspection section 4 divided from the previous reflected light map 3 (same as inspection section 4 in Figure 4). The right image shows one inspection section 4' divided from the reflected light map 3 after the passage of time dt. The positions (coordinates) of inspection section 4 and inspection section 4' in the reflected light map 3 are the same.

ステップS207:仮想面上の気流の方向及び速度を計測 Step S207: Measure the direction and speed of airflow on the virtual surface

演算装置132は、前回の検査区画4の反射光のパターンから今回の検査区画4'の反射光のパターンへの変化に基づき、仮想面6上の気流の方向及び速度を計測する。演算装置132が計測する仮想面6上の気流の方向及び速度は、仮想面6の面方向の気流の方向及び速度である。上述のように、仮想面6は、照射される全ての光の光軸に交差し、中でも、仮想面6上の1点(例えば、中心点)に照射される光の光軸に直交する。また、例えば、航空機1が通過する予定である飛行経路が、仮想面6上の中心点又は中心点付近に直交する。この場合、演算装置132が計測する気流の方向及び速度は、予定される飛行経路に直交する方向(即ち、上下方向及び左右方向の成分を含む方向)の気流の速度である。 The arithmetic device 132 measures the direction and speed of the airflow on the virtual surface 6 based on the change from the pattern of reflected light of the previous inspection section 4 to the pattern of reflected light of the current inspection section 4'. The direction and speed of the airflow on the virtual surface 6 measured by the arithmetic device 132 are the direction and speed of the airflow in the surface direction of the virtual surface 6. As described above, the virtual surface 6 intersects with the optical axes of all the irradiated light, and is perpendicular to the optical axis of the light irradiated to one point (e.g., the center point) on the virtual surface 6. Also, for example, the flight path that the aircraft 1 is scheduled to pass through is perpendicular to the center point or near the center point on the virtual surface 6. In this case, the direction and speed of the airflow measured by the arithmetic device 132 are the speed of the airflow in a direction perpendicular to the planned flight path (i.e., a direction including vertical and horizontal components).

具体的には、演算装置132は、前回の検査区画4の反射光のパターン及び今回の検査区画4'の反射光のパターンから相互相関関数を算出し、その最大値を与える移動ベクトルを算出し、速度を算出する。図5の例の場合、前回の検査区画4の反射光のパターンを右方向に1個及び下方向に2個移動して、dt時間経過後の検査区画4'の反射光のパターンとの相互相関を取ると、最大値が得られる。よって、演算装置132は、時間dtの間に、仮想面6上の右方向にLθ/N、且つ、仮想面6上の下方向に2Lθ/Nだけ大気が移動したと推定する。よって、演算装置132は、気流の速度が、仮想面6上の右方向にLθ/N/dt、且つ、仮想面6上の下方向に2Lθ/N/dtであると計測する。これにより、演算装置132は、気流の流れの方向及び速度を計測することができる。 Specifically, the arithmetic device 132 calculates a cross-correlation function from the pattern of reflected light of the previous inspection section 4 and the pattern of reflected light of the current inspection section 4', calculates a movement vector that gives the maximum value, and calculates the speed. In the example of FIG. 5, the pattern of reflected light of the previous inspection section 4 is moved one to the right and two to the downward direction, and the cross-correlation with the pattern of reflected light of the inspection section 4' after dt time has elapsed is obtained to obtain the maximum value. Therefore, the arithmetic device 132 estimates that the air has moved by Lθ 1 /N 1 to the right on the virtual surface 6 and by 2Lθ 2 /N 2 to the downward direction on the virtual surface 6 during the time dt. Therefore, the arithmetic device 132 measures the speed of the airflow to be Lθ 1 /N 1 /dt to the right on the virtual surface 6 and 2Lθ 2 /N 2 /dt to the downward direction on the virtual surface 6. This allows the arithmetic device 132 to measure the direction and speed of the airflow.

ステップS208:進行方向の気流の速度を計測 Step S208: Measure the airflow speed in the direction of travel

演算装置132は、記憶装置140から、仮想面6で反射した反射光のドップラーシフト量と、位置情報とを読み出す。演算装置132は、ドップラーシフト量に基づき、仮想面6上の光軸の方向(即ち、航空機1の進行方向)の気流の速度を計測する。本実施形態によれば、仮想面6の面方向の気流の速度を計測するのに用いられる反射光を利用して、航空機1の進行方向の気流の速度も計測することができる。 The arithmetic device 132 reads the Doppler shift amount of the reflected light reflected by the virtual surface 6 and the position information from the storage device 140. The arithmetic device 132 measures the airflow speed in the direction of the optical axis on the virtual surface 6 (i.e., the traveling direction of the aircraft 1) based on the Doppler shift amount. According to this embodiment, the reflected light used to measure the airflow speed in the surface direction of the virtual surface 6 can also be used to measure the airflow speed in the traveling direction of the aircraft 1.

演算装置132は、速度表示画面151を、表示装置150(例えば、コックピット内のディスプレイ)に表示する。速度表示画面151は、計測した仮想面6の面方向の気流の速度と、進行方向の気流の速度を表示するGUIである。 The computing device 132 displays a speed display screen 151 on the display device 150 (e.g., a display in the cockpit). The speed display screen 151 is a GUI that displays the measured airflow speed in the surface direction of the virtual surface 6 and the airflow speed in the traveling direction.

ステップS209:ループ Step S209: Loop

演算装置132は、スキャン回数Nが閾値を超えていなければ(ステップS209、NO)、スキャン回数Nをインクリメント(+1)し、照射装置110にN+1回目のレーザ光を照射させる(ステップS202)。一方、演算装置132は、スキャン回数Nが閾値を超えれば(ステップS209、YES)、新たに、N=1と、基準となる遅れ時間Tとを設定してから(ステップS201)、照射装置110にN回目のレーザ光を照射させる(ステップS202)。 If the number of scans N does not exceed the threshold (step S209, NO), the calculation device 132 increments the number of scans N by 1 and causes the irradiation device 110 to irradiate the laser light for the N+1th time (step S202). On the other hand, if the number of scans N exceeds the threshold (step S209, YES), the calculation device 132 sets a new N=1 and a reference delay time T (step S201), and then causes the irradiation device 110 to irradiate the laser light for the Nth time (step S202).

上述のように、仮想面6上に乱気流が存在する場合に、航空機1は乱気流を回避して仮想面6内の安全な(即ち、乱気流が存在しない)空域を飛行することが望まれる。従って、航空機1から仮想面6までの距離として、航空機1が仮想面6上の乱気流を回避できる程度の距離が必要である。このため、スキャン回数Nの閾値(上限値)を設定する。スキャン回数Nが閾値を超えると、演算装置132は、基準となる遅れ時間Tを新たに設定する。遅れ時間Tに依存して仮想面6までの距離Lが新たに定義される。このため、基準となる遅れ時間Tを新たに設定することで、新たな仮想面6上の乱気流を検出及び計測することができる。これにより、気流計測装置100は、常に、所定距離(即ち、乱気流を回避できる程度の距離)だけ離れた仮想面6上の乱気流を検出及び計測し続けることができる。 As described above, when turbulence exists on the virtual plane 6, it is desirable for the aircraft 1 to avoid the turbulence and fly in a safe airspace (i.e., no turbulence exists) within the virtual plane 6. Therefore, the distance from the aircraft 1 to the virtual plane 6 must be a distance that allows the aircraft 1 to avoid the turbulence on the virtual plane 6. For this reason, a threshold (upper limit) for the number of scans N is set. When the number of scans N exceeds the threshold, the computing device 132 sets a new reference delay time T. The distance L to the virtual plane 6 is newly defined depending on the delay time T. For this reason, by setting a new reference delay time T, it is possible to detect and measure turbulence on a new virtual plane 6. As a result, the airflow measuring device 100 can always continue to detect and measure turbulence on the virtual plane 6 that is a predetermined distance away (i.e., a distance that allows turbulence to be avoided).

4.粒子画像流速測定法(PIV)との相違点 4. Differences from particle image velocimetry (PIV)

流体の流れ方向及び速度を計測する方法として、粒子画像流速測定法(PIV:Particle Image Velocimetry)が知られている。仮想面6上の気流の方向及び速度を計測する方法(ステップS207)は、PIVとの相違点も多く、また、PIVに対して有利な点もある。それを説明する。 Particle Image Velocimetry (PIV) is known as a method for measuring the flow direction and speed of a fluid. The method for measuring the direction and speed of the airflow on the virtual surface 6 (step S207) has many differences from PIV, but also has advantages over PIV. These differences will be explained below.

PIVでは、アルミ微粉末等のトレーサ粒子を、流れ場(例えば、実験室内の風洞や回流水槽)にシーディングする。上方から流れ場に向けて、レーザライトシート(即ち、シート状のレーザ光)を照射する。このときライトシートの面方向が、流体の流れの方向(即ち、速度を検出すべき流れの方向)と一致するように、レーザライトシートを照射する。レーザライトシートに直交する方向から、CCD又はCMOSカメラを用いて、レーザライトシートに照射されたトレーサ粒子を経時的に複数回撮像する。経時的に撮像された複数の粒子画像に含まれるトレーサ粒子のパターンを解析し、トレーサ粒子の移動速度、即ち、流体の速度を計測する。これがPIVの概要である。 In PIV, tracer particles such as aluminum powder are seeded into a flow field (for example, a wind tunnel or circulating water tank in a laboratory). A laser light sheet (i.e., sheet-shaped laser light) is irradiated from above toward the flow field. The laser light sheet is irradiated so that the surface direction of the light sheet coincides with the direction of the fluid flow (i.e., the direction of the flow whose speed is to be detected). A CCD or CMOS camera is used to capture images of the tracer particles irradiated onto the laser light sheet multiple times over time from a direction perpendicular to the laser light sheet. The patterns of the tracer particles contained in the multiple particle images captured over time are analyzed to measure the movement speed of the tracer particles, i.e., the speed of the fluid. This is an overview of PIV.

PIVは、トレーサ粒子を流れ場にシーディングし、流れ場に向けて上方からレーザライトシートを照射する。これに対して、航空機1が遭遇する可能性のある気流を計測する場合は、航空機1から離れた空域にトレーサ粒子をシーディングし、空域の上方からレーザライトシートを照射することは、不可能である。これが本実施形態の前提条件である。 PIV seeds tracer particles into a flow field and irradiates the flow field with a laser light sheet from above. In contrast, when measuring air currents that aircraft 1 may encounter, it is not possible to seed tracer particles into an airspace away from aircraft 1 and irradiate the airspace with a laser light sheet from above. This is a prerequisite for this embodiment.

PIVは、画像解析に適した良好な粒子画像を撮像する必要がある。このため、トレーサ粒子の粒径や、流れ場中の粒子密度を適切にコントロールする必要がある。これに対して、航空機1が遭遇する可能性のある気流を計測する場合は、トレーサ粒子の代わりとなる大気中のエアロゾル(塵、微粒子)の粒径や粒子密度をコントロールすることは出来ない。また、航空機1から速度計測対象である大気までの距離が、PIVでのカメラからライトシートまでの距離に比べて、遥かに大きい。このため、航空機1から離れた空域の大気中の塵等の画像を、画像解析に適したレベルで、航空機1から撮像するのは極めて困難である。 PIV needs to capture good particle images suitable for image analysis. For this reason, it is necessary to appropriately control the particle size of the tracer particles and the particle density in the flow field. In contrast, when measuring airflows that aircraft 1 may encounter, it is not possible to control the particle size or particle density of atmospheric aerosols (dust, fine particles) that serve as substitutes for tracer particles. In addition, the distance from aircraft 1 to the atmosphere, whose speed is to be measured, is much greater than the distance from the camera to the light sheet in PIV. For this reason, it is extremely difficult to capture images of dust and other particles in the atmosphere in airspace away from aircraft 1 at a level suitable for image analysis from aircraft 1.

以上のような事情に鑑み、本実施形態によれば、気流計測装置100は、照射から受光までの反射光の時間遅れを算出し、時間遅れに基づき、予定される飛行経路に直交する仮想面6で反射した反射光を抽出し、仮想面6で反射した反射光の強度のパターンを示す反射光マップ3を作成する。これにより、速度を検出すべき流れの方向と一致する面方向を有する仮想面6で反射した反射光の反射光マップ3を作成することができる。言い換えれば、PIVの粒子画像と同様に流れの方向及び速度を解析することが可能な反射光マップ3を作成することができる。PIVでは、粒子密度に基づきトレーサ粒子の分布を解析するが、本実施形態では、粒子密度の代わりに反射光の強度に基づき塵の分布を解析することができる。 In view of the above circumstances, according to this embodiment, the airflow measurement device 100 calculates the time delay of the reflected light from irradiation to reception, extracts the reflected light reflected from the virtual surface 6 perpendicular to the planned flight path based on the time delay, and creates a reflected light map 3 showing the pattern of the intensity of the reflected light reflected from the virtual surface 6. This makes it possible to create a reflected light map 3 of the reflected light reflected from the virtual surface 6 having a surface direction that matches the direction of the flow whose speed is to be detected. In other words, it is possible to create a reflected light map 3 that can analyze the flow direction and speed in the same way as a particle image of PIV. In PIV, the distribution of tracer particles is analyzed based on particle density, but in this embodiment, the distribution of dust can be analyzed based on the intensity of the reflected light instead of particle density.

PIVは、レーザライトシートを照射し、レーザライトシートに照らされた(即ち、レーザライトシート内に分布する)トレーサ粒子だけを撮像する。レーザライトシートは、必ず物理的な厚みを持つ。このため、レーザライトシートの面方向だけでなく、レーザライトシートの厚み方向に分布するトレーサ粒子も撮像されることになる。言い換えれば、レーザライトシートの面方向だけでなくレーザライトシートの厚み方向に分布するトレーサ粒子が、粒子画像に含まれる。このため、PIVでは、レーザライトシートの面方向だけの気流速度を正確に計測できないおそれがある。 PIV shines a laser light sheet and captures only the tracer particles illuminated by the laser light sheet (i.e., distributed within the laser light sheet). Laser light sheets always have a physical thickness. For this reason, tracer particles distributed not only in the surface direction of the laser light sheet but also in the thickness direction of the laser light sheet are captured. In other words, tracer particles distributed not only in the surface direction of the laser light sheet but also in the thickness direction of the laser light sheet are included in the particle image. For this reason, PIV may not be able to accurately measure the airflow speed only in the surface direction of the laser light sheet.

これに対して、本実施形態の気流計測装置100は、特定の遅れ時間で受光する反射光だけを抽出することにより、物理的な厚みの無い仮想面6で反射した反射光だけを抽出することが可能である。このため、本実施形態の気流計測装置100は、仮想面6の面方向だけの気流速度を正確に計測することを図れる。 In contrast, the airflow measurement device 100 of this embodiment extracts only the reflected light received with a specific delay time, making it possible to extract only the reflected light reflected by the virtual surface 6, which has no physical thickness. As a result, the airflow measurement device 100 of this embodiment can accurately measure the airflow speed only in the surface direction of the virtual surface 6.

ところで、エアロゾルの大きさが0である場合、仮想面6上にエアロゾルが存在する確率は0である。エアロゾルに大きさ(代表長さd)があればdに比例して存在確率は大きくなる。時間遅れTに時間幅dt'を持たせれば、仮想面6の面積Sにdt'Cの厚さを加えた直方体の中にあるエアロゾルを検出することになり、エアロゾルの大きさdが0だとしても、存在確率は有限になる。高い高度に存在するエアロゾルの密度は非常に低いため、エアロゾルの検出のための時間幅dt'を設けてもよい。 When the size of the aerosol is 0, the probability that the aerosol exists on the virtual surface 6 is 0. If the aerosol has a size (characteristic length d), the probability of its existence increases in proportion to d. If the time delay T is given a time width dt', the aerosol will be detected within a rectangular parallelepiped whose area S of the virtual surface 6 is added to a thickness of dt'C, and even if the size d of the aerosol is 0, the probability of its existence will be finite. Since the density of aerosols at high altitudes is very low, a time width dt' can be set for detecting the aerosol.

5.結語 5. Conclusion

本実施形態によれば、気流計測装置100は、照射から受光までの反射光の時間遅れを算出し、時間遅れに基づき、予定される飛行経路に直交する仮想面6で反射した反射光を抽出し、仮想面6で反射した反射光の強度のパターンを示す反射光マップ3を作成する。これにより、速度を検出すべき流れの方向と一致する面方向を有する仮想面6で反射した反射光の反射光マップ3を作成することができる。時間的に異なる複数の反射光マップ3を作成することで、仮想面6の面方向の気流の速度を計測することができる。 According to this embodiment, the airflow measurement device 100 calculates the time delay of the reflected light from irradiation to reception, extracts the reflected light reflected from a virtual surface 6 perpendicular to the planned flight path based on the time delay, and creates a reflected light map 3 that shows the intensity pattern of the reflected light reflected from the virtual surface 6. This makes it possible to create a reflected light map 3 of the reflected light reflected from a virtual surface 6 having a surface direction that matches the direction of the flow whose speed is to be detected. By creating multiple reflected light maps 3 that are different in time, it is possible to measure the speed of the airflow in the surface direction of the virtual surface 6.

仮想面6の面方向を予定される飛行経路に直交する方向とすることで、予定される飛行経路に直交する方向(即ち、上下方向及び左右方向の成分を含む方向)の気流の速度を計測することができる。特に、晴天乱気流の上下方向の速度を計測することができる。また、仮想面6のサイズを適切に設定することで、航空機1が乱気流を回避するために必要な空域の気流速度を精度よく適切に算出することを図れる。 By orienting the surface direction of the virtual surface 6 perpendicular to the planned flight path, it is possible to measure the air current speed in a direction perpendicular to the planned flight path (i.e., a direction including vertical and horizontal components). In particular, it is possible to measure the vertical speed of clear air turbulence. Furthermore, by appropriately setting the size of the virtual surface 6, it is possible to accurately and appropriately calculate the air current speed in the airspace required for the aircraft 1 to avoid turbulence.

また、本実施形態では、航空機1の位置の変化に依存せず、大気中の同じ位置にある仮想面6で反射した反射光の反射光マップ3を作成し続けた(ステップS204、ステップS205)。変形例として、航空機1の位置の変化に依存して、進行方向に対して異なる位置の仮想面6で反射した反射光の反射光マップ3を作成すれば、進行方向の成分を含む方向の気流の速度を計測することも可能である。 In addition, in this embodiment, the reflected light map 3 of the light reflected by the virtual surface 6 at the same position in the atmosphere is continuously created, regardless of the change in the position of the aircraft 1 (steps S204 and S205). As a modified example, if the reflected light map 3 of the light reflected by the virtual surface 6 at a different position relative to the traveling direction is created depending on the change in the position of the aircraft 1, it is also possible to measure the speed of the air current in a direction that includes a component in the traveling direction.

本実施形態では、最小限のハードウェアであるレーザ発振装置111及び受光装置120を航空機1に搭載するだけでよく、処理装置130は受光装置120が受光した反射光に基づき、ソフトウェア処理により仮想面6の面方向の気流の速度を計測する。また、同時にドップラーシフト量に基づく進行方向の気流速度を検出することも可能である。 In this embodiment, it is only necessary to mount the minimum amount of hardware, the laser oscillator 111 and the light receiving device 120, on the aircraft 1, and the processing device 130 measures the airflow speed in the surface direction of the virtual surface 6 through software processing based on the reflected light received by the light receiving device 120. It is also possible to simultaneously detect the airflow speed in the direction of travel based on the amount of Doppler shift.

気流計測装置100は、航空機1に搭載するだけではなく、地上装置にも適用可能である。また、航空機1に搭載された気流計測装置100が計測した気流を、気流計測装置100が地上装置や他の航空機に送信してもよい。また、気流計測装置100は、船舶や車両に搭載することも可能である。 The airflow measuring device 100 can be applied not only to the aircraft 1 but also to ground equipment. The airflow measuring device 100 mounted on the aircraft 1 may transmit the airflow measured by the airflow measuring device 100 to ground equipment or other aircraft. The airflow measuring device 100 can also be mounted on ships and vehicles.

本技術の各実施形態及び各変形例について上に説明したが、本技術は上述の実施形態にのみ限定されるものではなく、本技術の要旨を逸脱しない範囲内において種々変更を加え得ることは勿論である。 Although the embodiments and variations of the present technology have been described above, the present technology is not limited to the above-described embodiments, and various modifications can of course be made without departing from the spirit of the present technology.

1 航空機
2 照射及び受光装置
3 反射光マップ
4 検査区画
5 グリッドエリア
6 仮想面
100 気流計測装置
110 照射装置
111 レーザ発振装置
112 スキャン装置
120 受光装置
130 処理装置
131 算出装置
132 演算装置
140 記憶装置
150 表示装置
151 速度表示画面
REFERENCE SIGNS LIST 1 Aircraft 2 Illumination and light receiving device 3 Reflected light map 4 Inspection section 5 Grid area 6 Virtual surface 100 Airflow measurement device 110 Illumination device 111 Laser oscillation device 112 Scanning device 120 Light receiving device 130 Processing device 131 Calculation device 132 Arithmetic device 140 Storage device 150 Display device 151 Speed display screen

Claims (15)

照射装置が光を前記光の光軸に交差する大気中の仮想面に照射してから、受光装置が前記光の反射光を受光するまでの時間遅れを算出し、
前記時間遅れに基づき、前記仮想面で反射した反射光を抽出する
算出装置と、
前記仮想面に対応し複数のグリッドエリアを含むマップの各グリッドエリアに、前記反射光の強度をプロットすることで、反射光の強度のパターンを示す反射光マップを作成し、
前記反射光マップと、所定時間経過後に作成された反射光マップとの前記パターンの変化を演算し、
前記パターンの変化に基づき、前記仮想面上の前記光軸に交差する方向の気流の速度を計測する
演算装置と
を具備する気流計測装置。
Calculating a time delay from when an irradiating device irradiates light onto a virtual surface in the atmosphere that intersects with an optical axis of the light until when a light receiving device receives reflected light of the light;
A calculation device that extracts the reflected light reflected by the virtual surface based on the time delay;
creating a reflected light map showing a pattern of the intensity of the reflected light by plotting the intensity of the reflected light in each grid area of a map corresponding to the virtual surface and including a plurality of grid areas;
Calculating a change in the pattern between the reflected light map and a reflected light map created after a predetermined time has elapsed;
and a calculation device that measures the speed of the airflow on the virtual surface in a direction intersecting the optical axis based on the change in the pattern.
請求項1に記載の気流計測装置であって、
前記複数の反射光マップに対応する前記仮想面は、大気中の同じ位置にあり、
前記演算装置は、前記光軸に交差する方向として、前記仮想面の面方向の気流の速度を計測する
気流計測装置。
The airflow measurement device according to claim 1,
the virtual surfaces corresponding to the plurality of reflected light maps are at the same position in the atmosphere;
The airflow measuring device, wherein the calculation device measures the speed of the airflow in a planar direction of the virtual surface as a direction intersecting the optical axis.
請求項1又は2に記載の気流計測装置であって、
前記仮想面は、前記仮想面上の1点に照射される光の光軸に直交する
気流計測装置。
The airflow measuring device according to claim 1 or 2,
The virtual plane is perpendicular to an optical axis of light irradiated to one point on the virtual plane.
請求項1乃至3の何れか一項に記載の気流計測装置であって、
前記気流計測装置は、航空機に搭載される
気流計測装置。
The airflow measuring device according to any one of claims 1 to 3,
The airflow measuring device is an airflow measuring device mounted on an aircraft.
請求項4に記載の気流計測装置であって、
前記仮想面は、飛行中である前記航空機が通過する予定である飛行経路を含む
気流計測装置。
The airflow measuring device according to claim 4,
The virtual surface includes a flight path along which the aircraft is scheduled to fly.
請求項4又は5に記載の気流計測装置であって、
前記仮想面は、前記仮想面上に乱気流が存在する場合に前記航空機が前記乱気流を回避するために飛行する飛行経路を含む
気流計測装置。
The airflow measuring device according to claim 4 or 5,
An airflow measurement device, wherein the virtual plane includes a flight path along which the aircraft will fly to avoid turbulence when turbulence exists on the virtual plane.
請求項4乃至6の何れか一項に記載の気流計測装置であって、
前記算出装置は、前記航空機が飛行中に飛行経路上の異なる位置から光を照射する毎に、抽出する反射光の時間遅れを異ならせる
気流計測装置。
The airflow measuring device according to any one of claims 4 to 6,
The calculation device varies the time delay of the extracted reflected light each time the light is irradiated from a different position on a flight path while the aircraft is flying.
請求項1乃至7の何れか一項に記載の気流計測装置であって、
前記演算装置は、
前記複数の反射光マップから、前記仮想面上の同じエリアに対応する検査区画をそれぞれ抽出し、
前記複数の反射光マップからそれぞれ抽出した複数の前記検査区画の前記パターンの変化を演算することにより、
前記複数の反射光マップの前記パターンの変化を演算する
気流計測装置。
The airflow measuring device according to any one of claims 1 to 7,
The computing device includes:
Extracting inspection sections from the plurality of reflected light maps, each of which corresponds to a same area on the virtual surface;
Calculating changes in the patterns of the plurality of inspection sections extracted from the plurality of reflected light maps, respectively;
An airflow measurement device that calculates changes in the patterns of the plurality of reflected light maps.
請求項8に記載の気流計測装置であって、
前記検査区画に含まれる前記グリッドエリアの数は可変であり、
前記検査区画に含まれる、反射光の強度がプロットされたグリッドエリアの数は、数個以上数十個以下である
気流計測装置。
The airflow measurement device according to claim 8,
The number of the grid areas included in the inspection section is variable;
The number of grid areas in which the intensity of reflected light is plotted and included in the inspection section is between several and several tens.
請求項1乃至9の何れか一項に記載の気流計測装置であって、
前記照射装置は、
前記仮想面を分割したエリアと同形の断面形状を有するレーザ光を発振するレーザ発振装置と、
前記レーザ光の照射方向を制御して、前記レーザ光で前記仮想面の全てのエリアをスキャンするスキャン装置と
を有する
気流計測装置。
The airflow measuring device according to any one of claims 1 to 9,
The irradiation device includes:
a laser oscillation device that oscillates a laser beam having a cross-sectional shape identical to that of an area obtained by dividing the imaginary surface;
and a scanning device that controls an irradiation direction of the laser light and scans all areas of the virtual surface with the laser light.
請求項10に記載の気流計測装置であって、
前記算出装置又は前記演算装置は、前記レーザ光で前記仮想面をスキャンすることにより生じる照射タイミングのずれに応じて時間遅れのずれを補正する
気流計測装置。
The airflow measurement device according to claim 10 ,
The calculation device or the arithmetic device corrects a time delay deviation in accordance with a deviation in irradiation timing caused by scanning the virtual surface with the laser light.
請求項1乃至11の何れか一項に記載の気流計測装置であって、
前記算出装置は、さらに、照射された光と、前記仮想面で反射した反射光との間の周波数のドップラーシフト量を算出し、
前記演算装置は、さらに、前記ドップラーシフト量に基づき、前記仮想面上の前記光軸の方向の気流の速度を計測する
気流計測装置。
The airflow measuring device according to any one of claims 1 to 11,
The calculation device further calculates a Doppler shift amount of a frequency between the irradiated light and the reflected light reflected by the virtual surface,
The calculation device further measures a speed of the airflow in the direction of the optical axis on the virtual plane based on the amount of Doppler shift.
請求項1乃至12の何れか一項に記載の気流計測装置であって、
前記照射装置と、
前記受光装置と
をさらに具備する気流計測装置。
The airflow measuring device according to any one of claims 1 to 12,
The irradiation device;
The airflow measuring device further comprises the light receiving device.
照射装置が光を前記光の光軸に交差する大気中の仮想面に照射してから、受光装置が前記光の反射光を受光するまでの時間遅れに基づき、前記仮想面で反射した反射光を抽出し、
前記仮想面に対応し複数のグリッドエリアを含むマップの各グリッドエリアに、前記反射光の強度をプロットすることで、反射光の強度のパターンを示す反射光マップを作成し、
前記反射光マップと、所定時間経過後に作成された反射光マップとの前記パターンの変化を演算し、
前記パターンの変化に基づき、前記仮想面上の前記光軸に交差する方向の気流の速度を計測する
気流計測方法。
extracting light reflected by the virtual surface based on a time delay between when an irradiating device irradiates light onto a virtual surface in the atmosphere that intersects with an optical axis of the light and when a light receiving device receives reflected light of the light;
creating a reflected light map showing a pattern of the intensity of the reflected light by plotting the intensity of the reflected light in each grid area of a map corresponding to the virtual surface and including a plurality of grid areas;
Calculating a change in the pattern between the reflected light map and a reflected light map created after a predetermined time has elapsed;
based on the change in the pattern, measuring a speed of the airflow on the virtual surface in a direction intersecting the optical axis.
気流計測装置の演算装置を、
照射装置が光を前記光の光軸に交差する大気中の仮想面に照射してから、受光装置が前記光の反射光を受光するまでの時間遅れに基づき、前記仮想面で反射した反射光を抽出し、
前記仮想面に対応し複数のグリッドエリアを含むマップの各グリッドエリアに、前記反射光の強度をプロットすることで、反射光の強度のパターンを示す反射光マップを作成し、
前記反射光マップと、所定時間経過後に作成された反射光マップとの前記パターンの変化を演算し、
前記パターンの変化に基づき、前記仮想面上の前記光軸に交差する方向の気流の速度を計測する
よう動作させる気流計測プログラム。
A calculation unit of the airflow measuring device,
extracting light reflected by the virtual surface based on a time delay between when an irradiating device irradiates light onto a virtual surface in the atmosphere that intersects with an optical axis of the light and when a light receiving device receives reflected light of the light;
creating a reflected light map showing a pattern of the intensity of the reflected light by plotting the intensity of the reflected light in each grid area of a map corresponding to the virtual surface and including a plurality of grid areas;
Calculating a change in the pattern between the reflected light map and a reflected light map created after a predetermined time has elapsed;
an airflow measurement program that operates to measure the speed of airflow on the virtual surface in a direction intersecting the optical axis based on the change in the pattern.
JP2021043041A 2021-03-17 2021-03-17 Airflow measurement device, airflow measurement method, and airflow measurement program Active JP7597367B2 (en)

Priority Applications (1)

Application Number Priority Date Filing Date Title
JP2021043041A JP7597367B2 (en) 2021-03-17 2021-03-17 Airflow measurement device, airflow measurement method, and airflow measurement program

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
JP2021043041A JP7597367B2 (en) 2021-03-17 2021-03-17 Airflow measurement device, airflow measurement method, and airflow measurement program

Publications (2)

Publication Number Publication Date
JP2022142821A JP2022142821A (en) 2022-10-03
JP7597367B2 true JP7597367B2 (en) 2024-12-10

Family

ID=83454093

Family Applications (1)

Application Number Title Priority Date Filing Date
JP2021043041A Active JP7597367B2 (en) 2021-03-17 2021-03-17 Airflow measurement device, airflow measurement method, and airflow measurement program

Country Status (1)

Country Link
JP (1) JP7597367B2 (en)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002236174A (en) 2001-02-07 2002-08-23 Mitsubishi Electric Corp Fog observation device and fog observation method
JP2002286840A (en) 2001-03-23 2002-10-03 Mitsubishi Electric Corp Wake turbulence detection device and wake turbulence detection method
JP2009501915A (en) 2005-07-21 2009-01-22 エーアブス・ドイチュラント・ゲーエムベーハー Turbulence measurements on aircraft and methods and lidar systems for airports and wind farms
JP2010241412A (en) 2009-03-31 2010-10-28 Honeywell Internatl Inc System and method for assessing weather approaching an aircraft in flight
WO2017208375A1 (en) 2016-05-31 2017-12-07 淳一 古本 Doppler shift analysis device

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
JP2002236174A (en) 2001-02-07 2002-08-23 Mitsubishi Electric Corp Fog observation device and fog observation method
JP2002286840A (en) 2001-03-23 2002-10-03 Mitsubishi Electric Corp Wake turbulence detection device and wake turbulence detection method
JP2009501915A (en) 2005-07-21 2009-01-22 エーアブス・ドイチュラント・ゲーエムベーハー Turbulence measurements on aircraft and methods and lidar systems for airports and wind farms
JP2010241412A (en) 2009-03-31 2010-10-28 Honeywell Internatl Inc System and method for assessing weather approaching an aircraft in flight
WO2017208375A1 (en) 2016-05-31 2017-12-07 淳一 古本 Doppler shift analysis device

Also Published As

Publication number Publication date
JP2022142821A (en) 2022-10-03

Similar Documents

Publication Publication Date Title
US8554395B2 (en) Method and system for facilitating autonomous landing of aerial vehicles on a surface
US10043404B2 (en) Method and system for aircraft taxi strike alerting
Sabatini et al. LIDAR obstacle warning and avoidance system for unmanned aircraft
US7683928B2 (en) Lidar with streak-tube imaging, including hazard detection in marine applications; related optics
EP3372508B1 (en) Method and system for aircraft taxi strike alerting
US8666570B1 (en) Volcanic ash detection by optical backscatter using standard aircraft lights
Hann UAV icing: Ice accretion experiments and validation
CN108387885A (en) Inversion Method of Aircraft Wake Characteristic Parameters in Clear Sky Conditions Based on Lidar Detection
US20180321273A1 (en) Particle image velocimetry and control method therefor
Chan Application of LIDAR-based F-factor in windshear alerting
CN110162927A (en) Takeoff and landing method for early warning based on flight simulation platform and anemometry laser radar
CN106840598B (en) Condition of raining based on side looking radar is got off the plane wake flow circular rector estimation method
Gong et al. Laser scanning applied for ice shape measurements
EP3372509B1 (en) Method and system for aircraft taxi strike alerting
Shen et al. Aircraft wake recognition and strength classification based on deep learning
CN111679285A (en) Optical detection method and device for aircraft wake vortex
JP7597367B2 (en) Airflow measurement device, airflow measurement method, and airflow measurement program
Harris et al. Aircraft wake vortices: a comparison of wind-tunnel data with field trial measurements by laser radar
JP5511196B2 (en) Back turbulence detector
US10723466B2 (en) Device and method for determining the water content in the atmosphere, device and method for detecting icing conditions and computer programs
Douxchamps et al. On-board axial detection of wake vortices using a 2 µm LiDAR
Yoshikawa et al. Wake vortex observation campaign by ultra fast-scanning lidar in narita airport, japan
CN117630970A (en) An intelligent wind shear monitoring and early warning method based on shipborne lidar
Mutuel et al. ATM decision support tool for wake vortex hazard management combining sensors and modeling
US11519997B1 (en) Broad-area laser awareness sensor

Legal Events

Date Code Title Description
A621 Written request for application examination

Free format text: JAPANESE INTERMEDIATE CODE: A621

Effective date: 20231212

RD02 Notification of acceptance of power of attorney

Free format text: JAPANESE INTERMEDIATE CODE: A7422

Effective date: 20231212

A977 Report on retrieval

Free format text: JAPANESE INTERMEDIATE CODE: A971007

Effective date: 20240926

A131 Notification of reasons for refusal

Free format text: JAPANESE INTERMEDIATE CODE: A131

Effective date: 20241001

A521 Request for written amendment filed

Free format text: JAPANESE INTERMEDIATE CODE: A523

Effective date: 20241023

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

A61 First payment of annual fees (during grant procedure)

Free format text: JAPANESE INTERMEDIATE CODE: A61

Effective date: 20241121

R150 Certificate of patent or registration of utility model

Ref document number: 7597367

Country of ref document: JP

Free format text: JAPANESE INTERMEDIATE CODE: R150