中国激光, 2017, 44 (11): 1105001, 网络出版: 2017-11-17   

海洋湍流中扩展物体漫反射光成像的数值模拟 下载: 926次

Numerical Simulation of Imaging by Diffuse Reflection Light from Extended Objects in Oceanic Turbulence
作者单位
四川师范大学物理与电子工程学院, 四川 成都 610068
摘要
采用空域滤波方法和海洋湍流功率谱, 编制了海洋湍流中扩展物体漫反射光成像的模拟仿真程序, 研究了成像质量与海洋湍流参数和接收孔径之间的关系。研究发现, 随着海水湍流功率谱中由温度引起的海水折射率变化与由盐度变化引起的海水折射率变化的比率的增大, 海水温度方差耗散率的增大, 以及单位质量海水动能耗散率的减小, 海洋湍流的可见参数减小, 成像质量变差。该数值模拟方法为研究海洋湍流中扩展物体的漫反射光成像提供了一种可行、正确的方法。
Abstract
By using the spatial filtering method and the oceanic turbulence power spectra, the numerical simulation program of imaging by diffuse reflection light from extended objects in oceanic turbulence is designed. The dependence of the imaging quality on the oceanic turbulence parameters and the receiver size is studied. It is found that, with the increases of the ratio of seawater refractive index change caused by temperature to that caused by salinity in the oceanic turbulence power spectra and the dissipation rate of seawater temperature variance, and with the decrease of the dissipation rate of kinetic energy per unit fluid mass, the seeing parameter in oceanic turbulence decreases and the imaging quality degrades. The proposed numerical simulation method provides a feasible and correct approach for the study of imaging by diffuse reflection light from extended objects in oceanic turbulence.

1 引言

激光照明目标漫反射回光特性和成像技术已被广泛应用于诸多激光领域[1-8]。湍流会降低光学成像质量,引起光闪烁和漂移等[9-17]。大气湍流中漫反射光的特性和成像已有研究[18-20]。近年来,随着水下光通信、传感和激光雷达等应用的兴起,深入研究海洋湍流对光成像的影响具有重要意义[21-23]。大气湍流折射率起伏的空间功率谱仅含有一个峰值,而海洋湍流折射率起伏的空间功率谱具有复杂的双峰值结构[24]。本课题组推导出了光波在海洋湍流中的结构函数和空间相干长度的解析公式,并得到Kolmogorov结构函数的5/3幂律在海洋湍流的惯性区域中仍然成立的结论[25];推导出了海洋湍流中可见参数的解析表达式,并修正了已有的用光学传递函数定量描述水下光学成像的模型[26]。然而,关于海洋湍流中漫反射光的传输及成像的数值模拟研究鲜有报道。本文利用空域滤波和海洋湍流功率谱,对扩展物体漫反射光在海洋湍流中的成像进行了数值模拟,研究了光学成像质量与海洋湍流参数和接收孔径之间的关系。

2 海洋湍流中漫反射光传输理论

在傍轴近似条件下,光在随机介质中沿z轴方向的传输满足标量波动方程[27],即

2ikUz+2U+2k2(n2-1)U=0,(1)

式中U为初始光波复振幅; 2=∂2/∂2x+∂2/∂2y为拉普拉斯算子,xy是与传输方向垂直的横向坐标;n为海水折射率;光波波数k=2π,其中λ为波长。通常,采用多层相位屏法求解(1)式,即将光束传输路径分成若干段,每一段距离中湍流引起的相位扰动可近似用薄的随机相位屏来处理,光在湍流介质中的传输等价为光在真空和薄随机相位屏中的交替传输[27]

对于均匀各向同性的海水介质,采用Nikishov等[24]提出的海洋湍流功率谱Φn(κ)来描述由温度和盐度共同变化引起的海水折射率起伏,即

Φn(κ)=0.388×10-8ε-1/3κ-11/3[1+2.35(κη)23]×χT[ω2exp(-ATγ)+exp(-ASγ)-2ωexp(-ATSγ)]/ω2,(2)

式中常数AT=1.863×10-2,AS=1.9×10-4,ATS=9.41×10-3,γ=8.284(κη)4/3+12.978(κη)2;κ为湍流功率谱空间波数;η为Kolmogorov内尺度;χT为海水的温度方差耗散率,其取值范围是为10-4~10-10 K2·s-1;ε为单位质量的海水湍流动能耗散率,其取值范围为10-1~10-10 m2·s-3;w为海水湍流功率谱中由温度引起的海水折射率变化与由盐度变化引起的海水折射率变化的比率,其取值范围为[-5, 0],-5和0分别表示仅由温度或盐度变化引起海水折射率变化的两种情况。

由Fried[28]首次提出的可见参数r0是大气光学中一个重要物理量,它可表征光传输路径上的大气湍流强度。本课题组推导出的海洋湍流中球面波可见参数r0的表达式[26]

r0=[5.237×10-8k2zε-1/3(χT/ω2)×(0.419ω2-0.838ω+1.119)]-3/5,(3)

式中z为光传输距离。r0值越大,可视条件越好。

由于漫反射光场的发散角非常大,漫反射光传输过程不满足近轴近似光传输波动方程。利用激光照明物体漫反射光成像时,仅关心接收孔径内的光波信息。如果在接受孔径内可得到精确的光波数值,就能够研究其成像特性[18]。因此,在数值模拟计算时,可采用辅助模板进行空域滤波,辅助模板函数[29]

M(x,y)=1,(x,y)S0,otherwise,(4)

式中S为辅助模板的通光区域,通常辅助模板的通光直径需略大于接收孔径。数值计算时,将传输路径分成若干段,传输每段路径后,将光场复振幅与辅助模板函数相乘,这样可使得S区域内的光场保持不变,而区域外的光场强度设置为0。

湍流介质中的长曝光调制传递函数[28]可表示为

fMT=exp-12D(λΩ),(5)

式中Ω为空间角频率,波结构函数D(·)为在Rytov微扰近似中获得的对数振幅和相位起伏的结构函数之和[9],可表示为

D(ρ)=8π2k2z0101-J0(κξρ)]Φn(κ)κdκdξ,(6)

式中J0(·)为零阶贝塞尔函数,ρ为垂直于光波传输方向截面上两点的距离,归一化距离变量ξ=1-L/z(z为物面到像面的距离;L为光传输的实际距离,其变化范围为[0,z])。将(2)式代入(6)式,ρη时,可得到海洋湍流的波结构函数[25]

D(ρ,z)=3.603×10-7k2zε-1/3(χT/ω2)ρ53×(0.419ω2-0.838ω+0.419)(7)

研究调制传递函数时,将波结构函数中的ρλΩ代替。

真空中,成像系统的调制传递函数[28]可以表示为

fMT0=2π[arccos(Ω/Ω0)-(Ω/Ω0)1-(Ω/Ω0)2],Ω<Ω00,ΩΩ0,(8)

式中Ω0=D0为截止频率,D0为接收孔径。

对一个光学系统分辨本领的描述有多种方法,例如可通过到达角起伏和瑞利判据等简单分析分辨本领[30],可使用调制传递函数的全频谱空间积分作为光学系统分辨本领的度量参数[30]。在海洋湍流中,长曝光成像分辨率R[28]表示为

R=2π0fMT0fMTΩdΩ(9)

将(5)~(8)式代入(9)式,即可得到长曝光条件下海洋湍流的成像分辨率。

在真空中,系统的最小可分辨尺寸R主要由接收孔径D0决定,即R=1.22λz/D0。在海洋湍流中,当接收孔径D0远大于可见参数r0时,系统的最小可分辨尺寸R主要由海洋湍流决定,即R=1.22λz/r0

3 漫反射光成像的数值模拟结果及分析

研究图1所示的漫反射体对相干照明光的反射光在海洋湍流中的传输及成像特性,该漫反射体是由一系列尺寸为0.23 mm×0.23 mm的小方格组成,圆环的半径由中心向外依次分别为0.82,1.64,3.05,4.69 mm。采用在光场上施加高斯随机振幅和相位的方法模拟反射物体对相干照明光的漫反射。在数值模拟中,假设物体被均匀照明,且物体反射率均匀分布,即从物体散射出的漫反射光场的振幅均匀分布,相位在[0,2π]范围内均匀随机变化。取光波波长λ=0.6328 μm,网格数为512×512,像面距离z=50 m,模拟50次随机相位散射光的独立传输,海洋湍流中像面光强为100次数值模拟的平均结果。

图 1. 漫反射物体

Fig. 1. Diffuse reflection object

下载图片 查看所有图片

3.1 接收孔径对成像质量的影响

r0=17.7 mm时,海洋湍流中不同接收孔径D0对应的成像最小可分辨尺寸R以及可分辨物体细节情况见表1,其中D1D2D3D4分别表示中心与第一环、第二环、第三环、第四环之间的距离,○表示可以分辨细节,×表示不能分辨细节。与表1对应的不同D0情况下的成像及其二维(2D)光强分布如图2所示。可以看出,图2的结果与表1的理论分析结果是一致的。D0越大,能够分辨的细节越多,成像质量越高。但是,当D0增大到一定值后,可分辨细节不再增多,这时成像质量由海洋湍流强度决定。

表 1. 海洋湍流中,不同D0下的R以及可分辨物体细节情况

Table 1. R and distinguishable details under different D0 in oceanic turbulence

D0 /mmR /mmDistinguishable detail
D1=0.82 mmD2=1.64 mmD3=3.05 mmD4=4.69 mm
66.43××××
103.86×××
202.18××
302.18××

查看所有表

图 2. 海洋湍流中,不同D0下的(a1)(b1)(c1)(d1)像及(a2)(b2)(c2)(d2)其二维光强分布

Fig. 2. (a1)(b1)(c1)(d1) Images and (a2)(b2)(c2)(d2) their 2D intensity distributions under different D0 in oceanic turbulence

下载图片 查看所有图片

图 3. 海洋湍流中,不同r0下的R随D0的变化

Fig. 3. R versus D0 under different r0 in oceanic turbulence

下载图片 查看所有图片

在海洋湍流中,不同可见参数r0下分辨率R随接收孔径D0的变化曲线如图3所示,可以看出,R随着D0的增大而增大。当D0足够大时,R趋于稳定值。因此,表1图2的结果与图3的结果是一致的,这表明该数值模拟方法是可行和正确的。此外,由图3还可知,不同r0下的R的差值随D0的增大而增大。

图 4. 在海洋湍流中,不同r0下的(a1)(b1)(c1)(d1)(e1)像及(a2)(b2)(c2)(d2)(e2)其二维光强分布

Fig. 4. (a1)(b1)(c1)(d1)(e1) Images and (a2)(b2)(c2)(d2)(e2) their 2D intensity distributions under different r0 in oceanic turbulence

下载图片 查看所有图片

3.2 海洋湍流对成像质量的影响

D0=40 mm和不同可见参数r0条件下,海洋湍流中的成像最小可分辨尺寸R以及可分辨物体细节情况见表2,与表2对应的像及其2D光强分布如图4所示。可以看出,图4表2的理论分析结果是一致的,r0越大,可分辨细节越多,成像质量越高。

表 2. 海洋湍流中,不同r0下的R以及可分辨物体细节情况

Table 2. R and distinguishable details under different r0 in oceanic turbulence

r0 /mmR /mmDistinguishable detail
D1=0.82 mmD2=1.64 mmD3=3.05 mmD4=4.69 mm
4.269.06××××
10.433.70×××
16.972.27××
29.061.33×
72.990.53

查看所有表

在海洋湍流中,不同接收孔径D0下的分辨率R随不同海洋湍流参数的变化如图5所示。图5(a)表明,随着ω的增大,R逐渐减小,成像质量变差。如图5(b)、(c)所示,随着海水温度方差耗散率χT的增大和单位质量海水湍流动能耗散率ε的减小,R减小,成像质量下降。因此,表2图4的结果与图5的结果是一致的,进一步表明该数值模拟方法是可行的,结果是正确的。此外,由图5可知,随着ωχT的减小以及ε的增大,不同D0下的R的差别增大,即接收孔径对成像质量的影响增大。反之,当ωχT足够大且ε足够小时,不同D0对应的R值接近,即此时接收孔径对成像质量的影响可以忽略,而海洋湍流的强度是决定成像质量的主要因素。海洋湍流的可见参数r0随不同海洋湍流参数的变化如图6所示。

图 5. 海洋湍流中,不同D0下的R随不同海洋湍流参数的变化。(a) ω;(b) lg χT;(c) lgε

Fig. 5. R in oceanic turbulence versus different oceanic turbulence parameters under different D0. (a) ω; (b) lg χT; (c) lg ε

下载图片 查看所有图片

图 6. 海洋湍流中,r0随不同海洋湍流参数的变化。(a) ω;(b) lg χT;(c) lg ε

Fig. 6. r0 in oceanic turbulence versus different oceanic turbulence parameters. (a) ω; (b) lg χT; (c) lg ε

下载图片 查看所有图片

4 结论

采用空域滤波的方法,基于海洋湍流功率谱,编制了海洋湍流中扩展物体漫反射光成像的模拟仿真程序,研究了海洋湍流中扩展物体漫反射光成像的规律,并研究了成像质量与海洋湍流参数和接收孔径之间的关系。研究发现,随着海水湍流功率谱中温度与盐度起伏的比率的增大,海水温度方差耗散率的增大,以及单位质量海水湍流动能耗散率的减小,海洋湍流的可见参数r0减小,成像质量变差。研究结果表明,这种数值模拟方法是可行的,结果是正确的,为研究海洋湍流中扩展物体漫反射光成像提供了一种方法。研究结论对工作于水下湍流环境中的扩展物体漫反射光成像的相关应用具有重要意义。

参考文献

[1] Borah D K, Voelz D G. Cramer-Rao lower bounds on estimation of laser system pointing parameters by use of the return photon signal[J]. Opt Lett, 2006, 31(8): 1029-1031.

[2] Borah D K, Voelz D G. Estimation of laser beam pointing parameters in the presence of atmospheric turbulence[J]. Appl Opt, 2007, 46(23): 6010-6018.

[3] Lukesh G, Chandler S, Voelz D G. Estimation of laser system pointing performance by use of statistics of return photons[J]. Appl Opt, 2000, 39(9): 1359-1371.

[4] Andrejczuk M, Moszkowic S, Haman K E . Radar-echo tracking by use of invariant moments[J]. Appl Opt, 2003, 42(29): 5891-5896.

[5] 董志伟, 张瑞环, 张伟斌, 等. 条纹管激光雷达回波噪声特性研究[J]. 光学学报, 2016, 36(9): 0928001.

    Dong Zhiwei, Zhang Ruihuan, Zhang Weibin, et al. Noise features in streak tube lidar echo signal[J]. Acta Optica Sinica, 2016, 36(9): 0928001.

[6] 张金生, 王仕成, 徐萍, 等. 激光制导武器仿真系统漫反射屏研制及标定[J]. 光电子·激光, 2005, 16(4): 421-424.

    Zhang Jinsheng, Wang Shicheng, Xu Ping, et al. Development and demarcating of diffuse reflector of laser guided weapon simulation system[J]. Journal of Optoelectronics·Laser, 2005, 16(4): 421-424.

[7] 周鑫, 孙剑锋, 姜鹏, 等. 基于532 nm激光目标回波脉冲宽度的成像实验研究[J]. 光学学报, 2015, 35(s1): s111002.

    Zhou Xin, Sun Jiangfeng, Jiang Peng, et al. Imaging experimental study based on the width of 532 nm-laser target echo pulse[J]. Acta Optica Sinica, 2015, 35(s1): s111002.

[8] 赵鹏, 张艳, 张小亚, 等. 激光漫反射测距回波信号在线提取方法[J]. 光学学报, 2015, 35(10): 1028002.

    Zhao Peng, Zhang Yan, Zhang Xiaoya, et al. Echo online extraction for diffuse reflection laser ranging[J]. Acta Optica Sinica, 2013, 35(10): 1028002.

[9] Andrews LC, Phillips RL. Laser beam propagation through random media[M]. Bellingham: SPIE Press, 2005.

[10] Gu Y, Gbur G. Scintillation of nonuniformly correlated beams in atmospheric turbulence[J]. Opt Lett, 2013, 38(9): 1395-1397.

[11] Ji X L, Deng J P. Effect of spherical aberration on scintillations of Gaussian beams in atmospheric turbulence[J]. Phys Lett A, 2014, 378(36): 2729-2735.

[12] Liu X H, Pu J X. Investigation on the scintillation reduction of elliptical vortex beams propagating in atmospheric turbulence[J]. Opt Express, 2011, 19(27): 26444-26450.

[13] Liu X L, Shen Y, Liu L, et al. Experimental demonstration of vortex phase-induced reduction in scintillation of a partially coherent beam[J]. Opt Lett, 2013, 38(24): 5323-5326.

[14] 贾锐, 韦宏艳, 张洪建, 等. 斜程大气湍流中点目标回波的闪烁研究[J]. 中国激光, 2015, 42(11): 1113001.

    Jia Rui, Wei Hongyan, Zhang Hongjian, et al. Scintillation index of echo wave in slant atmospheric turbulence[J]. Chinese J Lasers, 2015, 42(11): 1113001.

[15] 陈振, 楚兴春, 赵尚弘, 等. 艾里涡旋光束在大气湍流中的漂移特性研究[J]. 中国激光, 2015, 42(12): 1213002.

    Chen Zhen, Chu Xingchun, Zhao Shanghong, et al. Study of the drift characteristics of Airy vortex beams in atmospheric turbulence[J]. Chinese J Lasers, 2015, 42(12): 1213002.

[16] 王加安, 王馨兰, 郭林炀, 等. 艾里光束的光强闪烁研究[J]. 光学学报, 2017, 37(6): 0626001.

    Wang Jiaan, Wang Xinlan, Guo Linyang, et al. Light intensity scintillation of Airy beam[J]. Acta Optica Sinica, 2017, 37(6): 0626001.

[17] 柯熙政, 张宇. 部分相干光在大气湍流中的光强闪烁效应[J]. 光学学报, 2015, 35(1): 0106001.

    Ke Xizheng, Zhang Yu. Scintillation of partially coherent beam in atmospheric turbulence[J]. Acta Optica Sinica, 2015, 35(1): 0106001.

[18] 王英俭, 吴毅. 扩展物体漫反射光传输及成像的数值模拟研究[J]. 光学学报, 1998, 18(10): 1470-1472.

    Wang Yingjian, Wu Yi. Numerical simulation of propagation of diffuse reflection light of extended object[J]. Acta Optica Sinica, 1998, 18(10): 1470-1472.

[19] 张宇, 李新阳, 饶长辉. 基于空间频率域滤波的漫反射光大气湍流传输数值模拟仿真方法[J]. 光学与光电技术, 2011, 9(5): 52-57.

    Zhang Yu, Li Xinyang, Rao Changhui. Numerical simulation method of propagation of diffuse reflection optics based on frequency filtering[J]. Optics and Optoelectronic Technology, 2011, 9(5): 52-57.

[20] Wu Z S, Li Y Q. Scattering of a partially coherent Gaussian-Schell beam from a diffuse target in slant atmospheric turbulence[J]. JOSA A, 2011, 28(7): 1531-1539.

[21] Hou W, Woods S. Optical turbulence on underwater image degradation in natural environments[J]. Appl Phys, 2012, 51(14): 2678-2686.

[22] Hou W, Jarosz E. Impacts of underwater turbulence on acoustical and optical signals and their linkage[J]. Opt Express, 2013, 21(4): 4367-4375.

[23] Hou W. A simple underwater imaging model[J]. Opt Lett, 2009, 34(17): 2688-2690.

[24] Nikishov V V, Nikishov V I. Spectrum of turbulent fluctuations of the sea-water refraction index[J]. International Journal of Fluid Research, 2000, 27(1): 82-98.

[25] Lu L, Ji X L, Baykal Y. Wave structure function and spatial coherence radius of plane and spherical waves propagating through oceanic turbulence[J]. Opt Express, 2014, 22(22): 27112-27122.

[26] Pu H, Ji X L. Oceanic turbulence effects on long-exposure and short-exposure imaging[J]. J Opt, 2016, 18(10): 105704.

[27] Fleck J A, Morris J R, Feit M D. Time-dependent propagation of high energy laser beams through the atmosphere[J]. Appl Phys A, 1976, 10(2): 129-160.

[28] Fried D L. Optical resolution through a randomly inhomogeneous medium for very long and very short exposures[J]. JOSA, 1966, 56(10): 1372-1379.

[29] Vadim V D, Mikhail A V, Valeriy V K. Speckle-field propagation in "frozen" turbulence: Brightness function approach[J]. JOSA A, 2006, 23(8): 1924-1936.

[30] 饶瑞中. 现代大气光学[M]. 北京: 科学出版社, 2012.

    RaoRuizhong. Modern atmospheric optics[M]. Beijing: Science Press, 2012.

罗燏娟, 季小玲, 李晓庆. 海洋湍流中扩展物体漫反射光成像的数值模拟[J]. 中国激光, 2017, 44(11): 1105001. Luo Yujuan, Ji Xiaoling, Li Xiaoqing. Numerical Simulation of Imaging by Diffuse Reflection Light from Extended Objects in Oceanic Turbulence[J]. Chinese Journal of Lasers, 2017, 44(11): 1105001.

本文已被 2 篇论文引用
被引统计数据来源于中国光学期刊网
引用该论文: TXT   |   EndNote

相关论文

加载中...

关于本站 Cookie 的使用提示

中国光学期刊网使用基于 cookie 的技术来更好地为您提供各项服务,点击此处了解我们的隐私策略。 如您需继续使用本网站,请您授权我们使用本地 cookie 来保存部分信息。
全站搜索
您最值得信赖的光电行业旗舰网络服务平台!