蔣大鋼，鄧 科，黃 健，朱 彬，秦開宇

(電子科技大學(xué)航空航天學(xué)院，成都611731)

1 引言

相對射頻通信而言，激光通信具有體積小、質(zhì)量輕、功耗低、抗干擾、抗截獲等優(yōu)點，是一種極具發(fā)展?jié)摿Φ奶旎ㄐ偶夹g(shù)。近年來，衛(wèi)星激光通信技術(shù)發(fā)展迅猛，中國、德國、美國、日本等多個國家都已經(jīng)成功開展了星間、星地激光通信試驗［1-3］。

激光通信與激光測距具有許多技術(shù)共性，如:激光收發(fā)單元，光束捕獲、跟蹤和瞄準單元，時鐘模塊等。通過應(yīng)答式測距方式可以實現(xiàn)通信與測距功能的最大程度復(fù)用，這不僅能夠提高有效載荷集成化性能，還能為航天測控通信提供一種新興的技術(shù)體制［4］。

國外已經(jīng)在空間探測任務(wù)中使用了激光應(yīng)答通信測距技術(shù)，測距精度約半米，如:2013年NASA計劃開展第一個月球激光通信演示系統(tǒng)(LLC其D計劃)［5］;2010年由 ESA和 NASA聯(lián)合開展的引力波探測計劃(LISA計劃)［6］;2005年NASA 開展的用于哥達德飛行中心與水星飛行器Messenger測控的激光雙向異步應(yīng)答測距試驗［7］等。

近年來，國內(nèi)在激光通信測距領(lǐng)域的研究已經(jīng)從最初的原理樣機實現(xiàn)［8-10］發(fā)展到了地面演示驗證開展［11-13］，新研究成果逐年產(chǎn)生，研究水平也在迅速提升。

可以預(yù)見，在載波頻段統(tǒng)一的天地一體化新測控體制的發(fā)展過程中，深空、星間、星地激光通信測距技術(shù)的發(fā)展需求和潛力都是十分巨大的，但星地激光通信測距技術(shù)實現(xiàn)上還諸多需要探索并解決諸如大氣光傳輸?shù)葐栴}，在諸多大氣光傳輸效應(yīng)中，大氣折射對于激光通信測距鏈路的影響是最基本問題之一。

文獻［14-15］從星地激光通信技術(shù)出發(fā)，研究了大氣折射星地鏈路的影響，并建議在天頂角60°范圍內(nèi)開始捕獲，而且在這樣的天頂角范圍內(nèi)可以進一步忽略大氣折射對光束瞄準的影響［15］。

文獻［16-17］從星地激光測距技術(shù)出發(fā)，基于經(jīng)典大氣折射修正模型Marini-Murray模型，研究了大氣折射對激光光程差的影響。Marini-Murray模型的修正精度隨天頂角的增大而提高，在天頂角10°時，修正殘差約0.9 mm;當天頂角為80°時，修正后殘差約1.6 mm。

星地激光通信測距鏈路的應(yīng)用特點是從航天測控需求出發(fā)，滿足航天器長時段觀測和定軌需求，因此在應(yīng)用層面，需要能夠在大天頂角條件下建立鏈路;在技術(shù)層面，需要進一步探索一種適用于激光通信測距的大氣折射角度和距離修正模型。本文將以從大氣層同心球殼模型出發(fā)，采用光束追跡法分析大氣折射對光束傳輸?shù)挠绊?，進一步梳理星地激光通信測距鏈路的大氣折射補償方案。

2 工作原理

應(yīng)答式激光通信與測距系統(tǒng)的工作原理是檢測測距碼在主測端和被測端之間的傳遞時間，并扣除兩端的電路處理時間，得到測距碼在兩個終端間的飛行時間，從而可以估計兩個終端距離［8-13］。具體實現(xiàn)分為2個步驟:一是通信鏈路建立;二是通過通信鏈路進行測距碼轉(zhuǎn)發(fā)。

星地激光通信鏈路，其建立過程示意圖如圖1所示。

圖1 星地激光通信鏈路建立過程示意圖Fig.1 Illustration of link establishing process for LEOGround laser communication

地面站在衛(wèi)星可能出現(xiàn)的不確定區(qū)域FOU內(nèi)進行上行信標光掃描。當衛(wèi)星接收到地面信標光后，向相應(yīng)方位返回一束下行信標光。通過信標光握手，雙方獲取對方位置信息，完成捕獲過程，進而通過跟蹤和瞄準過程建立星地激光通信鏈路。

基于星地激光通信鏈路的測距碼轉(zhuǎn)發(fā)原理(圖2)，測距碼轉(zhuǎn)發(fā)過程包括:(1)根據(jù)任務(wù)需求和通信協(xié)議規(guī)范，若需要開展測距任務(wù)，主測端立即在數(shù)據(jù)幀中的特定位置插入測距碼并準備將數(shù)據(jù)發(fā)送給被測端。當主測端測距碼處理器檢測到待發(fā)數(shù)據(jù)中存在測距碼時，便在T0時刻開啟測距計數(shù)器;(2)被測端接收數(shù)據(jù)后，如果當T1時刻被測端檢測到數(shù)據(jù)幀的測距碼，便立即在T2時刻給主測端傳輸?shù)臄?shù)據(jù)幀中插入測距碼;(3)T3時刻，主測端檢測到被測端數(shù)據(jù)流中的測距碼后，便關(guān)閉測距計數(shù)器，通過傳輸時延T3-T0，并扣除電路延遲T，獲得兩個終端距離估計T=0.5c(T3-T0-T)。

在星地激光通信測距鏈路中，大氣折射一方面使得光束傳播方向發(fā)生變化，影響光束捕獲和跟瞄;另一方面會光程變得更長，影響實際距離估算，下面將通過光束追跡法量化分析，并提出相應(yīng)補償措施。

圖2 測距碼轉(zhuǎn)發(fā)原理Fig.2 Retransmission principle of ranging codes

3 大氣折射影響分析

地球大氣層分布可認為是在重力場條件下形成的同心球殼分布(圖3)。假設(shè)與地心距離相等的位置的大氣折射率相同，則大氣層可分為多個薄層。當光束穿越每個薄層時，會因為大氣折射率不同而產(chǎn)生光束折射。光束折射一方面會造成角度對準誤差;另一方面會增加光傳輸路程，影響光程估計。

假設(shè)距地心Ri處大氣厚度為dh的薄層內(nèi)折射率為ni，光束從前一層進入該層的折射角為θi，由該層進入折射率為nn=1的后一層的入射角為，折射角為 θi+1，則根據(jù) Snell定律和正弦定理有:

圖3 大氣折射模型Fig.3 Model of atmospheric refraction

其中，Ri=Re+h，Ri=1=Ri+dh，地球半徑Re，h為折射率ni的薄層的海拔高度。

光束在經(jīng)過折射率為ni的薄層進入折射率為ni=1的薄層時產(chǎn)生的折射角為:

則當光束穿越整個大氣層后，總折射角可表示為:

同樣當光束穿過折射率為ni的薄層進入折射率為ni+1的薄層時，行進的光程di可表達為:

則光束穿越整個大氣層后，總的行進光程與同等距離的真空光程的光程差Δd為:

則可根據(jù)式(4)即可分析大氣折射角度，并在通信鏈路建立過程中進行折射角度補償，根據(jù)式(6)即可分析由于大氣折射產(chǎn)生的附加光程差，并在測距鏈路中進行光程差補償。

4 光束追跡法仿真分析

光束追跡法仿真需要獲取大氣折射率分布模型。

大氣折射率是大氣溫度、濕度、氣壓和波長的函數(shù)。濕度對折射率的影響很小，可以忽略。根據(jù)文獻［14，18］，可以根據(jù)大氣溫度和壓強隨海拔的分布模型間接推算大氣折射率隨海拔分布模型。

折射率n與波長λ(μm)、氣壓P(Pa)及溫度T(K)的關(guān)系為:

其中，T及P與海拔高度有關(guān)，可以換算出n隨海拔的分布。

T(h)定義為:

其中，T0為海平面溫度。取各個溫度轉(zhuǎn)折點海拔高度ha=11 km，hb，=22 km，hc=32 km，hd=50 km，溫度變化率 Δ1=-6.5 K/km，Δ2=1.0 K/km，Δ3=3.0 K/km。

P(h)定義為:

其中，P0為海平面處標準大氣壓。

仿真時，設(shè)置光載波波長λ=1550 nm，對于任意給定的初始發(fā)射光束天頂角θ0，通過光束追跡法可以得到各層的入射角θi+1、折射角θ'i和光程di，進而逐層累積得到大氣折射角和光程差。

仿真得到大氣折射角隨天頂角變化關(guān)系(圖4)。當天頂角大于80°時，大氣折射迅速增加到數(shù)個毫弧度。當天頂角小于75°時，大氣折射處于亞毫弧度量級。

圖4 大氣折射角隨天頂角變化關(guān)系Fig.4 Relationship between atmospheric refraction and zenith angle

這表明星地激光通信鏈路在大天頂角條件下建立鏈路時，必須首先補償大氣折射對捕獲過程的影響，其次還需要進一步補償大氣折射對光束跟蹤瞄準的影響，否則難以建立通信鏈路，測距功能更無從談起。因此，在光束捕獲、跟蹤和瞄準過程中，可以結(jié)合不同天頂角下的大氣折射理論仿真結(jié)果，在光軸對準算法中直接進行大氣折射角補償;也可借鑒文獻［14］提出的方案，采用增加捕獲不確定域的方式進行補償，但這樣會延長鏈路建立時間，因為典型不確定域是毫弧度量級，與大氣折射角量級相當。

仿真得到的大氣光程差修正值隨天頂角變化關(guān)系，如圖5所示。當天頂角大于80°時，大氣光程差修正值迅速增加，達到數(shù)十厘米。當天頂角小于75°時，大氣光程差修正值為1～4 cm。

圖5 大氣光程差修正值隨天頂角變化Fig.5 Relationship between optical path difference and zenith angle

這就意味著星地激光通信鏈路建立后，進入測距、測控環(huán)節(jié)時，必須進行大氣光程差補償，將光束的大氣傳輸距離折算成真空距離，從而實現(xiàn)精密測距。典型的應(yīng)答式通信測距精度為半米量級，而大氣附加的額外光程差在大天頂角條件下就達到了0.1 m量級，如果不進行大氣光程差補償，將大大降低測距精度，進而影響測控精度。

5 結(jié)論

通過光束追跡法仿真表明:大氣折射對星地激光通信測距鏈路的影響主要體現(xiàn)在附加折射角和附加光程兩個方面。在大天頂角(＞75°)情況下，附加折射角對通信測距鏈路建立的影響尤其明顯，可結(jié)合大氣折射理論仿真結(jié)果，在光軸對準算法中對大氣折射進行直接補償;附加光程影響主要體現(xiàn)在測距階段，典型的應(yīng)答式通信測距精度為半米量級，而大氣附加的額外光程差在大天頂角條件下就達到了0.1 m量級，必須根據(jù)大氣光程差理論仿真結(jié)果進行補償，以保障測控精度。

綜合以上分析，大氣折射對星地激光通信測距鏈路的影響貫穿在整個運行過程中，在大天頂角情況下表現(xiàn)尤為明顯。目前尚無針對性的大氣折射修正模型。但是通過本文分析可以明確的是:獲取實際大氣折射率分布廓線是實現(xiàn)星地激光通信測距鏈路大氣折射補償?shù)年P(guān)鍵。在星地激光測控技術(shù)發(fā)展過程中，有必要加強大氣折射率分布廓線實時監(jiān)測技術(shù)研究。

［1］ Toyoshima M.Trends in satellite communications and the role of optical free-space communications［J］.Opt Networking，2005，4(6):300-311.

［2］ Toyoshima M，Jono T，Suzuki K，etc.Laser beam propagation in ground-to-OICETS laser communication experiments［C］.//Proc.SPIE.Orlando，F(xiàn)londa，USA，2007，(6551):65510A1-A12.

［3］ Fields RA，Kozlowski DA，Yura HT，etc.5.625 Gbps bidirectional laser communications measurements between the NFIRE satellite and an ptical ground station［C］.//Proc.SPIE.Proc.SPIE.Preague，Czech Republic，2011，(8184):81840D1-D17.

［4］ 錢衛(wèi)平.航天測控通信新體制研究［J］.載人航天，2009，15(1):1-8.

［5］ Xiaoli Sun，Skillman DR，Hoffman ED，etc.Free space laser communication experiments from Earth to the Lunar Reconnaissance Orbiter in lunar orbit［J］.OPTICS EXPRESS:2013，21(2):1865-1871.

［6］ Juan Jos'e Esteban，Antonio F.Garc'a，Simon Barke，etc.Experimental demonstration of weak-light laser ranging and data communication for LISA［J］.OPTICS EXPRESS:2011，19(17):15937-15946.

［7］ Neumann GA，Cavanaughl JF.Laser ranging at interplanetary distances［C］//Proceedings of the 15th ，International Workshop on Laser Ranging.Carberra，Austulia，Ctober 2006:451-456.

［8］ 楊紅宇，鄧科，幺周石等.星間衛(wèi)星通信與測距新技術(shù)研究［J］.光學(xué)技術(shù)，2006，30(11):157-162.

［9］ 胡瑋，蔣大鋼，鄧科等.無線激光通信與測距一體機的測距精度分析［J］.紅外與激光工程，2008 S3:1231-1235.

［10］ Ke Deng，Dagang Jiang，Wei Hu，etc.A novel technology combined with free space optics communication and laser ranging［C］.//Proc.SPIE.Hangzhou，China，2008，(7136):7136301-5.

［11］ 李艦艇，邢強林，王治中等.衛(wèi)星激光應(yīng)答測距關(guān)鍵技術(shù)實驗研究［J］.紅外與激光工程，2009(S38):110-113.

［12］ 唐嘉，高昕，邢強林等.異步應(yīng)答激光測距技術(shù)測量精度實驗［J］.紅外與激光工程，2011，40(5):939-943.

［13］ 唐嘉，高昕，邢強林等.異步應(yīng)答激光測距技術(shù)［J］.光電工程，2010 37(5):25-31.

［14］ 向勁松.地-星上行激光通信中大氣折射及色散引起的對準誤差［J］.重慶郵電大學(xué)學(xué)報 (自然科學(xué)版)，2008，20(4):423-426.

［15］ 榮健，鄧代竹，鐘曉春.大氣折射對空地衛(wèi)星光通信鏈路影響的研究［J］.激光雜志，2004，25(5):49-50.

［16］ 袁宏武，梅海平，黃印博等.星載激光測距大氣校正算法與模型研究［J］.光學(xué)學(xué)報，2011 31(4):1-7.

［17］ 祝立，吳斌.利用衛(wèi)星激光觀測對大氣折射模型的綜合比較［J］.天文研究與技術(shù)，2011，8(3):262-267.

［18］ Clifford SF.The classical theory of wave propagation in a turbulent medium［J］.Topics in Applied Physics，1978(25):9-43.