西安地区大气相干长度的实验测量及分析 下载: 1051次
1 引言
激光通信具有低成本、高速率、部署灵活和强抗干扰能力等优势,因此备受关注。激光在传输过程中受到大气湍流的影响,导致激光波阵面发生畸变[1]、光束漂移、扩展、闪烁[2]等现象,从而影响通信质量,因此需要准确分析激光传输过程中大气湍流对通信系统的影响。一般采用大气相干长度r0对大气湍流介质进行定量描述,由于大气湍流效应受温度、风速等因素的影响,其值会不断发生变化,故快速测量大气相干长度对大气激光通信具有重要意义[3]。测量大气相干长度r0的方法包括光闪烁、激光雷达等,与这些方法相比,差分像运动法(DIMM)因成本低、操作简便得到了广泛应用。
O'Donovan等[4-5]设计研发了一种基于卡塞格林望远镜的DIMM测量仪,该仪器便携性较好,通过与风廓线雷达法测量结果进行对比,验证了该仪器稳定性较好。Okita等[6]使用DIMM测量仪,测量了南极科考站上方的大气湍流强度,同时在白天进行了观星实验。2012年,于龙昆等[7]利用DIMM测量仪得到了海边整层大气湍流参数的全天变化数据。数据分析表明,海边整层大气相干长度日变化明显,夜晚起伏无规则,且日变化没有陆地上的剧烈。2017年,董振龙等[8]分别在北京和成都两个城市地区搭建了光波传输特性测试链路,分析了不同地理区域内大气相干长度值的变化趋势及误差来源。2018年,王钰茹等[9]设计了折反路径上的激光探测成像系统,通过拟合的方法推算出湍流折射率结构常数与大气相干长度的值。
r0可表示为光波传输路径上的总湍流积分,因此可通过测量r0提高通信系统中的目标识别跟踪效率。本文在现有理论基础上,通过数值仿真模拟了不同大气相干长度下的湍流情况,定量分析了大气湍流强度对通信系统产生的影响,为实际中进行激光通信提供了理论支撑。最后搭建了光传输系统实验平台,对西安地区大气相干长度进行测量,分析了大气湍流在白天和夜间、晴天和多云条件下的变化趋势。
2 差分像运动法原理分析
反演法测量大气相干长度的方法,主要通过测量大气折射率结构常数进行反演,两者的关系可表示为[10]
式中:k为波数;
因反演法得到的大气相干长度值并不准确,所以本实验使用差分成像运动法测量r0,主要通过测量光波在大气中传输时引起的到达角起伏,即两个光斑相对位置的随机抖动计算大气相干长度。
结合CCD单像素的尺寸,可知到达角与光斑实际质心的位置关系为
式中:α为到达角;f为光学系统等效焦长。
差分像运动法是通过计算图像上两个光斑质心位置的相对运动求
式中:λ为光波波长;D为子瞳口径;d为子瞳中心间距,要求d≥2D;
3 不同湍流条件下远场光斑仿真分析
3.1 湍流相位屏模拟
使用功率谱反演法仿真湍流相位屏的原理为:首先生成一个频域内零均值、单位方差的Hermitian复高斯随机数矩阵,然后用湍流相位谱的功率谱密度函数Φ(kx,ky,kz)对其进行滤波处理,最后进行傅里叶逆变换得到大气湍流随机相位屏[12]。
假定大气湍流是各向同性且局地均匀的,则Kolmogorov谱的功率谱密度函数可表示为
与z轴方向垂直,通过任意相位板调制后的大气湍流相位功率谱可表示为
式中:Δz为湍流板的切面厚度;kr=(
引入大气湍流的相干长度r0,可将(5)式简化为
利用大气湍流相位谱与r0的对应关系,使用均匀大气条件下的功率谱结构模型,对不同湍流强度的湍流相位屏进行仿真,结果如
图 2. r0=0.01 m时模拟湍流相位屏的三维空间分布
Fig. 2. Three-dimensional space distribution of simulated turbulent phase screen with r0=0.01 m
图 3. r0=0.1 m时模拟湍流相位屏的三维空间分布
Fig. 3. Three-dimensional space distribution of simulated turbulent phase screen with r0=0.1 m
对比
图 4. r0=0.2 m时模拟湍流相位屏的三维空间分布
Fig. 4. Three-dimensional space distribution of simulated turbulent phase screen with r0=0.2 m
3.2 不同湍流强度下的激光传输仿真
1) 真空条件下的数值仿真
图 5. 真空中传输时的远场光斑三维分布图
Fig. 5. Three-dimensional map of far-field light spot during vacuum transmission
2) 不同强度湍流下的数值仿真
图 6. r0=0.01 m的光强分布。(a)三维分布图;(b) 二维分布图
Fig. 6. Light intensity distributions with r0=0.01 m. (a) Three-dimensional distribution; (b) two-dimensional distribution
图 7. r0=0.1 m的光强分布。(a)三维分布图;(b) 二维分布图
Fig. 7. Light intensity distributions with r0=0.1 m. (a) Three-dimensional distribution; (b) two-dimensional distribution
图 8. r0=0.2 m的光强分布。(a)三维分布图;(b) 二维分布图
Fig. 8. Light intensity distributions with r0=0.2 m. (a) Three-dimensional distribution; (b) two-dimensional distribution
对比真空中传输后的光斑分布可以发现,经不同强度大气湍流传输后,接收像面上的光斑发生了明显的畸变,随着湍流强度的增加,受到湍流影响也越严重,特别是
4 实验组成
4.1 实验系统搭建
采用差分像运动法对西安地区大气相干长度进行测量,采集全天接收到的信号光。选定发射地点与接收地点,在接收端设置一块白色幕布可明显观察光斑的实时位置。用CCD相机实时拍摄光斑图像,利用计算机软件对每一帧图像进行预处理,提取出光斑的质心位置进行分析。
4.2 测量地点选取
4.3 激光发射系统
4.4 接收系统
接收端望远镜镜筒前有一块4个直径为2 cm小孔的圆形挡板,用一块白色幕布作为接收屏与激光传输方向垂直放置,CCD相机拍摄尽量正对接收屏上的光斑,且在固定位置采集差分运动图像。
表 1. CCD相机性能指标
Table 1. Performance indexes of CCD camera
|
5 实验结果分析
由于西安冬季雾霾比较严重,因此选择在2018年8月(夏季)晴朗无云的天气进行测试,对采集到的光斑,使用文献[ 14]中的算法进行处理。
1) 测量时间:2018-08-17T19:00,晴天。使用DIMM测量仪进行实验观测,采样频率为20 Hz,采样数目为55个,对采集得到的连续多帧差分图像进行处理和计算,以
图 16. 连续多帧图像大气相干长度值
Fig. 16. Atmospheric coherence length of continuous multi-frame images
从
2) 测量时间:2018-08-15T08:00—23:30。使用DIMM测量仪平均每15 min测量一次,采样方法与上述相同,得到该时刻各个方向上r0的统计平均值如
图 17. 全天内4个方向大气相干长度测量值
Fig. 17. Measured atmospheric coherence length in four directions in a whole day
从
3) 测量时间:2018-08-15T08:00—23:30、2018-08-16T08:00—23:30。使用DIMM测量仪平均每15 min测量一次,采样方法与上述相同,得到该时刻内所有方向上r0的平均值,该时段内大气相干长度变化趋势如
图 18. 白天大气相干长度测量数值对比
Fig. 18. Numerical comparison of measured atmospheric coherence lengths at daytime
从
4) 测量时间:2018-08-16T00:00—08:00。使用DIMM测量仪平均每15 min测量一次,采样方法与上述相同,得到该时段内所有方向上r0的平均值,该时段内大气相干长度变化趋势如
对比
通过长时间的观测,得出光斑在晴天不同时段内的变化情况为:9:00前,光斑抖动范围较小,此时大气相干长度值总体变化范围也较小,且变化较为平稳;11:00左右,太阳逐渐升高,湍流开始变强,接收面上的光斑有小部分出现破碎现象;12:00—14:00左右,光斑抖动范围较大,光斑分布为很多破碎的小块,闪烁比较严重,12:30左右湍流强度达到极大值,此时大气相干长度最小,分别为18.5 cm和19.0 cm;15∶00以后光斑闪烁逐渐减小,光斑抖动情况开始变弱,大气相干长度值逐渐增大,变化比较稳定,变化范围在30.0~38.0 cm。
6 结论
因为大气湍流是随机运动的,现如今并没有一个衡量大气相干长度测量结果是否准确的标准,常用的判断方式有两种:一种是通过与已经开发并投入使用的设备进行对比,观察两种设备测量结果的一致性程度;另一种是分析测量结果与理论研究的湍流强度变化趋势是否一致。由于没有同类型的测量仪器进行对比,因此通过在不同时间段进行测量实验,证实了该系统测量得到的西安地区大气相干长度变化趋势基本与近地面层的湍流强度随时间变化的趋势相吻合,稳定性较好,与已有的大气相干长度测量系统相比较,该设备可进行不间断测量,在外场具有便于操作等优势,采取的多光斑检测算法可快速准确地提取所需数据,因而在实际工程中应用是可行的。
[1] 郭振华, 许德胜, 王世鹏, 等. 大气湍流对激光传输的影响[J]. 激光技术, 1992, 16(2): 65-72.
Guo Z H, Xu D S, Wang S P, et al. Effect of atomospheric turbulence on laser beam propagation[J]. Laser Technology, 1992, 16(2): 65-72.
[2] 姚海峰, 倪小龙, 陈纯毅, 等. 基于脉冲激光在大气中传输的信道补偿[J]. 光学学报, 2018, 38(1): 0101003.
[3] 宋卢军, 刘智, 倪小龙, 等. 大气相干长度测量实验与分析[J]. 长春理工大学学报(自然科学版), 2015, 38(2): 125-127, 132.
Song L J, Liu Z, Ni X L, et al. Measurement experiment and analysis of atmospheric coherence length[J]. Journal of Changchun University of Science and Technology (Natural Science Edition), 2015, 38(2): 125-127, 132.
[4] O'DonovanB, Seneta E B, John S, et al. DIMMWIT measurements of the spatial and temporal scale of atmospheric turbulence at COAST and MROI[C]∥SPIE Astronomical Telescopes + Instrumentation. Proc SPIE 5491, New Frontiers in Stellar Interferometry, Glasgow, United Kingdom, 2004: 1310- 1319.
[5] O'DonovanB. Differential image motion monitor, which is transportable[D]. Cambridge, 2006: 20- 35.
[6] Okita H, Ichikawa T. Ashley M C B, et al. Excellent daytime seeing at Dome Fuji on the Antarctic plateau[J]. Astronomy & Astrophysics, 2013, 554: L5.
[7] 于龙昆, 吴毅, 侯再红, 等. 南海海边整层大气相干长度的测量与分析[J]. 红外与激光工程, 2012, 41(11): 3046-3049.
Yu L K, Wu Y, Hou Z H, et al. Measurement and analysis of whole layer atmospheric coherence length at seaside of South China Sea[J]. Infrared and Laser Engineering, 2012, 41(11): 3046-3049.
[8] 董振龙, 高超, 李黎, 等. 城市信道中大气相干长度测试与分析[J]. 光通信技术, 2017, 41(5): 52-54.
Dong Z L, Gao C, Li L, et al. Test and analysis of atmospheric coherence length in optical wireless communication channel for urban area[J]. Optical Communication Technology, 2017, 41(5): 52-54.
[9] 王钰茹, 梅海平, 康丽, 等. 湍流大气中折反路径激光成像探测实验[J]. 中国激光, 2018, 45(4): 0401008.
[10] Fried D L. Optical resolution through a randomly inhomogeneous medium for very long and very short exposures[J]. Journal of the Optical Society of America, 1966, 56(10): 1372-1379.
[11] AoJ, Liu JQ, Ma CB. A new measurement of atmospheric coherent length using differential image motion[C]∥2016 11th International Symposium on Antennas, Propagation and EM Theory (ISAPE), October 18-21, 2016. Guilin, China. IEEE, 2016: 264- 267.
[12] 韩冬, 刘云清, 赵馨, 等. 光束腰半径对空间相干光通信外差效率的影响[J]. 长春理工大学学报(自然科学版), 2016, 39(3): 36-40.
Han D, Liu Y Q, Zhao X, et al. Effect of beam radius on heterodyne efficiency for space coherent optical communication[J]. Journal of Changchun University of Science and Technology(Natural Science Edition), 2016, 39(3): 36-40.
[13] 胡明勇, 封志伟, 范二荣, 等. 一种可变远距离大功率激光聚焦系统的设计[J]. 激光与光电子学进展, 2018, 55(6): 062204.
[14] 吴鹏飞, 贾璐瑶. 一种多目标光斑自适应快速检测算法[J]. 激光杂志, 2019( 6): 27- 31.
Wu PF, Jia LY. A fast adaptive detection algorithm for multi-objective spot[J]. Laser Journal, 2019( 6): 27- 31.
Article Outline
吴鹏飞, 贾璐瑶. 西安地区大气相干长度的实验测量及分析[J]. 激光与光电子学进展, 2020, 57(9): 090101. Pengfei Wu, Luyao Jia. Experimental Measurement and Analysis of Atmospheric Coherence Length in Xi'an Area[J]. Laser & Optoelectronics Progress, 2020, 57(9): 090101.