EP1613980A1 - Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system - Google Patents

Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system

Info

Publication number
EP1613980A1
EP1613980A1 EP04725596A EP04725596A EP1613980A1 EP 1613980 A1 EP1613980 A1 EP 1613980A1 EP 04725596 A EP04725596 A EP 04725596A EP 04725596 A EP04725596 A EP 04725596A EP 1613980 A1 EP1613980 A1 EP 1613980A1
Authority
EP
European Patent Office
Prior art keywords
doppler
interferogram
lidar system
laser beam
interferometer
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.)
Granted
Application number
EP04725596A
Other languages
German (de)
French (fr)
Other versions
EP1613980B1 (en
Inventor
Thorsteinn Halldorsson
Michael Hofmann
Gerald Sobotta
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.)
Airbus Defence and Space GmbH
Original Assignee
EADS Deutschland GmbH
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 EADS Deutschland GmbH filed Critical EADS Deutschland GmbH
Publication of EP1613980A1 publication Critical patent/EP1613980A1/en
Application granted granted Critical
Publication of EP1613980B1 publication Critical patent/EP1613980B1/en
Anticipated expiration legal-status Critical
Expired - Lifetime legal-status Critical Current

Links

Classifications

    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/88Lidar systems specially adapted for specific applications
    • G01S17/95Lidar systems specially adapted for specific applications for meteorological use
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S17/00Systems using the reflection or reradiation of electromagnetic waves other than radio waves, e.g. lidar systems
    • G01S17/02Systems using the reflection of electromagnetic waves other than radio waves
    • G01S17/50Systems of measurement based on relative movement of target
    • G01S17/58Velocity or trajectory determination systems; Sense-of-movement determination systems
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01SRADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
    • G01S7/00Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
    • G01S7/48Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S17/00
    • G01S7/497Means for monitoring or calibrating
    • GPHYSICS
    • G01MEASURING; TESTING
    • G01WMETEOROLOGY
    • G01W1/00Meteorology
    • G01W2001/003Clear air turbulence detection or forecasting, e.g. for aircrafts
    • 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

Definitions

  • the present invention relates to a method for detecting wind speeds with a Doppler lidar system according to the preamble of claim 1 and a Doppler lidar system for detecting wind speeds, in particular on board aircraft, according to the preamble of claim 13.
  • a resolution of the speed of the air particles is of the order of 1 m / sec. and their location in front of the aircraft with a resolution in the area of the wing extension of about 5 - 10 m depth necessary.
  • Weather radar, laser anemometer and lidar are examples of such remote measurement methods.
  • weather radar systems have large, high-density suspended particles, such as, for their backscatter measurement
  • Wavelength of a few centimeters compared to the size of the
  • Air particles is relatively long.
  • pulse length and low pulse frequency can only be used for long ranges in the range of 5 - 20 km and for large measurement volumes with edge lengths of 50 - 500 m. Therefore, such systems do not allow control by rudder deflections, but only timely flying around large-scale turbulence zones.
  • the received photons are distributed over an interference pattern with an imaging, spatially resolving detector over several reception channels, namely in the Fabry-Perot interferometer on concentric interference rings and in the Fiz ⁇ au interferometer on interference fringes.
  • an imaging, spatially resolving detector over several reception channels, namely in the Fabry-Perot interferometer on concentric interference rings and in the Fiz ⁇ au interferometer on interference fringes.
  • Doppler lidar measurements therefore have the general problem of the low intensity of the backscattering of both molecules and aerosols.
  • a photodetector that is as sensitive as possible must be used and at the same time the influence of noise from background radiation from the
  • Atmosphere and can be suppressed efficiently by electronic noise in the detector and amplifier.
  • the weak backscatter signals can be integrated over several seconds, since the winds are stable over longer periods.
  • such systems and methods are not suitable for measuring air movements from an aircraft, since for this purpose much more compact systems with lasers of lower power and smaller ones Reception optics are necessary.
  • the signal evaluation must take place in just a few tens of milliseconds.
  • a photomultiplier with a conventional photocathode but with a special anode structure is used, concentric anode rings of different widths being analog If the concentric, ring-shaped interferogram is imaged on the photocathode of the photomultiplier, the photoelectrons are amplified and reach the concentrically arranged photoanode rings, and a stepped charge imprint of the interference pattern then arises on the anode rings By shifting the charge levels on the anode rings, the position of the interference ring of the Fabry-Perot interferometer can then be approximately determined.
  • this method has the disadvantage that only a rough scanning of the interferogram is possible. Measurement errors also occur if the interferogram is not imaged concentrically with the anode rings.
  • a laser beam of a predetermined frequency is emitted to a spatial area and the light scattered back from the spatial area is received, an interferogram being generated with an interferometer for determining a Doppler shift and with an
  • the intensity distribution of the interferogram is measured by the photodetector, and the intensity distribution of the interferogram is compared with one or more reference patterns, which were previously determined for defined parameters, and the Doppler shift is determined from the comparison as a measure of the wind speed.
  • the method according to the invention ensures a sufficiently high level of measurement reliability for flight control at all flight heights and at all Weather conditions, because in addition to backscattering the light on aerosols, the Doppler shifted backscatter on air molecules is determined with high accuracy and used for wind measurement. Since the density of the air molecules only varies by a maximum of a factor of 5 in the typical flight altitudes up to 40,000 feet, the reference patterns and expected measurement signals are very reliable and can be determined stably for every flight altitude in all regions of the world.
  • the invention improves the reception and evaluation method in Doppler measurements so that with a reduced transmission power of the laser, a sufficiently high level of measurement certainty is possible for flight control in a wide variety of atmospheric conditions. Shear winds and gusts can thus be determined with a compact system at measurement distances of, for example, 50-200 m from the aircraft using direct, incoherent reception in order to enable regulation by aircraft control.
  • the disadvantages of the prior art described above are eliminated by the invention, since the light which is transmitted through the interferometer can be fully used and additionally all the information which is contained in the geometric distribution of the light intensity in the imaging plane of the interferometer, can be used to register the Doppler shift.
  • the method used is particularly designed so that it is as insensitive as possible to signal noise and optical aberrations and interference in the interferometer.
  • the parts used for the selected photo receiver are available conventionally.
  • short-wave laser light in particular in the UV range, is preferably used. All states of the atmosphere can thus be recorded. It is taken into account that the increase in the intensity of the backscatter signals on air molecules is proportional to ⁇ "4 , where ⁇ denotes the wavelength of the laser beam.
  • Laser wavelengths that are as short as possible are therefore used, which are preferably in the UV range
  • a short pulse duration for example ⁇ 10 ns
  • a high pulse repetition frequency for example> 100 Hz
  • YAG lasers at 0.355 ⁇ m and 0.266 ⁇ m can, so that a sufficiently high spatial resolution is achieved.
  • a sufficient atmospheric transmission is up to a lower limit in the UV of about 0.230 ⁇ m for
  • the invention also solves the problem of spectral broadening of laser backscattering on molecules due to collisions and thermal movement of the molecules.
  • This broadening can be greater, for example, by a factor of 10-100 than the smallest Doppler frequency shift in air turbulence that is to be recorded.
  • This broadening of the received signal to a larger spectral range has so far made detection of the Doppler shift in the spectrally uniformly distributed noise so difficult that sufficiently precise measurements were not possible.
  • the invention allows the Doppler shift to be determined precisely despite the spectral broadening.
  • the optical beam path of the lidar is designed such that the entire interference pattern or interferogram is imaged directly, that is to say without reshaping, scaling or locating on a two-dimensional detector array.
  • all previous knowledge about the expected local distribution of the optical received signals on the detector surface is used for the signal evaluation to determine the Doppler shift.
  • the invention is based on the knowledge that the shape and broadening of the molecular backscatter spectrum at a narrow-band laser frequency that is single-mode represents a well-defined spectral line, the intensity curve of which depends on the atmospheric density and temperature, ie the measurement height of the aircraft above sea level, is predictable.
  • the expectation profiles of the intensity distribution in the photodetector plane of the interferometer can be determined for each density and air temperature and as historical values in, for example, a two-dimensional look-up
  • the Doppler line shift can be determined in real time for each flight altitude depending on the air speed.
  • the ring-shaped interferogram of a Fabry-Perot interferometer is advantageously imaged directly on the 2-dimensional photodetector. This results in an image on the flat photo receiver without optical conversion and without
  • the scattered, Doppler-shifted laser signal of an aircraft from the atmosphere can be fully used in the Doppler-Lidar system with the Fabry-Perot interferometer, so that there is no loss of information due to the optical conversion.
  • the reference patterns preferably contain different densities and / or temperatures of the atmosphere as parameters.
  • reference patterns are available for different atmospheric conditions, for example, at different heights, from which the reference pattern is selected which has the least deviation from the recorded interferogram.
  • the complete theoretically expected can be used for parameter calculation Intensity distribution can be used, which is previously obtained from a model calculation.
  • the model calculation includes the modeling of the laser transmission beam, the modeling of the molecular velocity distribution and the modeling of the interferometer transfer function.
  • an intensified CCD camera which is available, for example, in the form of 2-dimensional gray-scale images, a determination of the air speed in the area to be measured as well as a determination of various atmospheric parameters such as can be made by comparison with the family of reference patterns for example, the measurement volume temperature or the pressure of the measurement volume.
  • the at least one reference pattern thus contains, for example, the speed of the atmosphere relative to the Doppler-Lidar system as a parameter.
  • the variation in the speed of the atmosphere relative to the Doppler-Lidar system is advantageously determined from several successive measurements. From the fast variations compared to the airspeed caused by turbulence or shear winds, the
  • the laser beam is preferably pulsed and in each case a part of a laser pulse is used to determine a temporal reference point.
  • the distance of the backscattering area of the room can be determined by means of the transit time of the remaining part of the laser pulse.
  • Part of the laser beam is preferably received and registered directly, ie without backscattering in the atmosphere, it being possible to determine a transfer function of the optical components of the system from the intensity distribution or to carry out a calibration.
  • the density and / or the temperature of the spatial area is preferably determined from the reference pattern with the smallest deviation from the measured interferogram.
  • the method according to the invention is advantageously carried out on board a moving system, for example on board an aircraft, helicopter or another aircraft, just as it can also be used on board a ship or on the ground.
  • the expected intensity distribution of the reference pattern can also be calculated from measured atmospheric parameters and / or flight parameters of an aircraft, or the stored reference pattern with the smallest deviation from the registered interferogram can be further adapted using measured parameters.
  • the laser beam is advantageously sent in different directions in order to determine the wind speed vector, that is to say the amount and direction of the wind speed in the area to be measured.
  • a Doppler lidar system for recording wind speeds, in particular on board aircraft, with a transmitter device for emitting a laser beam, a receiver device for receiving the laser beam scattered back into the atmosphere, and an interferometer for generating a Interferogram from the backscattered laser beam, a photodetector for determining the intensity distribution of the interferogram, and an evaluation unit for determining the Doppler shift as a measure of the wind speed of the atmosphere, the interferogram being imaged directly on the photodetector, and the evaluation unit having a memory with or more
  • the invention contributes to flight safety on the one hand, on the other hand, flight comfort is noticeably improved, and there is also a reduction in fuel costs.
  • the photodetector is advantageously a 2-dimensional photodetector which comprises an image intensifier and a CCD or CMOS array. This results in low costs, since the photodetector is available commercially and conventionally.
  • a 'transfer path for a part of the laser beam is provided between the transmitting device and the receiving device to the generated laser beam to register directly in the receiving device.
  • a glass fiber cable is provided which couples the laser beam of the transmitting device branched off by means of a beam splitter to a receiving telescope of the receiving device. This makes it possible to register the emitted laser beam directly and thus to determine the transfer function of the receiving device with the components present therein.
  • the receiving device can be calibrated from time to time.
  • the laser beam guided directly to the receiving device can be used as a time reference in order to define the transit time of the emitted and backscattered laser beam. This allows the measurement distance, ie the distance of the area in which the air speed is determined can be determined.
  • the interferometer is preferably a Fabry-Perot interferometer.
  • the interferogram is preferably ring-shaped and comprises concentric rings of the intensity distribution.
  • the laser advantageously generates pulsed laser beams in the UV range. This makes it possible to measure the backscatter on air molecules very precisely and to determine the Doppler shift.
  • the Doppler lidar system preferably comprises field-programmable gate arrays for calculating the reference pattern. This enables a quick calculation in real time.
  • the evaluation unit advantageously comprises a module for determining the transfer function of the components of the receiving device.
  • the Doppler lidar system according to the invention is preferably designed with corresponding components for carrying out the method according to the invention.
  • FIG. 1 shows a schematic representation of a Doppler lidar system according to a preferred embodiment of the invention
  • Figure 2 shows schematically an image of an interferogram on the photoreceptor of the system shown in Figure 1;
  • a laser 11 is used to generate a pulsed laser beam and is optically coupled to a transmission telescope 12, which is used to transmit a laser beam into the atmosphere in a predetermined direction.
  • the laser 11 and the transmitting telescope 12 thus form a transmitting device for emitting the laser beam.
  • a receiving device in the form of a receiving telescope 13 is used to receive the laser beam scattered back into the atmosphere.
  • the received laser beam is fed via an optical fiber 14 and a filter unit 15 to an interferometer 16, which in this preferred embodiment is a Fabry-Perot interferometer.
  • a photodetector 17 is used to determine the intensity distribution of an interferogram, which is generated by the interferometer 16 from the laser beam supplied.
  • the photodetector 17 is arranged in relation to the interferometer 16 so that e.g. annular interferogram is imaged directly on the photodetector 17.
  • the output of the photodetector 17 is electrically coupled to an evaluation unit 18, which comprises a memory 18a and a comparison unit 18b in the form of a microprocessor.
  • an evaluation unit 18 which comprises a memory 18a and a comparison unit 18b in the form of a microprocessor.
  • One or more reference patterns for interferograms, those for defined ones, are stored in the memory 18a atmospheric and possibly other parameters apply.
  • the comparison unit 17b serves to compare the interferogram imaged on the photoreceiver 17 with the one or more reference patterns and to determine the wind speed from the comparison. Additional parameters can be supplied to the comparison unit 18b as input data 22 in order to adapt, calculate or change the reference pattern.
  • a beam splitter 19 is connected between the laser 11 and the transmitting telescope 12 and is used to branch off part of the laser beam generated by the laser 11 and to feed it to the receiving telescope 13 via a transfer path 20 in the form of an optical fiber or reference measuring fiber.
  • a small part of the radiation emitted by the laser 11 reaches the receiving telescope 13 via the transfer path 20 and is used to determine the starting time of the transit time of the laser beam emitted and backscattered in the atmosphere at a defined distance. This defines the distance at which the wind speed is to be measured by specifying the time of reception.
  • the transfer function of the optical receiving system with its various components can be determined.
  • the directly supplied laser beam is fed to the interferometer 16 and the resulting interferogram is then imaged on the photodetector 17.
  • the receiving telescope 13 serves to receive the light scattered back from the measurement volume, which is a defined spatial area. Via the fiber 14 in the form of a glass fiber, this light is parallelized by means of an optical lens 21 and passed on to the filter 15.
  • the filter 15 serves to filter out the laser beam scattered back from the atmosphere from the background radiation of the sun.
  • the filter 15 is followed by the interferometer 16, which forms an interferogram from the light supplied.
  • the interferometer 16 can also be a Fizeau interferometer, for example.
  • the pattern generated by the interferometer consists of concentric circles, the intensity distribution in the pattern and the distance of the circles from the center from the wavelength of the light supplied, from the spectrum of the laser, from the Doppler shift, the double Widening in the atmosphere and the transfer function of the components of the receiving channel depends.
  • the generated interferogram contains the entire measured information about the atmosphere and is mapped as a whole without optical shaping or fading out or centering onto the detector plane of the two-dimensional photodetector 17 with its 2-dimensional pixel structure.
  • FIG. 2 shows a typical intensity distribution in an interferogram 34 reproduced by the photoreceptor or the photodetector 17. Because of the discrete pixel structure and the typically weak optical reception signals, a step-like 2-dimensional intensity image is produced which is fed to the comparison unit 18b implemented by a microprocessor.
  • this intensity profile or intensity pattern is compared with reference patterns which are available in the storage unit 18a in the form of look-up tables.
  • the microprocessor looks for one from the one or more look-up tables Intensity distribution pattern that comes closest to the measured one.
  • the value of the Doppler shift, which is a parameter of the stored pattern, is derived from this comparison, and the speed of the air particles in the measurement volume is derived from this.
  • the reference patterns or the set of reference patterns can be generated in various ways. For example, different reference patterns can be generated for different atmospheric temperatures and different air velocities relative to the measuring system and stored in the look-up tables. Furthermore, the print and the serve
  • Instrument function as a parameter for determining the reference pattern. It is also possible to determine one or more of these parameters during the measuring operation by sensors and to adapt the stored reference patterns on the basis of the current values.
  • the current parameters determined by the further sensors can then be in the form of
  • Input data 22 are fed to the comparison unit 18b (see FIG. 1).
  • the interference patterns contain the various parameters that influence a measured interferogram, so that the reference patterns apply to different heights, temperatures, speeds of the aircraft and also different instrument functions.
  • the instrument function can either be measured or calculated.
  • part of the radiation from the laser is directed directly into the receiving telescope 13.
  • the transfer function can then be determined in the meantime between the measurement of the backscattered laser light from the atmosphere and taken into account in the further evaluation.
  • the intensity distributions of the reference patterns are essentially compared with the intensity distribution of the measured interferogram compared.
  • the spectrum of the light scattered back from the atmosphere contains a collision broadening of the molecules as well as a thermal broadening which is dependent on the respective temperature.
  • the center of gravity of the spectrum is shifted by the Doppler shift.
  • the instrument function ie the transfer function of the optical components of the receiver, also has an effect on the spectrum determined. All these factors are taken into account in the reference samples.
  • the Doppler shift is measured with a certain measuring frequency, which can be freely determined.
  • the intensity distribution in this interferogram is the convolution of the spectrum of the laser used and the transfer function of all components of the receiving system, such as receiving telescope 13, glass fiber or fiber 14, lens 21, filter 15 and interferometer 16. This intensity distribution is used in the measurement to complete the look -Up tables and / or used to monitor the settings of the entire system, especially for calibration.
  • the laser 11 When measuring the wind speed within a spatial area spaced from the measuring system with the Doppler-Lidar system 0 shown in FIG. 1, the laser 11 is operated with a typical pulse duration of 10 ns, which corresponds to a pulse length of 1.5 m.
  • the receiving device is set so that the backscatter signal is received after a fixed pulse transit time from a fixed distance, for example 100 m.
  • the point in time of the passage of the pulsed laser radiation through the aperture of the transmitting telescope 12 is used as the time reference. This point in time is determined by directing part of the radiation via the transfer path 20 designed as glass fiber
  • Receiving telescope 13 is supplied, so without the way through the atmosphere.
  • the spectral intensity distribution of the directly supplied laser beam can also be recorded and stored in the imaging plane of the interferometer 16 on the photodetector 1. This results in the transfer function F 0 of all optical components as:
  • this function is the inherent instrument function and only rarely, i.e. changed only when the system is misaligned and degraded, does not need to be saved with every pulse of the laser, but can be measured again and again for initial self-calibration and for monitoring the system. If the measuring system is stable enough, there is also the possibility to determine the instrument function purely by calculation and to take it into account in further evaluations. The evaluation of these tests can be used in the calculation to support the data.
  • the light of the laser pulse scattered back from the measurement volume in the atmosphere at a later point in time is picked up by the receiving telescope 13 and the directly supplied one is carried out in the same way as before
  • Laser light is fed via the receiving telescope 13, the lens 21, the filter 14, the fiber 15 and the interferometer 16 to the photo receiver or photodetector 17.
  • the laser light scattered back from the atmosphere is broadened spectrally and due to the Doppler effect from the relative speed of the aircraft Air and the natural movement of the air mass shifted spectrally.
  • the measured intensity curve is compared with the stored 2-dimensional field of possible intensity curves depending on the atmospheric parameters density and temperature.
  • the intensity curve with the least difference between measurement and calculation is used in a model-based evaluation, or it is improved by an appropriate interpolation or parameter adjustment.
  • the Doppler function F D can now be obtained from this intensity curve by computational unfolding. From the width and the course of the Doppler function, the density and temperature of the air at corresponding flight altitudes can now be obtained.
  • the Doppler shift and the alignment of the measuring axis from the aircraft can then be used to calculate the speed of the aircraft in relation to the air mass and the internal movement of the air mass.
  • the reference values of the interferogram function are stored for different parameter sets, in particular laser wavelength ⁇ , atmospheric density p and air temperature T, in order to reduce the computational effort.
  • FIG. 4 schematically shows a wind speed measurement with the Doppler-Lidar system according to the invention on board an aircraft 50.
  • Different measuring beams 51 are sent in different directions in order to scan the air space in front of the aircraft 50.
  • the laser beams scattered back from different spatial areas 52a, 52b, 52c, 52d on aerosols and molecules are recorded and analyzed as described above, so that the wind speeds in the respective measuring direction are determined in the respective spatial area.
  • the distance D of the spatial areas 52a-d from the aircraft 50 is 100 m, for example.
  • the division of the measuring beam into the plurality of measuring beams 51 can e.g. with a superior scanner, the air space in front of the aircraft 50 being scanned simultaneously or successively in different directions.
  • the spatial areas 52a-d have, for example, an extension in the measuring direction of approximately 10 m.
  • the wind speed vector can be determined from the individual measurement results by projection in the different directions.
  • the Doppler speed components measured in the respective laser transmission beam direction are processed into an estimated speed vector, for example, by an over-determined system of equations with at least four individual measurements in four measuring directions, weighting factors for reducing the amount of past measurements being taken into account.
  • the invention enables an accurate measurement of the wind speed vector far in front of an aircraft in real time with the help of the interferometric measurement of the Doppler line shift of pulsed laser beams scattered back on air molecules and aerosols. Because of the arrangement according to the invention and the evaluation method according to the invention, a large range is achieved, i.e. only reduced backscatter intensities are required for the evaluation.
  • the gust measurement in real time far in front of the aircraft allows targeted control and trimming maneuvers to reduce the harmful effects of gusts on the aircraft.
  • a distance that corresponds to a flight time of 0.5 seconds allows targeted control and trimming maneuvers to reduce the harmful effects of gusts on the aircraft.
  • the measurement can be carried out from moving measurement platforms, for example from airplanes, helicopters, balloons, satellites or ships, or also from floor-mounted measurement stations.
  • a time window is set for the measurement data acquisition, which is synchronized with the laser beam pulsing.
  • the measurement data can be processed in real time, the reference values or patterns being used in parallel by using field-programmable arrays be calculated.
  • the optimal parameters for the wavelength ⁇ and the temperature T of the measurement volume can be determined from this.
  • the reference patterns which represent reference values of the interferogram function, can be stored for different parameter sets ⁇ and T. In addition to the memory accesses, only addition and formation of amounts for positive, whole numbers are required for calculating the reference values, as a result of which the area requirements of the local computer units in the field programmable gate array are kept very small.
  • the wavelength ⁇ and the measurement volume temperature T are determined during the evaluation.
  • the complete, theoretically expected intensity distribution which is previously obtained from a model calculation, is used to calculate the parameters.
  • the model calculation includes the modeling of the laser beam, the modeling of the molecular velocity distribution and the modeling of the interferometer transfer function.
  • the Doppler lidar according to the invention and the method according to the invention for detecting wind speeds can also be used for applications and operations on the ground, for example for the detection of shear winds at airports or for the measurement of wind profiles in the vertical direction through the atmosphere.

Landscapes

  • Engineering & Computer Science (AREA)
  • Physics & Mathematics (AREA)
  • Computer Networks & Wireless Communication (AREA)
  • General Physics & Mathematics (AREA)
  • Radar, Positioning & Navigation (AREA)
  • Remote Sensing (AREA)
  • Electromagnetism (AREA)
  • Optical Radar Systems And Details Thereof (AREA)
  • Radar Systems Or Details Thereof (AREA)

Abstract

In the case of a method of detecting wind velocities by means of a Doppler-lidar system (10), a laser beam of a defined frequency generated by means of a laser (11) is emitted by a transmitting device (12) toward a space area and the light backscattered from the space area is received by means of a receiving telescope (13). For determining a Doppler shift, an interferogram is generated by means of an interferometer (16), the intensity distribution of the interferogram being directly measured by means of a photodetector (17). The measured intensity distribution is compared with one or more reference patterns which had previously been determined for defined parameters and are filed in a memory device (18a). From the comparison, the Doppler shift is determined as a measurement for the wind velocity. The Doppler-lidar system (10) comprises an analyzing unit (18a, 18b) for implementing the method, having a comparison unit (18b) for the comparison of reference patterns with the measured interferogram.

Description

Verfahren zur Erfassung von Windgeschwindigkeiten mit einem Doppler-Lidar-System, insbesondere an Bord von Flugzeugen, und Doppler-Lidar-System Method for recording wind speeds with a Doppler lidar system, in particular on board aircraft, and Doppler lidar system
Die vorliegende Erfindung betrifft ein Verfahren zur Erfassung von Windgeschwindigkeiten mit einem Doppler-Lidar-System gemäß dem Oberbegriff von Patentanspruch 1 und ein Doppler-Lidar-System zur Erfassung von Windgeschwindigkeiten, insbesondere an Bord von Flugzeugen, gemäß dem Oberbegriff von Patentanspruch 13.The present invention relates to a method for detecting wind speeds with a Doppler lidar system according to the preamble of claim 1 and a Doppler lidar system for detecting wind speeds, in particular on board aircraft, according to the preamble of claim 13.
Die Erfassung von Windgeschwindigkeiten und Turbulenzen in der Atmosphäre ist insbesondere für den Flugverkehr von erheblicher Bedeutung. So können z.B. auf ein Flugzeug schädliche Böen einwirken, welche zu hohen dynamischen Belastungen der Flugzeugstruktur führen. In vielen Fällen kommt es sogar zu kritischen Flugzuständen aufgrund von besonders starken Luftturbulenzen in der Atmosphäre, die in einigen extremen Situationen in der Vergangenheit sogar zu Unfällen geführt haben. Darüber hinaus wirken sich Scherwinde und Luftturbulenzen negativ auf den Flugkomfort aus und führen zu einem erhöhten Treibstoffverbrauch.The detection of wind speeds and turbulence in the atmosphere is of particular importance for air traffic. For example, Inflicting gusts on an aircraft, which lead to high dynamic loads on the aircraft structure. In many cases there are even critical flight conditions due to particularly strong air turbulence in the atmosphere, which has even resulted in accidents in some extreme situations in the past. In addition, shear winds and air turbulence have a negative impact on flight comfort and lead to increased fuel consumption.
Deshalb wird in der Luftfahrttechnik versucht, Scherwinde und Luftturbulenzen zu vermessen und anschließend mit entsprechenden Ruderausschlägen zu kompensieren, d.h. aktiv auszuregeln. Jedoch ist eine direkte Messung der Luftbewegung am Flugzeugrumpf oder an den Tragflächen mit mechanischen Sonden, wie beispielsweise Beschleunigungsaufnehmern oder Drucksensoren, nicht geeignet, da keine ausreichende Zeit zur Betätigung der Flügelaktorik verbleibt. Es ist daher notwendig, die Luftbewegung mit Sensoren vor dem Flugzeug auf seiner Flugstrecke im Vorfeld zu messen, also einige Sekundenbruchteile, bevor sie durchflogen wird. Zur aktiven Ausregelung der Turbulenzen und Scherwinde ist z.B. bei Verkehrsflugzeugen eine Mindestentfernung von etwa 50 - 200 m vor dem Flugzeug erforderlich. Zusätzlich ist eine Auflösung der Geschwindigkeit der Luftteilchen in der Größenordnung von 1 m/Sek. und ihres Ortes vor dem Flugzeug mit einer Auflösung im Bereich der Flügelausdehnung von etwa 5 - 10 m Tiefe notwendig. Für derartige Fernmessverfahren kommen beispielsweise Wetter-Radar, Laser-Anemometer und Lidar in Frage.For this reason, aeronautical engineering tries to measure shear winds and air turbulence and then to compensate for them with appropriate rudder deflections, ie to actively correct them. However, a direct measurement of the air movement on the fuselage or on the wings with mechanical probes, such as, for example, acceleration sensors or pressure sensors, is not suitable, since there is no sufficient time to actuate the wing actuators. It is therefore necessary to measure the air movement in advance with sensors in front of the aircraft on its flight path, i.e. a few fractions of a second before it is flown through. For active control of turbulence and shear winds, for example in the case of commercial aircraft, a minimum distance of approximately 50-200 m is required in front of the aircraft. In addition, a resolution of the speed of the air particles is of the order of 1 m / sec. and their location in front of the aircraft with a resolution in the area of the wing extension of about 5 - 10 m depth necessary. Weather radar, laser anemometer and lidar are examples of such remote measurement methods.
Bei Wetter-Radar-Systemen besteht das Problem, dass sie für ihre Rückstreumessung große Schwebeteilchen hoher Dichte, wie beispielsweiseThe problem with weather radar systems is that they have large, high-density suspended particles, such as, for their backscatter measurement
Wassertropfen in Wolken, benötigen und deshalb nur zur Frühwarnung vonDrops of water in clouds, and therefore only for early warning of
Wetterwolken in größerer Entfernung bis zu einigen 10 Kilometern geeignet sind.Weather clouds at a greater distance up to a few 10 kilometers are suitable.
Die häufigsten Luftstörungen bei klarer Sicht, wie beispielsweise Scherwinde und sogenannte „clear air turbulances" in der Troposphäre und Jetstreams in der Stratosphäre können dagegen vom Wetter-Radar nicht erfasst werden, da dieThe most common air disturbances with clear visibility, such as shear winds and so-called "clear air turbulances" in the troposphere and jet streams in the stratosphere, on the other hand, cannot be detected by the weather radar, since the
Wellenlänge von einigen Zentimetern (X-Band) im Vergleich zur Größe derWavelength of a few centimeters (X-band) compared to the size of the
Luftteilchen relativ lang ist.Air particles is relatively long.
Deshalb wurde versucht, Laser-Doppler-Anemometer zu verwenden, die in der Strömungsmechanik zur Diagnose von Luftströmungen eingesetzt werden. Sie benötigen jedoch eine scharfe Fokussierung des Laser-Messstrahls, um ausreichend starke Rückstreusignale von einzelnen Aerosolen zu erzeugen und können daher nur für Messungen im Nahbereich bis höchstens einige Meter Entfernung verwendet werden. Ein derartiges System und Verfahren ist z.B. in der Druckschrift DE 40 13 702 C und in S. Damp, „Battery-driven Miniature LDA-Therefore, attempts have been made to use laser Doppler anemometers, which are used in fluid mechanics to diagnose air flows. However, you need a sharp focus of the laser measuring beam in order to generate sufficiently strong backscatter signals from individual aerosols and can therefore only be used for measurements in the close range up to a few meters away. Such a system and method is e.g. in the publication DE 40 13 702 C and in S. Damp, "Battery-driven Miniature LDA-
System with Semiconductor Laser Diode", in Application of Laser Anemometry to Fluid Mechanics, 4th International Symposium, Lissabon, Portugal, 1989, beschrieben.System with Semiconductor Laser Diode ", in Application of Laser Anemometry to Fluid Mechanics, 4th International Symposium, Lisbon, Portugal, 1989.
Um Bewegungen kleinster Luftteilchen eines ausgedehnten Volumens in größerer Entfernung zu messen, wurde deshalb versucht, sogenannte Doppler-Lidar- Verfahren, die in kohärente und inkohärente Lidars unterteilt werden, einzusetzen. Beim Aerosol-Lidar-Verfahrβn können jedoch nur Schwebeteilchen bzw. Aerosole bei einer entsprechend geringen Frequenzverbreiterung des Doppler- Rückstreusignals detektiert werden.In order to measure movements of the smallest air particles of an extended volume at a greater distance, attempts were made to use so-called Doppler lidar Use methods that are divided into coherent and incoherent lidars. With the aerosol lidar method, however, only suspended particles or aerosols can be detected with a correspondingly low frequency broadening of the Doppler backscatter signal.
In dem Artikel "Tracking weathers flight path", in IEEE Spectrum, Sept. 2000, Seiten 38 - 45, wird ein Aerosol-Lidar-System zur groben Anzeige von Wind- oder Turbulenzzonen beschrieben. Derartige Aerosol-Lidarsysteme mit Lasern im Infrarot-Wellenlängenbereich von 2 μm oder 10 μm bringen den Vorteil, dass Messungen an kleineren Aerosolen möglich sind. Auf Grund ihrer langenIn the article "Tracking weathers flight path", in IEEE Spectrum, Sept. 2000, pages 38-45, an aerosol lidar system for the rough display of wind or turbulence zones is described. Such aerosol lidar systems with lasers in the infrared wavelength range of 2 μm or 10 μm have the advantage that measurements on smaller aerosols are possible. Because of their long
Pulslänge und niedriger Pulsfrequenz sind sie jedoch nur für große Reichweiten im Bereich von 5 - 20 km und für große Messvolumen mit Kantenlängen von 50 - 500 m verwendbar. Deshalb ist mit derartigen Systemen keine Ausregelung durch Ruderausschläge möglich, sondern nur ein rechtzeitiges Umfliegen großräumiger Turbulenzzonen.However, pulse length and low pulse frequency can only be used for long ranges in the range of 5 - 20 km and for large measurement volumes with edge lengths of 50 - 500 m. Therefore, such systems do not allow control by rudder deflections, but only timely flying around large-scale turbulence zones.
Deshalb wurde versucht, Aerosol-Lidar-Systeme mit kürzeren Wellenlängen und höherer Pulsfrequenz zu entwickeln, wie z.B. in „Aerosol and cloud backscatter at 1.064, 1.54, and 0.53 μm by airborne hard target calibrated Nd: YAG/methane Raman lidar", Appl. Opt. Vol 36, pp. 3475 - 3490 (1997) beschrieben. Dabei besteht jedoch das weitere Problem, dass die Dichte der Aerosole in vielen Regionen der Erde und insbesondere in Luftstraßen in großen Höhen über den Ozeanen zu niedrig ist, um zuverlässige Messungen zur Ausregelung von Luftturbulenzen durchzuführen.Therefore, attempts have been made to develop aerosol lidar systems with shorter wavelengths and higher pulse frequency, e.g. in "Aerosol and cloud backscatter at 1,064, 1.54, and 0.53 μm by airborne hard target calibrated Nd: YAG / methane Raman lidar", Appl. Opt. Vol 36, pp. 3475 - 3490 (1997). However, there is further Problem that the density of aerosols in many regions of the world and especially in airways at high altitudes above the oceans is too low to carry out reliable measurements to control air turbulence.
Aus diesem Grunde wurde versucht, zusätzlich zu der Aerosolstreuung auch die Doppler-verschobene Rückstreuung an Luftmolekülen zu verwenden, um ausreichend hohe Messsicherheit für Flugregelungen in allen Flughöhen und bei allen Wetterbedingungen zu erreichen. In dem Artikel von D. Rees und I.S. McDermid, „Doppler lidar atmospheric wind sensor: reevaluation of a 355-nm iήcoherent Doppler lidar", Appl. Opt. Vol. 29, No. 28, pp 4133 - 4144 (1990), wird die Messung der molekularen Doppler-Verschiebung in Doppler-Lidar-Systemen mittels Fizeau- bzw. Fabry-Perot-Interferometern mit inkohärentem Empfang beschrieben. Bei der spektralen Messung werden die empfangenen Photonen über ein Interferenzmuster mit einem abbildenden, ortsauflösenden Detektor auf mehrere Empfangskanäle verteilt, und zwar im Fabry-Perot-Interferometer auf konzentrische Interferenzringe und beim Fizβau-Interferometer auf Interferenzstreifen. Dabei besteht jedoch das Problem, dass die Erkennung von schwachen Signalen gegenüber dem Rauschen mehrerer einzelner Empfangskanäle erschwert ist.For this reason, attempts have been made to use the Doppler-shifted backscatter on air molecules in addition to the aerosol scattering in order to achieve sufficiently high measurement reliability for flight controls at all flight altitudes and in all weather conditions. In the article by D. Rees and IS McDermid, "Doppler lidar atmospheric wind sensor: reevaluation of a 355-nm iήcoherent Doppler lidar", Appl. Opt. Vol. 29, No. 28, pp 4133-4144 (1990) described the measurement of the molecular Doppler shift in Doppler-Lidar systems using Fizeau or Fabry-Perot interferometers with incoherent reception. In the spectral measurement, the received photons are distributed over an interference pattern with an imaging, spatially resolving detector over several reception channels, namely in the Fabry-Perot interferometer on concentric interference rings and in the Fizβau interferometer on interference fringes. However, there is the problem that the detection of weak signals is more difficult compared to the noise of several individual reception channels.
Bei Doppler-Lidarmessungen besteht deshalb das generelle Problem der geringen Intensität der Rückstreuung sowohl von Molekülen als auch von Aerosolen. Um die Laserleistung, die für zuverlässige Messungen notwendig ist, trotzdem gering zu halten, muss ein möglichst empfindlicher Photodetektor verwendet werden und gleichzeitig der Einfluss von Rauschen durch Hintergrundstrahlung derDoppler lidar measurements therefore have the general problem of the low intensity of the backscattering of both molecules and aerosols. In order to keep the laser power, which is necessary for reliable measurements, low, a photodetector that is as sensitive as possible must be used and at the same time the influence of noise from background radiation from the
Atmosphäre und durch elektronisches Rauschen im Detektor und Verstärker entsprechend effizient unterdrückt werden.Atmosphere and can be suppressed efficiently by electronic noise in the detector and amplifier.
In dem Artikel von C. Nardell et al.,"GroundWinds New Hampshire and the LIDARFest 2000 Campaign", SPI Internat. Symp. on Opt. Science and Technology, San Diego 2001, SPIE 4484-05 wird ein Doppler-Lidar mit Direktempfang und UV-Lasern beschrieben, das zur Messung von Winden vom Erdboden aus in verschiedenen Höhen geeignet ist. Dabei wird die Luftbewegung von einer Bodenstation mit Hilfe von Lasern im UV-Bereich mit einer Ausgangsleistung von mehreren Watt und Empfangsteleskopen mit einemIn the article by C. Nardell et al., "GroundWinds New Hampshire and the LIDARFest 2000 Campaign", SPI Internat. Symp. On Opt. Science and Technology, San Diego 2001, SPIE 4484-05 describes a Doppler lidar with direct reception and UV lasers, which is suitable for measuring winds from the ground at different heights. The air movement from a ground station with the help of lasers in the UV range with an output power of several watts and receiving telescopes with one
Durchmesser von einigen zehn Zentimern bis zu einigen Kilometern Höhe in der Atmosphäre gemessen. Eine Integration der schwachen Rückstreusignale über mehrere Sekunden ist hierbei möglich, da die Winde über längere Zeiten stabil sind. Derartige Systeme und Verfahren sind jedoch zur Messung von Luftbewegungen von einem Flugzeug aus nicht geeignet, da für diesen Zweck wesentlich kompaktere Systeme mit Lasern geringerer Leistung und kleineren Empfangsoptiken notwendig sind. Weiterhin muss die Signalauswertung in nur einigen zehn Millisekunden erfolgen.Diameters measured from a few tens of centimeters to a few kilometers in height in the atmosphere. The weak backscatter signals can be integrated over several seconds, since the winds are stable over longer periods. However, such systems and methods are not suitable for measuring air movements from an aircraft, since for this purpose much more compact systems with lasers of lower power and smaller ones Reception optics are necessary. Furthermore, the signal evaluation must take place in just a few tens of milliseconds.
In dem Artikel von P.B. Hays und J. Wang „Image plane detector for Fabry-Perot interferometers: physical model and improvement with anticoincidence detection", Appl. Opt. Vol. 30, No. 22, pp 3100 - 3107 (1991 ) wird ein Verfahren zur Registrierung der Bewegung von Interferenzringen eines Fabry-Perot- Interferometers mit einem speziell dafür ausgelegten Detektor gezeigt. Dabei ist die Geometrie des Detektors dem Intensitätsverlauf im Interferenzmuster angepasst. Es wird ein Photomultiplier mit einer konventionellen Photokathode, jedoch mit spezieller Anodenstruktur, verwendet, wobei konzentrische Anodenringe unterschiedlicher Breite analog dem radialen Intensitätsverlauf im Interferogramm verwendet werden. Wird das konzentrische, ringförmige Interferogramm auf die Photokathode des Photomultipliers abgebildet, werden die Photoelektronen verstärkt und gelangen auf die konzentrisch angeordneten Photoanodenringe. Auf den Anodenringen entsteht dann ein abgestufter Ladungsabdruck des Interferenzmusters. Aus der Verschiebung der Ladungsstufen auf den Anodenringen kann die Lage des Interferenzringes des Fabry-Perot-Interferometers dann näherungsweise bestimmt werden. Dieses Verfahren hat jedoch den Nachteil, dass nur eine grobe Abtastung des Interferogramms möglich ist. Weiterhin entstehen Messfehler, wenn das Interferogramm nicht konzentrisch zu den Anodenringen abgebildet wird.In the article by P.B. Hays and J. Wang "Image plane detector for Fabry-Perot interferometers: physical model and improvement with anticoincidence detection", Appl. Opt. Vol. 30, No. 22, pp 3100-3107 (1991) becomes a method for registering the movement of interference rings of a Fabry-Perot interferometer with a specially designed detector. The geometry of the detector is adapted to the intensity curve in the interference pattern. A photomultiplier with a conventional photocathode but with a special anode structure is used, concentric anode rings of different widths being analog If the concentric, ring-shaped interferogram is imaged on the photocathode of the photomultiplier, the photoelectrons are amplified and reach the concentrically arranged photoanode rings, and a stepped charge imprint of the interference pattern then arises on the anode rings By shifting the charge levels on the anode rings, the position of the interference ring of the Fabry-Perot interferometer can then be approximately determined. However, this method has the disadvantage that only a rough scanning of the interferogram is possible. Measurement errors also occur if the interferogram is not imaged concentrically with the anode rings.
Um diese Nachteile zu überwinden, beschreiben J. Wu et al. in „Performance of a circle-to-line optical System for a Fabry-Perot-Interferometer: a laboratory study", Appl. Opt. Vol. 33 pp. 7823 - 7828 (1994) einen Empfänger, bei dem mittels einer verspiegelten Innenfläche eines Konus ein Viertel des kreisförmigen Interferogramms des Fabry-Perot-Interferometers in Zeilenstrukturen umgeformt wird und auf einen linearen Detektor oder ein CCD-Array abgebildet wird. Dieses Verfahren hat jedoch den Nachteil, dass nur ein Viertel von dem Interferogramm genutzt wird und optische Verzerrungen bei der Abbildung entstehen, die zu zusätzlichen Messfehlern führen.To overcome these disadvantages, J. Wu et al. in "Performance of a circle-to-line optical system for a Fabry-Perot interferometer: a laboratory study", Appl. Opt. Vol. 33 pp. 7823 - 7828 (1994) a receiver, in which a A quarter of the circular interferogram of the Fabry-Perot interferometer is converted into line structures and imaged on a linear detector or a CCD array, but this method has the disadvantage that only a quarter of the interferogram is used and optical distortions occur in the image, which lead to additional measurement errors.
In C.L. Korb et al., „Edge technique Doppler lidar wind measurements with high vertical resolution", Appl. Opt. Vol. 36, 5976 - 5983 (1997) wird ein Verfahren zur Bestimmung der Dopplerverschiebung in Molekül-Lidars durch Flankendetektion beschrieben. Dabei wird mit einem schmalbandigen spektralen Filter wie dem Fabry-Perot-Interferometer der spektrale Durchlassbereich des Filters gegenüber der Doppler-verbreiterten Linie der Rückstreuung des Lasers aus der Atmosphäre so eingestellt, dass nur ein Teil des Spektrums an seiner Intensitätsflanke durchgelassen wird. Geringste spektrale Verschiebungen der Linie werden dann als Intensitätsänderung entlang der Flanke mit einem einkanaiigen Detektor registriert, die dann ein Maß der Doppler-Verschiebung des gesamten Spektrums darstellen.In C.L. Korb et al., "Edge technique Doppler lidar wind measurements with high vertical resolution", Appl. Opt. Vol. 36, 5976-5983 (1997) describes a method for determining the Doppler shift in molecular lidars by edge detection a narrow-band spectral filter such as the Fabry-Perot interferometer, the spectral passband of the filter compared to the Doppler-broadened line of backscattering the laser from the atmosphere is set in such a way that only a part of the spectrum is transmitted on its intensity flank then registered as a change in intensity along the flank with a single-channel detector, which then represent a measure of the Doppler shift of the entire spectrum.
Bei diesem Verfahren der Flankendetektion in Molekulardoppler-Lidars besteht ebenso wie bei der sogenannten „Fringe-Detektion" der Nachteil, dass nur ein Teil des empfangenen Lichts genutzt wird und weiterhin nur ein geringer Teil der verfügbaren spektralen Information im Empfangslicht zur Ermittlung der Doppler- Verschiebung zum Tragen kommt. Darüber hinaus entstehen hohe Kosten, da spezielle und sehr teure Photoempfänger eingesetzt werden müssen.With this method of edge detection in molecular Doppler lidars, just like with the so-called “fringe detection”, there is the disadvantage that only part of the received light is used and furthermore only a small part of the available spectral information in the received light for determining the Doppler shift In addition, there are high costs because special and very expensive photo receivers have to be used.
In Todd D. Irgang et al., „Two-channel direct-detection Doppler lidar employing a charge-coupled device as a detector", Appl. Optics, Vol. 41, No. 6, S. 1145 - 1155 (2002) wird ein Doppler-Lidar-System beschrieben, das einen CCD-Detektor umfasst. Durch Verwendung einer Konvertierungseinrichtung, die kreisförmige Interferenzmuster eines Fabry-Perot-Etalons in streifenförmige Muster umwandelt, kann diese lineare Vorrichtung verwendet werden. Dadurch sollen durch Rückstreuung an Aerosolen und Molekülen Winde gemessen werden. Die kreisförmigen Ringe, die vom Fabry-Perot-Etalon ausgegeben werden, werden mittels eines reflektierenden Konus in eine lineare Streifenserie umgewandelt. Jedoch hat auch dieses System und Verfahren den Nachteil, dass aufgrund des Konverters optische Fehler entstehen und ein großer Aufwand notwendig ist.In Todd D. Irgang et al., "Two-channel direct-detection Doppler lidar employing a charge-coupled device as a detector", Appl. Optics, Vol. 41, No. 6, pp. 1145-1155 (2002) described a Doppler-Lidar system which comprises a CCD detector. By using a conversion device which converts circular interference patterns of a Fabry-Perot etalon into strip-like patterns, this linear device can be used. This should result in winds due to backscattering on aerosols and molecules The circular rings emitted by the Fabry-Perot etalon are converted into a linear series of strips using a reflective cone. However, this system and method also has the disadvantage that optical errors occur due to the converter and a great deal of effort is required.
Es ist die Aufgabe der vorliegenden Erfindung, ein Verfahren zur Messung von Windgeschwindigkeiten mit einem Doppler-Lidar-System anzugeben, das für eine nachfolgende Flugregelung zur Ausregelung von Turbulenzen oder Scherwinden im Bereich der Luftfahrt geeignet ist und eine hohe Messsicherheit gewährleistet. Weiterhin soll ein Doppier-Lidar-System angegeben werden, mit dem Scherwinde und Turbulenzen mit hoher Genauigkeit detektiert werden können, so dass eine Ausregelung von Turbulenzen oder Scherwinden im Flugzeug möglich ist.It is the object of the present invention to provide a method for measuring wind speeds with a Doppler-Lidar system which is suitable for a subsequent flight control for regulating turbulence or shear winds in the field of aviation and which ensures a high level of measurement reliability. Furthermore, a double lidar system is to be specified, with which shear winds and turbulence can be detected with high accuracy, so that regulation of turbulence or shear winds in the aircraft is possible.
Diese Aufgabe wird gelöst durch das Verfahren zur Erfassung von Windgeschwindigkeiten gemäß Patentanspruch 1 und durch das Doppler-Lidar- System zur Erfassung von Windgeschwindigkeiten, insbesondere an Bord von Flugzeugen gemäß Patentanspruch 13. Weitere vorteilhafte Merkmale, Aspekte und Details der Erfindung ergeben sich aus den abhängigen Ansprüchen, der Beschreibung und den Zeichnungen.This object is achieved by the method for detecting wind speeds according to claim 1 and by the Doppler-Lidar system for detecting wind speeds, in particular on board aircraft according to claim 13. Further advantageous features, aspects and details of the invention result from the dependent Claims, the description and the drawings.
Bei dem erfindungsgemäßen Verfahren zur Erfassung von Windgeschwindigkeiten mit einem Doppler-Lidar-System wird ein Laserstrahl vorgegebener Frequenz zu einem Raumbereich hin ausgesandt und das aus dem Raumbereich zurückgestreute Licht empfangen, wobei zur Bestimmung einer Doppler-Verschiebung mit einem Interferometer ein Interferogramm erzeugt und mit einem Photodetektor die Intensitätsverteilung des Interferogramms gemessen wird, und wobei die Intensitätsverteilung des Interferogramms mit ein oder mehreren Referenzmustern verglichen wird, die zuvor für definierte Parameter bestimmt wurden, und aus dem Vergleich die Doppler-Verschiebung als Maß für die Windgeschwindigkeit ermittelt wird.In the method according to the invention for detecting wind speeds with a Doppler lidar system, a laser beam of a predetermined frequency is emitted to a spatial area and the light scattered back from the spatial area is received, an interferogram being generated with an interferometer for determining a Doppler shift and with an The intensity distribution of the interferogram is measured by the photodetector, and the intensity distribution of the interferogram is compared with one or more reference patterns, which were previously determined for defined parameters, and the Doppler shift is determined from the comparison as a measure of the wind speed.
Das erfindungsgemäße Verfahren gewährleistet eine ausreichend hohe Messsicherheit für Flugregelung in allen Flughöhen und bei allen Wetterbedingungen, da zusätzlich zur Rückstreuung des Lichts an Aerosolen auch die Doppler-verschobene Rückstreuung an Luftmolekülen mit hoher Genauigkeit ermittelt und zur Windmessung verwendet wird. Da die Dichte der Luftmoleküle nur um höchstens einen Faktor 5 in den typischen Flughöhen bis zu 40.000 Fuß variiert, sind die Referenzmuster bzw. erwarteten Messsignale sehr zuverlässig und lassen sich stabil für jede Flughöhe in allen Regionen der Erde ermitteln. Durch die Erfindung wird das Empfangs- und Auswerteverfahren bei Doppler-Messungen so verbessert, dass bei verringerter Sendeleistung des Lasers eine ausreichend hohe Messsicherheit für die Flugregelung bei unterschiedlichsten atmosphärischen Bedingungen möglich ist. Somit können Scherwinde und Böen mit einem kompakten System bei Messentfernungen von beispielsweise 50 - 200 m aus dem Flugzeug unter Verwendung des direkten, inkohärenten Empfangs ermittelt werden, um eine Ausregelung durch Flugzeugsteuerung zu ermöglichen.The method according to the invention ensures a sufficiently high level of measurement reliability for flight control at all flight heights and at all Weather conditions, because in addition to backscattering the light on aerosols, the Doppler shifted backscatter on air molecules is determined with high accuracy and used for wind measurement. Since the density of the air molecules only varies by a maximum of a factor of 5 in the typical flight altitudes up to 40,000 feet, the reference patterns and expected measurement signals are very reliable and can be determined stably for every flight altitude in all regions of the world. The invention improves the reception and evaluation method in Doppler measurements so that with a reduced transmission power of the laser, a sufficiently high level of measurement certainty is possible for flight control in a wide variety of atmospheric conditions. Shear winds and gusts can thus be determined with a compact system at measurement distances of, for example, 50-200 m from the aircraft using direct, incoherent reception in order to enable regulation by aircraft control.
Durch die Erfindung werden die oben beschriebenen Nachteile des Standes der Technik beseitigt, da das Licht, das durch das Interferometer hindurchgelassen wird, vollständig genutzt werden kann und zusätzlich die gesamte Information, die in der geometrischen Verteilung der Lichtintensität in der Abbildungsebene des Interferometers enthalten ist, zur Registrierung der Doppler-Verschiebung verwendet werden kann. Dabei wird das verwendete Verfahren insbesondere so ausgelegt, dass es gegenüber Signalrauschen und optischen Aberationen und Störungen im Interferometer möglichst unempfindlich ist. Weiterhin sind die verwendeten Teile des ausgewählten Photoempfängers konventionell erhältlich.The disadvantages of the prior art described above are eliminated by the invention, since the light which is transmitted through the interferometer can be fully used and additionally all the information which is contained in the geometric distribution of the light intensity in the imaging plane of the interferometer, can be used to register the Doppler shift. The method used is particularly designed so that it is as insensitive as possible to signal noise and optical aberrations and interference in the interferometer. Furthermore, the parts used for the selected photo receiver are available conventionally.
Um die Rückstreuung des Laserlichts in dem Gemisch von Molekülen und Aerosolen und zusätzlich auch an kondensierten Wassertröpfchen in Wolken zu ermöglichen, wird vorzugsweise kurzwelliges Laserlicht, insbesondere im UV- Bereich, verwendet. Somit können alle Zustände der Atmosphäre erfasst werden. Dabei wird berücksichtigt, dass die Zunahme der Intensität der Rückstreusignale an Luftmolekülen proportional zu λ"4 erfolgt, wobei λ die Wellenlänge des Laserstrahls bezeichnet. Es werden deshalb möglichst kurze Laserwellenlängen verwendet, die vorzugsweise im UV-Bereich liegen. Es ergeben sich dadurch noch die weiteren Vorteile, dass auch eine kurze Pulsdauer, beispielsweise <10 ns, und eine hohe Pulswiederholfrequenz, beispielsweise > 100 Hz, mit marktgängigen Lasern, wie z.B. Frequenz-verdreifachten oder -vervierfachten Nd: YAG-Lasern bei 0,355 μm bzw. 0,266 μm erreicht werden kann, so dass eine ausreichend hohe Ortsauf lösung erzielt wird. Eine ausreichende atmosphärische Transmission ist bis zu einer unteren Grenze im UV von etwa 0,230 μm fürIn order to enable the backscattering of the laser light in the mixture of molecules and aerosols and additionally also on condensed water droplets in clouds, short-wave laser light, in particular in the UV range, is preferably used. All states of the atmosphere can thus be recorded. It is taken into account that the increase in the intensity of the backscatter signals on air molecules is proportional to λ "4 , where λ denotes the wavelength of the laser beam. Laser wavelengths that are as short as possible are therefore used, which are preferably in the UV range Further advantages that a short pulse duration, for example <10 ns, and a high pulse repetition frequency, for example> 100 Hz, can be achieved with commercially available lasers, such as frequency tripled or quadrupled Nd: YAG lasers at 0.355 μm and 0.266 μm can, so that a sufficiently high spatial resolution is achieved. A sufficient atmospheric transmission is up to a lower limit in the UV of about 0.230 μm for
Fernmessungen bis zu einigen Hundert Metern bei allen Flughöhen vorhanden.Remote measurements up to a few hundred meters available at all flight altitudes.
Durch die Erfindung wird insbesondere auch das Problem der spektralen Verbreiterung der Laserrückstreuung an Molekülen auf Grund von Stößen und thermischer Bewegung der Moleküle gelöst. Diese Verbreiterung kann beispielsweise um einen Faktor von 10 - 100 größer sein als die geringste Doppler-Frequenzverschiebung an Luftturbulenzen, die erfasst werden sollen. Diese Verbreiterung des Empfangssignals auf einen größeren spektralen Bereich hat bisher die Detektion der Doppler-Verschiebung in dem spektral gleichverteilten Rauschen so erschwert, dass ausreichend genaue Messungen nicht möglich waren. Durch die Erfindung kann die Doppler-Verschiebung trotz der spektralen Verbreiterung genau ermittelt werden.In particular, the invention also solves the problem of spectral broadening of laser backscattering on molecules due to collisions and thermal movement of the molecules. This broadening can be greater, for example, by a factor of 10-100 than the smallest Doppler frequency shift in air turbulence that is to be recorded. This broadening of the received signal to a larger spectral range has so far made detection of the Doppler shift in the spectrally uniformly distributed noise so difficult that sufficiently precise measurements were not possible. The invention allows the Doppler shift to be determined precisely despite the spectral broadening.
Insbesondere wird bei der Erfindung der optische Strahlengang des Lidars so ausgelegt, dass das gesamte Interferenzmuster bzw. Interferogramm direkt, also ohne Umformung, Skalierung oder Ortseinstellung auf ein zweidimensionales Detektor-Array, abgebildet wird. Gleichzeitig werden bei der Signalauswertung z.B. alle Vorkenntnisse über die zu erwartende örtliche Verteilung der optischen Empfangssignale auf der Detektorfläche zur Ermittlung der Doppler-Verschiebung verwendet. Die Erfindung basiert auf der Erkenntnis, dass die Gestalt und Verbreiterung des molekularen Rückstreuspektrums bei einer schmalbandigen Laserfrequenz, die monomodig ist, eine wohldefinierte spektrale Linie darstellt, deren Intensitätsverlauf in Abhängigkeit von der atmosphärischen Dichte und Temperatur, d.h. der Messhöhe des Flugzeugs über dem Meeresspiegel, vorausberechenbar ist. Insbesondere können mit der zusätzlichen Kenntnis der Transferfunktion des Interferomβters, die auch mit hinreichender Genauigkeit berechenbar bzw. messbar ist, die Erwartungsprofile der Intensitätsverteilung in der Photodetektorebene des Interferometers für jede Dichte und Lufttemperatur bestimmt werden und als Verlaufswerte in einer beispielsweise zweidimensionalen „look-up-table" abgespeichert werden. Durch den Vergleich dieser modellierten spektralen Profile mit den tatsächlich gemessenen Spektren und den dabei festgestellten Abweichungen, kann die Doppler-Linienverschiebung in Abhängigkeit der Luftgeschwindigkeit für jede Flughöhe in Echtzeit bestimmt werden.In particular, in the invention the optical beam path of the lidar is designed such that the entire interference pattern or interferogram is imaged directly, that is to say without reshaping, scaling or locating on a two-dimensional detector array. At the same time, all previous knowledge about the expected local distribution of the optical received signals on the detector surface is used for the signal evaluation to determine the Doppler shift. The invention is based on the knowledge that the shape and broadening of the molecular backscatter spectrum at a narrow-band laser frequency that is single-mode represents a well-defined spectral line, the intensity curve of which depends on the atmospheric density and temperature, ie the measurement height of the aircraft above sea level, is predictable. In particular, with the additional knowledge of the transfer function of the interferometer, which can also be calculated or measured with sufficient accuracy, the expectation profiles of the intensity distribution in the photodetector plane of the interferometer can be determined for each density and air temperature and as historical values in, for example, a two-dimensional look-up By comparing these modeled spectral profiles with the actually measured spectra and the deviations found, the Doppler line shift can be determined in real time for each flight altitude depending on the air speed.
Vorteilhafterweise wird bei dem erfindungsgemäßen Verfahren das ringförmige Interferogramm eines Fabry-Perot-Interferometers direkt auf den 2-dimensionalen Photodetektor abgebildet. Dadurch ergibt sich eine Abbildung auf dem flächenhaften Photoempfänger ohne optische Umwandlung und ohneIn the method according to the invention, the ring-shaped interferogram of a Fabry-Perot interferometer is advantageously imaged directly on the 2-dimensional photodetector. This results in an image on the flat photo receiver without optical conversion and without
Abschattung. Somit kann beispielsweise das gestreute, Doppler-verschobene Lasersignal eines Flugzeugs aus der Atmosphäre in dem Doppler-Lidar-System mit dem Fabry-Perot-Interferometer vollständig genutzt werden, so dass kein Informationsverlust durch die optische Umwandlung erfolgt.Shadowing. Thus, for example, the scattered, Doppler-shifted laser signal of an aircraft from the atmosphere can be fully used in the Doppler-Lidar system with the Fabry-Perot interferometer, so that there is no loss of information due to the optical conversion.
Bevorzugt enthalten die Referenzmuster als Parameter unterschiedliche Dichten und/oder Temperaturen der Atmosphäre. Dadurch sind für unterschiedliche atmosphärische Bedingungen, die z.B. in verschiedenen Höhen herrschen, Referenzmuster vorhanden, aus denen dasjenige Referenzmuster ausgewählt wird, das die geringst Abweichung zum aufgenommenen Interferogramm hat. Zur Parameterberechnung kann die vollständige theoretisch zu erwartende Intensitätsverteilung genutzt werden, die vorher aus einer Modellrechnung gewonnen wird. Dabei schließt die Modellrechnung die Modellierung des Lasersendestrahls, die Modellierung der Molekülgeschwindigkeitsverteilung und die Modellierung der Interferometertransferfunktion ein. Aus dem mit einer intensivierten CCD-Kamera aufgenommenen Interferogramm, das beispielsweise in Form von 2-dimensionalen Grauwertbildern vorliegt, kann dann durch Vergleich mit der Schar von Referenzmustern sowohl eine Bestimmung der Luftgeschwindigkeit in dem zu messenden Raumbereich als auch eine Bestimmung verschiedener atmosphärischer Parameter, wie beispielsweise die Messvolumentemperatur oder der Druck des Messvolumens, durchgeführt werden.The reference patterns preferably contain different densities and / or temperatures of the atmosphere as parameters. As a result, reference patterns are available for different atmospheric conditions, for example, at different heights, from which the reference pattern is selected which has the least deviation from the recorded interferogram. The complete theoretically expected can be used for parameter calculation Intensity distribution can be used, which is previously obtained from a model calculation. The model calculation includes the modeling of the laser transmission beam, the modeling of the molecular velocity distribution and the modeling of the interferometer transfer function. From the interferogram recorded with an intensified CCD camera, which is available, for example, in the form of 2-dimensional gray-scale images, a determination of the air speed in the area to be measured as well as a determination of various atmospheric parameters such as can be made by comparison with the family of reference patterns for example, the measurement volume temperature or the pressure of the measurement volume.
Das mindestens eine Referenzmuster enthält somit beispielsweise die Geschwindigkeit der Atmosphäre relativ zum Doppler-Lidar-System als Parameter.The at least one reference pattern thus contains, for example, the speed of the atmosphere relative to the Doppler-Lidar system as a parameter.
Vorteilhafterweise wird aus mehreren aufeinanderfolgenden Messungen die Variation der Geschwindigkeit der Atmosphäre relativ zum Doppler-Lidar-System ermittelt. Aus den im Vergleich zur Fluggeschwindigkeit schnellen Variationen, die durch Turbulenzen oder Scherwinde verursacht werden, lässt sich dann dieThe variation in the speed of the atmosphere relative to the Doppler-Lidar system is advantageously determined from several successive measurements. From the fast variations compared to the airspeed caused by turbulence or shear winds, the
Luftgeschwindigkeit in dem vor dem Flugzeug liegenden Raumbereich ermitteln.Determine the air speed in the area in front of the aircraft.
Bevorzugt ist der Laserstrahl gepulst und jeweils ein Teil eines Laserpulses wird zur Festlegung eines zeitlichen Referenzpunktes verwendet. Mittels der Laufzeit des restlichen Teils des Laserpulses kann die Entfernung des rückstreuenden Raumbereichs festgelegt werden.The laser beam is preferably pulsed and in each case a part of a laser pulse is used to determine a temporal reference point. The distance of the backscattering area of the room can be determined by means of the transit time of the remaining part of the laser pulse.
Vorzugsweise wird ein Teil des Laserstrahls direkt, d.h. ohne Rückstreuung in der Atmosphäre, empfangen und registriert, wobei aus der Intensitätsverteilung eine Transferfunktion der optischen Komponenten des Systems bestimmt oder auch eine Kalibrierung durchgeführt werden kann. Bevorzugt wird aus dem Referenzmuster mit der geringsten Abweichung zum gemessenen Interferogramm die Dichte und/oder die Temperatur des Raumbereichs ermittelt.Part of the laser beam is preferably received and registered directly, ie without backscattering in the atmosphere, it being possible to determine a transfer function of the optical components of the system from the intensity distribution or to carry out a calibration. The density and / or the temperature of the spatial area is preferably determined from the reference pattern with the smallest deviation from the measured interferogram.
Vorteilhafterweise wird das erfindungsgemäße Verfahren an Bord eines bewegten Systems, beispielsweise an Bord eines Flugzeugs, Hubschraubers oder eines anderen Luftfahrzeugs durchgeführt, ebenso wie es auch an Bord eines Schiffes oder am Boden eingesetzt werden kann.The method according to the invention is advantageously carried out on board a moving system, for example on board an aircraft, helicopter or another aircraft, just as it can also be used on board a ship or on the ground.
Weiterhin kann die erwartete Intensitätsverteilung des Referenzmusters auch aus gemessenen atmosphärischen Parametern und/oder Flugparametern eines Luftfahrzeugs berechnet werden, oder das gespeicherte Referenzmuster mit der geringsten Abweichung zum registrierten Interferogramm kann durch gemessene Parameter weiter angepasst werden .Furthermore, the expected intensity distribution of the reference pattern can also be calculated from measured atmospheric parameters and / or flight parameters of an aircraft, or the stored reference pattern with the smallest deviation from the registered interferogram can be further adapted using measured parameters.
Vorteilhafterweise wird der Laserstrahl in verschiedene Richtungen gesandt, um den Windgeschwindigkeitsvektor, also Betrag und Richtung der Windgeschwindigkeit in dem zu messenden Raumbereich, zu bestimmen.The laser beam is advantageously sent in different directions in order to determine the wind speed vector, that is to say the amount and direction of the wind speed in the area to be measured.
Gemäß einem anderen Aspekt der Erfindung wird ein Doppler-Lidar-System zur Erfassung von Windgeschwindigkeiten, insbesondere an Bord von Flugzeugen, geschaffen, mit einer Sendeeinrichtung zum Aussenden eines Laserstrahls, einer Empfangseinrichtung zum Empfang des in der Atmosphäre zurückgestreuten Laserstrahls, einem Interferometer zur Erzeugung eines Interferogramms aus dem zurückgestreuten Laserstrahl, einem Photodetektor zur Bestimmung der Intensitätsverteilung des Interferogramms, und einer Auswerteeinheit zur Bestimmung der Doppler-Verschiebung als Maß für die Windgeschwindigkeit der Atmosphäre, wobei das Interferogramm direkt auf den Photodetektor abgebildet wird, und die Auswerteeinheit einen Speicher mit ein oder mehrerenAccording to another aspect of the invention, a Doppler lidar system for recording wind speeds, in particular on board aircraft, is provided, with a transmitter device for emitting a laser beam, a receiver device for receiving the laser beam scattered back into the atmosphere, and an interferometer for generating a Interferogram from the backscattered laser beam, a photodetector for determining the intensity distribution of the interferogram, and an evaluation unit for determining the Doppler shift as a measure of the wind speed of the atmosphere, the interferogram being imaged directly on the photodetector, and the evaluation unit having a memory with or more
Referenzmustern umfasst, die für zuvor definierte atmosphärische Parameter gelten, und eine Vergleichseinheit zum Vergleichen des abgebildeten Interferogramms mit den Referenzmustern, um aus dem Vergleich die Windgeschwindigkeit zu ermitteln.Includes reference patterns for previously defined atmospheric parameters apply, and a comparison unit for comparing the depicted interferogram with the reference patterns in order to determine the wind speed from the comparison.
Mit dem erfindungsgemäßen Doppler-Lidar-System ist es möglich, mit hoherWith the Doppler-Lidar system according to the invention, it is possible with high
Genauigkeit die Windgeschwindigkeit bzw. Turbulenzen vor einem Flugzeug so zu bestimmen, dass beim Durchfliegen des zu messenden Raumbereichs eine Ausregelung durch entsprechende Klappensteuerung erfolgen kann. Dadurch trägt die Erfindung zum einen zur Flugsicherheit bei, zum anderen wird der Flugkomfort spürbar verbessert, und es erfolgt darüber hinaus eine Reduktion der Treibstoffkosten.Accuracy to determine the wind speed or turbulence in front of an aircraft so that when flown through the area to be measured can be regulated by appropriate flap control. As a result, the invention contributes to flight safety on the one hand, on the other hand, flight comfort is noticeably improved, and there is also a reduction in fuel costs.
Vorteilhafterweise ist der Photodetektor ein 2-dimensionaler Photodetektor, der einen Bildverstärker und ein CCD- oder CMOS-Array umfasst. Dadurch ergeben sich geringe Kosten, da der Photodetektor marktüblich und konventionell erhältlich ist.The photodetector is advantageously a 2-dimensional photodetector which comprises an image intensifier and a CCD or CMOS array. This results in low costs, since the photodetector is available commercially and conventionally.
Bevorzugt ist zwischen der Sendeeinrichtung und der Empfangseinrichtung ein 'Transferpfad für einen Teil des Laserstrahls vorgesehen, um den erzeugten Laserstrahl direkt in der Empfangseinrichtung zu registrieren. Beispielsweise ist ein Glasfaserkabel vorgesehen, das den mittels eines Strahlteilers abgezweigten Laserstrahl der Sendeeinrichtung an ein Empfangsteleskop der Empfangseinrichtung koppelt. Dadurch ist es möglich, den ausgesandten Laserstrahl direkt zu registrieren und somit die Übertragungsfunktion der Empfangseinrichtung mit den darin vorhandenen Komponenten zu bestimmen. Weiterhin kann von Zeit zu Zeit eine Kalibrierung der Empfangseinrichtung durchgeführt werden.Preferably a 'transfer path for a part of the laser beam is provided between the transmitting device and the receiving device to the generated laser beam to register directly in the receiving device. For example, a glass fiber cable is provided which couples the laser beam of the transmitting device branched off by means of a beam splitter to a receiving telescope of the receiving device. This makes it possible to register the emitted laser beam directly and thus to determine the transfer function of the receiving device with the components present therein. Furthermore, the receiving device can be calibrated from time to time.
Weiterhin kann der direkt zur Empfangseinrichtung geführte Laserstrahl als zeitliche Referenz verwendet werden, um die Laufzeit des ausgesandten und zurückgestreuten Laserstrahls zu definieren. Dadurch kann die Messentfernung, d.h. die Entfernung des Raumbereichs, in dem die Luftgeschwindigkeit ermittelt wird, festgelegt werden.Furthermore, the laser beam guided directly to the receiving device can be used as a time reference in order to define the transit time of the emitted and backscattered laser beam. This allows the measurement distance, ie the distance of the area in which the air speed is determined can be determined.
Bevorzugt ist das Interferometer ein Fabry-Perot-Interferometer. Vorzugsweise ist das Interferogramm ringförmig und umfasst konzentrische Ringe der Intensitätsverteilung.The interferometer is preferably a Fabry-Perot interferometer. The interferogram is preferably ring-shaped and comprises concentric rings of the intensity distribution.
Vorteilhafterweise erzeugt der Laser gepulste Laserstrahlen im UV-Bereich. Dadurch ist es möglich, die Rückstreuung an Luftmolekülen sehr genau zu messen und die Doppler-Verschiebung zu ermitteln.The laser advantageously generates pulsed laser beams in the UV range. This makes it possible to measure the backscatter on air molecules very precisely and to determine the Doppler shift.
Bevorzugt umfasst das Doppler-Lidar-System Feld-programmierbare-Gate-Arrays zur Berechnung der Referenzmuster. Dadurch kann eine schnelle Berechnung in Echtzeit erfolgen.The Doppler lidar system preferably comprises field-programmable gate arrays for calculating the reference pattern. This enables a quick calculation in real time.
Vorteilhafterweise umfasst die Auswerteeinheit ein Modul zur Bestimmung der Transferfunktion der Bauteile der Empfangseinrichtung.The evaluation unit advantageously comprises a module for determining the transfer function of the components of the receiving device.
Bevorzugt ist das erfindungsgemäße Doppler-Lidar-System mit entsprechenden Komponenten zur Durchführung des erfindungsgemäßen Verfahrens ausgestaltet.The Doppler lidar system according to the invention is preferably designed with corresponding components for carrying out the method according to the invention.
Vorteile und Merkmale, die im Zusammenhang mit dem erfindungsgemäßen Verfahren beschrieben werden, gelten auch für das erfindungsgemäße Doppler- Lidar-System, ebenso wie Vorteile und Merkmale, die für das Doppler-Lidar- System gelten, auch für das erfindungsgemäße Verfahren gelten.Advantages and features that are described in connection with the method according to the invention also apply to the Doppler-Lidar system according to the invention, as well as advantages and features that apply to the Doppler-Lidar system also apply to the method according to the invention.
Im Folgenden wird die Erfindung beispielhaft anhand der Zeichnungen beschrieben, in denenThe invention is described below by way of example with reference to the drawings, in which:
Fig. 1 eine schematische Darstellung eines Doppler-Lidar-Systems gemäß einer bevorzugten Ausführungsform der Erfindung zeigt; Fig. 2 schematisch eine Abbildung eines Interferogramms auf dem Photoempfänger des in Fig. 1 dargestellten Systems zeigt;1 shows a schematic representation of a Doppler lidar system according to a preferred embodiment of the invention; Figure 2 shows schematically an image of an interferogram on the photoreceptor of the system shown in Figure 1;
Fig. 3 zeigt ein typisches 2-dimensionales Interferogramm bei relativ niedriger Photonenzahl; und3 shows a typical 2-dimensional interferogram with a relatively low number of photons; and
Fig. 4 zeigt schematisch eine Windgeschwindigkeitsmessung mit dem erfindungsgemäßen Doppler-Lidar-System.4 schematically shows a wind speed measurement with the Doppler-Lidar system according to the invention.
Fig. 1 zeigt ein Doppler-Lidar-System 10 als bevorzugte Ausführungsform der Erfindung. Ein Laser 11 dient zur Erzeugung eines gepulsten Laserstrahls und ist optisch an ein Sendeteleskop 12 gekoppelt, das dazu dient, einen Laserstrahl in einer vorgegebenen Richtung in die Atmosphäre zu senden. Der Laser 11 und das Sendeteleskop 12 bilden somit eine Sendeeinrichtung zum Aussenden des Laserstrahls. Zum Empfang des in der Atmosphäre zurückgestreuten Laserstrahls dient eine Empfangseinrichtung in Form eines Empfangsteleskops 13. Über eine Lichtleitfaser 14 und eine Filtereinheit 15 wird der empfangene Laserstrahl einem Interferometer 16 zugeführt, das in dieser bevorzugten Ausführungsform ein Fabry-Perot-Interferometer ist.1 shows a Doppler lidar system 10 as a preferred embodiment of the invention. A laser 11 is used to generate a pulsed laser beam and is optically coupled to a transmission telescope 12, which is used to transmit a laser beam into the atmosphere in a predetermined direction. The laser 11 and the transmitting telescope 12 thus form a transmitting device for emitting the laser beam. A receiving device in the form of a receiving telescope 13 is used to receive the laser beam scattered back into the atmosphere. The received laser beam is fed via an optical fiber 14 and a filter unit 15 to an interferometer 16, which in this preferred embodiment is a Fabry-Perot interferometer.
Ein Photodetektor 17 dient zur Bestimmung der Intensitätsverteilung eines Interferogramms, das vom Interferometer 16 aus dem zugeführten Laserstrahl erzeugt wird. Dabei ist der Photodetektor 17 in Bezug auf das Interferometer 16 so angeordnet, dass das z.B. ringförmige Interferogramm direkt auf den Photodetektor 17 abgebildet wird.A photodetector 17 is used to determine the intensity distribution of an interferogram, which is generated by the interferometer 16 from the laser beam supplied. The photodetector 17 is arranged in relation to the interferometer 16 so that e.g. annular interferogram is imaged directly on the photodetector 17.
Der Ausgang des Photodetektors 17 ist elektrisch an eine Auswerteeinheit 18 gekoppelt, die einen Speicher 18a und eine Vergleichseinheit 18b in Form eines Mikroprozessors umfasst. In dem Speicher 18a sind ein oder mehrere Referenzmuster für Interferogramme gespeichert, die für definierte atmosphärische und gegebenenfalls weitere Parameter gelten. Die Vergleichseinheit 17b dient dazu, das auf dem Photoempfänger 17 abgebildete Interferogramm mit den ein oder mehreren Referenzmustern zu vergleichen und aus dem Vergleich die Windgeschwindigkeit zu ermitteln. Zur Anpassung, Berechnung oder Veränderung der Referenzmuster können der Vergleichseinheit 18b zusätzliche Parameter als Eingangsdaten 22 zugeführt werden.The output of the photodetector 17 is electrically coupled to an evaluation unit 18, which comprises a memory 18a and a comparison unit 18b in the form of a microprocessor. One or more reference patterns for interferograms, those for defined ones, are stored in the memory 18a atmospheric and possibly other parameters apply. The comparison unit 17b serves to compare the interferogram imaged on the photoreceiver 17 with the one or more reference patterns and to determine the wind speed from the comparison. Additional parameters can be supplied to the comparison unit 18b as input data 22 in order to adapt, calculate or change the reference pattern.
Zwischen dem Laser 11 und dem Sendeteleskop 12 ist ein Strahlteiler 19 geschaltet, der dazu dient, ein Teil des vom Laser 11 erzeugten Laserstrahls abzuzweigen und über einen Transferpfad 20 in Form einer Lichtleitfaser oder Referenzmessfaser dem Empfangsteleskop 13 zuzuführen. Dadurch gelangt im Betrieb ein kleiner Teil der vom Laser 11 abgegebenen Strahlung über den Transferpfad 20 zum Empfangsteleskop 13 und dient zur Bestimmung des Anfangszeitpunktes der Laufzeit des ausgesandten und in der Atmosphäre in definierter Entfernung zurückgestreuten Laserstrahls. Dadurch wird durch Festlegung des Zeitpunktes des Empfangs die Entfernung definiert, in der die Windgeschwindigkeit gemessen werden soll.A beam splitter 19 is connected between the laser 11 and the transmitting telescope 12 and is used to branch off part of the laser beam generated by the laser 11 and to feed it to the receiving telescope 13 via a transfer path 20 in the form of an optical fiber or reference measuring fiber. As a result, a small part of the radiation emitted by the laser 11 reaches the receiving telescope 13 via the transfer path 20 and is used to determine the starting time of the transit time of the laser beam emitted and backscattered in the atmosphere at a defined distance. This defines the distance at which the wind speed is to be measured by specifying the time of reception.
Weiterhin kann über den direkt von der Sendeeinrichtung zugeführten Laserstrahl die Transferfunktion des optischen Empfangssystems mit seinen verschiedenen Komponenten, wie z.B. Faser 14, Filter 15 und Interferometer 16 bestimmt werden. Dabei wird der direkt zugeführte Laserstrahl dem Interferometer 16 zugeführt und anschließend das dabei entstehende Interferogramm auf dem Photodetektor 17 abgebildet.Furthermore, the transfer function of the optical receiving system with its various components, such as e.g. Fiber 14, filter 15 and interferometer 16 can be determined. The directly supplied laser beam is fed to the interferometer 16 and the resulting interferogram is then imaged on the photodetector 17.
Das Empfangsteleskop 13 dient zur Aufnahme des aus dem Messvolumen, das ein definierter Raumbereich ist, zurückgestreuten Lichts. Über die Faser 14 in Form einer Glasfaser wird dieses Licht mittels einer optischen Linse 21 parallelisiert und zum Filter 15 weitergeleitet. Der Filter 15 dient dazu, den aus der Atmosphäre zurückgestreuten Laserstrahl aus der Hintergrundstrahlung der Sonne herauszufiltern. Im Strahlengang schließt sich nach dem Filter 15 das Interferometer 16 an, das aus dem zugeführten Licht ein Interferogramm bildet. Das Interferometer 16 kann neben einem Fabry-Perot-Interferometer beispielsweise auch ein Fizeau- Interferometer sein.The receiving telescope 13 serves to receive the light scattered back from the measurement volume, which is a defined spatial area. Via the fiber 14 in the form of a glass fiber, this light is parallelized by means of an optical lens 21 and passed on to the filter 15. The filter 15 serves to filter out the laser beam scattered back from the atmosphere from the background radiation of the sun. The filter 15 is followed by the interferometer 16, which forms an interferogram from the light supplied. In addition to a Fabry-Perot interferometer, the interferometer 16 can also be a Fizeau interferometer, for example.
Das vom Interferometer erzeugte Muster besteht beim Fabry-Perot-Interferometer aus konzentrischen Kreisen, wobei die Intensitätsverteilung in dem Muster und der Abstand der Kreise vom Zentrum von der Wellenlänge des zugeführten Lichts, vom Spektrum des Lasers, von der Doppler-Verschiebung, der Doppier- Verbreiterung in der Atmosphäre und der Transferfunktion der Komponenten des Empfangskanals abhängt.In the Fabry-Perot interferometer, the pattern generated by the interferometer consists of concentric circles, the intensity distribution in the pattern and the distance of the circles from the center from the wavelength of the light supplied, from the spectrum of the laser, from the Doppler shift, the double Widening in the atmosphere and the transfer function of the components of the receiving channel depends.
Das erzeugte Interferogramm enthält die gesamten gemessenen Informationen über die Atmosphäre und wird ohne optische Umformung oder Ausblendung bzw. Zentrierung als Ganzes auf die Detektorebene des flächenhaften Photodetektors 17 mit seiner 2-dimensionalen Pixelstruktur abgebildet.The generated interferogram contains the entire measured information about the atmosphere and is mapped as a whole without optical shaping or fading out or centering onto the detector plane of the two-dimensional photodetector 17 with its 2-dimensional pixel structure.
Figur 2 zeigt eine typische von dem Photoempfänger oder dem Photodetektor 17 wiedergebene Intensitätsverteilung in einem Interferogramm 34. Aufgrund der diskreten Pixelstruktur und der typischerweise schwachen optischen Empfangssignale entsteht ein stufenartiges 2-dimensionales Intensitätsbild, das der durch einen Mikroprozessor realisierten Vergleichseinheit 18b zugeführt wird.FIG. 2 shows a typical intensity distribution in an interferogram 34 reproduced by the photoreceptor or the photodetector 17. Because of the discrete pixel structure and the typically weak optical reception signals, a step-like 2-dimensional intensity image is produced which is fed to the comparison unit 18b implemented by a microprocessor.
Fig. 3 dient zur besseren Veranschaulichung des typischen 2-dimensionalen Interferogramms 34 bei einer relativ niedrigen Photonenzahl. Dort ist die Intensitätsverteilung der konzentrischen Ringe in Form von Grauwertabstufungen gezeigt. In der Vergleichseinheit 18b (siehe Fig. 1) wird dieses Intensitätsprofil oder Intensitätsmuster mit Referenzmustern verglichen, die in der Speichereinheit 18a in Form von Look-Up-Tables vorliegen. Bei der Auswertung sucht der Mikroprozessor aus den ein oder mehreren Look-Up-Tables dasjenige Intensitätsverteilungsmuster aus, das dem gemessenen am nächsten kommt. Aus diesem Vergleich ergibt sich der gesuchte Wert der Doppler-Verschiebung, die ein Parameter des gespeicherten Musters ist, und daraus abgeleitet die Geschwindigkeit der Luftteilchen im Messvolumen.3 serves to better illustrate the typical 2-dimensional interferogram 34 with a relatively low number of photons. There the intensity distribution of the concentric rings is shown in the form of gray value gradations. In the comparison unit 18b (see FIG. 1), this intensity profile or intensity pattern is compared with reference patterns which are available in the storage unit 18a in the form of look-up tables. During the evaluation, the microprocessor looks for one from the one or more look-up tables Intensity distribution pattern that comes closest to the measured one. The value of the Doppler shift, which is a parameter of the stored pattern, is derived from this comparison, and the speed of the air particles in the measurement volume is derived from this.
Die Referenzmuster oder die Schar von Referenzmustern können auf verschiedene Arten erzeugt werden. Beispielsweise können für verschiedene Temperaturen der Atmosphäre und verschiedene Luftgeschwindigkeiten relativ zum Messsystem unterschiedliche Referenzmuster erzeugt und in den Look-Up- Tables abgelegt werden. Weiterhin dienen der Druck und dieThe reference patterns or the set of reference patterns can be generated in various ways. For example, different reference patterns can be generated for different atmospheric temperatures and different air velocities relative to the measuring system and stored in the look-up tables. Furthermore, the print and the serve
Instrumentenfunktion als Parameter für die Bestimmung der Referenzmuster. Es ist ebenso möglich, ein oder mehrere dieser Parameter während des Messbetriebs durch Sensoren festzustellen und die gespeicherten Referenzmuster aufgrund der aktuellen Werte anzupassen. Die von den weiteren Sensoren ermittelten aktuellen Parameter können dann in Form vonInstrument function as a parameter for determining the reference pattern. It is also possible to determine one or more of these parameters during the measuring operation by sensors and to adapt the stored reference patterns on the basis of the current values. The current parameters determined by the further sensors can then be in the form of
Eingangsdaten 22 der Vergleichseinheit 18b zugeführt werden (s. Fig. 1).Input data 22 are fed to the comparison unit 18b (see FIG. 1).
D.h., in den Interferenzmustern sind die verschiedenen Parameter, die auf ein gemessenes Interferogramm Einfluss haben, enthalten, so dass die Referenzmuster für verschiedene Höhen, Temperaturen, Geschwindigkeiten des Flugzeugs und auch unterschiedliche Instrumentenfunktionen gelten. Dabei kann die Instrumentenfunktion entweder gemessen oder berechnet sein.This means that the interference patterns contain the various parameters that influence a measured interferogram, so that the reference patterns apply to different heights, temperatures, speeds of the aircraft and also different instrument functions. The instrument function can either be measured or calculated.
Zur Vermessung der Transferfunktion des Empfangssystems wird ein Teil der Strahlung des Lasers direkt in das Empfangsteleskop 13 geleitet. DieTo measure the transfer function of the receiving system, part of the radiation from the laser is directed directly into the receiving telescope 13. The
Transferfunktion kann dann in den Zwischenzeiten zwischen der Vermessung des rückgestreuten Laserlichts der Atmosphäre ermittelt und bei der weiteren Auswertung berücksichtigt werden.The transfer function can then be determined in the meantime between the measurement of the backscattered laser light from the atmosphere and taken into account in the further evaluation.
Bei der Auswertung werden im Wesentlichen die Intensitätsverteilungen der Referenzmuster mit der Intensitätsverteilung des gemessenen interferogramms verglichen. Das Spektrum des aus der Atmosphäre zurückgestreuten Lichts enthält eine Stoßverbreiterung der Moleküle sowie eine thermische Verbreiterung, die von der jeweiligen Temperatur abhängig ist. Durch die Doppler-Verschiebung ist der Schwerpunkt des Spektrums verschoben. Wie oben erwähnt, wirkt sich auch die Instrumentenfunktion, d.h. die Transferfunktion der optischen Komponenten des Empfängers, auf das ermittelte Spektrum aus. In den Referenzmustern sind all diese Faktoren berücksichtigt.During the evaluation, the intensity distributions of the reference patterns are essentially compared with the intensity distribution of the measured interferogram compared. The spectrum of the light scattered back from the atmosphere contains a collision broadening of the molecules as well as a thermal broadening which is dependent on the respective temperature. The center of gravity of the spectrum is shifted by the Doppler shift. As mentioned above, the instrument function, ie the transfer function of the optical components of the receiver, also has an effect on the spectrum determined. All these factors are taken into account in the reference samples.
Die Messung der Doppler-Verschiebung geschieht mit einer bestimmten Messfrequenz, die frei festgelegt werden kann. Zwischen dem Abstrahlen des Sendepulses bis zum Empfang des rückgestreuten Laserlichts aus dem Messvolumen in der Atmosphäre besteht ein Zeitintervall, in dem der Sendepuls, der direkt über den Transferpfad 20 dem Empfangsteleskop 13 zugeführt wird, vom Photodetektor 17 gemessen wird. Die Intensitätsverteilung in diesem Interferogramm ist die Faltung des Spektrums des verwendeten Lasers und der Transferfunktion aller Bauteile des Empfangssystems, wie beispielsweise Empfangsteleskop 13, Glasfaser oder Faser 14, Linse 21, Filter 15 und Interferometer 16. Diese Intensitätsverteilung wird bei der Messung zur Vervollständigung der Look-Up-Tables und/oder zur Überwachung der Einstellungen des gesamten Systems, insbesondere zur Kalibrierung, verwendet.The Doppler shift is measured with a certain measuring frequency, which can be freely determined. Between the emission of the transmission pulse until reception of the backscattered laser light from the measurement volume in the atmosphere, there is a time interval in which the transmission pulse, which is fed directly to the reception telescope 13 via the transfer path 20, is measured by the photodetector 17. The intensity distribution in this interferogram is the convolution of the spectrum of the laser used and the transfer function of all components of the receiving system, such as receiving telescope 13, glass fiber or fiber 14, lens 21, filter 15 and interferometer 16. This intensity distribution is used in the measurement to complete the look -Up tables and / or used to monitor the settings of the entire system, especially for calibration.
Bei der Messung der Windgeschwindigkeit innerhalb eines vom Messsystem beabstandeten Raumbereichs mit dem in Fig. 1 gezeigten Doppler-Lidar-System 0 wird der Laser 11 mit einer typischen Pulsdauer von 10 ns betrieben, was einer Pulslänge von 1 ,5 m entspricht. Die Empfangseinrichtung wird dabei so eingestellt, dass das Rückstreusignal nach einer festen Pulslaufzeit aus einer festen Entfernung, beispielsweise 100 m, empfangen wird. Als Zeitreferenz dient dabei der Zeitpunkt des Durchgangs der gepulsten Laserstrahlung durch die Apertur des Sendeteleskops 12. Dieser Zeitpunkt wird ermittelt, indem ein Teil der Strahlung über den als Glasfaser ausgestalteten Transferpfad 20 direkt demWhen measuring the wind speed within a spatial area spaced from the measuring system with the Doppler-Lidar system 0 shown in FIG. 1, the laser 11 is operated with a typical pulse duration of 10 ns, which corresponds to a pulse length of 1.5 m. The receiving device is set so that the backscatter signal is received after a fixed pulse transit time from a fixed distance, for example 100 m. The point in time of the passage of the pulsed laser radiation through the aperture of the transmitting telescope 12 is used as the time reference. This point in time is determined by directing part of the radiation via the transfer path 20 designed as glass fiber
Empfangsteleskop 13 zugeführt wird, also ohne den Weg durch die Atmosphäre. Außer diesem Zeitpunkt bzw. Trigger, der von dem Photodetektor 17 registriert wird, kann auch die spektrale Intensitätsverteilung des direkt zugeführten Laserstrahls in der Abbildungsebene des Interferometers 16 auf dem Photodetektor 1 aufgezeichnet und abgespeichert werden. Dadurch ergibt sich die Transferfunktion F0 aller optischen Komponenten als :Receiving telescope 13 is supplied, so without the way through the atmosphere. In addition to this point in time or trigger, which is registered by the photodetector 17, the spectral intensity distribution of the directly supplied laser beam can also be recorded and stored in the imaging plane of the interferometer 16 on the photodetector 1. This results in the transfer function F 0 of all optical components as:
Fo = fLa0 fFi0 fFa0 f|F, (1)Fo = f La 0 fFi0 fFa0 f | F, (1)
wobei 0 den Faltungsoperator bezeichnet und f|_a, fn. Fa und f|F diewhere 0 denotes the convolution operator and f | _ a , fn. F a and f | F the
Transferfunktionen der optischen Komponenten Laser, Filter, Faser und Interferometer darstellen.Represent transfer functions of the optical components laser, filter, fiber and interferometer.
Da diese Funktion die inhärente Instrumentenfunktion ist und nur selten, d.h. nur bei Dejustage und Degradation des Systems verändert wird, braucht sie nicht bei jedem Puls des Lasers abgespeichert werden, sondern kann zur anfänglichen Selbstkalibrierung und zur Überwachung des Systems immer wieder neu gemessen werden. Wenn das Messsystem stabil genug ist, besteht auch die Möglichkeit, die Instrumentenfunktion rein rechnerisch zu ermitteln und bei weiteren Auswertungen zu berücksichtigen. Die Auswertung dieser Tests kann bei der Berechnung zur Stützung der Daten verwendet werden.Since this function is the inherent instrument function and only rarely, i.e. changed only when the system is misaligned and degraded, does not need to be saved with every pulse of the laser, but can be measured again and again for initial self-calibration and for monitoring the system. If the measuring system is stable enough, there is also the possibility to determine the instrument function purely by calculation and to take it into account in further evaluations. The evaluation of these tests can be used in the calculation to support the data.
Das zum späteren Zeitpunkt von dem Messvolumen in der Atmosphäre zurückgestreute Licht des Laserpulses wird vom Empfangsteleskop 13 aufgenommen und über denselben Weg wie vorher das direkt zugeführteThe light of the laser pulse scattered back from the measurement volume in the atmosphere at a later point in time is picked up by the receiving telescope 13 and the directly supplied one is carried out in the same way as before
Laserlicht über das Empfangsteleskop 13, die Linse 21 , den Filter 14, die Faser 15 und das Interferometer 16 dem Photoempfänger bzw. Photodetektor 17 zugeführt.Laser light is fed via the receiving telescope 13, the lens 21, the filter 14, the fiber 15 and the interferometer 16 to the photo receiver or photodetector 17.
Das aus der Atmosphäre zurückgestreute Laserlicht ist spektral verbreitert und durch den Doppler-Effekt aus der Relativgeschwindigkeit des Flugzeugs durch die Luft und die Eigenbewegung der Luftmasse spektral verschoben. Nach dem Durchgang durch das Empfangsteleskop 13 wird es noch mit der Transferfunktion der Optik F0 gefaltet und stellt dann die gesamte gemessene Funktion F dar:The laser light scattered back from the atmosphere is broadened spectrally and due to the Doppler effect from the relative speed of the aircraft Air and the natural movement of the air mass shifted spectrally. After passing through the receiving telescope 13, it is folded with the transfer function of the optics F 0 and then represents the entire measured function F:
F = FD 0 Fo (2)F = F D 0 F o (2)
Der gemessene Intensitätsverlauf wird mit dem gespeicherten 2-dimensionaIen Feld von möglichen Intensitätsverläufen in Abhängigkeit der atmosphärischen Parameter Dichte und Temperatur verglichen. Der Intensitätsverlauf mit dem geringsten Unterschied zwischen Messung und Berechnung wird in einer modellbasierten Auswertung verwendet, oder er wird durch eine entsprechende Interpolation oder Parameteranpassung verbessert. Aus diesem Intensitätsverlauf kann nun die Doppler-Funktion FD durch rechnerische Entfaltung gewonnen werden. Aus der Breite und dem Verlauf der Doppler-Funktion können nun Dichte und Temperatur der Luft in entsprechenden Flughöhen zusätzlich gewonnen werden.The measured intensity curve is compared with the stored 2-dimensional field of possible intensity curves depending on the atmospheric parameters density and temperature. The intensity curve with the least difference between measurement and calculation is used in a model-based evaluation, or it is improved by an appropriate interpolation or parameter adjustment. The Doppler function F D can now be obtained from this intensity curve by computational unfolding. From the width and the course of the Doppler function, the density and temperature of the air at corresponding flight altitudes can now be obtained.
Andererseits besteht auch die Möglichkeit, aus den bekannten Werten Luftdichte und Lufttemperatur, die mit anderen Messinstrumenten des Flugzeugs gemessen werden, den erwarteten Verlauf der Doppler-Funktion im Voraus zu berechnen und mit der zusätzlichen Kenntnis der Instrumentenfunktion die Doppler- Verschiebung rechnerisch aus der gemessenen Verschiebung der Interferogramme 34 zu gewinnen. Aus der Doppler-Verschiebung und der Ausrichtung der Messachse aus dem Flugzeug lässt sich dann die Geschwindigkeit des Flugzeugs gegenüber der Luftmasse und die interne Bewegung der Luftmasse berechnen.On the other hand, it is also possible to calculate the expected course of the Doppler function from the known values of air density and air temperature, which are measured with other measuring instruments of the aircraft, and to calculate the Doppler shift from the measured shift with the additional knowledge of the instrument function to win the interferograms 34. The Doppler shift and the alignment of the measuring axis from the aircraft can then be used to calculate the speed of the aircraft in relation to the air mass and the internal movement of the air mass.
Nach der direkten Abbildung des Interferogramms als Ganzes und ohne optische Umformung oder Ausblendung auf dem flächenhaften Detektor, z.B. einem CCD- Array oder CMOS-Array, erfolgt eine Lichtverstärkung unter Verwendung einer Multi-Channel-Plate, die einen Bildverstärker 17a bildet und zwischen der Abbildungsebene des Interferometers 16 und dem flächenartigen Photodetektor 17 vorgesehen ist.After direct imaging of the interferogram as a whole and without optical shaping or masking on the area detector, for example a CCD array or CMOS array, light amplification is carried out using a multi-channel plate, which forms an image intensifier 17a and between Imaging plane of the interferometer 16 and the planar photodetector 17 is provided.
Anschließend w rden die Interferogrammβ mit den Parameterberechnungen und/oder abgespeicherten Informationen über den erwarteten Intensitätsverlauf bei dem verwendeten Interferometer, der Lufttemperatur und Dichte mit Hilfe von parallelrechnenden Mikroprozessoren, wie z.B. Feld-Programmierbare-Gate- Arrays (FPGA), verglichen. Die Referenzwerte der Interferogrammfunktion werden für unterschiedliche Parametersätze, insbesondere Laserwellenlänge λ, Atmosphärendichte p und Lufttemperatur T gespeichert, um den Rechenaufwand zu reduzieren. Dadurch können alle im System selbst vorhandenen Informationen, sowie die in jedem Flugzeug vorhandenen Informationen über Luftdichte, Temperatur und Eigengeschwindigkeit gegenüber der Luftmasse genutzt werden.Then the interferogram with the parameter calculations and / or stored information about the expected intensity curve for the interferometer used, the air temperature and density with the help of parallel computing microprocessors, e.g. Field Programmable Gate Arrays (FPGA) compared. The reference values of the interferogram function are stored for different parameter sets, in particular laser wavelength λ, atmospheric density p and air temperature T, in order to reduce the computational effort. As a result, all information available in the system itself, as well as the information about air density, temperature and airspeed in relation to the air mass, can be used in every aircraft.
Fig. 4 zeigt schematisch eine Windgeschwindigkeitsmessung mit dem erfindungsgemäßen Doppler-Lidar-System an Bord eines Flugzeugs 50. Dabei werden verschiedene Messstrahlen 51 in unterschiedliche Richtungen gesendet, um den Luftraum vor dem Flugzeug 50 abzutasten. Die aus verschiedenen Raumbereichen 52a, 52b, 52c, 52d an Aerosolen und Molekülen zurückgestreuten Laserstrahlen werden wie oben beschrieben aufgezeichnet und analysiert, so dass die Windgeschwindigkeiten in der jeweiligen Messrichtung in dem jeweiligen Raumbereich ermittelt werden. Dabei beträgt die Entfernung D der Raumbereiche 52a - d vom Flugzeug 50 beispielsweise 100 m. Die Aufteilung des Messstrahls in die Vielzahl von Messstrahlen 51 kann z.B. mit einem vorgesetzten Scanner erfolgen, wobei gleichzeitig oder nacheinander in verschiedenen Richtungen der Luftraum vor dem Flugzeug 50 abgetastet wird. Die Raumbereiche 52a - d haben beispielsweise eine Ausdehnung in Messrichtung von ca. 10 m.4 schematically shows a wind speed measurement with the Doppler-Lidar system according to the invention on board an aircraft 50. Different measuring beams 51 are sent in different directions in order to scan the air space in front of the aircraft 50. The laser beams scattered back from different spatial areas 52a, 52b, 52c, 52d on aerosols and molecules are recorded and analyzed as described above, so that the wind speeds in the respective measuring direction are determined in the respective spatial area. The distance D of the spatial areas 52a-d from the aircraft 50 is 100 m, for example. The division of the measuring beam into the plurality of measuring beams 51 can e.g. with a superior scanner, the air space in front of the aircraft 50 being scanned simultaneously or successively in different directions. The spatial areas 52a-d have, for example, an extension in the measuring direction of approximately 10 m.
Durch Projektion in die verschiedenen Richtungen kann aus den einzelnen Messergebnissen der Windgeschwindigkeitsvektor bestimmt werden. Die Verarbeitung der in der jeweiligen Laser-Sendestrahlrichtung gemessenen Doppler-Geschwindigkeitskomponenten zu einem geschätzten Geschwindigkeitsvektor erfolgt beispielsweise durch ein überbestimmtes Gleichungssystem mit mindestens vier Einzelmessungen bei vier Messrichtungen, wobei Wichtungsfaktoren zur Reduzierung des Betrags vergangener Messungen berücksichtigt werden. Der Geschwindigkeitsberechnung von vn aus derThe wind speed vector can be determined from the individual measurement results by projection in the different directions. The The Doppler speed components measured in the respective laser transmission beam direction are processed into an estimated speed vector, for example, by an over-determined system of equations with at least four individual measurements in four measuring directions, weighting factors for reducing the amount of past measurements being taken into account. The speed calculation of v n from the
Messfolge m-t, m2, m3, m mk folgt dann die Geschwindigkeitsberechnung vn+1 aus der Messfolge m2, m3, m4, m5 rrik+i.Measurement sequence m- t , m 2 , m 3 , mm k then the speed calculation v n + 1 follows from the measurement sequence m 2 , m 3 , m 4 , m 5 rrik + i.
Die Erfindung ermöglicht eine genaue Messung des Windgeschwindigkeitsvektors weit vor einem Flugzeug in Echtzeit mit Hilfe der interferometrischen Messung der Doppler-Linienverschiebung von an Luftmolekülen und Aerosolen rückgestreuten gepulsten Laserstrahlen. Aufgrund der erfindungsgemäßen Anordnung und dem erfindungsgemäßen Auswerteverfahren wird eine große Reichweite erreicht, d.h. bei der Auswertung werden nur verringerte Rückstreuintensitäten benötigt.The invention enables an accurate measurement of the wind speed vector far in front of an aircraft in real time with the help of the interferometric measurement of the Doppler line shift of pulsed laser beams scattered back on air molecules and aerosols. Because of the arrangement according to the invention and the evaluation method according to the invention, a large range is achieved, i.e. only reduced backscatter intensities are required for the evaluation.
Die Böenvermessung in Echtzeit weit vor dem Flugzeug, beispielsweise in einer Entfernung, die 0,5 Sek. Flugzeit entspricht, erlaubt gezielte Steuerungs- und Trimmmanöver zur Reduktion schädlicher Böeneinwirkungen auf das Flugzeug. Für Systeme mit kürzerer Reichweite, beispielsweise zur Messung desThe gust measurement in real time far in front of the aircraft, for example at a distance that corresponds to a flight time of 0.5 seconds, allows targeted control and trimming maneuvers to reduce the harmful effects of gusts on the aircraft. For systems with a shorter range, for example for measuring the
Anströmgeschwindigkeitsvektors, sind noch weiter reduzierte Einbaumaß möglich.Inlet velocity vector, further reduced installation dimensions are possible.
Die Messung kann von bewegten Messplattformen, beispielsweise von Flugzeugen, Hubschraubern, Ballonen, Satelliten oder Schiffen erfolgen, oder auch von bodenfesten Messstationen aus durchgeführt werden.The measurement can be carried out from moving measurement platforms, for example from airplanes, helicopters, balloons, satellites or ships, or also from floor-mounted measurement stations.
Zur Einstellung der Messentfernung wird bei der Messdatenaufnahme ein Zeitfenster festgelegt, das mit der Laser-Strahlpulsung synchronisiert wird.In order to set the measurement distance, a time window is set for the measurement data acquisition, which is synchronized with the laser beam pulsing.
Die Messdatenverarbeitung kann in Echtzeit erfolgen, wobei die Referenzwerte bzw. -muster durch Verwendung von Feld-Programmierbaren-Arrays parallel berechnet werden. Daraus können die optimalen Parameter für die Wellenlänge λ und die Temperatur T des Messvolumens bestimmt werden. Zur Reduzierung von Rechenaufwand können die Referenzmuster, die Referenzwerte der Interferogrammfunktion darstellen, für unterschiedliche Parametersätze λ und T gespeichert werden. Neben den Speicherzugriffen werden nur Addition und Betragsbildung für positive, ganze Zahlen zur Berechnung der Referen∑werte benötigt, wodurch der Flächenbedarf der lokalen Rechnereinheiten im Feld- Programmierbaren-Gate-Array sehr klein gehalten wird.The measurement data can be processed in real time, the reference values or patterns being used in parallel by using field-programmable arrays be calculated. The optimal parameters for the wavelength λ and the temperature T of the measurement volume can be determined from this. To reduce computational effort, the reference patterns, which represent reference values of the interferogram function, can be stored for different parameter sets λ and T. In addition to the memory accesses, only addition and formation of amounts for positive, whole numbers are required for calculating the reference values, as a result of which the area requirements of the local computer units in the field programmable gate array are kept very small.
Bei dem Vergleich des Interferogramms der aus der Atmosphäre zurückgestreuten Laserstrahlung mit der Schar von Referenzmustern wird während der Auswertung eine Bestimmung der Wellenlänge λ und der Messvolumentemperatur T durchgeführt. Zur Parameterberechnung wird die vollständige, theoretisch zu erwartende Intensitätsverteilung genutzt, die vorher aus einer Modellrechnung gewonnen wird. Dabei schließt die Modellrechnung die Modellierung des Laserstrahls, die Modellierung der Molekülgeschwindigkeitsverteilung und die Modellierung der Interferometertransferfunktion ein.When the interferogram of the laser radiation backscattered from the atmosphere is compared with the host of reference patterns, the wavelength λ and the measurement volume temperature T are determined during the evaluation. The complete, theoretically expected intensity distribution, which is previously obtained from a model calculation, is used to calculate the parameters. The model calculation includes the modeling of the laser beam, the modeling of the molecular velocity distribution and the modeling of the interferometer transfer function.
Neben der Detektion von Böen oder Luftturbulenzen kann das erfindungsgemäße Doppler-Lidar und das erfindungsgemäße Verfahren zur Erfassung von Windgeschwindigkeiten auch für Anwendungen und Einsätze am Boden verwendet werden, beispielsweise zur Detektion von Scherwinden an Flughäfen oder zur Messung von Windprofilen in vertikaler Richtung durch die Atmosphäre. In addition to the detection of gusts or air turbulence, the Doppler lidar according to the invention and the method according to the invention for detecting wind speeds can also be used for applications and operations on the ground, for example for the detection of shear winds at airports or for the measurement of wind profiles in the vertical direction through the atmosphere.

Claims

Patentansprüche claims
1. Verfahren zur Erfassung von Windgeschwindigkeiten mit einem Doppler-Lidar-System (10), bei dem ein Laserstrahl vorgegebener Wellenlänge zu einem Raumbereich (52a, 52b, 52c, 52d) hin ausgesandt wird und das aus dem Raumbereich zurückgestreute Licht empfangen wird, wobei zur Bestimmung einer Dopplerverschiebung mit einem Interferometer (16) ein Interferogramm (34) erzeugt wird und mit einem Photodetektor (17) die Intensitätsverteilung des Interferogramms (34) gemessen wird, dadurch gekennzeichnet, dass die Intensitätsverteilung mit einer Schar von Referenzmustern verglichen wird, die zuvor für definierte atmosphärische Parameter bestimmt wurden, die unterschiedliche Dichten und/oder Temperaturen der Atmosphäre umfassen, wobei aus dem Vergleich mit der Schar von unterschiedlichen Referenzmustern die Dopplerverschiebung als Maß für die Windgeschwindigkeit ermittelt wird.1. A method for detecting wind speeds with a Doppler lidar system (10), in which a laser beam of a predetermined wavelength is emitted to a spatial area (52a, 52b, 52c, 52d) and the light scattered back from the spatial area is received, whereby for determining a Doppler shift, an interferogram (34) is generated with an interferometer (16) and the intensity distribution of the interferogram (34) is measured with a photodetector (17), characterized in that the intensity distribution is compared with a family of reference patterns that were previously were determined for defined atmospheric parameters, which include different densities and / or temperatures of the atmosphere, the Doppler shift being determined as a measure of the wind speed from the comparison with the family of different reference patterns.
2. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass das Interferogramm (34) ringförmig ist und direkt auf den 2-dimensionalen2. The method according to claim 1, characterized in that the interferogram (34) is annular and directly on the 2-dimensional
Photodetektor (17) abgebildet wird.Photodetector (17) is imaged.
3. Verfahren nach Anspruch 1 , dadurch gekennzeichnet, dass das Interferogramm streifenförmig ist und direkt auf den 2-dimensionalen Photodetektor abgebildet wird.3. The method according to claim 1, characterized in that the interferogram is strip-shaped and is imaged directly on the 2-dimensional photodetector.
4. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass das Referenzmuster mit der geringsten Abweichung zum gemessenen Interferogramm (34) zur Bestimmung der Dopplerverschiebung verwendet wird. 4. The method according to any one of the preceding claims, characterized in that the reference pattern with the smallest deviation from the measured interferogram (34) is used to determine the Doppler shift.
Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass das Referenzmuster die Geschwindigkeit der Atmosphäre relativ zum Doppler-Lidar-Systern (10) als Parameter enthält. Method according to one of the preceding claims, characterized in that the reference pattern contains the speed of the atmosphere relative to the Doppler-Lidar system (10) as a parameter.
Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass aus mehreren aufeinanderfolgenden Messungen die Variation der Geschwindigkeit der Atmosphäre relativ zum Doppler- Lidar-System (10) ermittelt wird.Method according to one of the preceding claims, characterized in that the variation in the speed of the atmosphere relative to the Doppler-Lidar system (10) is determined from a plurality of successive measurements.
7. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass der Laserstrahl gepulst ist und jeweils ein Teil eines Laserpulses zur Festlegung eines zeitlichen Referenzpunktes verwendet wird, um mittels der Laufzeit des restlichen Teils des Laserpulses die Entfernung des rückstreuenden Raumbereichs (52a,7. The method according to any one of the preceding claims, characterized in that the laser beam is pulsed and in each case part of a laser pulse is used to determine a temporal reference point in order to use the transit time of the remaining part of the laser pulse to remove the distance from the backscattering spatial region (52a,
52b, 52c, 52d) festzulegen.52b, 52c, 52d).
8. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass ein Teil des Laserstrahls direkt und ohne Rückstreuung empfangen und registriert wird, wobei aus der8. The method according to any one of the preceding claims, characterized in that a part of the laser beam is received and registered directly and without backscattering, from which
Intensitätsverteilung eine Transferfunktion optischer Komponenten bestimmt und/oder ein Kalibrierung durchgeführt wird.Intensity distribution determines a transfer function of optical components and / or a calibration is carried out.
9. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass aus dem Referenzmuster mit der geringsten9. The method according to any one of the preceding claims, characterized in that from the reference pattern with the lowest
Abweichung zum gemessenen Interferogramm (34) Dichte und/oder Temperatur des Raumbereichs (52a, 52b, 52c, 52d) ermittelt werden.Deviation from the measured interferogram (34) density and / or temperature of the spatial area (52a, 52b, 52c, 52d) can be determined.
10. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass es an Bord eines bewegten Systems, bevorzugt an Bord eines Luftfahrzeugs (50), durchgeführt wird.10. The method according to any one of the preceding claims, characterized in that it is preferred on board a moving system on board an aircraft (50).
11. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass aus gemessenen atmosphärischen Parametern und/oder Flugparametern eines Luftfahrzeugs (50) die erwartete11. The method according to any one of the preceding claims, characterized in that the expected from measured atmospheric parameters and / or flight parameters of an aircraft (50)
Intensitätsverteilung des Referenzmusters berechnet wird.Intensity distribution of the reference pattern is calculated.
12. Verfahren nach einem der vorhergehenden Ansprüche, dadurch gekennzeichnet, dass der Laserstrahl in verschiedene Richtungen gesandt wird, um durch Messung der Dopplerverschiebung in diesen12. The method according to any one of the preceding claims, characterized in that the laser beam is sent in different directions by measuring the Doppler shift in this
Richtungen den Windgeschwindigkeitsvektor zu bestimmen.Directions to determine the wind speed vector.
13. Doppler-Lidar-System (10) zur Erfassung von Windgeschwindigkeiten, insbesondere an Bord von Flugzeugen (50), mit einer Sendeeinrichtung (12) zum Aussenden eines Laserstrahls, einer Empfangseinrichtung (13) zum Empfang des in der Atmosphäre zurückgestreuten Laserstrahls, einem Interferometer (16) zur Erzeugung eines Interferogramms (34) aus dem zurückgestreuten Laserstrahl, einem Photodetektor (17) zur Bestimmung der Intensitätsverteilung des13. Doppler lidar system (10) for detecting wind speeds, in particular on board aircraft (50), with a transmitting device (12) for emitting a laser beam, a receiving device (13) for receiving the laser beam scattered back into the atmosphere, one Interferometer (16) for generating an interferogram (34) from the backscattered laser beam, a photodetector (17) for determining the intensity distribution of the
Interferogramms (34), wobei das Interferogramm (34) direkt auf den Photodetektor (17) abgebildet wird, und einer Auswerteeinheit (18a, 18b) zur Bestimmung der Dopplerverschiebung als Maß für die Windgeschwindigkeit der Atmosphäre, dadurch gekennzeichnet, dass die Auswerteeinheit (18a, 18b) einen Speicher (18a) mit einer Schar von Referenzmustern umfasst, die für zuvor definierte atmosphärische Parameter gelten, welche unterschiedliche Dichten und/oder Temperaturen der Atmosphäre umfassen, und eine Vergleichseinheit (18b) vorgesehen ist, die aus einem Vergleich des abgebildeten Interferogramms (34) mit der Schar von Referenzmustern die Windgeschwindigkeit ermittelt.Interferogram (34), the interferogram (34) being imaged directly on the photodetector (17), and an evaluation unit (18a, 18b) for determining the Doppler shift as a measure of the wind speed of the atmosphere, characterized in that the evaluation unit (18a, 18b) comprises a memory (18a) with a family of reference patterns which apply to previously defined atmospheric parameters, which comprise different densities and / or temperatures of the atmosphere, and a comparison unit (18b) is provided which is based on a comparison of the interferogram (34) shown, the wind speed is determined with the family of reference patterns.
14. Doppler-Lidar-System nach Anspruch 13 , dadurch gekennzeichnet, dass der Photodetektor (17) ein 2-dimensionaler Photodetektor ist, der einen Bildverstärker (17a) und ein CCD- oder CMOS-Array umfasst.14. Doppler lidar system according to claim 13, characterized in that the photodetector (17) is a 2-dimensional photodetector which comprises an image intensifier (17a) and a CCD or CMOS array.
15. Doppler-Lidar-System nach Anspruch 13 oder 14 , dadurch gekennzeichnet, dass zwischen der Sendeeinrichtung (12) und der Empfangseinrichtung (13) ein Transferpfad (20) für einen Teil des15. Doppler lidar system according to claim 13 or 14, characterized in that between the transmitting device (12) and the receiving device (13) a transfer path (20) for part of the
Laserstrahls vorgesehen ist, um den erzeugten Laserstrahl direkt in der Empfangseinrichtung (13) zu registrieren.Laser beam is provided to register the laser beam generated directly in the receiving device (13).
16. Doppler-Lidar-System nach einem der Ansprüche 13 bis 15, dadurch gekennzeichnet, dass das Interferometer (34) ein Fabry-Perot- Interferometer ist, das ringförmige Interferenzmuster erzeugt.16. Doppler-lidar system according to one of claims 13 to 15, characterized in that the interferometer (34) is a Fabry-Perot interferometer that generates ring-shaped interference patterns.
17. Doppler-Lidar-System nach einem der Ansprüche 13 bis 16, dadurch gekennzeichnet, dass das Interferometer (34) ein Fizeau-Interferometer ist, das streifenförmige Interferenzmuster erzeugt.17. Doppler lidar system according to one of claims 13 to 16, characterized in that the interferometer (34) is a Fizeau interferometer which generates stripe-shaped interference patterns.
18. Doppler-Lidar-System nach einem der Ansprüche 13 bis 17, gekennzeichnet durch einen Laser (11), der gepulste Laserstrahlen im UV-Bereich erzeugt.18. Doppler lidar system according to one of claims 13 to 17, characterized by a laser (11) which generates pulsed laser beams in the UV range.
19. Doppler-Lidar-System nach einem der Ansprüche 13 bis 18, gekennzeichnet durch Feld-programmierbare Gate-Arrays zur Berechnung der Referenzmuster.19. Doppler lidar system according to one of claims 13 to 18, characterized by field-programmable gate arrays for calculating the reference pattern.
20. Doppler-Lidar-System nach einem der Ansprüche 13 bis 19 , dadurch gekennzeichnet, dass die Auswerteeinheit (18a, 18b) ein Modul zur Bestimmung der Transferfunktion von Bauteilen auf der Empfangsseite des Doppler-Lidar-Systems (10) umfasst. 20. Doppler lidar system according to one of claims 13 to 19, characterized in that the evaluation unit (18a, 18b) is a module for Determination of the transfer function of components on the receiving side of the Doppler lidar system (10) includes.
EP04725596A 2003-04-10 2004-04-03 Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system Expired - Lifetime EP1613980B1 (en)

Applications Claiming Priority (2)

Application Number Priority Date Filing Date Title
DE10316762A DE10316762B4 (en) 2003-04-10 2003-04-10 Method for detecting wind speeds with a Doppler-Lidar system, in particular on board aircraft, and Doppler Lidar system
PCT/DE2004/000694 WO2004092767A1 (en) 2003-04-10 2004-04-03 Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system

Publications (2)

Publication Number Publication Date
EP1613980A1 true EP1613980A1 (en) 2006-01-11
EP1613980B1 EP1613980B1 (en) 2008-01-23

Family

ID=33154165

Family Applications (1)

Application Number Title Priority Date Filing Date
EP04725596A Expired - Lifetime EP1613980B1 (en) 2003-04-10 2004-04-03 Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system

Country Status (6)

Country Link
US (1) US7463341B2 (en)
EP (1) EP1613980B1 (en)
AT (1) ATE384962T1 (en)
DE (2) DE10316762B4 (en)
ES (1) ES2298741T3 (en)
WO (1) WO2004092767A1 (en)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106546766A (en) * 2016-11-02 2017-03-29 中国海洋大学 Clinoplain scan method based on two anemometry laser radars
CN111077338A (en) * 2019-12-29 2020-04-28 中国科学院西安光学精密机械研究所 High-time-resolution medium-high atmosphere wind field measurement interferometer system for foundation

Families Citing this family (58)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US7495774B2 (en) 2002-03-01 2009-02-24 Michigan Aerospace Corporation Optical air data system
FR2870942B1 (en) 2004-05-25 2006-08-25 Airbus France Sas ANTICIPATED MEASUREMENT SYSTEM FOR TURBULENCE UPSTREAM OF AN AIRCRAFT
CN100360957C (en) * 2005-03-10 2008-01-09 中国科学院合肥物质科学研究院 Unattended full time operating portable polarization-metre laser radar and detecting method thereof
DE102005034729B3 (en) * 2005-07-21 2007-02-08 Eads Deutschland Gmbh Method and lidar system for measuring air turbulence on board aircraft, airports and wind farms
EP1814006B1 (en) * 2006-01-25 2016-09-21 Airbus Opérations SAS Minimizing dynamic structural loads of an aircraft
DE102008005129B4 (en) 2007-11-09 2017-11-23 Deutsches Zentrum für Luft- und Raumfahrt e.V. Non-linear optical frequency converter, uses thereof and method of generating pulsed tunable laser radiation
FR2926894B1 (en) * 2008-01-25 2010-03-12 Thales Sa METHOD OF ESTIMATING AT ANY POINT OF AN AIRCRAFT AIR DATA PATH
US9733392B2 (en) * 2008-06-27 2017-08-15 Deep Sciences, LLC Methods of using environmental conditions in sports applications
US8864606B2 (en) * 2008-06-27 2014-10-21 The Invention Science Fund I, Llc Sports applications for wind profile systems
US20090326894A1 (en) * 2008-06-27 2009-12-31 Chan Alistair K Methods of processing wind profile information in sports applications
DE102008031682A1 (en) * 2008-07-04 2010-03-11 Eads Deutschland Gmbh Direct Receive Doppler LIDAR Method and Direct Receive Doppler LIDAR Device
DE102008031681A1 (en) * 2008-07-04 2010-01-14 Eads Deutschland Gmbh LIDAR method for measuring velocities and LIDAR device with timed detection
FR2938075B1 (en) * 2008-11-05 2015-09-25 Airbus France DEVICE AND METHOD FOR DETECTING AND MEASURING WIND FOR AIRCRAFT
WO2010052385A1 (en) * 2008-11-05 2010-05-14 Airburs Operations Device and method for detecting and measuring wind for an aircraft
DE102008056869B3 (en) 2008-11-12 2009-11-12 Eads Deutschland Gmbh Method for image processing of intensity images having annular structure for determining position, involves defining sectors of intensity images with sector vertex on sector origin point which lies in annular structure
US20100195089A1 (en) * 2009-01-30 2010-08-05 General Electric Company Wind velocity measurement system and method
US8797550B2 (en) 2009-04-21 2014-08-05 Michigan Aerospace Corporation Atmospheric measurement system
EP2430392B1 (en) * 2009-05-15 2015-07-22 Michigan Aerospace Corporation Range imaging lidar
US8144325B2 (en) * 2009-07-23 2012-03-27 Rosemount Aerospace, Inc. In-flight multiple field of view detector for supercooled airborne water droplets
DE102009039016B4 (en) 2009-08-28 2012-05-03 Deutsches Zentrum für Luft- und Raumfahrt e.V. flow investigation
FR2951275B1 (en) * 2009-10-09 2012-11-02 Epsiline DEVICE FOR MEASURING WIND SPEED
US8390791B2 (en) * 2009-11-30 2013-03-05 General Electric Company Light detection and ranging system
US8355120B2 (en) * 2010-03-10 2013-01-15 The United States Of America, As Represented By The Secretary Of The Navy Doppler asymmetric spatial heterodyne spectroscopy light detection and ranging receiver
US8842995B2 (en) 2010-05-11 2014-09-23 The Invention Science Fund I, Llc Optical power transmission systems and methods
FR2962553B1 (en) * 2010-07-08 2012-08-03 Centre Nat Etd Spatiales LASER REMOTE SENSING DEVICE AND INTERFEROMETRY METHOD
US8542999B2 (en) * 2011-02-01 2013-09-24 Vello Systems, Inc. Minimizing bandwidth narrowing penalties in a wavelength selective switch optical network
US8913124B2 (en) * 2011-02-03 2014-12-16 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Lock-in imaging system for detecting disturbances in fluid
DE102011010411B4 (en) 2011-02-04 2013-02-21 Deutsches Zentrum für Luft- und Raumfahrt e.V. Laser system and control of a laser system
RU2468387C1 (en) * 2011-05-06 2012-11-27 Государственное Учреждение "Арктический и антарктический научно-исследовательский институт" Device for measuring vertical component of wind velocity for detecting wind shift
US9026278B2 (en) 2011-12-23 2015-05-05 Optical Air Data Systems, Llc LDV system for measuring wind at high altitude
DE102012103755A1 (en) 2012-04-27 2013-10-31 Centrotherm Sitec Gmbh Synthesizing trichlorosilane, comprises converting metallic silicon and hydrogen chloride in reactor using aluminum as a catalyst, recycling aluminum chloride in the reactor leaked out from reactor, and reusing recycled aluminum chloride
US20130311013A1 (en) * 2012-05-16 2013-11-21 Optical Air Data Systems, Llc Measurement Assisted Aerodynamic State Estimator
US9546954B2 (en) 2014-09-18 2017-01-17 Vision Engineering Solutions, LLC Atmosphere profiling systems
WO2016057171A1 (en) * 2014-09-18 2016-04-14 Vision Engineering Solutions, LLC Atmosphere profiling systems
IL237717A (en) 2015-03-12 2017-04-30 Elbit Systems Ltd Detection, characterization and presentation of adverse airborne phenomena
WO2016181493A1 (en) * 2015-05-12 2016-11-17 三菱電機株式会社 Laser radar device and wind speed observation method
US10444367B2 (en) * 2016-02-26 2019-10-15 Honeywell International Inc. Enhanced LiDAR air data using supplementary sensor outputs
US10451518B2 (en) * 2016-05-10 2019-10-22 Rd2, Llc All fiber temperature and air density sensor
KR102657118B1 (en) 2016-10-26 2024-04-15 엘지이노텍 주식회사 Apparatus for Light Detection and Ranging
DE102017115710A1 (en) 2017-07-12 2019-02-07 Airbus Defence and Space GmbH LIDAR arrangement and lidar method
CN108663537B (en) * 2018-04-02 2021-02-09 福建省新能海上风电研发中心有限公司 Offshore wind measuring method and system
US10613229B2 (en) 2018-08-24 2020-04-07 Ball Aerospace & Technologies Corp. Compact quadrature mach-zehnder interferometer
US20210063429A1 (en) * 2018-08-28 2021-03-04 Ball Aerospace & Technologies Corp. Optical Wind Lidar-Based Multifunctional Instrument for Enhanced Measurements and Prediction of Clear Air Turbulence and Other Wind-Based Aviation Related Phenomena
US11131685B2 (en) 2019-05-15 2021-09-28 The Boeing Company Air data system using magnetically induced voltage
CN112711031B (en) * 2019-07-08 2023-10-10 中国科学院合肥物质科学研究院 Improved quasi-blind area-free Doppler coherent laser radar wind speed measurement system and method
US11556000B1 (en) 2019-08-22 2023-01-17 Red Creamery Llc Distally-actuated scanning mirror
US11120699B2 (en) * 2019-09-27 2021-09-14 Lockheed Martin Corporation Predictive wind gust alleviation using optical velocity sensing
US11332260B2 (en) 2019-11-18 2022-05-17 The Boeing Company Electrode-arc sensor air data system for an aircraft
CN110927702B (en) * 2019-11-22 2022-11-01 成都凯天电子股份有限公司 Laser wind measuring radar flow field environment testing system
DE102020003570A1 (en) 2020-06-16 2020-08-06 FEV Group GmbH Assistance system for a vehicle for determining a speed of an air flow flowing transversely to the direction of travel
CN112698357A (en) * 2020-12-16 2021-04-23 中国科学院合肥物质科学研究院 Dual-beam laser radar wind field detection device and method
DE102021002239A1 (en) 2021-04-28 2022-11-03 Leonardo Germany Gmbh Doppler lidar for detecting wind and/or vortex situations
CN113340560B (en) * 2021-07-05 2022-05-17 中国空气动力研究与发展中心低速空气动力研究所 Doppler interferometer calibration and debugging system, calibration method and debugging method
CN113447999B (en) * 2021-07-08 2022-05-13 天津大学 Atmospheric parameter measuring method and device based on laser data
US20240183876A1 (en) * 2022-12-02 2024-06-06 The Boeing Company Aircraft Laser System with a Corrective Optical Device
CN115902834B (en) * 2022-12-02 2024-05-17 中国科学技术大学 Helium laser radar double-frequency temperature and wind measuring system and method based on Fizeau interferometer
DE102023110581B3 (en) 2023-04-25 2024-07-04 Deutsches Zentrum für Luft- und Raumfahrt e.V. Optical measuring system for recording environmental parameters in an environment
CN117092662B (en) * 2023-10-17 2024-02-23 中国科学技术大学 Quantum interference laser radar system and method for wind field detection

Family Cites Families (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4195931A (en) * 1978-05-18 1980-04-01 The United States Of America As Represented By The Secretary Of The Army Clear air turbulence detector
FR2648915B1 (en) * 1989-06-26 1991-09-27 Centre Nat Rech Scient MEDIUM ALTITUDE WIND SPEED MEASURING DEVICE
DE4013702C2 (en) * 1990-04-28 1996-01-11 Christian Dipl Phys Heupts Process for detecting turbulence in the atmosphere and device for carrying out the process
FR2663752B1 (en) * 1990-06-25 1993-01-22 Seso METEOROLOGICAL PARAMETER MEASURING DEVICE.
US5216477A (en) * 1991-05-20 1993-06-01 Korb Charles L Edge technique for measurement of laser frequency shifts including the doppler shift
US6313908B1 (en) * 1998-07-02 2001-11-06 The United States Of America As Represented By The Administrator Of The National Aeronautics And Space Administration Apparatus and method using a holographic optical element for converting a spectral distribution to image points
CA2494458C (en) * 2002-08-02 2009-06-30 Ophir Corporation Optical air data measurement systems and methods

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
See references of WO2004092767A1 *

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN106546766A (en) * 2016-11-02 2017-03-29 中国海洋大学 Clinoplain scan method based on two anemometry laser radars
CN106546766B (en) * 2016-11-02 2018-12-25 中国海洋大学 Clinoplain scan method based on two anemometry laser radars
CN111077338A (en) * 2019-12-29 2020-04-28 中国科学院西安光学精密机械研究所 High-time-resolution medium-high atmosphere wind field measurement interferometer system for foundation
CN111077338B (en) * 2019-12-29 2021-04-20 中国科学院西安光学精密机械研究所 High-time-resolution medium-high atmosphere wind field measurement interferometer system for foundation

Also Published As

Publication number Publication date
ES2298741T3 (en) 2008-05-16
EP1613980B1 (en) 2008-01-23
WO2004092767A1 (en) 2004-10-28
US7463341B2 (en) 2008-12-09
US20070171397A1 (en) 2007-07-26
DE10316762B4 (en) 2007-01-25
ATE384962T1 (en) 2008-02-15
DE10316762A1 (en) 2004-11-11
DE502004006046D1 (en) 2008-03-13

Similar Documents

Publication Publication Date Title
EP1613980B1 (en) Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system
EP1910868B1 (en) Method and lidar system for measuring air turbulences on board aircraft and for airports and wind farms
US9977045B2 (en) Atmospheric measurement system
DE60313616T2 (en) SYSTEM AND METHOD FOR THE OPTICAL MEASUREMENT OF AIR DATA
EP2297594B1 (en) Lidar method for measuring speeds and lidar device with time-controlled detection
EP0204295B1 (en) Measuring device for determining wind direction and wind speed in the atmosphere
Eloranta et al. The determination of wind speeds in the boundary layer by monostatic lidar
DE102008031682A1 (en) Direct Receive Doppler LIDAR Method and Direct Receive Doppler LIDAR Device
Schmitt et al. The AWIATOR airborne LIDAR turbulence sensor
EP2955545B1 (en) Atmospheric measurement system
US6813020B2 (en) Device for determining the values of at least one parameter of particles, especially of water droplets
EP0009533B1 (en) Doppler velocimeter for aircraft
EP4083660A1 (en) Doppler lidar for detecting wind and / or vortex flows
DE102014216368A1 (en) MULTICOLOR LIDAR SYSTEM
DE10359868B3 (en) Detecting condensation trails forming behind aircraft involves comparing parts of camera-acquired image that should show condensation trail according to imaging geometry with image parts that cannot contain trail using image processing
DE2607169A1 (en) Air pollution concentration profile measuring system - uses several radiation beams emitted from aircraft directed onto earth surface and reflected beam evaluation
CN110261874A (en) Real-time clear-air bumpiness detection method and system based on coherent laser
WO2019234047A1 (en) Method and assembly for calibrating flying-object-mounted sensors
Schmitt Research results, lessons learned and future perspective of forward-looking LIDAR for aircraft
Figueroa Cotorogea Development and Test of a Automatic Gain Controller for a Lidar-Detector Based on Airborne Measurements of Cirrus Clouds
DE3316600A1 (en) Eye-safe slant visual range meter
DE3314807C2 (en) Overflight and trail sensor
Schmitt et al. Forward looking clear air turbulence measurement with the AWIATOR LIDAR sensor
Series–Volume Optical Air Flow Measurements in Flight

Legal Events

Date Code Title Description
PUAI Public reference made under article 153(3) epc to a published international application that has entered the european phase

Free format text: ORIGINAL CODE: 0009012

17P Request for examination filed

Effective date: 20051105

AK Designated contracting states

Kind code of ref document: A1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PL PT RO SE SI SK TR

AX Request for extension of the european patent

Extension state: AL HR LT LV MK

DAX Request for extension of the european patent (deleted)
RIN1 Information on inventor provided before grant (corrected)

Inventor name: SOBOTTA, GERALD

Inventor name: HALLDORSSON, THORSTEINN

Inventor name: HOFMANN, MICHAEL

GRAP Despatch of communication of intention to grant a patent

Free format text: ORIGINAL CODE: EPIDOSNIGR1

GRAS Grant fee paid

Free format text: ORIGINAL CODE: EPIDOSNIGR3

GRAA (expected) grant

Free format text: ORIGINAL CODE: 0009210

AK Designated contracting states

Kind code of ref document: B1

Designated state(s): AT BE BG CH CY CZ DE DK EE ES FI FR GB GR HU IE IT LI LU MC NL PL PT RO SE SI SK TR

REG Reference to a national code

Ref country code: GB

Ref legal event code: FG4D

Free format text: NOT ENGLISH

REG Reference to a national code

Ref country code: CH

Ref legal event code: EP

REG Reference to a national code

Ref country code: IE

Ref legal event code: FG4D

Free format text: LANGUAGE OF EP DOCUMENT: GERMAN

REF Corresponds to:

Ref document number: 502004006046

Country of ref document: DE

Date of ref document: 20080313

Kind code of ref document: P

REG Reference to a national code

Ref country code: SE

Ref legal event code: TRGR

GBT Gb: translation of ep patent filed (gb section 77(6)(a)/1977)

Effective date: 20080420

REG Reference to a national code

Ref country code: ES

Ref legal event code: FG2A

Ref document number: 2298741

Country of ref document: ES

Kind code of ref document: T3

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: FI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BG

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080423

ET Fr: translation filed
PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: PT

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080623

Ref country code: PL

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

Ref country code: SI

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

BERE Be: lapsed

Owner name: EADS DEUTSCHLAND G.M.B.H.

Effective date: 20080430

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CZ

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

Ref country code: DK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

Ref country code: SK

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: RO

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

Ref country code: MC

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20080430

PLBE No opposition filed within time limit

Free format text: ORIGINAL CODE: 0009261

STAA Information on the status of an ep patent application or granted ep patent

Free format text: STATUS: NO OPPOSITION FILED WITHIN TIME LIMIT

26N No opposition filed

Effective date: 20081024

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: EE

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: BE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20080430

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: CY

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LU

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20080403

Ref country code: HU

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080724

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: TR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080123

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GR

Free format text: LAPSE BECAUSE OF FAILURE TO SUBMIT A TRANSLATION OF THE DESCRIPTION OR TO PAY THE FEE WITHIN THE PRESCRIBED TIME-LIMIT

Effective date: 20080424

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: ES

Payment date: 20110426

Year of fee payment: 8

Ref country code: CH

Payment date: 20110428

Year of fee payment: 8

Ref country code: IE

Payment date: 20110426

Year of fee payment: 8

Ref country code: SE

Payment date: 20110414

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: NL

Payment date: 20110426

Year of fee payment: 8

Ref country code: AT

Payment date: 20110414

Year of fee payment: 8

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: IT

Payment date: 20110420

Year of fee payment: 8

REG Reference to a national code

Ref country code: NL

Ref legal event code: V1

Effective date: 20121101

REG Reference to a national code

Ref country code: CH

Ref legal event code: PL

REG Reference to a national code

Ref country code: SE

Ref legal event code: EUG

REG Reference to a national code

Ref country code: AT

Ref legal event code: MM01

Ref document number: 384962

Country of ref document: AT

Kind code of ref document: T

Effective date: 20120403

REG Reference to a national code

Ref country code: IE

Ref legal event code: MM4A

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: LI

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120430

Ref country code: AT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120403

Ref country code: CH

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120430

Ref country code: IE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120403

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: IT

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120403

Ref country code: SE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120404

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: NL

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20121101

REG Reference to a national code

Ref country code: ES

Ref legal event code: FD2A

Effective date: 20130716

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: ES

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20120404

REG Reference to a national code

Ref country code: DE

Ref legal event code: R081

Ref document number: 502004006046

Country of ref document: DE

Owner name: AIRBUS DEFENCE AND SPACE GMBH, DE

Free format text: FORMER OWNER: EADS DEUTSCHLAND GMBH, 85521 OTTOBRUNN, DE

Effective date: 20140819

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 13

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 14

REG Reference to a national code

Ref country code: FR

Ref legal event code: PLFP

Year of fee payment: 15

PGFP Annual fee paid to national office [announced via postgrant information from national office to epo]

Ref country code: GB

Payment date: 20220420

Year of fee payment: 19

Ref country code: FR

Payment date: 20220420

Year of fee payment: 19

Ref country code: DE

Payment date: 20220420

Year of fee payment: 19

REG Reference to a national code

Ref country code: DE

Ref legal event code: R119

Ref document number: 502004006046

Country of ref document: DE

GBPC Gb: european patent ceased through non-payment of renewal fee

Effective date: 20230403

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20230403

PG25 Lapsed in a contracting state [announced via postgrant information from national office to epo]

Ref country code: GB

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20230403

Ref country code: FR

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20230430

Ref country code: DE

Free format text: LAPSE BECAUSE OF NON-PAYMENT OF DUE FEES

Effective date: 20231103