一种基于多层复数相位屏表征的远场激光大气传输仿真方法技术

本发明专利技术公开了一种基于多层复数相位屏表征的远场激光大气传输仿真方法，利用复数相位屏综合表征不同传输距离处的大气衰减、大气湍流和大气热晕三种大气光学效应；基于求得的复数相位屏，通过构建离散传输路径下的多层相位屏数值仿真方法，得到远场激光的功率密度分布，从而实现远场激光大气传输仿真过程。本发明专利技术通过引入复数相位屏，其实质是将大气光学中复杂的光学效应等效表征为复折射率的扰动效应，可以综合表征激光传输过程中大气衰减、大气湍流和大气热晕等光学效应，为激光大气传输仿真方法提供了理论和算法基础。

【技术实现步骤摘要】
一种基于多层复数相位屏表征的远场激光大气传输仿真方法
本专利技术属于大气光学
，特别是一种基于多层复数相位屏表征的远场激光大气传输仿真方法。
技术介绍
大气对激光测距、激光通信、激光雷达、激光武器等应用会造成诸多影响。这些影响包括：大气气体分子、大气中悬浮的气溶胶粒子以及各种降水粒子(如：霾、云、雾、雨、雪、沙尘、烟尘等)对光波的散射和吸收引起的衰减效应；大气温度、气压和湿度等随机起伏造成的大气折射率随机波动，导致光束强度的起伏、光束的扩展和漂移等大气湍流效应；高能激光大气传输时产生的非线性效应(如热晕和气体击穿效应)。大气对激光传输的影响不仅会造成激光能量的衰减，还会引起激光光强的起伏，时间、空间相干性下降，从而导致激光器的效能降低，甚至无法正常工作。目前，人们已经对陆地大气环境中的激光传输特性进行了大量的研究，提出了不少理论模型与数值计算模型，试验经验模型等，但对于海洋大气环境激光传输模型的研究较少，模型中考虑的影响因素通常只有单个或少数，无法全面分析各种大气光学效应对激光传输特性的影响，如气溶胶的指数衰减模型、大气湍流的光强起伏模型、光斑漂移模型等。为了弄清激光在大气中的传输特性，学者们提出了许多的激光传输模型，根据适用条件大致分为三个方面：大气衰减、大气湍流、大气热晕。大气衰减研究混合介质中的光传播问题，即对辐射传输方程求解，将光强作为求解对象；大气湍流研究湍流介质中的光传播问题，即对波传播方程求解，将光场作为求解对象；大气热晕研究高能激光在大气介质中的光传播问题，即对流体力学的质量守恒、动能守恒、动量守恒方程和波传播方程求解，将大气密度作为求解对象。其中大气衰减的激光传输模型存在一定的缺陷，只能得到激光传输特性中的光强信息，光场空间分布以及相位信息无法得到，不适用于激光通信等方面的应用；大气湍流的激光传输模型无法体现衰减的信息，故而无法体现大气环境中温度、湿度、能见度等大气参数对激光传输特性的影响；大气热晕的激光传输模型研究主要以数值仿真为主，由于瞬时热晕的复杂性和难可控性，稳态热晕效应是目前的主要研究手段。以上三种大气光学效应的产生机理不同，因此模型构建时通常是将每种光学效应单独处理后再综合评估传输效能，但缺点是不能对激光大气传输过程中的功率密度分布进行实时跟踪和图像可视化，三种大气光学效应的综合影响不能同时评估。
技术实现思路
本专利技术的目的在于针对现有各种大气光学效应不能综合表征的缺点，提供一种基于多层复数相位屏表征的远场激光大气传输仿真方法，实现远场激光传输特性的评估。实现本专利技术目的的技术方案为：一种基于多层复数相位屏表征的远场激光大气传输仿真方法，包括以下步骤：利用复数相位屏综合表征不同传输距离处的大气衰减、大气湍流和大气热晕三种大气光学效应；其中大气衰减由复数相位屏的虚部表征，其虚部值由Beer定律等效计算得到；大气湍流和大气热晕由复数相位屏的实部表征，其实部值由流体力学方程计算得到湍流和热晕引起的大气密度起伏并根据Gladstone-Dale关系式间接计算得到；基于求得的复数相位屏，通过构建离散传输路径下的多层相位屏数值仿真方法，得到远场激光的功率密度分布，从而实现远场激光大气传输仿真过程。与现有技术相比，本专利技术的显著优点为：(1)本专利技术通过引入复数相位屏，其实质是将大气光学中复杂的光学效应等效表征为复折射率的扰动效应，可以综合表征激光传输过程中大气衰减、大气湍流和大气热晕等光学效应，为激光大气传输仿真方法提供了理论和算法基础；(2)本专利技术基于复数相位屏拓展了传统的多层相位屏法，提出了一种新的激光大气传输数值仿真方法，该方法可以解决大气衰减、大气湍流和大气热晕等光学效应不能综合表征的困难；(3)本专利技术可以方便地对激光大气传输的远场光斑功率密度分布进行仿真计算并可视化显示，预计为海洋大气环境中远场激光特性评估及其军事应用提供重要的技术支持。附图说明图1是单层复数相位屏的构建示意图。图2是多层相位屏法的算法流程图。图3是基于多层复数相位屏的激光大气传输仿真过程的信息流程图。图4是仿真计算的到靶激光光斑的功率密度分布灰度图。具体实施方式多层相位屏法作为激光大气传输的基本数值仿真方法，主要适用于描述大气湍流扰动下的激光传输过程，其模型中主要以大气折射率功率谱结构函数体现大气湍流对光学扰动的影响，其仿真过程具有全路径跟踪、易参数控制的优点。基于多层相位屏法，研究将大气衰减、大气湍流、大气热晕效应综合表征为多层复数相位屏上的相位扰动，完善激光大气传输仿真理论及实现远场激光传输特性的评估是本专利技术的主要目的。本专利技术的一种基于复数相位屏表征的远场激光大气传输仿真方法，包括以下步骤：第一步，利用复数相位屏综合表征不同传输距离处的大气衰减、大气湍流和大气热晕三种大气光学效应；其中大气衰减由复数相位屏的虚部表征，其虚部值由Beer定律等效计算得到；大气湍流和大气热晕由复数相位屏的实部表征，其实部值由流体力学方程计算得到湍流和热晕引起的大气密度起伏并根据Gladstone-Dale关系式间接计算得到；如图1所示。复数相位屏表示为：式(1)中，r表示垂直于传输方向横截面上两点间的距离，λ为激光波长，Δz为复数相位屏的厚度，n1(r,zi)表征传输距离zi处的大气湍流和大气热晕综合引起的折射率屏起伏，κ(λ,zi)等效表征传输距离zi处的光强衰减系数。其中，复数相位屏S(r,zi)的实部等效于由大气密度扰动间接导致折射率起伏对应的附加相位扰动，因此大气湍流与大气热晕的综合效果(即大气密度的扰动)需要利用流体力学方程进行求解。大气湍流的影响因素较为复杂，其扰动主要与全空间域的温度脉动有关(可近似为大气分子内能的涨落)，因此求解流体力学方程并不现实，故大气湍流引起的光学效应常常表征为空间分布的大气折射率结构分布函数；大气热晕是由激光束能量与大气相互作用导致的局部大气密度起伏，故可直接由流体力学方程求解得到大气密度起伏值。综上分析，大气湍流与大气热晕可综合表征为温度脉动引起的大气分子内能涨落和激光束局部能量的叠加而导致的大气密度起伏，从而间接引起大气折射率起伏变化，可表征为复数相位屏的实部值；大气衰减等效为光强的衰减，可等效表征为复数相位屏的虚部值。(1)κ(λ,zi)的求解激光在实际大气中传输时，大气分子与气溶胶对激光的散射与吸收综合表现为对激光光强的衰减。相邻复数相位屏间的激光光强由Beer定律给出：dI(r,zi)＝-I(r,zi-1)μt(λ,zi)dz(2)其中，zi-1、zi分别为第i-1、i个复数相位屏与激光光源之间的传输距离；μt(λ,zi)为第i个复数相位屏的大气衰减系数；I(r,zi-1)、I(r,zi)分别为激光入射到距离zi-1、zi处的垂直于传输方向上的光强分布。式(2)积分得：I(r,zi)＝I(r,zi-1)exp[-Δz·μt(λ,zi)](3)其中，Δz为复数相位屏的厚度。式(3)中的大气衰减系数可由半经验公式计算得到：其中，V(λ)为大气能见度；联立式(1)、(3)、(4)可得到大气衰减系数为：(2)n1(r,zi)的求解大气湍流和热晕均能引起大气密度的起伏，而大气密度起伏的描述是从大气流体特性的质量守恒、动量守恒和能量守恒的流体力学方程组出发的。其中，ρ为大气介质的本文档来自技高网...

【技术保护点】
1.一种基于多层复数相位屏表征的远场激光大气传输仿真方法，其特征在于，包括以下步骤：利用复数相位屏综合表征不同传输距离处的大气衰减、大气湍流和大气热晕三种大气光学效应；其中大气衰减由复数相位屏的虚部表征，其虚部值由Beer定律等效计算得到；大气湍流和大气热晕由复数相位屏的实部表征，其实部值由流体力学方程计算得到湍流和热晕引起的大气密度起伏并根据Gladstone‑Dale关系式间接计算得到；基于求得的复数相位屏，通过构建离散传输路径下的多层相位屏数值仿真方法，得到远场激光的功率密度分布，从而实现远场激光大气传输仿真过程。

【技术特征摘要】
1.一种基于多层复数相位屏表征的远场激光大气传输仿真方法，其特征在于，包括以下步骤：利用复数相位屏综合表征不同传输距离处的大气衰减、大气湍流和大气热晕三种大气光学效应；其中大气衰减由复数相位屏的虚部表征，其虚部值由Beer定律等效计算得到；大气湍流和大气热晕由复数相位屏的实部表征，其实部值由流体力学方程计算得到湍流和热晕引起的大气密度起伏并根据Gladstone-Dale关系式间接计算得到；基于求得的复数相位屏，通过构建离散传输路径下的多层相位屏数值仿真方法，得到远场激光的功率密度分布，从而实现远场激光大气传输仿真过程。2.根据权利要求1所述的基于多层复数相位屏表征的远场激光大气传输仿真方法，其特征在于，复数相位屏表示为：式(1)中，r表示垂直于传输方向横截面上两点间的距离，λ为激光波长，Δz为复数相位屏的厚度，n1(r,zi)表征传输距离zi处的大气湍流和大气热晕综合引起的折射率屏起伏，κ(λ，,zi)等效表征传输距离zi处的光强衰减系数；(1)κ(λ，,zi)的求解相邻复数相位屏间的激光光强由Beer定律给出：dI(r,zi)＝-I(r,zi-1)μt(λ,zi)dz(2)其中，zi-1、zi分别为第i-1、i个复数相位屏与激光光源之间的传输距离；μt(λ,zi)为第i个复数相位屏的大气衰减系数；I(r,zi-1)、I(r,zi)分别为激光入射到zi-1、zi对应复数相位屏上的光强分布；式(2)积分得：I(r，zi)＝I(r，zi-1)exp[-Δz·μt(λ，，zi)](3)式(3)中的大气衰减系数可由半经验公式计算得到：其中，V(λ)为大气能见度；联立式(1)、(3)、(4)可得到大气衰减系数为：(2)n1(r，zi)的求解大气密度起伏与大气流体特性的质量守恒、动量守恒和能量守恒的流体力学方程组相关；其中，ρ为大气介质的密度；P为大气压强；v为大气介质速度；γ＝CP/CV为流动的定压比热和定容比热之比；μabs表示大气吸收系数；T为大气温度；I为光强；f(T)为与大气温度脉动有关的扰动能量；将大气介质的密度，压力和速度分别表示为平均量和一阶起伏量的和，ρ＝ρ0+ρ1，P＝P0+P1，v＝v0+v1，其中ρ0、ρ1分别为大气介质密度的平均量和一阶起伏量，P0、P1分别为大气介质的压力平均量和一阶起伏量，v0、v1分别为大气介质的速度平均量和一阶起伏量；线性化的方程组并以真微分的形式写为：式(11)中，cS为定义的风速在等压近似下，略去P1则可得：式(12)的解为：其中，F(T)为f(T)关于时间t的积分，表示温度T下的大气分子内能扰动量；依据流体介质折射率与密度的...

【专利技术属性】
技术研发人员：杨云涛，冷坤，武文远，黄雁华，龚艳春，章曦，
申请(专利权)人：中国人民解放军陆军工程大学，
类型：发明
国别省市：江苏,32