CN115902834B - Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer - Google Patents

Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer Download PDF

Info

Publication number
CN115902834B
CN115902834B CN202211544782.6A CN202211544782A CN115902834B CN 115902834 B CN115902834 B CN 115902834B CN 202211544782 A CN202211544782 A CN 202211544782A CN 115902834 B CN115902834 B CN 115902834B
Authority
CN
China
Prior art keywords
fizeau interferometer
laser
target
echo signal
helium
Prior art date
Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
Active
Application number
CN202211544782.6A
Other languages
Chinese (zh)
Other versions
CN115902834A (en
Inventor
赵若灿
刘映妤
薛向辉
兰家欣
刘振威
陈廷娣
窦贤康
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Science and Technology of China USTC
Original Assignee
University of Science and Technology of China USTC
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 University of Science and Technology of China USTC filed Critical University of Science and Technology of China USTC
Priority to CN202211544782.6A priority Critical patent/CN115902834B/en
Priority to PCT/CN2023/072058 priority patent/WO2024113472A1/en
Publication of CN115902834A publication Critical patent/CN115902834A/en
Application granted granted Critical
Publication of CN115902834B publication Critical patent/CN115902834B/en
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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
    • 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/483Details of pulse systems
    • YGENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
    • Y02TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
    • Y02ATECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
    • Y02A90/00Technologies having an indirect contribution to adaptation to climate change
    • Y02A90/10Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation

Landscapes

  • 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)

Abstract

The invention provides a helium laser radar double-frequency temperature and wind measuring system based on a Fizeau interferometer, which comprises: a transmission system, the transmission system comprising: at least two lasers, which are suitable for emitting pulse lasers with different frequencies to the region to be measured; a receiving system, the receiving system comprising: a telescope adapted to receive echo signals from the region under test having metastable helium scattering information; the filter is suitable for carrying out first-stage filtering processing on the echo signal with metastable helium scattering information received by the telescope; the Fizeau interferometer is suitable for sequentially carrying out second-stage filtering treatment on the echo signals with metastable helium scattering information after the filter processing; the detector is suitable for receiving the scattered echo signals processed by the Fizeau interferometer; the data processing system is suitable for processing the echo signals received by the detector to obtain the wind speed and the temperature corresponding to the region to be detected. The invention also provides a helium laser radar double-frequency temperature and wind measuring method based on the Fizeau interferometer.

Description

基于斐索干涉仪的氦激光雷达双频测温测风系统及方法Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer

技术领域Technical Field

本发明涉及激光雷达技术领域,具体涉及一种基于斐索干涉仪的氦激光雷达双频测温测风系统及方法。The present invention relates to the field of laser radar technology, and in particular to a helium laser radar dual-frequency temperature and wind measurement system and method based on a Fizeau interferometer.

背景技术Background Art

热层大气温度和风速是研究与热层有关的热力学和动力学过程的重要参数。然而,此高度范围内的中性大气颗粒密度过低,使得从地基测量热层的密度、温度和风场廓线变得极为困难。过去,科学家大多使用探空火箭、卫星和其他航天器进行热层测量。然而,探空火箭成本高昂,无法进行长时间连续探空;卫星空间分辨率低,无法获取特定区域的长期不间断观测。Thermosphere temperature and wind speed are important parameters for studying thermodynamic and kinetic processes related to the thermosphere. However, the density of neutral atmospheric particles in this altitude range is too low, making it extremely difficult to measure the density, temperature and wind field profile of the thermosphere from the ground. In the past, scientists mostly used sounding rockets, satellites and other spacecraft to measure the thermosphere. However, sounding rockets are expensive and cannot conduct long-term continuous sounding; satellites have low spatial resolution and cannot obtain long-term uninterrupted observations of specific areas.

随着激光雷达技术的发展,共振荧光激光雷达逐渐成为热层探测的常用手段之一。比起传统的航天器测量手段,激光雷达可以提供热层连续的、高时空分辨率的数据。传统的三频测温测风手段已经在钠共振荧光激光雷达中得到了比较成熟和稳定的应用,但比起钠激光雷达,氦激光雷达的观测高度高、信号弱,并且采用一台发射机交替发射三个频率的激光实现测温测风的方法,需要将激光频率进行快速切换并交替发射,相当于将能量在频域上和时域上都减弱了三分之一,因此不利于获得好的信噪比。基于此,开发出更优化的测温测风方法是目前丞待解决的技术问题。With the development of lidar technology, resonant fluorescence lidar has gradually become one of the common means of thermosphere detection. Compared with traditional spacecraft measurement methods, lidar can provide continuous data of the thermosphere with high temporal and spatial resolution. The traditional three-frequency temperature and wind measurement methods have been relatively mature and stably applied in sodium resonance fluorescence lidar, but compared with sodium lidar, helium lidar has a high observation altitude and weak signal, and uses a transmitter to alternately emit three-frequency lasers to achieve temperature and wind measurement. The method requires rapid switching of the laser frequency and alternate emission, which is equivalent to reducing the energy by one third in both the frequency domain and the time domain, so it is not conducive to obtaining a good signal-to-noise ratio. Based on this, the development of more optimized temperature and wind measurement methods is a technical problem that needs to be solved at present.

发明内容Summary of the invention

针对上述问题,本申请提供了一种基于斐索干涉仪的氦激光雷达双频测温测风系统,通过在发射系统引入至少两台激光器、接收系统引入新型斐索干涉仪以及采集待测区域的亚稳态氦散射信息可以实现了双频测温测风。不仅克服了现有技术中存在的技术缺陷,也为后续研究测温测风提供了一个新的研究方向。In view of the above problems, this application provides a dual-frequency temperature and wind measurement system of helium laser radar based on Fizeau interferometer, which can realize dual-frequency temperature and wind measurement by introducing at least two lasers in the transmitting system, introducing a new Fizeau interferometer in the receiving system, and collecting metastable helium scattering information in the area to be measured. It not only overcomes the technical defects existing in the prior art, but also provides a new research direction for subsequent research on temperature and wind measurement.

为达到上述目的,本发明的第一方面提供了一种基于斐索干涉仪的氦激光雷达双频测温测风系统,包括:To achieve the above object, the first aspect of the present invention provides a helium laser radar dual-frequency temperature and wind measurement system based on a Fizeau interferometer, comprising:

发射系统,包括:Launch system, including:

至少两台激光器,适用于向待测区域发射频率不同的脉冲激光;At least two lasers, adapted to emit pulsed lasers of different frequencies to the area to be measured;

接收系统,包括:Receiving system, including:

望远镜,适用于接收来自上述待测区域的具有亚稳态氦散射信息的回波信号;A telescope, adapted to receive an echo signal having metastable helium scattering information from the above-mentioned area to be measured;

滤波片,适用于将上述望远镜接收到的上述具有亚稳态氦散射信息的回波信号进行第一级滤波处理;A filter, adapted to perform a first-stage filtering process on the echo signal having metastable helium scattering information received by the telescope;

斐索干涉仪,适用于将上述滤波片处理后的上述具有亚稳态氦散射信息的回波信号进行第二级滤波处理;A Fizeau interferometer, adapted to perform a second-stage filtering process on the echo signal having the metastable helium scattering information after the filter has been processed;

探测器,适用于接收来自上述斐索干涉仪处理后的散射回波信号,其中,上述斐索干涉仪处理后的回波信号通过测量光谱表示,上述测量光谱为散射光谱和上述斐索干涉仪的透过率函数的卷积;A detector adapted to receive a scattered echo signal processed by the Fizeau interferometer, wherein the echo signal processed by the Fizeau interferometer is represented by a measurement spectrum, and the measurement spectrum is a convolution of a scattered spectrum and a transmittance function of the Fizeau interferometer;

数据处理系统,适用于对上述探测器接收的回波信号进行处理,得到上述待测区域对应的风速和温度,其中,根据上述激光的回波信号谱和仿真得到回波信号光谱库的线性插值比对结果,确定上述待测区域对应的风速和温度,上述激光的回波信号谱为上述测量光谱和上述斐索干涉仪的透过率函数的反卷积。The data processing system is suitable for processing the echo signal received by the above-mentioned detector to obtain the wind speed and temperature corresponding to the above-mentioned area to be measured, wherein the wind speed and temperature corresponding to the above-mentioned area to be measured are determined according to the linear interpolation comparison result of the echo signal spectrum library obtained by simulation based on the echo signal spectrum of the above-mentioned laser, and the echo signal spectrum of the above-mentioned laser is the deconvolution of the above-mentioned measurement spectrum and the transmittance function of the above-mentioned Fizeau interferometer.

根据本发明的实施例,上述发射系统包括上述至少两台激光器,包括:According to an embodiment of the present invention, the transmitting system comprises the at least two lasers, including:

第一可调谐连续激光器和第二可调谐连续激光器,适用于产生稳定且连续的种子激光;A first tunable continuous laser and a second tunable continuous laser are adapted to generate stable and continuous seed lasers;

第一脉冲激光器和第二脉冲激光器,适用于将上述种子激光进行光放大,并发射到上述待测区域;A first pulse laser and a second pulse laser, adapted to amplify the seed laser and emit the amplified seed laser to the region to be measured;

其中,上述第一脉冲激光器和上述第二脉冲激光器将光放大后的一部分发射激光传输到第一波长计和第二波长计中实时监测上述发射激光的频率,并反馈给上述第一可调谐连续激光器和上述第二可调谐连续激光器,以实现上述发射激光频率锁定。Among them, the above-mentioned first pulse laser and the above-mentioned second pulse laser transmit a part of the emitted laser after light amplification to the first wavelength meter and the second wavelength meter to monitor the frequency of the above-mentioned emitted laser in real time, and feed back to the above-mentioned first tunable continuous laser and the above-mentioned second tunable continuous laser to achieve the above-mentioned emitted laser frequency locking.

根据本发明的实施例,上述第一可调谐连续激光器和第二可调谐连续激光器的结构相同,上述第一脉冲激光器和第二脉冲激光器的结构相同,上述第一波长计和第二波长计的结构相同。According to an embodiment of the present invention, the first tunable continuous laser and the second tunable continuous laser have the same structure, the first pulse laser and the second pulse laser have the same structure, and the first wavelength meter and the second wavelength meter have the same structure.

根据本发明的实施例,上述至少两台激光器垂直向上发射上述激光。According to an embodiment of the present invention, the at least two lasers emit the lasers vertically upwards.

本发明的第二方面提供了一种基于斐索干涉仪的氦激光雷达双频测温测风方法,应用于上述实施例中任一项所述的基于斐索干涉仪的氦激光雷达双频测温测风系统,上述方法包括:A second aspect of the present invention provides a dual-frequency temperature and wind measurement method of a helium laser radar based on a Fizeau interferometer, which is applied to a dual-frequency temperature and wind measurement system of a helium laser radar based on a Fizeau interferometer as described in any one of the above embodiments. The method comprises:

根据预设发射激光信息以及预设一定温度和风速范围内对应的多个温度和风速信息,仿真得到回波信号光谱数据库,其中,上述预设发射激光信息为上述至少两台激光器对应的待发射激光的频率信息;According to the preset laser emission information and the corresponding multiple temperature and wind speed information within a preset certain temperature and wind speed range, the echo signal spectrum database is simulated, wherein the preset laser emission information is the frequency information of the lasers to be emitted corresponding to the at least two lasers;

利用上述至少两台激光器向待测区域发射频率不同的脉冲激光;Using the at least two lasers mentioned above to emit pulsed lasers with different frequencies to the area to be tested;

利用上述望远镜接收来自上述待测区域的上述具有亚稳态氦散射信息的回波信号;Using the telescope to receive the echo signal with the metastable helium scattering information from the area to be measured;

利用上述滤波片和上述斐索干涉仪对上述望远镜接收到的上述具有亚稳态氦散射信息的回波信号依次进行两级滤波处理,并将处理后的回波信号输送到接收系统对应的光路中;The echo signal with metastable helium scattering information received by the telescope is subjected to two-stage filtering in sequence by using the filter and the Fizeau interferometer, and the processed echo signal is transmitted to the optical path corresponding to the receiving system;

探测器接收来自接收系统光路中的回波信号,其中,上述接收系统光路中的回波信号通过测量光谱表示,上述测量光谱为散射光谱和斐索干涉仪的透过率函数的卷积;The detector receives an echo signal from an optical path of a receiving system, wherein the echo signal in the optical path of the receiving system is represented by a measurement spectrum, and the measurement spectrum is a convolution of a scattering spectrum and a transmittance function of a Fizeau interferometer;

将上述测量光谱和上述斐索干涉仪的透过率函数进行反卷积处理,得到激光的回波信号谱;Deconvolution is performed on the measurement spectrum and the transmittance function of the Fizeau interferometer to obtain an echo signal spectrum of the laser;

将上述激光的回波信号谱与上述回波信号光谱数据库进行线性插值比对,以确定上述激光的回波信号谱中最接近的一组数据为待测区域对应的温度和风速。The echo signal spectrum of the laser is linearly interpolated and compared with the echo signal spectrum database to determine that the closest set of data in the echo signal spectrum of the laser is the temperature and wind speed corresponding to the area to be measured.

根据本发明的实施例,利用上述至少两台激光器向待测区域发射频率不同的脉冲激光,包括:According to an embodiment of the present invention, using the at least two lasers to emit pulse lasers with different frequencies to the area to be measured includes:

利用上述第一可调谐连续激光器和上述第二可调谐连续激光器产生稳定且连续的上述种子激光;Using the first tunable continuous laser and the second tunable continuous laser to generate the stable and continuous seed laser;

利用上述第一脉冲激光器和上述第二脉冲激光器将上述种子激光进行光放大,并发射到上述待测区域;The seed laser is optically amplified by the first pulse laser and the second pulse laser, and emitted to the area to be measured;

其中,在上述第一脉冲激光器和上述第二脉冲激光器将上述种子激光进行光放大后,利用上述第一波长计和上述第二波长计实时监测接收到的一部分上述发射激光的频率,并反馈给上述第一可调谐连续激光器和上述第二可调谐连续激光器,以实现上述发射激光频率锁定。Among them, after the above-mentioned first pulse laser and the above-mentioned second pulse laser optically amplify the above-mentioned seed laser, the above-mentioned first wavelength meter and the above-mentioned second wavelength meter are used to monitor the frequency of a part of the above-mentioned emission laser received in real time, and feed back to the above-mentioned first tunable continuous laser and the above-mentioned second tunable continuous laser to achieve the above-mentioned emission laser frequency locking.

根据本发明的实施例,上述至少两台激光器向待测区域发射两束频率不同的脉冲激光,可以得到四个不同频率的上述具有亚稳态氦散射信息的回波信号。According to an embodiment of the present invention, the at least two lasers emit two pulsed lasers with different frequencies to the area to be measured, so that four echo signals with the metastable helium scattering information at different frequencies can be obtained.

本发明的第三方面提供了一种斐索干涉仪的参数确定方法,应用于上述实施例中任一项所述的基于斐索干涉仪的氦激光雷达双频测温测风系统,上述方法包括:A third aspect of the present invention provides a method for determining parameters of a Fizeau interferometer, which is applied to a helium laser radar dual-frequency temperature and wind measurement system based on a Fizeau interferometer as described in any one of the above embodiments, and the method comprises:

获取初始斐索干涉仪的多个属性信息,其中,上述多个属性信息包括:上述初始斐索干涉仪的多个半高全宽组成的集合、多个自由谱间距组成的集合以及上述滤波片处理后的上述具有亚稳态氦散射信息的回波信号到达上述初始斐索干涉仪的入射倾角;Acquire multiple attribute information of the initial Fizeau interferometer, wherein the multiple attribute information includes: a set of multiple half-widths of the initial Fizeau interferometer, a set of multiple free spectrum spacings, and an incident inclination angle of the echo signal with metastable helium scattering information after being processed by the filter to the initial Fizeau interferometer;

基于上述仿真得到回波信号光谱数据库,遍历上述多个半高全宽组成的集合以及上述多个自由谱间距组成的集合,以确定目标半高全宽和目标自由谱间距;Based on the above simulation, an echo signal spectrum database is obtained, and the set consisting of the above multiple full widths at half maximum and the set consisting of the above multiple free spectrum intervals are traversed to determine the target full width at half maximum and the target free spectrum interval;

根据目标半高全宽和目标自由谱间距,确定斐索干涉仪的目标精细因子;Determine the target fine factor of the Fizeau interferometer based on the target full width at half maximum and the target free spectrum spacing;

根据上述目标精细因子,确定上述斐索干涉仪的目标反射率;Determining the target reflectivity of the Fizeau interferometer according to the target fine factor;

根据目标自由谱间距对应的波长,确定上述斐索干涉仪的目标平均腔长;Determining a target average cavity length of the Fizeau interferometer according to a wavelength corresponding to a target free spectrum spacing;

根据上述目标平均腔长,确定上述斐索干涉仪的目标楔角;Determining a target wedge angle of the Fizeau interferometer according to the target average cavity length;

其中,上述目标楔角为上述斐索干涉仪的两块平板对应的夹角,上述目标楔角用于调控上述斐索干涉仪的干涉条纹的宽度。The target wedge angle is the included angle between the two plates of the Fizeau interferometer, and the target wedge angle is used to adjust the width of the interference fringes of the Fizeau interferometer.

根据本发明的实施例,上述斐索干涉仪的目标精细因子和斐索干涉仪的目标反射率通过公式(1)计算得到;According to an embodiment of the present invention, the target fineness factor of the Fizeau interferometer and the target reflectivity of the Fizeau interferometer are calculated by formula (1);

公式(1)中,F表示上述斐索干涉仪的目标精细因子,FSR表示上述目标自由谱间距,FWHM表示上述目标半高全宽,R表示上述斐索干涉仪的目标反射率。In formula (1), F represents the target fineness factor of the Fizeau interferometer, FSR represents the target free spectrum spacing, FWHM represents the target full width at half maximum, and R represents the target reflectivity of the Fizeau interferometer.

根据本发明的实施例,上述斐索干涉仪的目标平均腔长和目标楔角通过公式(2)计算得到;According to an embodiment of the present invention, the target average cavity length and target wedge angle of the above-mentioned Fizeau interferometer are calculated by formula (2);

公式(2)中,L0表示上述斐索干涉仪的目标平均腔长,n表示介质的折射率,λ表示上述发射激光的波长,λFSR表示上述目标自由谱间距对应的波长,α表示上述斐索干涉仪的目标楔角。In formula (2), L0 represents the target average cavity length of the Fizeau interferometer, n represents the refractive index of the medium, λ represents the wavelength of the emitted laser, λ FSR represents the wavelength corresponding to the target free spectrum spacing, and α represents the target wedge angle of the Fizeau interferometer.

根据本发明的实施例,通过在发射系统引入至少两台激光器,确保同时发射两束频率不同的脉冲激光,克服了现有技术中利用三台激光器实现三频轮流切换并发射三个频率激光进行测温测风方法时引入的缺陷;该双频发射系统不仅可以节省传统激光轮流发射的时间,还可以减少激光雷达的复杂度和维护难度。According to an embodiment of the present invention, at least two lasers are introduced into the transmitting system to ensure that two pulsed lasers with different frequencies are emitted simultaneously, thereby overcoming the defects introduced by the prior art method of using three lasers to achieve three-frequency switching and emit three-frequency lasers for temperature and wind measurement. The dual-frequency transmission system can not only save the time of traditional laser emission in turn, but also reduce the complexity and maintenance difficulty of the laser radar.

根据本发明的实施例,通过改造部分模块,即在接收系统引入新型斐索干涉仪代替传统的F-P标准具,可以产生能够替代环状条纹的线形条纹,不需要进行频率扫描,直接配合线列探测器即可完成对线形条纹位置和形状变化的探测;该新型接收系统不仅可以节省扫频时间,持续累积光子信号,同时还极大地简化了数据处理的方法。According to an embodiment of the present invention, by modifying some modules, that is, introducing a new Fizeau interferometer in the receiving system to replace the traditional F-P standard, linear fringes that can replace the annular fringes can be generated. There is no need to perform frequency scanning, and the position and shape changes of the linear fringes can be detected directly in conjunction with a linear detector. This new receiving system can not only save scanning time and continuously accumulate photon signals, but also greatly simplify the data processing method.

根据本发明的实施例,氦原子是200km以上的主要大气成分之一,经过光电子激发等物理过程,可产生比其更稳定、荧光更强的亚稳态氦[He(23S)];因此与现有技术中的采用钠共振荧光激光雷达测温测风相比,通过采集待测区域的亚稳态氦散射信息进行测温测风,可以为后续研究测温测风提供一个新的研究方向,也填补了氦激光雷达在测温测风方面的理论空白。从而能够对热层大气物理过程和多圈层耦合机制进行更为全方位的观测和研究。According to the embodiments of the present invention, helium atoms are one of the main atmospheric components above 200 km. After physical processes such as photoelectron excitation, metastable helium [He(2 3 S)] which is more stable and has stronger fluorescence can be produced. Therefore, compared with the temperature and wind measurement using sodium resonance fluorescence laser radar in the prior art, the temperature and wind measurement by collecting metastable helium scattering information in the area to be measured can provide a new research direction for subsequent research on temperature and wind measurement, and also fill the theoretical gap in the temperature and wind measurement of helium laser radar. This enables a more comprehensive observation and research on the physical processes of the thermosphere and the multi-sphere coupling mechanism.

附图说明BRIEF DESCRIPTION OF THE DRAWINGS

通过以下参照附图对本公开实施例的描述,本公开的上述内容以及其他目的、特征和优点将更为清楚,在附图中:The above contents and other purposes, features and advantages of the present disclosure will become more apparent through the following description of the embodiments of the present disclosure with reference to the accompanying drawings, in which:

图1示意性示出了根据本发明实施例的基于斐索干涉仪的氦激光雷达双频测温测风系统框图;FIG1 schematically shows a block diagram of a helium laser radar dual-frequency temperature and wind measurement system based on a Fizeau interferometer according to an embodiment of the present invention;

图2示意性示出了根据本发明实施例的基于斐索干涉仪的氦激光雷达双频测温测风系统示意图;FIG2 schematically shows a schematic diagram of a helium laser radar dual-frequency temperature and wind measurement system based on a Fizeau interferometer according to an embodiment of the present invention;

图3示意性示出了根据本发明实施例的亚稳态氦在1083nm附近的三条辐射线对应的后向散射截面随波长变化曲线图;FIG3 schematically shows a graph showing backscattering cross sections of metastable helium corresponding to three radiation lines around 1083 nm as a function of wavelength according to an embodiment of the present invention;

图4示意性示出了根据本发明实施例的基于斐索干涉仪的氦激光雷达双频测温测风方法流程图;FIG4 schematically shows a flow chart of a dual-frequency temperature and wind measurement method using a helium laser radar based on a Fizeau interferometer according to an embodiment of the present invention;

图5示意性示出了根据本发明实施例的后向散射截面随温度和风速的变化曲线图;FIG5 schematically shows a curve diagram of backscattering cross section changes with temperature and wind speed according to an embodiment of the present invention;

图6示意性示出了根据本发明实施例的回波信号光谱数据库中的两组数据示例图;FIG6 schematically shows two sets of data examples in an echo signal spectrum database according to an embodiment of the present invention;

图7示意性示出了根据本发明实施例的亚稳态氦峰的后向散射截面和回波信号的归一化强度示意图;FIG7 schematically shows a schematic diagram of a backscattering cross section of a metastable helium peak and a normalized intensity of an echo signal according to an embodiment of the present invention;

图8示意性示出了根据本发明实施例的斐索干涉仪的参数确定方法流程图;FIG8 schematically shows a flow chart of a method for determining parameters of a Fizeau interferometer according to an embodiment of the present invention;

图9示意性示出了根据本发明实施例的不同形状因子与斐索干涉仪的透过率曲线关系图;FIG9 schematically shows a relationship diagram of transmittance curves of a Fizeau interferometer and different shape factors according to an embodiment of the present invention;

图10示意性示出了根据本发明实施例的数据处理流程图;FIG10 schematically shows a data processing flow chart according to an embodiment of the present invention;

图11示意性示出了根据本发明实施例的探测器探测到的激光线型、总的回波信号以及斐索干涉仪的卷积结果示意图。FIG11 schematically shows a schematic diagram of the laser line shape, the total echo signal and the convolution result of the Fizeau interferometer detected by the detector according to an embodiment of the present invention.

附图标记说明:Description of reference numerals:

1:第一可调谐连续激光器;1: first tunable continuous laser;

2:第一脉冲激光器;2: first pulse laser;

3:第一波长计;3: First wavelength meter;

4:第二可调谐连续激光器;4: Second tunable continuous laser;

5:第二脉冲激光器;5: second pulse laser;

6:第二波长计;6: Second wavelength meter;

7:望远镜;7: Telescope;

8:光开关;8: Optical switch;

9:滤波片;9: filter;

10:斐索干涉仪;10: Fizeau interferometer;

11:探测器。11: Detector.

具体实施方式DETAILED DESCRIPTION

以下,将参照附图来描述本公开的实施例。但是应该理解,这些描述只是示例性的,而并非要限制本公开的范围。在下面的详细描述中,为便于解释,阐述了许多具体的细节以提供对本公开实施例的全面理解。然而,明显地,一个或多个实施例在没有这些具体细节的情况下也可以被实施。此外,在以下说明中,省略了对公知结构和技术的描述,以避免不必要地混淆本公开的概念。Hereinafter, embodiments of the present disclosure will be described with reference to the accompanying drawings. However, it should be understood that these descriptions are exemplary only and are not intended to limit the scope of the present disclosure. In the following detailed description, for ease of explanation, many specific details are set forth to provide a comprehensive understanding of the embodiments of the present disclosure. However, it is apparent that one or more embodiments may also be implemented without these specific details. In addition, in the following description, descriptions of known structures and technologies are omitted to avoid unnecessary confusion of the concepts of the present disclosure.

在此使用的术语仅仅是为了描述具体实施例,而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在,但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。The terms used herein are only for describing specific embodiments and are not intended to limit the present disclosure. The terms "include", "comprising", etc. used herein indicate the existence of the features, steps, operations and/or components, but do not exclude the existence or addition of one or more other features, steps, operations or components.

在此使用的所有术语(包括技术和科学术语)具有本领域技术人员通常所理解的含义,除非另外定义。应注意,这里使用的术语应解释为具有与本说明书的上下文相一致的含义,而不应以理想化或过于刻板的方式来解释。All terms (including technical and scientific terms) used herein have the meanings commonly understood by those skilled in the art, unless otherwise defined. It should be noted that the terms used herein should be interpreted as having a meaning consistent with the context of this specification, and should not be interpreted in an idealized or overly rigid manner.

在使用类似于“A、B和C等中至少一个”这样的表述的情况下,一般来说应该按照本领域技术人员通常理解该表述的含义来予以解释(例如,“具有A、B和C中至少一个的系统”应包括但不限于单独具有A、单独具有B、单独具有C、具有A和B、具有A和C、具有B和C、和/或具有A、B、C的系统等)。When using expressions such as "at least one of A, B, and C", they should generally be interpreted according to the meaning of the expression commonly understood by technical personnel in this field (for example, "a system having at least one of A, B, and C" should include but is not limited to a system having A alone, B alone, C alone, A and B, A and C, B and C, and/or A, B, C, etc.).

为使本发明的目的、技术方案和优点更加清楚明白,以下结合具体实施例,并参照附图,对本发明作进一步的详细说明。In order to make the objectives, technical solutions and advantages of the present invention more clearly understood, the present invention is further described in detail below in conjunction with specific embodiments and with reference to the accompanying drawings.

随着激光雷达技术的发展,共振荧光激光雷达逐渐成为热层探测的常用手段之一。氦原子是200km以上的主要大气成分之一,经过光电子激发等物理过程,可产生比其更稳定、荧光更强的亚稳态氦[He(23S)],因此亚稳态氦可以作为地基激光雷达观测该高度范围中性大气首选的示踪物之一。亚稳态氦共振荧光激光雷达探测大气参数的原理为:激光雷达的发射端竖直向上空发射波长为1083nm的窄线宽激光束,激发在热层中的氦原子,使其核外电子跃迁到高能级并退激发到较低能级,即发生共振荧光反应。在地面,激光雷达的接收端望远镜接收回波信号,利用光电探测器探测回波光子数,并进行一系列的信号处理和数据反演,从而得到热层亚稳态氦密度。与此同时,还可以通过让雷达激光发射机交替发射三个频率的激光,通过三个频率激发的共振荧光强度之间的相对关系,进一步反演得到亚稳态氦在探测高度范围内的温度和风场。再结合大气模式,并基于一些假设和边界条件,可以推导出背景大气的密度、温度和风场。With the development of lidar technology, resonant fluorescence lidar has gradually become one of the common means of thermosphere detection. Helium atoms are one of the main atmospheric components above 200 km. After physical processes such as photoelectron excitation, they can produce metastable helium [He(2 3 S)] that is more stable and has stronger fluorescence than helium. Therefore, metastable helium can be used as one of the preferred tracers for ground-based lidar to observe the neutral atmosphere in this altitude range. The principle of metastable helium resonance fluorescence lidar to detect atmospheric parameters is as follows: the transmitting end of the lidar emits a narrow linewidth laser beam with a wavelength of 1083nm vertically upward to excite helium atoms in the thermosphere, causing their extranuclear electrons to transition to high energy levels and de-excite to lower energy levels, that is, a resonant fluorescence reaction occurs. On the ground, the receiving end telescope of the lidar receives the echo signal, uses a photodetector to detect the number of echo photons, and performs a series of signal processing and data inversion to obtain the metastable helium density in the thermosphere. At the same time, the temperature and wind field of metastable helium within the detection altitude range can be further inverted by letting the radar laser transmitter emit three frequencies of laser alternately and using the relative relationship between the resonance fluorescence intensities excited by the three frequencies. Combined with the atmospheric model and based on some assumptions and boundary conditions, the density, temperature and wind field of the background atmosphere can be derived.

然而,相关技术中不能较好地实现亚稳态氦测温测风激光雷达探测技术。相关技术中的测温测风手段,一般也仅局限在使用钠激光雷达的三频测温测风。与之相比,氦激光雷达探测的高度更高、信号更弱,如果要实现氦激光雷达测温测风功能,需要开发出更优化的办法,才能实现较高的信噪比,进而得到更为准确的温度和风场测量结果。However, the metastable helium temperature and wind measurement lidar detection technology cannot be well implemented in the related technology. The temperature and wind measurement methods in the related technology are generally limited to three-frequency temperature and wind measurement using sodium lidar. In comparison, the helium lidar detects at a higher altitude and has a weaker signal. If the helium lidar temperature and wind measurement function is to be realized, a more optimized method needs to be developed to achieve a higher signal-to-noise ratio, thereby obtaining more accurate temperature and wind field measurement results.

基于此,本申请提供了一种基于斐索干涉仪的氦激光雷达双频测温测风系统,通过在发射系统引入至少两台激光器、接收系统引入新型斐索干涉仪以及采集待测区域的亚稳态氦散射信息实现了双频测温测风。不仅克服了现有技术中存在的技术缺陷,也为后续研究测温测风提供了一个新的研究方向。Based on this, the present application provides a dual-frequency temperature and wind measurement system of helium laser radar based on Fizeau interferometer, which realizes dual-frequency temperature and wind measurement by introducing at least two lasers in the transmitting system, introducing a new Fizeau interferometer in the receiving system, and collecting metastable helium scattering information in the area to be measured. It not only overcomes the technical defects existing in the prior art, but also provides a new research direction for subsequent research on temperature and wind measurement.

图1示意性示出了根据本发明实施例的基于斐索干涉仪的氦激光雷达双频测温测风系统框图。FIG1 schematically shows a block diagram of a dual-frequency temperature and wind measurement system of a helium laser radar based on a Fizeau interferometer according to an embodiment of the present invention.

如图1所示,基于斐索干涉仪的氦激光雷达双频测温测风系统包括发射系统、接收系统以及数据处理系统。As shown in Figure 1, the helium lidar dual-frequency temperature and wind measurement system based on Fizeau interferometer includes a transmitting system, a receiving system and a data processing system.

图2示意性示出了根据本发明实施例的基于斐索干涉仪的氦激光雷达双频测温测风系统示意图。FIG2 schematically shows a schematic diagram of a helium laser radar dual-frequency temperature and wind measurement system based on a Fizeau interferometer according to an embodiment of the present invention.

根据本发明的示例性实施例,本发明公开了一种基于斐索干涉仪的氦激光雷达双频测温测风系统,参考图2所示,发射系统包括:至少两台激光器,适用于向待测区域发射频率不同的脉冲激光;接收系统包括:望远镜7、光开关8、滤波片9、斐索干涉仪10以及探测器11,其中,望远镜7,适用于接收来自待测区域的具有亚稳态氦散射信息的回波信号;滤波片9,适用于将望远镜7接收到的具有亚稳态氦散射信息的回波信号进行第一级滤波处理;斐索干涉仪10,适用于将滤波片9处理后的具有亚稳态氦散射信息的回波信号进行第二级滤波处理;探测器11,适用于接收来自斐索干涉仪10处理后的散射回波信号,其中,斐索干涉仪10处理后的回波信号通过测量光谱表示,测量光谱为散射光谱和斐索干涉仪10的透过率函数的卷积;以及数据处理系统,适用于对探测器11接收的回波信号进行处理,得到待测区域对应的风速和温度,其中,根据激光的回波信号谱和仿真得到回波信号光谱库的线性插值比对结果,确定待测区域对应的风速和温度,激光的回波信号谱为测量光谱和斐索干涉仪10的透过率函数的反卷积。According to an exemplary embodiment of the present invention, the present invention discloses a dual-frequency temperature and wind measurement system of a helium laser radar based on a Fizeau interferometer. As shown in FIG2 , the transmitting system includes: at least two lasers, which are suitable for transmitting pulse lasers with different frequencies to the area to be measured; the receiving system includes: a telescope 7, an optical switch 8, a filter 9, a Fizeau interferometer 10 and a detector 11, wherein the telescope 7 is suitable for receiving an echo signal with metastable helium scattering information from the area to be measured; the filter 9 is suitable for performing a first-stage filtering process on the echo signal with metastable helium scattering information received by the telescope 7; the Fizeau interferometer 10 is suitable for filtering the echo signal with metastable helium scattering information processed by the filter 9. The echo signal of the information is subjected to a second-stage filtering process; the detector 11 is adapted to receive the scattered echo signal processed by the Fizeau interferometer 10, wherein the echo signal processed by the Fizeau interferometer 10 is represented by a measurement spectrum, and the measurement spectrum is the convolution of the scattered spectrum and the transmittance function of the Fizeau interferometer 10; and the data processing system is adapted to process the echo signal received by the detector 11 to obtain the wind speed and temperature corresponding to the area to be measured, wherein the wind speed and temperature corresponding to the area to be measured are determined according to the linear interpolation comparison result of the echo signal spectrum library obtained by simulation and the echo signal spectrum spectrum spectrum of the laser, and the echo signal spectrum of the laser is the deconvolution of the measurement spectrum and the transmittance function of the Fizeau interferometer 10.

根据本发明的实施例,通过在发射系统引入至少两台激光器,确保同时发射两束频率不同的脉冲激光,克服了现有技术中利用三台激光器实现三频轮流切换并发射三个频率激光进行测温测风方法时引入的缺陷;该双频发射系统不仅可以节省传统激光轮流发射的时间,还可以减少激光雷达的复杂度和维护难度。According to an embodiment of the present invention, at least two lasers are introduced into the transmitting system to ensure that two pulsed lasers with different frequencies are emitted simultaneously, thereby overcoming the defects introduced by the prior art method of using three lasers to achieve three-frequency switching and emit three-frequency lasers for temperature and wind measurement. The dual-frequency transmission system can not only save the time of traditional laser emission in turn, but also reduce the complexity and maintenance difficulty of the laser radar.

根据本发明的实施例,通过改造部分模块,即在接收系统引入新型斐索干涉仪10代替传统的F-P标准具,可以产生能够替代环状条纹的线形条纹,不需要进行频率扫描,直接配合线列探测器11即可完成对线形条纹位置和形状变化的探测;该新型接收系统不仅可以节省扫频时间,持续累积光子信号,同时还极大地简化了数据处理的方法。According to an embodiment of the present invention, by modifying some modules, that is, introducing a new Fizeau interferometer 10 in the receiving system to replace the traditional F-P standard, linear stripes that can replace the annular stripes can be generated, and there is no need to perform frequency scanning. The linear stripe position and shape change can be detected directly in conjunction with the linear array detector 11; this new receiving system can not only save scanning time and continuously accumulate photon signals, but also greatly simplify the data processing method.

根据本发明的实施例,氦原子是200km以上的主要大气成分之一,经过光电子激发等物理过程,可产生比其更稳定、荧光更强的亚稳态氦[He(23S)];因此与现有技术中的采用钠共振荧光激光雷达测温测风相比,通过采集待测区域的亚稳态氦散射信息进行测温测风,可以为后续研究测温测风提供一个新的研究方向,也填补了氦激光雷达在测温测风方面的理论空白。从而能够对热层大气物理过程和多圈层耦合机制进行更为全方位的观测和研究。According to the embodiments of the present invention, helium atoms are one of the main atmospheric components above 200 km. After physical processes such as photoelectron excitation, metastable helium [He(2 3 S)] which is more stable and has stronger fluorescence can be produced. Therefore, compared with the temperature and wind measurement using sodium resonance fluorescence laser radar in the prior art, the temperature and wind measurement by collecting metastable helium scattering information in the area to be measured can provide a new research direction for subsequent research on temperature and wind measurement, and also fill the theoretical gap in helium laser radar in temperature and wind measurement. This enables a more comprehensive observation and research on the physical processes of the thermosphere and the multi-sphere coupling mechanism.

根据本发明的实施例,发射系统包括至少两台激光器,包括:第一可调谐连续激光器1、第一脉冲激光器2、第二可调谐连续激光器4和第二脉冲激光器5;According to an embodiment of the present invention, the transmitting system includes at least two lasers, including: a first tunable continuous laser 1, a first pulse laser 2, a second tunable continuous laser 4 and a second pulse laser 5;

第一可调谐连续激光器1和第二可调谐连续激光器4,适用于产生波长在1083nm附近的稳定且连续的种子激光;A first tunable continuous laser 1 and a second tunable continuous laser 4 are adapted to generate a stable and continuous seed laser having a wavelength of about 1083 nm;

第一脉冲激光器2和第二脉冲激光器5,适用于将种子激光进行光放大,并发射到待测区域;The first pulse laser 2 and the second pulse laser 5 are adapted to amplify the seed laser light and emit it to the area to be measured;

其中,第一脉冲激光器2和第二脉冲激光器5分别将光放大后的一部分发射激光传输到第一波长计3和第二波长计6中实时监测发射激光的频率,并反馈给第一可调谐连续激光器1和第二可调谐连续激光器4,以实现发射激光频率锁定。Among them, the first pulse laser 2 and the second pulse laser 5 respectively transmit a part of the emitted laser after light amplification to the first wavelength meter 3 and the second wavelength meter 6 to monitor the frequency of the emitted laser in real time, and feed back to the first tunable continuous laser 1 and the second tunable continuous laser 4 to achieve frequency locking of the emitted laser.

根据本发明的实施例,第一可调谐连续激光器1和第二可调谐连续激光器4的结构相同,第一脉冲激光器2和第二脉冲激光器5的结构相同,第一波长计3和第二波长计6的结构相同。According to the embodiment of the present invention, the first tunable continuous laser 1 and the second tunable continuous laser 4 have the same structure, the first pulse laser 2 and the second pulse laser 5 have the same structure, and the first wavelength meter 3 and the second wavelength meter 6 have the same structure.

根据本发明的实施例,第一脉冲激光器2和第二脉冲激光器5为固体掺镱钇铝石榴石晶体(Yb:YAG)脉冲激光器,可以输出1083nm附近的脉冲光,并垂直向上发射到待测区域。According to an embodiment of the present invention, the first pulse laser 2 and the second pulse laser 5 are solid ytterbium-doped yttrium aluminum garnet crystal (Yb:YAG) pulse lasers, which can output pulse light around 1083nm and emit it vertically upward to the area to be measured.

根据本发明的实施例,对于氦原子,有单态和三重态两套能级,相互之间无跃迁。在单态能级中,氦的第一激发态为21S0,寿命约为19.5ms;在三重态能级中,氦的第一激发态为23S1,包含1082.909nm、1083.025nm和1083.034nm三条辐射线,在理想情况下,氦(23S1)自然存在的寿命约为8000s,适合作为探测热层大气温度和风速的示踪物。According to an embodiment of the present invention, for helium atoms, there are two sets of energy levels, singlet and triplet, and there is no transition between them. In the singlet energy level, the first excited state of helium is 2 1 S 0 , with a lifetime of about 19.5 ms; in the triplet energy level, the first excited state of helium is 2 3 S 1 , including three radiation lines of 1082.909 nm, 1083.025 nm and 1083.034 nm. Under ideal conditions, the natural lifetime of helium (2 3 S 1 ) is about 8000 s, which is suitable as a tracer for detecting the temperature and wind speed of the thermosphere.

图3示意性示出了根据本发明实施例的亚稳态氦在1083nm附近的三条辐射线对应的后向散射截面随波长变化曲线图。FIG3 schematically shows a graph showing backscattering cross sections of metastable helium corresponding to three radiation lines around 1083 nm as a function of wavelength according to an embodiment of the present invention.

如图3所示,曲线1-909、1-025以及1-034分别表示亚稳态氦在1082.909nm、1082.025nm以及1082.034nm处的三条辐射线对应的后向散射截面。相比之下,亚稳态氦在1082.909nm处的辐射线受羟基干扰严重,故在本申请中选用波长1083.025nm和1083.034nm对应的两条辐射线进行观测。As shown in Figure 3, curves 1-909, 1-025 and 1-034 represent the backscattering cross sections corresponding to the three radiation lines of metastable helium at 1082.909nm, 1082.025nm and 1082.034nm, respectively. In contrast, the radiation line of metastable helium at 1082.909nm is seriously interfered by hydroxyl groups, so in this application, two radiation lines corresponding to wavelengths of 1083.025nm and 1083.034nm are selected for observation.

图4示意性示出了根据本发明实施例的基于斐索干涉仪的氦激光雷达双频测温测风方法流程图。FIG4 schematically shows a flow chart of a dual-frequency temperature and wind measurement method using a helium laser radar based on a Fizeau interferometer according to an embodiment of the present invention.

如图4所示,基于斐索干涉仪的氦激光雷达双频测温测风系统方法可以包括操作S410~操作S470。As shown in FIG. 4 , the helium laser radar dual-frequency temperature and wind measurement system method based on the Fizeau interferometer may include operations S410 to S470.

在操作S410,根据预设发射激光信息以及预设一定温度和风速范围内对应的多个温度和风速信息,仿真得到回波信号光谱数据库,其中,预设发射激光信息为至少两台激光器对应的待发射激光的频率信息。In operation S410, an echo signal spectrum database is simulated based on preset laser emission information and multiple temperature and wind speed information corresponding to a preset temperature and wind speed range, wherein the preset laser emission information is frequency information of lasers to be emitted corresponding to at least two lasers.

根据本发明的实施例,氦原子的速度分布服从麦克斯韦分布定律,每个氦原子的统计平均吸收截面通过公式(一)计算得到:According to an embodiment of the present invention, the velocity distribution of helium atoms obeys Maxwell's distribution law, and the statistical average absorption cross section of each helium atom is calculated by formula (1):

公式(一)中,多普勒吸收截面v表示激光频率,v0表示每个氦原子跃迁线的共振频率,w表示径向风速,λ表示激光波长,i表示较低的能级,k表示较高的能级,gi和gk分别是低能级和高能级的退化因子,Aki表示从k能级到i能级的跃迁概率,即爱因斯坦系数(s-1),kB表示玻尔兹曼常数,T表示原子温度,M表示原子质量,c表示光速。径向风速为负时,则表示风向为靠近雷达方向;径向风速为正时,则表示风向为远离雷达方向。In formula (1), the Doppler absorption cross section v represents the laser frequency, v0 represents the resonant frequency of each helium atomic transition line, w represents the radial wind speed, λ represents the laser wavelength, i represents the lower energy level, k represents the higher energy level, gi and gk represent the degradation factors of the lower and higher energy levels, Aki represents the transition probability from the k energy level to the i energy level, i.e., the Einstein coefficient (s -1 ), kB represents the Boltzmann constant, T represents the atomic temperature, M represents the atomic mass, and c represents the speed of light. When the radial wind speed is negative, it means that the wind direction is close to the radar; when the radial wind speed is positive, it means that the wind direction is away from the radar.

假设所有发射的光子数等于吸收的光子数,有效截面通过公式(二)计算得到:Assuming that the number of all emitted photons is equal to the number of absorbed photons, the effective cross section is calculated using formula (II):

σeff(v,v0)=σabs(v,v0);(二)σ eff (v, v 0 ) = σ abs (v, v 0 ); (2)

公式(二)中,σeff表示每个氦原子的有效截面,σabs表示每个氦原子的统计平均吸收截面,v0表示每个氦原子跃迁线的共振频率,v表示激光频率。In formula (ii), σ eff represents the effective cross section of each helium atom, σ abs represents the statistical average absorption cross section of each helium atom, v 0 represents the resonant frequency of each helium atom transition line, and v represents the laser frequency.

根据本发明的实施例,在不同的热层温度和风速下,亚稳态氦在1083nm峰的后向散射截面会发生变化,进而影响回波信号的相对强度。通过改变温度和风速,即可得到对应的有效散射截面谱线。According to an embodiment of the present invention, under different thermosphere temperatures and wind speeds, the backscattering cross section of metastable helium at the 1083nm peak will change, thereby affecting the relative intensity of the echo signal. By changing the temperature and wind speed, the corresponding effective scattering cross section spectrum can be obtained.

图5示意性示出了根据本发明实施例的后向散射截面随温度和风速的变化曲线图。FIG. 5 schematically shows a curve diagram of backscattering cross section changes with temperature and wind speed according to an embodiment of the present invention.

图5绘制出了在不同温度和不同风速下的后向散射截面。图5中(a)图表示在径向风速为0m/s时,在温度分别为800K、1000K和1200K时的后向散射截面。可以发现,随着温度的升高,由于多普勒展宽,半高全宽增大,峰值强度降低,1083nm处主峰的最大值位置逐渐向短波方向移动。图5中(b)图表示在温度为1000K时,在径向风速分别为-200m/s、0m/s和+200m/s时对应的后向散射截面。若风速为正,则风向为大气指向雷达方向,谱线总体向短波方向移动;若风速为负则风向为大气远离雷达方向,谱线总体向长波方向移动。Figure 5 plots the backscattering cross sections at different temperatures and wind speeds. Figure 5 (a) shows the backscattering cross sections at 800K, 1000K and 1200K when the radial wind speed is 0m/s. It can be found that with the increase of temperature, due to Doppler broadening, the full width at half maximum increases, the peak intensity decreases, and the maximum position of the main peak at 1083nm gradually moves toward the shortwave direction. Figure 5 (b) shows the backscattering cross sections corresponding to -200m/s, 0m/s and +200m/s when the temperature is 1000K. If the wind speed is positive, the wind direction is the direction of the atmosphere pointing to the radar, and the spectrum line generally moves toward the shortwave direction; if the wind speed is negative, the wind direction is the direction of the atmosphere away from the radar, and the spectrum line generally moves toward the longwave direction.

在温度范围700K至2000K内、风速-200m/s至+200m/s范围内,将温度和风速的刻度取0.1K和0.1m/s,得到不同温度和风速的亚稳态氦散射截面。结合待发射的两束激光对应的频率信息,将各组唯一对应于不同温度、不同风速的散射截面信息组成一个数据集合,即可以仿真得到对应于特定温度和风速的回波信号光谱数据库。In the temperature range of 700K to 2000K and the wind speed range of -200m/s to +200m/s, the temperature and wind speed scales are set to 0.1K and 0.1m/s, and the metastable helium scattering cross sections at different temperatures and wind speeds are obtained. Combined with the frequency information corresponding to the two laser beams to be emitted, each group of scattering cross-section information corresponding to different temperatures and wind speeds is combined into a data set, that is, the echo signal spectrum database corresponding to a specific temperature and wind speed can be simulated.

图6示意性示出了根据本发明实施例的回波信号光谱数据库中的两组数据示例图。FIG. 6 schematically shows two sets of data examples in an echo signal spectrum database according to an embodiment of the present invention.

如图6所示,虚线表示(800K,+200m/s)时的回波信号,实线表示(1200K,-200m/s)时的回波信号。可以发现,在不同温度和风速的情况下,两组回波线型不同,即每一组回波线型对应于特定温度和风速。As shown in Figure 6, the dotted line represents the echo signal at (800K, +200m/s), and the solid line represents the echo signal at (1200K, -200m/s). It can be found that under different temperatures and wind speeds, the two groups of echo line types are different, that is, each group of echo line types corresponds to a specific temperature and wind speed.

在操作S420,利用至少两台激光器向待测区域发射频率不同的脉冲激光。In operation S420, at least two lasers are used to emit pulse lasers with different frequencies to the area to be measured.

根据本发明的实施例,利用至少两台激光器向待测区域同时发射频率不同的脉冲激光,包括:According to an embodiment of the present invention, at least two lasers are used to simultaneously emit pulse lasers with different frequencies to the area to be measured, including:

利用第一可调谐连续激光器1和第二可调谐连续激光器4产生波长在1083nm附近的稳定且连续的种子激光;A first tunable continuous laser 1 and a second tunable continuous laser 4 are used to generate a stable and continuous seed laser with a wavelength of about 1083 nm;

利用第一脉冲激光器2和第二脉冲激光器5将种子激光进行光放大,并垂直向上发射到待测区域;The seed laser is optically amplified by the first pulse laser 2 and the second pulse laser 5, and is emitted vertically upward to the area to be measured;

其中,在第一脉冲激光器2和第二脉冲激光器5将种子激光进行光放大后,利用第一波长计3和第二波长计6实时监测接收到的一部分发射激光的频率,并反馈给第一可调谐连续激光器1和第二可调谐连续激光器4,以实现发射激光频率锁定。Among them, after the first pulse laser 2 and the second pulse laser 5 amplify the seed laser, the first wavelength meter 3 and the second wavelength meter 6 are used to monitor the frequency of a part of the received emission laser in real time, and feed back to the first tunable continuous laser 1 and the second tunable continuous laser 4 to achieve emission laser frequency locking.

在操作S430,利用望远镜7接收来自待测区域的具有亚稳态氦散射信息的回波信号。In operation S430, the telescope 7 is used to receive an echo signal having metastable helium scattering information from the region to be measured.

根据本发明的实施例,至少两台激光器向待测区域同时发射两束频率不同的脉冲激光,望远镜7可以接收到四个不同频率的具有亚稳态氦散射信息的回波信号。According to an embodiment of the present invention, at least two lasers simultaneously emit two pulse lasers with different frequencies to the area to be measured, and the telescope 7 can receive four echo signals with different frequencies containing metastable helium scattering information.

根据本发明的实施例,假设原子的速度分布服从麦克斯韦分布定律,则由于电磁波的多普勒效应,脉冲激光器发射激光,使热层大气中的氦原子发生共振荧光散射。发射一束激光可以得到两个频率的回波信号,并且两个回波信号在频率域上分别与对应1083.025nm和1083.034nm的频率对称,且强度大小与亚稳态氦1083nm总峰的后向散射截面有关。在不同的热层温度和风速下,亚稳态氦1083nm的后向散射截面会发生变化,从而使得回波信号的强度和频率不同。According to an embodiment of the present invention, assuming that the velocity distribution of atoms obeys Maxwell's distribution law, due to the Doppler effect of electromagnetic waves, a pulsed laser emits laser light to cause resonant fluorescence scattering of helium atoms in the thermosphere atmosphere. Emitting a laser beam can obtain echo signals of two frequencies, and the two echo signals are symmetrical with the frequencies corresponding to 1083.025nm and 1083.034nm in the frequency domain, and the intensity is related to the backscattering cross section of the total peak of metastable helium 1083nm. Under different thermosphere temperatures and wind speeds, the backscattering cross section of metastable helium 1083nm will change, resulting in different intensity and frequency of the echo signal.

需要说明的是,根据回波信号进行反演温度和风速时,若只发射一束激光,则包含信息量太少,误差过大,至少需要发射两束激光,获得四个频率上的回波信号。It should be noted that when inverting temperature and wind speed based on echo signals, if only one laser beam is emitted, the amount of information contained is too little and the error is too large. At least two laser beams need to be emitted to obtain echo signals at four frequencies.

图7示意性示出了根据本发明实施例的亚稳态氦峰的后向散射截面和回波信号的归一化强度示意图。FIG. 7 schematically shows a schematic diagram of the backscattering cross section of a metastable helium peak and the normalized intensity of an echo signal according to an embodiment of the present invention.

如图7所示,黑色点线和黑色点划线分别表示1083.025nm和1083.034nm对应的波长,曲线A表示亚稳态氦在1083nm处的合峰,曲线1-1-025表示第一脉冲激光器2发射的激光对应的第一个回波信号,曲线1-2-034表示第一脉冲激光器2发射的激光对应的第二个回波信号,曲线2-1-025表示第二脉冲激光器5发射的激光对应的第一个回波信号,曲线2-2-034表示第二脉冲激光器5发射的激光对应的第二个回波信号。可以发现,靠近亚稳态氦在1083nm处的合峰的峰值的回波信号对应的后向散射截面和强度较高,远离亚稳态氦在1083nm处的合峰的峰值的回波信号对应的后向散射截面和强度较低。As shown in Fig. 7, the black dotted line and the black dot-dash line represent the wavelengths corresponding to 1083.025nm and 1083.034nm, respectively. Curve A represents the combined peak of metastable helium at 1083nm. Curve 1-1-025 represents the first echo signal corresponding to the laser emitted by the first pulse laser 2. Curve 1-2-034 represents the second echo signal corresponding to the laser emitted by the first pulse laser 2. Curve 2-1-025 represents the first echo signal corresponding to the laser emitted by the second pulse laser 5. Curve 2-2-034 represents the second echo signal corresponding to the laser emitted by the second pulse laser 5. It can be found that the backscattering cross section and intensity corresponding to the echo signal close to the peak of the combined peak of metastable helium at 1083nm are higher, and the backscattering cross section and intensity corresponding to the echo signal far from the peak of the combined peak of metastable helium at 1083nm are lower.

根据本发明的实施例,经望远镜7接收到的激光雷达信号可以表达大气对激光光束的散射、吸收等效应的作用集合。利用热层氦原子的回波信号,可以建立起接收系统接收光子数和系统参数之间的联系,进行热层大气对应的温度和风速的反演。令P表示两台脉冲激光器的功率,A表示望远镜7的口径,η为接收系统的总光学效率,则共振荧光激光雷达方程可以通过公式(三)表示:According to an embodiment of the present invention, the laser radar signal received by the telescope 7 can express the effect set of the atmosphere on the scattering, absorption and other effects of the laser beam. Using the echo signal of the thermospheric helium atoms, the connection between the number of photons received by the receiving system and the system parameters can be established to invert the temperature and wind speed corresponding to the thermospheric atmosphere. Let P represent the power of the two pulsed lasers, A represent the aperture of the telescope 7, and η represent the total optical efficiency of the receiving system. Then the resonance fluorescence laser radar equation can be expressed by formula (III):

公式(三)中,z表示测量高度,Δz表示测量的空间分辨率,N(z)表示在z-Δz/2和z+Δz/2之间收集到的总光子数;τ为积分时间(s);h0为Planck常量,h≈6.62607015×10- 34J·s;c为光速,c≈299792458m/s;σeff表示每个氦原子的有效截面,λ0为发射激光的中心波长;ρ(z)为亚稳态氦的数密度;R为分支比;T为大气的单程透过率;E(z)为大气的单程消光比;NB×σ表示由于背景辐射光和暗电流造成的光子数,与天顶角SZA有关。In formula (III), z represents the measurement altitude, Δz represents the spatial resolution of the measurement, N(z) represents the total number of photons collected between z-Δz/2 and z+Δz/2; τ represents the integration time (s); h0 represents the Planck constant, h≈6.62607015× 10-34 J·s; c represents the speed of light, c≈299792458m/s; σeff represents the effective cross section of each helium atom, λ0 represents the central wavelength of the emitted laser; ρ(z) represents the number density of metastable helium; R represents the branching ratio; T represents the one-way transmittance of the atmosphere; E(z) represents the one-way extinction ratio of the atmosphere; NB ×σ represents the number of photons caused by background radiation and dark current, which is related to the zenith angle SZA.

令h表示底层亚稳态氦原子的高度,则大气的单程消光比E(z)可以通过公式(四)表示:Let h represent the height of the bottom metastable helium atom, then the one-way extinction ratio E(z) of the atmosphere can be expressed by formula (IV):

公式(四)中,z表示测量高度;ρ(z)为亚稳态氦的数密度;σeff表示每个氦原子的有效截面。可以发现,望远镜7的接收光子数和发射光子数、望远镜7口径和系统效率成正比。In formula (IV), z represents the measurement height; ρ(z) is the number density of metastable helium; σ eff represents the effective cross section of each helium atom. It can be found that the number of received photons and the number of emitted photons of telescope 7 are proportional to the aperture of telescope 7 and the system efficiency.

根据本发明的实施例,通过公式(一)和(二)确定待发射激光的两个频率与热层温度、风速的联系后,可以建立回波信号光谱数据库,每一组回波线型对应于特定温度和风速。根据模拟出的回波信号光谱数据库,可以确定斐索干涉仪10的多项参数。According to the embodiment of the present invention, after determining the relationship between the two frequencies of the laser to be emitted and the temperature and wind speed of the thermosphere through formulas (I) and (II), an echo signal spectrum database can be established, and each set of echo line types corresponds to a specific temperature and wind speed. Based on the simulated echo signal spectrum database, multiple parameters of the Fizeau interferometer 10 can be determined.

图8示意性示出了根据本发明实施例的斐索干涉仪的参数确定方法流程图。FIG8 schematically shows a flow chart of a method for determining parameters of a Fizeau interferometer according to an embodiment of the present invention.

如图8所示,斐索干涉仪的参数确定方法可以包括操作S810~S860。As shown in FIG. 8 , the method for determining parameters of a Fizeau interferometer may include operations S810 to S860 .

在操作S810,获取初始斐索干涉仪的多个属性信息,其中,多个属性信息包括:初始斐索干涉仪的多个半高全宽组成的集合、多个自由谱间距组成的集合以及滤波片9处理后的具有亚稳态氦散射信息的回波信号到达初始斐索干涉仪的入射倾角。In operation S810, multiple attribute information of the initial Fizeau interferometer is obtained, wherein the multiple attribute information includes: a set of multiple half-widths at half maximum of the initial Fizeau interferometer, a set of multiple free spectrum spacings, and an incident inclination angle of the echo signal with metastable helium scattering information processed by the filter 9 arriving at the initial Fizeau interferometer.

根据本发明的实施例,初始斐索干涉仪由分束器、准直物镜和标准平面所组成,假设在初始斐索干涉仪中心选取很窄的一小段Δy,此处的初始斐索干涉仪可以被看作传统F-P标准具。在中心位置上,初始斐索干涉仪透过率函数通过公式(五)表示:According to an embodiment of the present invention, the initial Fizeau interferometer is composed of a beam splitter, a collimating objective lens and a standard plane. Assuming that a very narrow segment Δy is selected at the center of the initial Fizeau interferometer, the initial Fizeau interferometer here can be regarded as a traditional F-P standard tool. At the center position, the transmittance function of the initial Fizeau interferometer is expressed by formula (V):

公式(五)中,TP表示斐索干涉仪的峰值透过率;v0表示斐索干涉仪透过率的中心频率;v表示激光频率;F表示斐索干涉仪的目标精细因子;FSR表示目标自由谱间距。In formula (5), TP represents the peak transmittance of the Fizeau interferometer; v0 represents the center frequency of the transmittance of the Fizeau interferometer; v represents the laser frequency; F represents the target fineness factor of the Fizeau interferometer; and FSR represents the target free spectrum spacing.

当使用初始斐索干涉仪获得后向散射信号光谱时,光谱频率范围的间隔应在初始斐索干涉仪的接收范围内,确保完整的光谱图案被线列探测器11所接收。为了减少光子损失造成的误差,每一个回波信号对应的散射光谱的大部分能量集中在一个半高全宽(FWHM)内,并且需要通过控制自由谱间距(FSR)使得在散射光谱中的透过率曲线尽可能包含所有回波信号。When the initial Fizeau interferometer is used to obtain the backscattered signal spectrum, the interval of the spectrum frequency range should be within the receiving range of the initial Fizeau interferometer to ensure that the complete spectrum pattern is received by the linear detector 11. In order to reduce the error caused by photon loss, most of the energy of the scattered spectrum corresponding to each echo signal is concentrated in a half-maximum full width (FWHM), and the free spectrum spacing (FSR) needs to be controlled so that the transmittance curve in the scattered spectrum contains all echo signals as much as possible.

在操作S820,基于仿真得到回波信号光谱数据库,遍历多个半高全宽组成的集合以及多个自由谱间距组成的集合,以确定目标半高全宽和目标自由谱间距。In operation S820, an echo signal spectrum database is obtained based on simulation, and a set consisting of a plurality of half-maximum full widths and a set consisting of a plurality of free spectrum intervals are traversed to determine a target half-maximum full width and a target free spectrum interval.

根据本发明的实施例,结合被动仿真得到的回波信号光谱数据库以及回波信号的散射信息,从多个半高全宽组成的集合以及多个自由谱间距组成的集合中选择出最符合回波信号的散射光谱曲线的一个子集,即为目标半高全宽和目标自由谱间距。优选地,目标半高全宽为0.005nm,目标自由谱间距为0.01nm。According to an embodiment of the present invention, in combination with the echo signal spectrum database obtained by passive simulation and the scattering information of the echo signal, a subset of the scattering spectrum curve that best matches the echo signal is selected from a set consisting of multiple half-height full widths and a set consisting of multiple free spectrum spacings, that is, the target half-height full width and the target free spectrum spacing. Preferably, the target half-height full width is 0.005nm and the target free spectrum spacing is 0.01nm.

在操作S830,根据目标半高全宽和目标自由谱间距,确定斐索干涉仪10的目标精细因子。In operation S830, a target fine factor of the Fizeau interferometer 10 is determined based on the target full width at half maximum and the target free spectral spacing.

在操作S840,根据目标精细因子,确定斐索干涉仪10的目标反射率。In operation S840, a target reflectivity of the Fizeau interferometer 10 is determined according to the target fine factor.

根据本发明的实施例,高的精细度会提高斐索干涉仪10的分辨率,但同时也会带来加工上的困难。因此选择合适的精细度不仅需要考虑提高斐索干涉仪10的分辨率,还要便于后续加工。According to the embodiment of the present invention, high fineness will improve the resolution of the Fizeau interferometer 10, but it will also bring difficulties in processing. Therefore, choosing a suitable fineness not only needs to consider improving the resolution of the Fizeau interferometer 10, but also facilitate subsequent processing.

斐索干涉仪10的目标精细因子和目标反射率通过公式(六)确定:The target fineness factor and target reflectivity of the Fizeau interferometer 10 are determined by formula (VI):

公式(六)中,F表示斐索干涉仪10的目标精细因子,FSR表示目标自由谱间距,FWHM表示目标半高全宽,R表示斐索干涉仪10的目标反射率。In formula (VI), F represents the target fineness factor of the Fizeau interferometer 10 , FSR represents the target free spectrum spacing, FWHM represents the target full width at half maximum, and R represents the target reflectivity of the Fizeau interferometer 10 .

根据本发明的实施例,在已知目标半高全宽和目标自由谱间距的情况下,可以得到斐索干涉仪10的目标精细因子F=2,目标反射率R=0.24。According to the embodiment of the present invention, when the target full width at half maximum and the target free spectrum spacing are known, the target fineness factor F=2 and the target reflectivity R=0.24 of the Fizeau interferometer 10 can be obtained.

在操作S850,根据目标自由谱间距对应的波长,确定斐索干涉仪10的目标平均腔长。In operation S850, a target average cavity length of the Fizeau interferometer 10 is determined according to the wavelength corresponding to the target free spectrum interval.

在操作S860,根据目标平均腔长,确定斐索干涉仪10的目标楔角;其中,目标楔角为斐索干涉仪10的两块平板对应的夹角,目标楔角用于调控斐索干涉仪10的干涉条纹的宽度。In operation S860 , a target wedge angle of the Fizeau interferometer 10 is determined according to the target average cavity length; wherein the target wedge angle is the angle between two plates of the Fizeau interferometer 10 , and the target wedge angle is used to adjust the width of the interference fringes of the Fizeau interferometer 10 .

根据本发明的实施例,假设单色平行光以θ角入射到两块夹角为α的反射率很高的平板上,一部分光直接透过两板,另一部分光在两板间经过多次反射后出射,两部分光的夹角为2nα(n为介质的折射率),二者相干并产生等厚条纹,其中,分束器、准直物镜结合两块反射率很高的平板构成初始斐索干涉仪。According to an embodiment of the present invention, it is assumed that monochromatic parallel light is incident at an angle θ onto two flat plates with a high reflectivity and an included angle α, a portion of the light directly passes through the two plates, and the other portion of the light is emitted after multiple reflections between the two plates, and the included angle between the two portions of light is 2nα (n is the refractive index of the medium), the two portions of light are coherent and produce fringes of equal thickness, wherein a beam splitter, a collimating objective lens and two flat plates with a high reflectivity constitute an initial Fizeau interferometer.

根据本发明的实施例,斐索干涉仪的透过率与斐索干涉仪楔子上入射光束的空间位置有关。因此,在入射激光的不同位置,会有不同的透过率曲线。假设望远镜7接收光束后垂直进入斐索干涉仪10中,在计算斐索干涉仪10的目标平均腔长时,可以将初始斐索干涉仪中心处的狭小范围看作一个传统F-P标准具,利用传统F-P标准具的自由光谱范围和平均腔长之间的关系可以计算出斐索干涉仪10的目标平均腔长。此外,目标楔角决定了干涉条纹的宽度。According to an embodiment of the present invention, the transmittance of the Fizeau interferometer is related to the spatial position of the incident light beam on the wedge of the Fizeau interferometer. Therefore, at different positions of the incident laser, there will be different transmittance curves. Assuming that the telescope 7 receives the light beam and enters the Fizeau interferometer 10 vertically, when calculating the target average cavity length of the Fizeau interferometer 10, the narrow range at the center of the initial Fizeau interferometer can be regarded as a traditional F-P standard tool. The target average cavity length of the Fizeau interferometer 10 can be calculated using the relationship between the free spectral range and the average cavity length of the traditional F-P standard tool. In addition, the target wedge angle determines the width of the interference fringes.

斐索干涉仪10的目标平均腔长和目标楔角通过公式(七)确定:The target average cavity length and target wedge angle of the Fizeau interferometer 10 are determined by formula (VII):

公式(七)中,L0表示斐索干涉仪10的目标平均腔长,n表示介质的折射率,λ表示发射激光的波长,λFSR表示目标自由谱间距对应的波长,α表示斐索干涉仪10的目标楔角。In formula (VII), L 0 represents the target average cavity length of the Fizeau interferometer 10 , n represents the refractive index of the medium, λ represents the wavelength of the emitted laser, λ FSR represents the wavelength corresponding to the target free spectrum spacing, and α represents the target wedge angle of the Fizeau interferometer 10 .

当发射激光的波长为1083nm,楔形空间内介质的折射率为1,根据目标自由谱间距对应的波长可以计算出目标平均腔长L0,约为0.06m;并且在已知回波信号垂直入射到初始斐索干涉仪的情况下,还可以计算出斐索干涉仪10的目标楔角为9.2μrad。When the wavelength of the emitted laser is 1083 nm, the refractive index of the medium in the wedge-shaped space is 1, and the target average cavity length L 0 can be calculated according to the wavelength corresponding to the target free spectrum spacing, which is about 0.06 m; and when it is known that the echo signal is vertically incident on the initial Fizeau interferometer, the target wedge angle of the Fizeau interferometer 10 can also be calculated to be 9.2 μrad.

根据本发明的实施例,形状因子会影响干涉条纹的形状,不同的形状因子对应不同的斐索干涉仪10的透过率曲线。因此,可以通过形状因子来验证斐索干涉仪10对应的各项参数的合理性。According to the embodiment of the present invention, the shape factor affects the shape of the interference fringes, and different shape factors correspond to different transmittance curves of the Fizeau interferometer 10. Therefore, the rationality of various parameters corresponding to the Fizeau interferometer 10 can be verified by the shape factor.

根据本发明的实施例,形状因子通过以下公式(八)确定:According to an embodiment of the present invention, the shape factor is determined by the following formula (VIII):

公式(八)中,S表示形状因子,F表示斐索干涉仪10的目标精细因子,L0表示斐索干涉仪10的目标平均腔长,α表示斐索干涉仪10的目标楔角,v表示激光频率,c表示光速。In formula (eight), S represents the shape factor, F represents the target fine factor of the Fizeau interferometer 10, L0 represents the target average cavity length of the Fizeau interferometer 10, α represents the target wedge angle of the Fizeau interferometer 10, v represents the laser frequency, and c represents the speed of light.

图9示意性示出了根据本发明实施例的不同形状因子与斐索干涉仪10的透过率曲线关系图。FIG. 9 schematically shows a relationship diagram between different shape factors and transmittance curves of the Fizeau interferometer 10 according to an embodiment of the present invention.

如图9所示,结合公式(八),当形状因子小于0.6时,斐索干涉仪10的透过率曲线可以近似为Airy函数,进而通过Airy函数形式来反演风速。若计算得到S≤0.6,则可以保持入射角θ=0的垂直入射状态;若计算得到S>0.6,则需要通过调整入射角来对干涉条纹进行调整,直到S≤0.6。将斐索干涉仪10的各项参数值带入公式,计算得到的S符合要求。As shown in FIG9 , combined with formula (VIII), when the shape factor is less than 0.6, the transmittance curve of the Fizeau interferometer 10 can be approximated as an Airy function, and then the wind speed can be inverted by the Airy function. If S≤0.6 is calculated, the vertical incidence state of the incident angle θ=0 can be maintained; if S>0.6 is calculated, the interference fringes need to be adjusted by adjusting the incident angle until S≤0.6. Substituting the various parameter values of the Fizeau interferometer 10 into the formula, the calculated S meets the requirements.

根据本发明的实施例,还可以通过斐索干涉仪10的分辨率来验证斐索干涉仪10对应的各项参数的合理性。According to the embodiment of the present invention, the rationality of various parameters corresponding to the Fizeau interferometer 10 can also be verified by the resolution of the Fizeau interferometer 10 .

使用泰勒判据,斐索干涉仪10的分辨率通过公式(九)确定:Using Taylor's criterion, the resolution of the Fizeau interferometer 10 is determined by formula (IX):

公式(九)进一步可以表示为:make Formula (IX) can be further expressed as:

公式(九)中,Δv表示斐索干涉仪10能分辨出的最小频率范围,m表示分辨级数,Δλ表示斐索干涉仪10能分辨出的最小波长范围。结合前文中计算出的平均腔长L0和目标精细因子F,可以进一步得到斐索干涉仪10的频率分辨率约1.28GHz。而亚稳态He回波信号两两之前频率间隔最小的为2.23GHz,由此可以得知前文的各项参数值是可以满足分辨出四个回波信号的要求。In formula (IX), Δv represents the minimum frequency range that the Fizeau interferometer 10 can distinguish, m represents the resolution level, and Δλ represents the minimum wavelength range that the Fizeau interferometer 10 can distinguish. Combining the average cavity length L0 and the target fine factor F calculated in the previous article, it can be further obtained that the frequency resolution of the Fizeau interferometer 10 is about 1.28 GHz. The minimum frequency interval between the metastable He echo signals is 2.23 GHz, from which it can be known that the various parameter values in the previous article can meet the requirements of distinguishing four echo signals.

在操作S440,利用滤波片9和斐索干涉仪10对望远镜7接收到的具有亚稳态氦散射信息的回波信号依次进行两级滤波处理,并将处理后的回波信号输送到接收系统对应的光路中。In operation S440, the filter 9 and the Fizeau interferometer 10 are used to perform two-stage filtering processing on the echo signal with metastable helium scattering information received by the telescope 7, and the processed echo signal is transmitted to the optical path corresponding to the receiving system.

根据本发明的实施例,滤波片9将望远镜7接收到的具有亚稳态氦散射信息的回波信号进行第一级滤波处理;超窄带斐索干涉仪10将滤波片9处理后的具有亚稳态氦散射信息的回波信号进行第二级滤波处理,其中,滤波片9处理后的具有亚稳态氦散射信息的回波信号在超窄带斐索干涉仪10中发生等厚干涉,并将发生等厚干涉后的回波信号平行光束按照不同频率进行分类以及输送到接收系统对应的光路中。According to an embodiment of the present invention, the filter 9 performs a first-stage filtering process on the echo signal with metastable helium scattering information received by the telescope 7; the ultra-narrowband Fizeau interferometer 10 performs a second-stage filtering process on the echo signal with metastable helium scattering information processed by the filter 9, wherein the echo signal with metastable helium scattering information processed by the filter 9 undergoes equal-thickness interference in the ultra-narrowband Fizeau interferometer 10, and the parallel light beams of the echo signal after the equal-thickness interference are classified according to different frequencies and transmitted to the optical path corresponding to the receiving system.

图10示意性示出了根据本发明实施例的数据处理流程图。FIG. 10 schematically shows a data processing flow chart according to an embodiment of the present invention.

如图10所示,数据处理过程包括步骤S450~步骤S470。As shown in FIG. 10 , the data processing process includes steps S450 to S470 .

在操作S450,探测器11接收来自接收系统光路中的回波信号,其中,接收系统光路中的回波信号通过测量光谱表示,测量光谱为散射光谱和斐索干涉仪10的透过率函数的卷积。In operation S450 , the detector 11 receives an echo signal from an optical path of the receiving system, wherein the echo signal in the optical path of the receiving system is represented by a measurement spectrum, and the measurement spectrum is a convolution of a scattering spectrum and a transmittance function of the Fizeau interferometer 10 .

测量光谱可以通过公式(十)确定:The measured spectrum can be determined by formula (10):

公式(十)中,Sconv表示测量光谱,Sscatter表示散射光谱,SFizeau表示斐索干涉仪的透过率函数。In formula (10), S conv represents the measured spectrum, S scatter represents the scattered spectrum, and S Fizeau represents the transmittance function of the Fizeau interferometer.

图11示意性示出了根据本发明实施例的探测器探测到的激光线型、总的回波信号以及斐索干涉仪的卷积结果示意图。FIG11 schematically shows a schematic diagram of the laser line shape, the total echo signal and the convolution result of the Fizeau interferometer detected by the detector according to an embodiment of the present invention.

如图11所示,假设激光雷达系统发射的激光为高斯线型,线宽为200MHz,温度为850K以及风速为0m/s。在激光线型已知的情况下,根据探测器11的各个通道接收到的光子数,可以得到散射光谱和斐索干涉仪10的透过率函数的卷积结果。As shown in Figure 11, it is assumed that the laser emitted by the laser radar system is a Gaussian line shape, the line width is 200MHz, the temperature is 850K and the wind speed is 0m/s. When the laser line shape is known, according to the number of photons received by each channel of the detector 11, the convolution result of the scattering spectrum and the transmittance function of the Fizeau interferometer 10 can be obtained.

在操作S460,将测量光谱和斐索干涉仪10的透过率函数进行反卷积处理,得到激光的回波信号谱。In operation S460, the measured spectrum and the transmittance function of the Fizeau interferometer 10 are deconvolved to obtain an echo signal spectrum of the laser.

在操作S470,将激光的回波信号谱与回波信号光谱数据库进行线性插值比对,以确定激光的回波信号谱中最接近的一组数据为待测区域对应的温度和风速,其中,待测区域对应为热层大气层。In operation S470, the laser echo signal spectrum is linearly interpolated and compared with the echo signal spectrum database to determine that the closest set of data in the laser echo signal spectrum is the temperature and wind speed corresponding to the area to be measured, wherein the area to be measured corresponds to the thermosphere atmosphere.

根据本发明的实施例,基于斐索干涉仪的氦激光雷达双频测温测风方法,通过在发射端引入至少两台激光器,并改造接收端的部分模块,一方面节省了传统激光轮流发射的时间,还可以减少激光雷达的复杂度和维护难度;另一方面也可以节省扫频时间,持续累积光子信号,同时还极大地简化了数据处理的方法。最后通过采集待测区域的亚稳态氦散射信息进行测温测风,可以为后续研究测温测风提供一个新的研究方向,也填补了氦激光雷达在测温测风方面的理论空白,从而能够对热层大气物理过程和多圈层耦合机制进行更为全方位的观测和研究。According to an embodiment of the present invention, the dual-frequency temperature and wind measurement method of helium laser radar based on Fizeau interferometer, by introducing at least two lasers at the transmitting end and modifying some modules at the receiving end, saves the time of traditional laser emission in turns, and can also reduce the complexity and maintenance difficulty of laser radar; on the other hand, it can also save the frequency sweep time, continuously accumulate photon signals, and greatly simplify the data processing method. Finally, by collecting metastable helium scattering information of the area to be measured for temperature and wind measurement, a new research direction can be provided for subsequent research on temperature and wind measurement, and the theoretical gap in helium laser radar in temperature and wind measurement can be filled, so that more comprehensive observation and research can be carried out on the physical process of the thermosphere atmosphere and the multi-layer coupling mechanism.

以上所述的具体实施例,对本发明的目的、技术方案和有益效果进行了进一步详细说明,应理解的是,以上所述仅为本发明的具体实施例而已,并不用于限制本发明,凡在本发明的精神和原则之内,所做的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。The specific embodiments described above further illustrate the objectives, technical solutions and beneficial effects of the present invention in detail. It should be understood that the above description is only a specific embodiment of the present invention and is not intended to limit the present invention. Any modifications, equivalent substitutions, improvements, etc. made within the spirit and principles of the present invention should be included in the protection scope of the present invention.

Claims (6)

1. The parameter determining method of the Fizeau interferometer is applied to a helium laser radar double-frequency temperature and wind measuring system based on the Fizeau interferometer, and the helium laser radar double-frequency temperature and wind measuring system comprises the following components:
a transmitting system, comprising:
At least two lasers, which are suitable for emitting pulse lasers with different frequencies to the region to be measured;
a receiving system, comprising:
a telescope (7) adapted to receive echo signals from the region under test having metastable helium scattering information;
the filter (9) is suitable for carrying out first-stage filtering processing on the echo signal with metastable helium scattering information received by the telescope (7);
A Fizeau interferometer (10) adapted to perform a second stage filtering process on the echo signal with metastable helium scattering information processed by the filter (9);
-a detector (11) adapted to receive a scattered echo signal from the fizeau interferometer (10), wherein the echo signal from the fizeau interferometer (10) is represented by a measured spectrum, the measured spectrum being a convolution of the scattered spectrum and a transmittance function of the fizeau interferometer (10);
The data processing system is suitable for processing echo signals received by the detector (11) to obtain wind speed and temperature corresponding to the region to be detected, wherein the linear interpolation comparison result of an echo signal spectrum library is obtained according to the echo signal spectrum and simulation of the laser, the wind speed and the temperature corresponding to the region to be detected are determined, and the echo signal spectrum of the laser is the deconvolution of the measurement spectrum and the transmittance function of the Fizeau interferometer (10);
The method comprises the following steps:
Acquiring a plurality of attribute information of an initial fizeau interferometer (10), wherein the plurality of attribute information comprises: the initial Fizeau interferometer (10) comprises a set of multiple full width at half maximum components and a set of multiple free spectrum spacing components, and the echo signal with metastable helium scattering information processed by the filter (9) reaches the incidence dip angle of the initial Fizeau interferometer (10);
traversing the set of multiple full width at half maximum and the set of multiple free spectrum spacing based on the simulation to obtain an echo signal spectrum database so as to determine a target full width at half maximum and a target free spectrum spacing;
Determining a target fine factor of the Fizeau interferometer (10) according to the target full width at half maximum and the target free spectrum spacing;
determining a target reflectivity of the fizeau interferometer (10) from the target fine factor;
Determining a target average cavity length of the Fizeau interferometer (10) according to the wavelength corresponding to the target free spectrum interval;
Determining a target wedge angle of the fizeau interferometer (10) according to the target average cavity length;
the target wedge angle is an included angle corresponding to two flat plates of the Fizeau interferometer (10), and is used for regulating and controlling the width of interference fringes of the Fizeau interferometer (10).
2. The determination method of claim 1, wherein the target fine factor of the fizeau interferometer (10) and the target reflectivity of the fizeau interferometer (10) are calculated by formula (1);
In equation (1), F represents a target fine factor of the fizeau interferometer (10), FSR represents the target free spectral separation, FWHM represents the target full width at half maximum, and R represents a target reflectivity of the fizeau interferometer (10).
3. The determination method of claim 1, wherein the target average cavity length and target wedge angle of the fizeau interferometer (10) are calculated by formula (2);
In the formula (2), L 0 represents a target average cavity length of the fizeau interferometer (10), n represents a refractive index of a medium, λ represents a wavelength of the emitted laser, λ FSR represents a wavelength corresponding to the target free spectrum interval, and α represents a target wedge angle of the fizeau interferometer (10).
4. The determination method according to claim 1, wherein the emission system includes the at least two lasers, including:
A first tunable continuous laser (1) and a second tunable continuous laser (4) adapted to generate a stable and continuous seed laser;
a first pulse laser (2) and a second pulse laser (5) adapted to optically amplify the seed laser and emit it to the area to be measured;
The first pulse laser (2) and the second pulse laser (5) transmit part of the amplified emission laser to a first wavelength meter (3) and a second wavelength meter (6) to monitor the frequency of the emission laser in real time and feed back the frequency to the first tunable continuous laser (1) and the second tunable continuous laser (4) so as to realize the emission laser frequency locking.
5. The determination method according to claim 4, wherein the first tunable continuous laser (1) and the second tunable continuous laser (4) are identical in structure, the first pulse laser (2) and the second pulse laser (5) are identical in structure, and the first wavelength meter (3) and the second wavelength meter (6) are identical in structure.
6. The determination method according to claim 1, wherein the at least two lasers emit the laser light vertically upward.
CN202211544782.6A 2022-12-02 2022-12-02 Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer Active CN115902834B (en)

Priority Applications (2)

Application Number Priority Date Filing Date Title
CN202211544782.6A CN115902834B (en) 2022-12-02 2022-12-02 Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer
PCT/CN2023/072058 WO2024113472A1 (en) 2022-12-02 2023-01-13 Helium laser radar dual-frequency temperature and wind measurement system and method based on fizeau interferometer

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202211544782.6A CN115902834B (en) 2022-12-02 2022-12-02 Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer

Publications (2)

Publication Number Publication Date
CN115902834A CN115902834A (en) 2023-04-04
CN115902834B true CN115902834B (en) 2024-05-17

Family

ID=86479332

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202211544782.6A Active CN115902834B (en) 2022-12-02 2022-12-02 Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer

Country Status (2)

Country Link
CN (1) CN115902834B (en)
WO (1) WO2024113472A1 (en)

Families Citing this family (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN118566942B (en) * 2024-07-30 2024-09-27 中国海洋大学 Horizontal wind field detection method, device, equipment and readable storage medium

Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2648915A1 (en) * 1989-06-26 1990-12-28 Centre Nat Rech Scient DEVICE FOR MEASURING THE SPEED OF THE WIND AT AVERAGE ALTITUDE
JPH06185997A (en) * 1992-10-19 1994-07-08 Nikon Corp Fizeau interferometer
WO2004092767A1 (en) * 2003-04-10 2004-10-28 Eads Deutschland Gmbh Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system
CN103713293A (en) * 2013-12-26 2014-04-09 西安理工大学 All-fiber Doppler lidar wind field detection system and method
CN107192380A (en) * 2016-03-15 2017-09-22 莱卡地球系统公开股份有限公司 Laser tracker with two measurement functions
EP3428684A1 (en) * 2017-07-12 2019-01-16 Airbus Defence and Space GmbH Lidar arrangement and lidar method
CN110018497A (en) * 2019-03-07 2019-07-16 中国科学技术大学 Based on the atmosphere thermometric laser radar and atmosphere temp measuring method converted under frequency
CN110187362A (en) * 2019-05-24 2019-08-30 中国科学技术大学 A dual-frequency wind-measuring lidar with simultaneous operation of ultraviolet and infrared
CN112269033A (en) * 2020-10-15 2021-01-26 南京恒高光电研究院有限公司 A distributed laser Doppler anemometer and distributed measurement method based on Fizeau's interference principle
CN114994711A (en) * 2022-08-05 2022-09-02 南京信息工程大学 Laser radar based on Fizeau interferometer
CN115308774A (en) * 2022-08-11 2022-11-08 盐城师范学院 Rayleigh scattering Doppler laser radar system based on four-channel double-pass F-P interferometer

Family Cites Families (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
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
CN210572755U (en) * 2019-07-30 2020-05-19 南京信息工程大学 A lidar system for measuring atmospheric temperature based on Doppler broadening
CN113885048B (en) * 2021-10-18 2024-10-22 南京信息工程大学 Double-edge Rayleigh laser radar system based on iodine molecule absorption tank

Patent Citations (11)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
FR2648915A1 (en) * 1989-06-26 1990-12-28 Centre Nat Rech Scient DEVICE FOR MEASURING THE SPEED OF THE WIND AT AVERAGE ALTITUDE
JPH06185997A (en) * 1992-10-19 1994-07-08 Nikon Corp Fizeau interferometer
WO2004092767A1 (en) * 2003-04-10 2004-10-28 Eads Deutschland Gmbh Method for detecting wind speeds using a doppler-lidar system, especially on aeroplanes, and doppler-lidar system
CN103713293A (en) * 2013-12-26 2014-04-09 西安理工大学 All-fiber Doppler lidar wind field detection system and method
CN107192380A (en) * 2016-03-15 2017-09-22 莱卡地球系统公开股份有限公司 Laser tracker with two measurement functions
EP3428684A1 (en) * 2017-07-12 2019-01-16 Airbus Defence and Space GmbH Lidar arrangement and lidar method
CN110018497A (en) * 2019-03-07 2019-07-16 中国科学技术大学 Based on the atmosphere thermometric laser radar and atmosphere temp measuring method converted under frequency
CN110187362A (en) * 2019-05-24 2019-08-30 中国科学技术大学 A dual-frequency wind-measuring lidar with simultaneous operation of ultraviolet and infrared
CN112269033A (en) * 2020-10-15 2021-01-26 南京恒高光电研究院有限公司 A distributed laser Doppler anemometer and distributed measurement method based on Fizeau's interference principle
CN114994711A (en) * 2022-08-05 2022-09-02 南京信息工程大学 Laser radar based on Fizeau interferometer
CN115308774A (en) * 2022-08-11 2022-11-08 盐城师范学院 Rayleigh scattering Doppler laser radar system based on four-channel double-pass F-P interferometer

Non-Patent Citations (1)

* Cited by examiner, † Cited by third party
Title
亚稳态氦激光雷达的研制进展;赵若灿 等;红外与激光工程;20201130;第49卷(第S2期);1-6 *

Also Published As

Publication number Publication date
CN115902834A (en) 2023-04-04
WO2024113472A1 (en) 2024-06-06

Similar Documents

Publication Publication Date Title
Eloranta High spectral resolution lidar
Rezus et al. Single-photon spectroscopy of a single molecule
CN105911020B (en) Method for simultaneously measuring multi-component gas based on cavity ring-down spectroscopy
Salzmann et al. The LIDAR Thomson scattering diagnostic on JET
CN108303706B (en) An aerosol optical parameter detection method and hyperspectral lidar detection system
CN110794387A (en) A radiometric calibration method for airborne hyperspectral imaging lidar system
CN110441792B (en) Rayleigh scattering laser radar system capable of measuring wind and temperature simultaneously and related calibration method
CN105334519B (en) More atmospheric parameters based on triple channel F-P etalons while detected with high accuracy laser radar system
CN207882443U (en) A kind of EO-1 hyperion Airborne Lidar examining system
CN109990843B (en) Method and device for monitoring flight speed and environment parameters of aircraft
WO2021103716A1 (en) Device and method for measuring airborne hyperspectral imaging laser radar spectrum in real time
CN104777487A (en) Atmospheric aerosol optical property measuring method and laser radar system
CN107193015A (en) Ultraviolet three frequencies high spectral resolution lidar system and its detection method based on F P etalons
CN115902834B (en) Helium laser radar dual-frequency temperature and wind measurement system and method based on Fizeau interferometer
CN111239072B (en) Method for accurately measuring temperature of combustion gas
CN112859112B (en) Wind temperature detection laser radar and method based on rotating Raman-Doppler mechanism
Sun et al. Accuracy analysis of the Fabry–Perot etalon based doppler wind lidar
CN106249252A (en) Detect airborne near-infrared laser radar system and the inversion method of subcooled water in cloud
CN1159598C (en) Dual-wavelength high-altitude detection lidar
Dehring et al. Performance and comparison of 532-nm and 355-nm groundwinds lidars
CN110031427B (en) A kind of scanning detection method of ambient gas and lidar
CN112904308B (en) Laser radar system and method for detecting cloud phase state and cloud water content
CN115267817A (en) High spectral resolution temperature measurement lidar system based on double cascade F-P interferometer and multi-longitudinal mode laser
CN115308774A (en) Rayleigh scattering Doppler laser radar system based on four-channel double-pass F-P interferometer
Bukreev et al. Excimer-laser-based lidar system for stratospheric and tropospheric ozone measurements

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
SE01 Entry into force of request for substantive examination
SE01 Entry into force of request for substantive examination
GR01 Patent grant
GR01 Patent grant