-
直接探测系统可利用F-P标准具和双边缘检测技术对大气回波信号进行频率鉴别。双边缘技术通过将入射激光频率锁定在F-P标准具两个信号通道透过率曲线的交叉点来准确响应激光频率变化。两个通道的透过率曲线中心频率分开,频谱分布相同,分别位于Rayleigh 散射信号两侧。当没有多普勒频移时,两个通道输出光能量相同;当回波信号产生多普勒频移,两个通道的光能量一个增加,一个减小,可从回波信号的能量变化中提取出多普勒频移信息,进而反演大气风速。当仅考虑Rayleigh散射时,可将两个信号通道光能量的变化定义为响应函数[12]:
$$ {{R}}(v) = \frac{{{T_{R1}}(v) - {T_{R2}}(v)}}{{{T_{R1}}(v) + {T_{R2}}(v)}} $$ (1) 式中:
${T_{R1}}(v) $ 和${T_{R2}}(v) $ 分别为两个通道Rayleigh散射的鉴频函数,可由F-P标准具透过率方程和Rayleigh散射谱卷积计算得到。多普勒频移可表示为:$$ v_{d}=\frac{R\left(v_{0}+v_{d}\right)-R\left(v_{0}\right)}{R\left(v_{0}\right) \varTheta\left(v_{0}\right)} $$ (2) 式中:v0为出射波长对应激光频率。根据
$v_{r}= \lambda v_{d} / 2 $ ($\lambda$ 为中心波长),可得径向风速vr为:$$v_{r}=\frac{\lambda}{2} \frac{R\left(v_{0}+v_{d}\right)-R\left(v_{0}\right)}{R\left(v_{0}\right) \varTheta\left(v_{0}\right)}$$ (3) 在风速测量之前,需要通过扫描F-P标准具的透过率曲线确定响应函数R,然后根据实测两个通道能量变化得到的R值反解得到多普勒频率。
-
该星载测风激光雷达系统设计搭建于工作在400 km轨道高度的太阳同步轨道卫星平台上,根据应用需求,其直接探测系统的参数指标见表1。
表 1 系统参数
Table 1. System parameters
Parameters Values Detecting range/km 2-30 Satellite altitude/km 400 Range resolution 1 km @ 2-16 km
2 km @ 16-30 kmWind speed rang/m·s−1 0-60 Time resolution/s 30 Wavelength/nm 354.7 Pulse energy/J 0.18 Repetition/Hz 100 Telescope diameter/mm 650 Quantum efficiency 0.4 -
所设计的星载多普勒测风激光雷达接收机采用收发合置的结构,光路原理如图1所示。1064.1 nm的单纵模连续种子光注入主激光器,经放大倍频后发出354.7 nm线偏振脉冲光。考虑到激光频率抖动会对风速测量造成误差,直接探测技术要求在测量大气回波信号之前,先测量零频参考信号,即由激光器直接分出的零多普勒频移信号。出射激光先经过分束片分束,小部分反射光进入积分球,用于零频校准。分束片的大部分透射光经扩束准直后,分别经过偏振分光棱镜、四分之一波片和扩束望远镜进入大气中,和大气分子相互作用。其中偏振分光棱镜和四分之一波片组成“光开关”,将线偏光转化为圆偏光射入大气中。回波信号再次通过四分之一波片后又变成线偏光,且偏振方向和出射时垂直,在偏振分光棱镜处发生反射,从而实现光学接收机的收发合置。反射光经过后续光学元件后,变成圆偏光,再由光纤耦合进入接收机的鉴频部分。进入鉴频部分的圆偏光先准直为平行光,再被非偏振分光棱镜等分成两束分别进入F-P标准具的两个信号通道,最终出射信号汇聚到光电探测器端面上进行能量检测。其中鉴频部分包括透镜,非偏振分光棱镜和F-P标准具。干涉滤光片用于滤除信号中的背景噪声。所用的耦合光纤为芯径100 μm的多模光纤,限制望远镜视场的同时,采用光纤传输可以根据结构设计需求对接收机内鉴频部分的位置进行灵活布置。
-
F-P标准具作为直接测风雷达接收机中的核心部件,其参数的选取反映了透过率曲线的陡峭程度,直接影响探测的信噪比和误差。
-
自由谱间距(FSR)为透过率曲线连续两个波峰之间的宽度。一方面,为了不引起测量误差,需将大部分Rayleigh散射信号能量集中在一个FSR内。在波长354.7 nm,大气温度250 K的条件下,Rayleigh散射的多普勒展宽可描述为高斯线型,其带宽
$\delta= 1.51 \;\mathrm{GHz}$ 。通常认为99.73%的Rayleigh信号都集中在$6 \delta=9.06 \;\mathrm{GHz}$ 的频率范围内,因此FSR不能小于该值[13]。但同时FSR也不能过大,因为在F-P标准具的有效精细度保持不变的条件下,随着FSR的增大,透过率曲线带宽增加,会导致测量灵敏度降低,从而影响测量精度。另一方面,该激光雷达接收机根据卫星运行引起的多普勒频移进行了针对性设计。卫星轨道高度为400 km,为了测量大气水平风速,设计雷达观测方向与雷达正面方向成45°,同时与垂直地球表面方向成45°。在最优的卫星轨道设计方案下,卫星与地球自转的相对速度变化为7.77 km/s,根据探测方向相对卫星平台的指向,相对速度带来的径向分量为3.88 km/s,对应频移。利用标准具透过率曲线的周期性,设计FSR的值为 (n取正整数),可消除卫星运行引起多普勒频移对测量的影响。综合考虑,设计标准具的FSR为10.925 GHz,取卫星运行带来的多普勒平移的1/2(如图2所示)。
-
系统采用固定腔长的F-P标准具,以种子注入式脉冲激光器作为激光光源,通过调谐种子波长实现F-P标准具透过率曲线的扫描。通常激光器可调谐范围有限,且连续调谐过程中经常出现激光器跳模的现象,会导致透过率曲线的测量不准确。设计两信号通道之间的峰值间距为FSR的一半,此时两信号通道的透过率曲线关于交叉点对称,只需扫描1/2个FSR即可获得完整的透过率曲线,降低了对连续可调谐脉冲激光器稳定性和波长调谐范围的要求。
根据双边缘技术的理论,风速测量误差可表示为:
$$\varepsilon=\frac{1}{\varTheta \cdot S N R}$$ (4) 式中:SNR是测量信噪比,表示为:
$${SNR}=\frac{1}{\sqrt{\left(S N R_{1}\right)^{-2}+\left(S N R_{2}\right)^{-2}}}$$ (5) 式中:SNR1和SNR2分别为两个信号通道的信噪比。若采用光子计数模式,可将SNR表示为:
$$S N R_{i}=\frac{N_{i}}{\sqrt{N_{i}+N_{b i}+N_{d i}}}$$ (6) 式中:SNRi为第i个通道的信噪比;Ni为第i个通道接收的Rayleigh散射信号光子数;Nbi为第i个通道接收的太阳背景光子数;Ndi为第i个通道的探测器暗计数。在标准具透过率曲线的FSR和峰值间隔确定的情况下,风速测量误差为谱宽v1/2的单值函数。结合激光雷达方程,经仿真计算得到系统风速误差和谱宽之间的关系如图3 (a)所示。根据误差最小原理,确定谱宽v1/2=1.64 GHz。此时系统信噪比随探测高度的变化如图3 (b)所示,图中信噪比曲线出现突变是由于16 km处距离分辨率变化导致的。综上,F-P标准具的频率参数设计值为:FSR=4.6 pm,∆v=2.3 pm,v1/2=0.7 pm。
-
口径作为标准具的一个关键参数,其值也反映了光学接收机中空间光路的接收口径大小,并需要和望远镜系统、透镜参数以及光纤参数相匹配。入射到F-P标准具的光束由光纤耦合得到,耦合光纤为芯径100 μm,数值孔径为0.22的多模光纤。标准具对输入光的发散角非常敏感,通常控制入射到标准具的光束发散角小于1 mrad[14]。为了保证信号的接收效率,空间光路的接收孔径应满足公式(7):
$$\varOmega_{r} {A}_{r} \geqslant \varOmega_{f} A_{f}$$ (7) 式中:
$\varOmega_{r},\varOmega_{f}$ 分别为接收光路和光纤的接收立体角;Ar,Af分别为接收光路和光纤的通光面积。进一步简化上式得到接收光路的通光口径dr应满足:$$d_{r} \geqslant \frac{2 \arcsin ({ N A}) d_{f}}{2 \theta_{0}}$$ (8) 式中:
$2 \theta_{0}$ 为标准具最大接收角,取$2 \theta_{0}=2 \;\mathrm{mrad}$ ;NA为光纤数值孔径;df为光纤芯径。经计算,接收光路的通光口径为22 mm,考虑安装余量,取标准具的口径为25 mm。腔长为F-P标准具两平板之间的距离,采用在同一个基板上沉积镀膜的方式形成两个不同长度的腔长作为标准具的两个通道。根据干涉极大产生条件2 nd=mc/vm,腔长d表示为:
$${d}=\frac{c}{2 n\left(v_{m+1}-v_{m}\right)}=\frac{c}{2 n \cdot {\rm FSR}}$$ (9) 式中:n为平板间介质折射率;m为干涉级次。根据两反射面之间的光程是否可调整将标准具分为固定式和可调谐式。考虑卫星平台震动会给系统带来测量误差,文中介绍的直接探测接收系统中采用固定腔长式的F-P标准具,通道腔长为13.71 mm,两通道腔长差为88.7 nm。
F-P标准具作为高精度光学鉴频器是直接探测激光雷达中的关键部件,它的加工、装调和调节精度要求都非常高,表面平整度要求为λ/100。温度变化会对标准具的腔长及介质折射率产生影响,从而引起频率漂移。温度变化1 ℃,会引起标准具透过率谱线平移80 MHz,对应于径向风速误差为15 m/s。平板间材料选用零膨胀微晶玻璃,其热膨胀系数在10−7/K的量级。工作时,需将标准具放置于温度控制精度为0.1 K的热控箱内以保证测量性能[15]。
-
为了验证接收机的性能,对搭建完成的接收机进行透过率曲线的扫描,测试试验装置如图6 (a)所示。1064.1 nm种子光注入主激光器经放大倍频后发出354.7 nm的脉冲光,通过调谐种子波长来改变激光器的输出波长,调谐步长为0.12 pm。从光路中分出小部分光接入波长计实时记录波长数据。试验测得波形如图6 (b)所示,此时标准具的一个通道能量达到最强,另一个通道能量最弱,当调节波长时,两个通道能量强弱会交替变化。进一步对测得波形数据进行采样、积分就得到对应当前波长下两个通道的能量。波长扫描一个FSR后就得到一组完整的透过率散点图。
F-P标准具对单一频率的光的透过率hv级数展开形式为:
$$\begin{split} {h}(v)=&T_{p}\left\{1+\sum_{n=1}^{\infty} R_{e}^{n} \cos \left[\frac{2 \pi n v}{\rm F S R} \frac{1+\cos \theta_{0}}{2}\right] \right.\cdot \\ &\left.{\sin\rm c}\left[\frac{2 n v}{\rm F S R} \frac{1-\cos \theta_{0}}{2}\right]\right\} \end{split}$$ (10) 式中:Tp为峰值透过率;θ0为入射光发散角;Re为标准具内表面有效反射率。将透过率连续扫描散点图和公式(10)进行拟合,拟合结果如图7所示。取F-P标准具参数为多次实测结果的平均值,得到FSR为4.6 pm,带宽为0.91 pm,实测的通道带宽大于设计值0.7 pm。由于装调和加工存在误差,入射光无法严格满足平行光正入射的条件,所以实际带宽不可能严格等于设计值,通常带宽超过理论值约0.2 pm是可以接受的。另外两个通道间隔分别为2.5 pm和2.1 pm,和设计值FSR的一半之间存在偏差,可能加工过程中两个通道镀膜工艺存在误差,导致通道腔长差不严格等于设计值。
在扫描过程中,激光器存在频率和能量的抖动现象,会导致测量不稳定,影响透过率曲线扫描的准确性。计算多次连续扫描的透过率曲线的均值和方差,结果如图8所示,可以看出:波长扫描得到的透过率曲线稳定性较好,测得两个通道的能量波动的最大方差分别为0.033和0.031。
根据实际测量得到的不同带宽,对光学接收机的测风性能进行了仿真模拟,测风误差随风速变化情况如图9所示。可以看出在设计带宽0.7 pm处,最大风速的测风误差为2.82 m/s。取平均测量带宽0.91 pm时,最大风速误差为2.94 m/s。图10进一步分析了带宽增宽对测风精度的影响,可以看出:带宽增宽越大,测风误差增长越快。在30 km高度处,当带宽从设计值0.7 pm增宽0.43 pm,即实际带宽为1.13 pm时,对应最大风速的测量误差增加1 m/s。所以为了保证测风精度,加工和装调完成的接收机的实际带宽不应超过1.13 pm。
Small scale optical receiver of spaceborne Doppler wind lidar
-
摘要: 星载测风激光雷达可以提供全球范围高实时性、高精度、高分辨率的大气风场信息,已经被认为是解决全球化风场观测的最佳手段。我国也在积极开展星载多普勒测风激光雷达相关研究工作。针对400 km高度的卫星轨道,设计并研制了一套多普勒直接测风激光雷达光学接收机,结合双边缘检测原理和星载平台相关技术参数,对Fabry-Perot标准具的主要参数进行设计。为了满足星载平台对稳定性和小型化的需求,接收机中主要光学元件之间采用分子粘接方式紧密连接。整个光学接收机集成在450 mm×300 mm×80 mm密闭箱体中,内部光学元件采用倒插方式沉入接收机壳体的凹槽内,整体结构稳定可靠,集成度高。通过改变激光波长的方式扫描Fabry-Perot标准具的透过率曲线,对所研制的接收机进行了性能测试。并由透过率曲线的实测参数对接收机的测风性能进行仿真,仿真结果显示在30 km处的最大风速误差为2.94 m/s。并进一步分析了接收机带宽增宽对测风精度的影响,分析结果显示带宽偏差为0.43 pm时,会引起1 m/s的风速误差增量。Abstract: Spaceborne wind lidar is considered to be the best means of wind field observation due to its real-time, high-precision and high resolution. Researches on spaceborne Doppler wind lidar have been actively carried out in China. The direct detection receiver of Doppler wind lidar developed for the satellite orbit with an altitude of 400 km is introduced. The receiver employs the dual Fabry-Perot etalons as the discriminating frequency unit and its parameters are specially designed for the spaceborne application. The main optical elements of the receiver are connected by molecular force and reversely inserted into the fillister of the receiver shell. The structure of the receiver is stable, reliable and highly integrated with the size of 450 mm×300 mm×80 mm, adapting to the spaceborne requirements of stability and small scale. The Fabry-Perot etalon transmission curves are obtained by tuning the wavelength of laser and the parameters of transmission curves are tested and analyzed. The wind measurement performance of the optical receiver is simulated according to the actual measured parameters. The simulation results show that the maximum wind velocity error is 2.94 m/s at the height of 30 km. Then the influence of the receiver bandwidth on measurement accuracy is further analyzed. The analysis results show that when the bandwidth deviation is 0.43 pm, the measurement error increment of 1m/s is caused.
-
Key words:
- spaceborne lidar /
- Fabry-Perot etalon /
- Doppler wind measurement /
- dual edge detection
-
表 1 系统参数
Table 1. System parameters
Parameters Values Detecting range/km 2-30 Satellite altitude/km 400 Range resolution 1 km @ 2-16 km
2 km @ 16-30 kmWind speed rang/m·s−1 0-60 Time resolution/s 30 Wavelength/nm 354.7 Pulse energy/J 0.18 Repetition/Hz 100 Telescope diameter/mm 650 Quantum efficiency 0.4 -
[1] Liu Z S, Chen Z, Yu C R, et al. Doppler wind lidar: From vehicle-mounted to space-borne [J]. Journal of Atmospheric and Environment Optics, 2015, 10(2): 126-138. (in Chinese) [2] Yan Z A, Hu X, Guo W J, et al. Near space Doppler lidar tech-niques and applications [J]. Infrared and Laser Engineering, 2021, 50(3): 20210100. (in Chinese) doi: 10.3788/IRLA20210100 [3] Kavaya M J, Spiers G D, Frehlich R G. Potential pitfalls related to space-based lidar remote sensing of the Earth with an emphasis on wind measurement[C]// Proceedings of SPIE, 2001, 4153: 385-393. [4] Guo J P, Liu B, Gong W, et al. Technical note: First comparison of wind observations from ESA ’s satellite mission Aeolus and ground-based radar wind profiler network of China [J]. Atmospheric Chemistry and Physics, 2021, 21(4): 2945-2958. doi: 10.5194/acp-21-2945-2021 [5] Reitebuch O, Lemmerz C, Paffrath U, et al. The airborne demonstrator for the direct detection Doppler lidar ALADIN on ADM-Aeolus. Part I: Instrument design and comparison to satellite instrument [J]. Journal of Atmospheric and Oceanic Technology, 2009, 26(12): 2516-2530. doi: 10.1175/2009JTECHA1314.1 [6] Marksteiner U, Lemmerz C, Lux O, et al. Calibrations and wind observations of an airborne direct detection wind lidar supporting ESA’s Aeolus mission [J]. Remote Sensing, 2018, 10(12): 2056. doi: 10.3390/rs10122056 [7] Benjamin W, Christian L, Alexander G, et al. First validation of Aeolus wind observations by airborne Doppler wind lidar measurements [J]. Atmospheric Measurement Techniques, 2020, 13(5): 2381-2396. doi: 10.5194/amt-13-2381-2020 [8] Wood S A, Emmitt G D, Greco S, et al. DLSM: a coherent and direct detection lidar simulation model for simulating space-based and aircraft-based lidar winds[C]//Proceedings of SPIE, 2000, 4035: 2-12. [9] Emmitt G D. Combining direct and coherent detection for Doppler wind lidar[C]//Proceedings of SPIE, 2004, 5575: 31-37. [10] Marx C T, Gentry B, Jordan P, et al. Lab demonstration of the hybrid Doppler wind lidar (HDWL) transceiver[C]//Proceedings of SPIE, 2013, 8872: 887207. [11] Chen B L, Yang Z D, Min M, et al. Application requirements and research progress of spaceborne Doppler wind lidar [J]. Laser & Optoelectronics Progress, 2020, 57(19): 190003. (in Chinese) doi: 10.3788/LOP57.190003 [12] Wang J J, Shu Z F, Dou X K, et al. Research of the frequency response function of the Rayleigh Doppler wind lidar [J]. Infrared and Laser Engineering, 2012, 41(9): 2364-2369. (in Chinese) doi: 10.3969/j.issn.1007-2276.2012.09.021 [13] Wang G C, Sun D S, Duan L F, et al. Analysis and design of Fabry-Perot etalon of Doppler wind lidar [J]. Acta Optica Sinica, 2011, 31(3): 0301001. (in Chinese) [14] Shen F H, Sun D S, Liu C L, et al. Single solid F-P etalon-based dual-frequency Doppler lidar [J]. Infrared and Laser Engineering, 2013, 42(11): 2944-2950. (in Chinese) doi: 10.3969/j.issn.1007-2276.2013.11.014 [15] Sun J, Feng Y T, Bai L Q, et al. Design of thermal stable Fabry-Perot etalon for measuring wind [J]. Optics and Precision Engineering, 2013, 21(5): 1167-1173. (in Chinese) doi: 10.3788/OPE.20132105.1167