劉 麗，何志偉

(商丘工學(xué)院，河南 商丘 476000)

從地面到60km高空激光傳輸?shù)钠D(zhuǎn)特性

劉 麗，何志偉

(商丘工學(xué)院，河南 商丘 476000)

激光在大氣中傳輸要考慮大氣折射的影響。文中給出了光學(xué)頻率范圍內(nèi)，對(duì)流層中地球大氣的空氣折射率模型，并給出近似計(jì)算橫向誤差的公式。利用MATLAB計(jì)算出在大氣層60km處、不同波長(zhǎng)的情況下，激光偏轉(zhuǎn)角和橫向誤差隨初始仰角的變化關(guān)系，認(rèn)為激光的偏轉(zhuǎn)主要依賴于初始仰角而不是激光的波長(zhǎng)。

激光傳輸；大氣折射；激光偏轉(zhuǎn)角；橫向誤差；仰角誤差

0 引言

近年來(lái)激光技術(shù)在軍事上的應(yīng)用發(fā)展很快[1-4]，尤其是利用激光對(duì)空中運(yùn)動(dòng)目標(biāo)進(jìn)行跟蹤、瞄準(zhǔn)、探測(cè)或打擊。許多激光應(yīng)用技術(shù)中涉及激光在大氣中的傳輸，因此激光在大氣中傳輸特性的研究吸引眾多學(xué)者對(duì)該領(lǐng)域進(jìn)行研究。

本文分析了激光在大氣中傳輸時(shí)，由于大氣折射等線性非線性問題[5-6]，對(duì)光束傳輸和光束質(zhì)量的影響，并編出程序得出激光在大氣中傳輸?shù)臄?shù)值模擬，指出了不同波長(zhǎng)的光束在大氣中傳輸時(shí)對(duì)偏轉(zhuǎn)角的影響，以及初始仰角對(duì)偏轉(zhuǎn)角的影響，從而得出較為系統(tǒng)的激光大氣傳輸受大氣折射影響的偏轉(zhuǎn)特性。

1 激光在大氣中傳輸時(shí)的偏轉(zhuǎn)

在涉及到大氣折射率的問題時(shí)，一般用修正的Edlén公式計(jì)算大氣折射率，在修正的Edlén公式[7]中，大氣折射率是由溫度、壓力、濕度和光波長(zhǎng)決定的。

(1)

n是在一個(gè)溫度Tt(℃)和大氣壓力P(Pa)下的折射率，(n-1)s是從(2)式的修正色散公式中得到的。

(n-1)s×10-8=8342.54+2406147(130-1/λ2)-1+15998(38.9-1/λ2)-1

(2)

λ是光波波長(zhǎng)，單位是μm。

對(duì)于潮濕空氣的折射率的差別，包含水汽壓和干燥空氣的壓力總和，差量被寫為：

nf=n-Pf(3.7345-0.00401/λ2)×10-10

(3)

前面只討論了折射率的垂直分布，而未涉及水平分布問題。一般而言，氣象要素在水平方向的不均勻性，除了在水陸交界面、山區(qū)以及存在一些天氣過(guò)程的地區(qū)外，同垂直不均勻性相比通常要小幾個(gè)數(shù)量級(jí)。即使在那些特殊區(qū)域，在一定高度以上(一般是離地面500m以上)水平方向也將趨向均勻化，因此許多實(shí)際問題中認(rèn)為大氣是水平或球面分層是合理的[8]。

(4)

從圖1中很容易得到下面的幾何關(guān)系：

(i=1,2,…,M)(5)

(i=1,2,…,M)(6)

α0=π/2-β0，在第i層和第i+1層大氣層分界處，激光光束的折射由下式?jīng)Q定：

ni(Δh)sinθi=ni+1(iΔh)sinαi

(7)

φi=αi-θi

(8)

(9)

(10)

其中，

(11)

(12)

考慮到熱效應(yīng)對(duì)大氣折射率的影響，最大的總的大氣偏斜距離大約等于兩個(gè)偏斜距離的和，即

ΔSM=ΔS+ΔSlaser

(13)

ΔSlaser是由于激光熱效應(yīng)引起的偏斜距離。

按照式(1)的大氣折射模型,根據(jù)理論計(jì)算分析式(9)，利用MATLAB仿真模擬得出結(jié)果如圖2所示。其條件是在大氣壓為1013mPa，溫度為293k的情況下進(jìn)行計(jì)算的(地球平均半徑為6371.004km，大氣層厚度為60km，大氣層分層數(shù)m為10000)。此圖說(shuō)明了偏斜角Φ隨初始仰角β0的變化。從圖中可以比較得出初始仰角β0越小，偏斜角Φ越大；并且偏斜角Φ在高度為20km、30km、50km、60km時(shí)幾乎相同，這是因?yàn)檎凵渎试谶@些高度近似真空的值。所以激光向空中發(fā)射時(shí)要選擇大的仰角，這樣命中率才會(huì)高。

根據(jù)公式(10)，可以得到在不同高度的偏斜距離隨初始仰角變化的函數(shù)關(guān)系，如圖3所示。根據(jù)圖3可以看出某一高度的光束偏斜距離與初始仰角之間有很大的關(guān)系，尤其是在仰角為0°～15°時(shí)的情況下，偏斜距離很大。例如，當(dāng)初始仰角β0是60°時(shí)，在60km高度的偏斜距離ΔS是15m，而在同樣的高度β0為10°時(shí)的ΔS是2000m。

利用公式(10)還可以得到偏斜距離ΔS隨著海拔高度H變化的函數(shù)關(guān)系，如圖4所示。從圖中可以看出，大氣折射率對(duì)光束傳播的影響不能忽視，特別是從地面到高空的激光定位和激光探測(cè)方面的應(yīng)用。

2 激光在大氣中傳輸時(shí)的橫向誤差

激光在大氣中傳輸時(shí)，由于射線發(fā)生彎曲而產(chǎn)生到達(dá)點(diǎn)偏離目標(biāo)點(diǎn)，到達(dá)點(diǎn)和目標(biāo)點(diǎn)之間的距離即為橫向誤差[10]，在距離60km的高空中，橫向誤差近似為：

ΔR=r0×(θ0-α0)

(14)

式中α0是發(fā)射點(diǎn)與出射點(diǎn)之間的視線仰角。

(15)

式中R0為地球平均半徑，r0為大氣層厚度約為60km，φ為出射點(diǎn)繞地心的偏轉(zhuǎn)角也叫地心張角。圖5，6，7分別是激光波長(zhǎng)為1.06μm，3.8μm，10.6μm在大氣中傳輸60km時(shí)的橫向誤差，其計(jì)算條件同樣也是在大氣壓為1013mPa、溫度為293K時(shí)，從這3個(gè)圖中可以看出激光的波長(zhǎng)對(duì)橫向誤差的影響非常小。

3 大氣層折射率對(duì)激光出射點(diǎn)仰角的影響

由Snell定律

n0Rsinθ0=C

(16)

可求得整個(gè)大氣層外邊緣處光束傳輸方向，即出射點(diǎn)仰角θm的值，

sinθm=C/(R0+r0)

(17)

由圖8，9可以看出隨著初始仰角的增大，出射點(diǎn)仰角不斷變大，仰角誤差也不斷增大，這些說(shuō)明射線是彎向地球的，但不可能折射回地面，這符合正常折射。特別是在低初始仰角情況下，激光束在低高度區(qū)域主要受到大氣的折射作用[11]。

4 結(jié)論

利用激光進(jìn)行地對(duì)空的跟蹤、瞄準(zhǔn)、探測(cè)時(shí)，首先要考慮大氣折射的影響，對(duì)大氣折射進(jìn)行修正，這樣才能更準(zhǔn)確地進(jìn)行作用。本文只考慮在近似理想的情況下，給出近似計(jì)算橫向誤差的公式，算出的結(jié)果與實(shí)驗(yàn)數(shù)值基本符合。但沒有把大氣湍流、大氣分子和氣溶膠的吸收散射等作用考慮進(jìn)去。折射效應(yīng)計(jì)算和預(yù)測(cè)修正的關(guān)鍵在于折射指數(shù)與時(shí)間變化規(guī)律的掌握。對(duì)于不同的應(yīng)用，預(yù)測(cè)修正的要求各不相同，所需采用的折射指數(shù)分布模式和修正方法也將不同。但在高精度修正要求的情況下，需要采用實(shí)時(shí)實(shí)地的實(shí)測(cè)分布剖面，甚至要求取得沿射線路經(jīng)的三維結(jié)構(gòu)，進(jìn)行精確的數(shù)值積分計(jì)算，才能滿足誤差修正的精度要求。

[1] 夏新仁,馮金平.激光制導(dǎo)武器的現(xiàn)在與將來(lái)[J].中國(guó)航天,2009,(12):22-26.

[2] 莊昕宇,陳兆兵.半主動(dòng)激光精確末制導(dǎo)武器的發(fā)展現(xiàn)狀與趨勢(shì)[J].艦船電子工程,2011,31(6)：6-10.

[3] 王狂飆.激光制導(dǎo)武器的現(xiàn)狀、關(guān)鍵技術(shù)與發(fā)展[J].紅外與激光工程,2007,36(5)：651-655.

[4] 謝莉,孫昊,賈玉林.激光制導(dǎo)武器的對(duì)抗技術(shù)研究[J].航天電子對(duì)抗,2010,26(2)：19-22.

[5] 梅海平,饒瑞中,吳曉慶,等.大氣折射率結(jié)構(gòu)常數(shù)和湍流特征譜的計(jì)算[J].強(qiáng)激光與粒子束,2003,(12)：1155-1158.

[6] 孫剛,翁寧泉,肖黎明,等.大氣溫度分布特性及對(duì)折射率結(jié)構(gòu)常數(shù)的影響[J].光學(xué)學(xué)報(bào),2004,24(5)：592-596.

[7] K.P.Birch.Correction to the Updated Edlén Equation for the Refractive Index of Air[J].Metrology,1994,31(4):315-316.

[8] 翁寧泉,曾宗濤,肖黎明,等.大氣折射率結(jié)構(gòu)常數(shù)垂直分布特征[J].強(qiáng)激光與粒子束,1999,12(11):673-676.

[9] 李雙剛,聶勁松,孫曉泉.大氣折射和色散對(duì)激光傳輸?shù)挠绊慬J].量子電子學(xué)報(bào),2004,21(5):679-682.

[10] 路遠(yuǎn),凌永順,聶勁松.激光穿越大氣層時(shí)的折射誤差研究[J].紅外與激光工程,2003,32(1):69-72.

[11] Yufeng Peng,Li Liu,V.S.Egorov, et al.Propagation offset characteristics of annular laser beams from confocal unstable resonators through the natural atmosphere[J].Optics Communication,2008,(281):705-717.

Deflection Characteristics of Laser Propagation from the Ground to 60km-High

LIULi，HEZhi-wei

(ShangqiuInstituteofTechnology,Shangqiu476000,China)

When laser is transmitted in the atmosphere, we should first consider the influence of atmospheric refraction on laser transmission. This paper gives the refraction rate model of air in the troposphere of the earth within optical frequency range. In the atmosphere at 60km high and under different wavelengths, the connection of laser deflection angle and lateral error changing with the initial elevation is computed with MATLAB. The deflection of the laser beam is mainly dependent on the initial elevation angle instead of the laser wavelength.

laser propagation; atmospheric refraction; laser deflection angle; lateral error; elevation error

2016-11-20

2015年河南省教育技術(shù)裝備和實(shí)踐教育研究立項(xiàng)課題(GZS134)

劉麗(1981-)，女，碩士，商丘工學(xué)院基礎(chǔ)教學(xué)部講師，研究方向：激光物理與器件。

TN248.1

A

1674-3229(2017)01-0064-05