刘大畅,曾 晶,孙石峰,谢寅亮,吴绍华(.哈尔滨工业大学(深圳),广东深圳 58055;.广东南方电信规划咨询设计院有限公司,广东深圳 58038)
空间光通信中的光波随机特性相较于射频通信和光纤通信,会受到大气信道中光湍流的影响,需要充分理解激光在大气信道中的传播理论特性,进而分析光湍流对空间光通信系统性能的影响机理。早在上世纪六十年代,国内外学者就展开了针对空间光通信系统大气信道理论的研究,经过数十年的研究积累和技术发展,逐渐形成了光湍流对空间光通信系统影响的分析理论[1-2]。
光波在大气湍流中传输的基础理论是由Chernov L A和Tatarskii V I通过理论与实验的一致性验证的,在1960年和1961年分别提出了大气湍流中平面波和球面波的传输波动方程[3-4]。直至2004年,Andrews L C提出的理论打破了经典理论的局限性,提出了在弱湍流和强湍流2种情况下,接收端强度起伏的二阶和四阶统计模型,从而得到光强起伏方差计算模型,进一步推导出平面波和球面的闪烁指数模型[5]。2008年Wu B等分析了在烟尘以及沙粒环境中传输的光信号的色散变换,并推导了接收光信号空间色散与角色散的关系[6]。
目前国内的科研院校等研究团队在激光大气传输方面也取得了一些研究进展。比如中科院安徽光机所在2008年研究了部分相干光在大气湍流中的传输特性。哈尔滨工业大学研究团队针对不同调制方案通信系统,分析并比较了光波在大气信道中的性能。2009年长春理工大学团队系统地总结了各种大气湍流效应对光通信的影响形式[7-8]。
本文首先阐述空间光通信的不同传输光类型经过大气湍流的传输特性区别,然后从光场表达式出发,推导了描述光束经过湍流的2个重要参数光束尺寸和相前曲率半径的表达式,并模拟了接收端光束传输扩散特性在不同相干性状态下随距离的变化趋势,最后针对空间光通信系统整体链路,根据接收端的信号功率以及数据速率的表达式,从多个方面推导了链路余量和链路可靠性的表达式。
在空间光通信系统中,经过大气湍流的信号到达探测器端时,信号强度会产生波动,从而影响整个通信系统的性能。为了减少空间光通信链路中大气湍流对光信号的影响,需要研究不同类型传输光的变化特征(见图1),如平面波、球面波、高斯光束、贝塞尔光束等,选择应用在空间光通信系统中的最佳光束类型[9]。
a)平面波是指其相位面是平行的电磁波,当平面波经过随机介质传播到达探测器表面时,其振幅和相位都发生了波动。
b)球面波是指相位面是球面的电磁波,当空间光通信的源被认为是点光源的时候可以假定为球面波。
c)高斯光束是一种横向电磁场和强度接近高斯方程的电磁波,最重要的特点是高斯光束经过透镜后依然是另一个高斯光束。大多数激光器的基横模和激光器光学谐振腔出射的基模都是一种近似高斯光束。
d)非衍射光束是指完全消除衍射的理想光场,这种光束通过大气湍流受到的影响比其他光束要小,适合于长远距离的空间光通信系统。在实际中,只能实现接近理想的非衍射光束,其光束强度取决于传播坐标,但在定义区域长度的传播区域内,光束轮廓几乎保持不变。
图1 空间光通信常见光束类型及其经过大气湍流的情况
e)贝赛尔光束是一种非衍射的自愈型光束,如果一部分光束在一点被遮挡了,沿着光轴方向的下一点会重新形成光束,这一点特性在空间光通信领域将非常有应用价值。在实际中,可以用轴棱锥或圆锥透镜聚焦高斯光束来产生近似的贝塞尔光束。
对于空间光通信系统来说,当激光束与大气湍流相互作用时,路径中大气的随机波动导致激光束功率在时间和空间上发生随机的起伏,从而增加了通信系统的误码率。一方面,光湍流模型可以估算特定地理位置和特定时间下的光湍流水平,另一方面,基于激光束在大气信道的传输理论框架来讨论激光束在光湍流传输过程中的真实情况,从而可以充分描述光湍流对空间光通信的影响,推算在各种应用场景下,激光束在大气信道的各方面性能变化,最优系统性能条件下,研发人员可以很好地权衡各种设计细节参数,并更好地选择通信系统的站点位置。
首先考虑激光束在大气空间的传输情况,设激光束是一个简单的抛物线波前相位包络的高斯光束模型,该模型广泛应用于激光束在随机介质中传播状况的研究,可以完全表征聚焦或发散的激光束[10]。那么沿z轴传播的单位振幅、低阶轴旁高斯光束在z=0处的大气空间电磁场表达式为:
式中:
ω0——发射光束半径
R0——相前曲率半径
k——光波数,k=2π/λ
r——发射面光束中心的横向距离,r=(x2+y2)1/2
光束传输了距离z之后的光场表达式变为:
式中:
ρ——接收面光束光轴的横向距离,ρ=(x'2+y'2)1/2
式(2)中描述发射光束的2个关键参数表达式分别为:
式中:
在这个表达式中,会聚光束用正R0表示,离焦光束用负R0表示。在束腰位置R0是无穷的,光束尺寸取最小值。根据聚焦参数的大小分类,<1时称为会聚光束,=1时称为准直光束,>1时称为离焦光束。另外和都是关于传播距离z的函数。那么大气空间激光束在接收端的特性表达式为:
根据波数尺寸ω(z)得到接收端的平均强度表达式为:
当光束经过存在光湍流的大气传输时,经过z距离传输后的光束尺寸表达式变为:
式中:
ξ——总相干参数,,其中是在光湍流下传播的球面波相干长度
ξs——携带信号的激光束在发射端的空间相干性,ξs=1代表相干光束,ξs>1代表部分相干光束
图2 接收端光束传输扩散特性随距离的变化趋势
由图2可知波长为1 550 nm的聚焦光束经过中湍流传输3 km后,发射光束尺寸从2 cm变成7.67 cm,但理想的相干光束条件很难实现,一旦ξs增加到100、200、300、400、500,对应的光束尺寸分别为74.07 cm,104.73 cm,128.26 cm,148.10 cm,165.58 cm。通过在发射孔径前直接放置一个相位扩散器,可以增加接收端的波数,但同时也会减少入射功率,所以不能一味的增加接收光尺寸,要通过适当的接收光束尺寸来平衡波数和功率[7][11]。
评价空间光通信系统的性能还可以参考很多参数,比如功耗、发射功率有效性、大气通道可用带宽、系统信噪比、数字通信系统的误码率等。空间光通信系统的整体性能还可以通过链路预算公式来量化表示,包括链路余量和链路可靠性,该公式利用各种通信系统的参数,包括发射功率、传输损耗、接收灵敏度等在不同通信场景下进行链路计算。在空间光通信大气信道中,总体的路径损失一部分是几何损失,另一部分是大气的衰减造成[12]。几何损失是由于激光光束的发散性,导致只有一部分光功率能够进入探测器,如图3所示。
下面研究计算光通信链路性能的几项必要参数,考虑两点之间的大气空间光通信,假设激光发射机天线的增益是GT=其中θT是发射机发散角,传输总功率是PT,发射机光效率是τT,在通信系统探测端接收到的信号功率表达式为[13]:
图3 通信链路几何损失示意图
式中:
τATM——在激光发射波长处的大气传输值,在波长λ处与大气衰减参数α的关系为α=-10ln(τATM)/L
S——大气空间损耗,S=(λ/4πL)2,其中L是传输距离
GR——接收天线增益,GR=(πD/λ)2,其中D是接收机直径
τR——接收机光效率
那么接收信号的表达式可以推导为:
链路余量也是光通信链路的一个重要参数,表示在一定速率下可接收功率与达到特定BER接收功率的比。已知实际接收到的光功率P'R与速率R和接收灵敏度Nb的关系为P'R=NbREp,数据速率表达式为:
式中:
τopt——发射机和接收机的光效效率
A——接收机孔径面积
Ep——光子能量
那么链路余量M的表达式为:
另外一个链路参数是光链路可靠性,它表示在一定BER下保证链路的功能,用可接受的功率范围来定义其动态范围。假设用传输损耗Lp代表所有平均损失,dBm为单位,可以得到光链路可靠性LM表达式:
设定一个饱和值PS,当接收功率PR>PS时接收机是饱和的,如果PR<Nb就是没有达到信噪比。另外链路可靠性可以通过数据速率超过需要值的时间百分比Tav(%)来量化表示,也可以用由大气效应引起的额外功率损耗LA<LM的概率来表示,假设衰减系数αΑ(dB/km)的概率密度用p(αΑ)来表示:
其中αΑlim是衰减系数的极限值,对于距离L,p(αΑ)的概率分布可以通过长时间测量实际链路接收到的信号电平或在同地点之前测量的数据来确定[14]。如果知道通信系统所处大气信道αΑ的概率分布,通过衰减系数极限就可以得到链路可靠性。由此可以根据链路余量和链路可靠性为空间光通信链路预算分析提供理论依据。
本文首先详细介绍了常应用于空间光通信的不同传输光类型,包括平面波、球面波、高斯光束、贝塞尔光束等经过大气湍流的传输特性区别,为应用在空间光通信系统中的最佳光束类型提供选择依据;紧接着,基于光场表达式推导了描述光束经过湍流的2个重要参数表达式,一个是光束尺寸,另一个是相前曲率半径,并在中湍流条件下,模拟了发射光束尺寸为2 cm的1 550 nm激光束,在不同相干性状态下接收端光束传输扩散特性经3 km传输的变化趋势,为今后设计接收端细节参数以及选择通信系统的站点位置提供理论分析方法;最后针对空间光通信系统整体链路,根据发射功率、传输损耗、接收灵敏度、接收端信号功率以及数据速率的表达式,从多个角度推导了链路余量和链路可靠性的表达式,为空间光通信链路预算分析提供理论依据。
另外,大气信道中空间光通信的性能研究还可以从其他几个方面入手,比如信道编码和信道建模,对空间光通信信道模型及信道容量进行分析,研究光在大气信道的衰落和随机变化特性,为进一步改善大气湍流对传输光的影响提供理论支持。