冰晶粒子对量子干涉雷达探测性能的影响 下载: 887次
1 引言
量子雷达是一种将量子信息光学和量子动力学引入遥感探测领域的新型技术[1-4],本质是利用光量子作为光频中的电磁波微观粒子,实现探测目标、检测信号、识别并分辨射频隐身平台和**系统等。此外,量子雷达还可以进行行星防御和空间探测。随着量子信息科学技术的不断发展,美国海军研究办公室(ONR)、哈佛大学、以及洛克希德·马丁公司等对量子雷达和远距离量子遥感技术都已经展开了研究[5-6]。Lanzagorta[7]在《量子雷达》中,引入了量子雷达的概念,采用量子动力学阐述并讨论了不同类型的雷达,包括量子干涉雷达、量子照射和量子激光雷达。其中量子干涉雷达属于纠缠光子雷达系统,主要使用光量子的纠缠相干原理提高量子雷达的性能[8]。冰晶粒子是一种在大气物理运动过程中产生的不规则粒子,主要存在于大气卷云层(高度为4.5~10 km)中,其形状大小和高度、温度、相对湿度以及在云中是否发生碰撞与合并有关[9],其复折射率、消光系数、散射系数、单次散射反照率、不对称因子以及单次散射相函数等散射特性与云型、光散射计算方法有关[10]。冰晶粒子具有的理化特性会使量子干涉雷达系统在大气空间中探测时,系统中的光量子信号产生散射、折射和吸收等光学作用。在量子干涉雷达发射的探测光子从目标反射回来时,反射光量子会经过包含冰晶粒子的大气环境,而冰晶粒子自身的散射特性会影响光量子的偏振态,不仅造成光量子信号强度的耗散,还会引发相位延迟,降低干涉雷达的探测性能。目前关于量子雷达的性能研究与方案设计已经比较完善。文献[ 11]对应用于非平面凸目标的量子雷达散射截面(QRCS)公式进行了拓展推导,理论计算表明,目标的QRCS主瓣基本不变时,量子旁瓣在原子稀疏排列的目标中更为突出。文献[ 12]研究了不同干扰背景下量子雷达的生存性,并根据不同干扰等级的损伤模型建立了量子雷达生存策略。文献[ 13]利用分光链路模拟大气介质,推导了衰减条件下的QRCS方程,并通过仿真得到单曲面QRCS的主瓣峰值与衰减系数和波长的关系。
关于大气环境对量子雷达的性能研究较少,目前只讨论了大气介质与大气损耗对量子雷达探测机理与性能的影响[14-15],并未探究冰晶粒子对量子干涉雷达性能的影响。本文根据冰晶粒子的光散射特性研究了其对光量子态的偏振影响,用偏振对比度描述雷达光量子在冰晶粒子背景下的偏振态;分析了冰晶粒子的有效尺度与雷达探测光子能量、传输距离之间的关系,并通过仿真分析,定量研究了冰晶粒子对量子干涉雷达分辨率和系统误码率的影响。
2 冰晶粒子对量子干涉雷达探测性能的影响
2.1 冰晶粒子背景下的量子干涉雷达系统
用量子雷达系统进行探测时,经目标反射后的探测光子,在传输时会经过大气卷云层,和其中的冰晶粒子群发生碰撞,导致量子信号强度衰减、能量降低。从而增加雷达系统的误码率,对量子干涉雷达的探测性能造成影响。量子干涉雷达属于纠缠光子量子雷达,采用量子计量技术,通过干涉测量相位的方法,提取纠缠量子态中包含的信息。量子干涉雷达的基本构架和工作原理如
图 1. 量子干涉雷达系统。(a)系统原理;(b) Mach-Zehnder干涉仪
Fig. 1. Quantum interference radar system. (a) Principles of system; (b) Mach-Zehnder interferometer
2.2 冰晶粒子的光散射特性
云物理学中将云中的冰晶粒子浓度(单位体积内冰晶粒子的个数)与尺度变化的关系称作冰晶谱,可用标准的伽马分布描述卷云中冰晶粒子的尺度分布(PSD)[16]
式中,D为粒子的特征尺度,N0为冰晶粒子浓度,μ为形状参数,b为斜率参数。冰晶粒子的高度Hc与有效半径vc的关系为Hc=vcDμ,其中vc和μ取决于冰晶粒子的类型和温度。
根据Van de Hulst近似理论和多项式近似方法[17],冰晶粒子的吸收效率因子为
式中,Dmax为粒子的最大特征尺度,mr、mi分别为折射率的实部和虚部形式,γe=4χemi为非球形粒子的有效尺度参数,χe=2πre/λ为球形粒子的有效尺度参数(半径为re),ε1为外部反射能量占总能量的比率,ε2为近似的不对称因子,η1、η2分别为球形与非球形粒子外部反射的影响系数,Q'abs(γe)为球形粒子的吸收效率因子,Q″abs(γe)为非球形粒子的吸收效率因子。当γe~¥时,Q'abs和Q″abs都趋近于1,系数ε1、ε2、η1、η2都可通过非线性拟合的方法得到。通过Van de Hulst近似理论得到球形粒子(半径为re)的吸收效率因子为
设非球形粒子平均光程为V/A,V和A分别为非球形粒子的等效球面积和体积。(3)式中存在一个可以移除的奇点(γe=0),若γe≤10-4,用泰勒级数的展开表达式可去除奇点。根据反常衍射理论(ADT)得到非球形粒子(平均光程为V/A)的吸收效率因子为
同理可得冰晶粒子的消光效率因子为
式中,ε为二次折射系数,β1为球形粒子的相位推移,β2为非球形粒子的相位推移,且β=arctan[mi/(mr-1)],ρe为有效相位延迟,ρe=4πre
实验仿真了四瓣子弹花(bullet rosetters-4)、六瓣子弹花(bullet rosetters-6)、聚合物(plate)和实心六棱柱(hollow column)四种类型冰晶粒子[18]的吸收效率因子和有效尺度的关系,结果如
图 2. 冰晶粒子的吸收效率因子随有效尺度的变化
Fig. 2. Absorption efficiency factor of ice crystal particles changes with effective scale
2.3 冰晶粒子对光量子偏振态的影响
量子通信系统中一般用量子偏振编码保存信息,对于量子雷达系统,发射端与接收端校准信息的过程中,大气中的冰晶粒子引起的光散射造成的相位延迟会使接收端处的入射光相对原参考系发生偏移,进而对量子误码率造成影响。散射光与入射光的传播方向发生偏离,导致光偏振特性出现退偏振效应。相位延迟是由光学元素和大气条件决定的,会造成偏振椭圆极化、降低信号的偏振对比程度(PCR),如
图 3. 相位延迟造成的椭圆偏振。(a)基态R、T下的椭圆极化光;(b)椭圆角率为ω的偏振光
Fig. 3. Elliptical polarization caused by phase delay. (a) Elliptic polarized light detected by the R、T base; (b) polarized light with ellipticity angle ω
偏振对比度可表示为
式中,a为椭圆的长半轴,b为椭圆的短半轴,a/b为横纵轴比,ω为椭圆率角,表示相位延迟导致的椭圆极化光。从(8)式中可以发现,偏振对比程度随相位延迟的增加而减小。当存在相位延迟时,在雷达接收端仍然使用基态R、T对已经发生偏振的圆极化光进行检测会得到两种结果:水平偏振光子和垂直偏振光子,对应的探测比率为XPCR∶1。由此可得到相位延迟引起的量子误码率为[19]
设雷达发射端光子偏振角为θ0,则发射端的量子态|ϕ>为
式中,|0>和|1>分别为通过水平和垂直端口的光子的状态,θ0与偏振椭圆投影相关,cosθ0|0>和sinθ0|1>分别为水平和垂直分量的偏振相位态。
设偏振角的变化量Δθ为随机值,经过链路传输后的量子态|
不失一般性,以量子态|0>为例,且只考虑基态相同,若Δθ是服从[0,π]的均匀分布,则被误判为|1>的概率为
由偏振态偏移造成的量子误码率就是将|0>解码为|1>的概率,可表示为
2.4 偏振态演化对量子干涉雷达性能的影响
冰晶粒子在大气中发生散射时,其散射角可由Henyey-Greenstein(H-G)相函数表示为[20]
式中,ψ为非对称因子,ψ∈[0,1]为冰晶粒子在大气中前后散射不对称的程度。
对冰晶粒子散射相位函数进行抽样,可得散射后的方位角
式中,c1、c2为[0,1]内的均匀分布的随机数。
假设雷达系统接收端接收到的水平偏振光偏离原坐标系的角度ϕ'=θ+σ+Ωτ。其中σ为雷达系统中发射端与接收端的参考系未进行补偿情况下的对准误差,即雷达发射的光量子与目标进行相对运动产生多普勒效应,进而产生的偏振误差角。τ为受到大气漩涡干扰后形成的偏差角度,Ωτ为τ的波动标准差,即线偏振光在大气空间中传输时湍流导致的光偏振变化[21],可表示为
式中,λ为入射光波长,L为线偏振光在大气湍流空间中的平均传输距离,l为折射率的三维谱密度高斯近似比系数,Δκ为折射率的归一化偏差,rm为空间均匀介质中任两点之间的相对位置。当偏振光在大气空间中传播10 km后,由(16)式计算出Ωτ约为10-7 rad。
当雷达接收端在基态R、T下探测接收信号时,水平与竖直方向能量比Er为
式中,XPP'为椭圆偏振光投影在R轴方向上的偏移轨迹,XQQ'为椭圆偏振光投影在T轴方向上的偏移轨迹。
由(9)式和(17)式可计算出整个量子干涉雷达系统的量子误码率QBER为
量子干涉雷达发射出的探测光子经目标反射,被本地雷达接收后,可利用探测光子与本地参考光子之间的纠缠相关性,得到相应的目标位置与区域信息。由于探测光子的能量受冰晶粒子等大气环境的影响会产生衰减,导致探测光子的传输距离减小,最终影响量子干涉雷达的探测性能。在冰晶粒子背景下,根据布格尔-朗勃特(Bouguer-Lambert)定律,可得到探测光子在传输距离为d时的光量子能量为
式中,Ei为探测光子的初始能量强度,K为衰减系数。取冰晶粒子的粒径范围为0~1000 μm,则衰减系数可表示为
式中,r为冰晶粒子的粒径,n(r)为冰晶粒子的分布函数,Kdot为衰减截面,表示散射截面和吸收截面之和[22]。由光子能量的定义可知
式中,E0为光子能量,普朗克常量h≈6.6×10-34 J·s,λp为光子波长,c为光速。
为了研究冰晶粒子对量子雷达分辨率的影响,设量子雷达的空间分辨率为δRQ=λpδRφ/2π[15]。其中,δRφ为脉冲信号的角度分辨率,与脉冲信号的半峰全宽有关,一般而言,脉冲信号的半峰全宽越小,信号的窄带宽度越小,得到的分辨率也就越高。根据衰减量子干涉测量法计算得到δRφ为
式中,Ci为媒介透射率,C1=exp(-nL1),C2=exp(-KL2),n为大气衰减系数,L1为探测光子在大气中的传输距离,K为冰晶粒子的衰减系数,L2为冰晶粒子群的几何光学厚度,N为发射波束光子数,ϕ0=
式中,TA=exp(-κA)为本地参考光束的能量透过率系数,κA为本地参考光束的大气衰减系数,TB=exp(-κB)+exp
由(18)式可得,当XPCR较大时对系统的量子误码率QBER的影响可以忽略不计。此外,雷达发射端的量子光源[23]和单光子探测器的物理特性,包括系统暗计数、发射光子数量、传输距离、门时间、接收端孔径面积以及探测器效率等,对量子干涉雷达的探测性能也会造成不同程度的影响。
3 仿真结果与分析
实验根据文献[ 24]选取仿真参数自,取探测器的探测效率为30%,普朗克常量h=6.62607015×10-34 J·s,真空光速c为3×105 km/s,接收孔径面积为30 cm2,接收器门时间为200 ps,波长窗口为450~550 nm,探测器视场立体角为π(1-cosXFOV),视场角XFOV=10°。
设入射光波波长为480 nm,冰晶粒子波长为1.315 μm,
图 4. 冰晶粒子的有效尺度对量子误码率影响
Fig. 4. Effect of effective scale of ice crystal particles on quantum bit error rate
图 5. 不同类型冰晶粒子的非对称因子对量子误码率的影响
Fig. 5. Effect of asymmetric factors of different types of ice crystal particles on quantum bit error rate
取探测光子的波长为1.45 μm,形状参数μ=0,冰晶粒子的粒子浓度N0=5.469×108 m-3,斜率参数b=0.018 μm-1,这些参数均由拟合观测数据得到[25]。冰晶粒子的复介电常数可通过文献[
26]中提出的LH模型确定。以λp=1.43 μm的探测光子为例,由(21)式计算出对应的光子能量为1.38×10-6 J,忽略大气损耗等其他因素的影响,将E0代入(19)式中,得到探测光子的传输距离d、冰晶粒子的有效尺度r与光量子信号能量Ed之间的关系,仿真结果如
图 6. 传输距离、冰晶粒子的有效尺度与光子能量的关系
Fig. 6. Relationship between transmission distance, effective scale of ice crystal particles and photon energy
从
取冰晶粒子有效尺度为200 μm,N0=5.469×108 m-3,将能见度为300 m时的大气衰减系数值0.0016 dB/km作为本地参考光束与探测光束的大气衰减系数值[10]。根据(23)式可以得到冰晶粒子浓度、冰晶光学厚度和量子雷达的系统分辨率之间的关系,如
图 7. 冰晶粒子的粒子浓度与量子雷达的分辨率的关系。(a) N=5;(b) N=10
Fig. 7. Relationship between particle concentration of ice crystal particles and resolution of quantum radar. (a) N=5; (b) N=10
由(18)式可知,量子干涉雷达发射的光量子经大气散射,并与冰晶粒子相互作用导致相位延迟,形成椭圆率角。对椭圆率角ω和偏振对比度XPCR与系统QBER的关系进行仿真,结果如
图 8. 椭圆率角与偏振对比度对量子误码率的影响。(a) ω=0°~90°;(b) ω=90°~180°
Fig. 8. Effect of ellipticity angle and polarization contrast on quantum bit error rate. (a) ω=0°-90°; (b) ω=90°-180°
4 结论
采用Van de Hulst近似理论,分析了量子干涉雷达的探测光量子在经过卷云时,冰晶粒子对量子干涉雷达探测性能的影响。引入了标准伽马分布描述冰晶粒子的尺度分布,探讨了不同类型冰晶粒子存在的散射特性差距,研究了冰晶粒子对光量子偏振作用造成的相位延迟与偏振对比度的变化,并结合冰晶粒子的散射相函数推导出系统量子误码率的公式;分析了探测光量子在冰晶粒子背景下的能量衰减对量子干涉雷达探测光子传输距离的影响;分析了冰晶粒子对量子干涉雷达分辨率的影响,并建立了两者之间的函数关系。仿真结果表明,冰晶粒子的有效尺度及非对称因子会对量子干涉雷达误码率产生不同的影响;当处于高偏振对比度时,系统量子误码率随着冰晶粒子非对称因子的增加而增大;冰晶粒子固有的消光效应,会造成探测光子能量的衰减,从而影响探测光子的传输距离;当发射波束的光子数目不变时,量子干涉雷达的系统分辨率随冰晶粒子浓度的增加而降低。因此,在冰晶粒子作为耗散通道时,可根据冰晶粒子的散射特性选取对光量子偏振影响最小的参数,同时适当提高发射波束的光子数,可达到减小量子干涉雷达误码率、增加量子干涉雷达探测光子的传输距离以及提高量子干涉雷达系统分辨率的目的。
[1] 王宏强, 刘康, 程永强, 等. 量子雷达及其研究进展[J]. 电子学报, 2017, 45(2): 492-500.
Wang H Q, Liu K, Cheng Y Q, et al. The advances in quantum radar[J]. Acta Electronica Sinica, 2017, 45(2): 492-500.
[2] 江涛, 孙俊. 量子雷达探测目标的基本原理与进展[J]. 中国电子科学研究院学报, 2014, 9(1): 10-16.
Jiang T, Sun J. The principle and development of quantum radar detection target[J]. Journal of China Academy of Electronics and Information Technology, 2014, 9(1): 10-16.
[3] Jiang K B, Lee H, Gerry C C, et al. Super-resolving quantum radar: coherent-state sources with homodyne detection suffice to beat the diffraction limit[J]. Journal of Applied Physics, 2013, 114(19): 193102.
[4] 吴琼, 柏业超, 张兴敢. 干涉量子雷达测量方法性能分析[J]. 南京大学学报(自然科学), 2016, 52(5): 939-945.
Wu Q, Bai Y C, Zhang X G. Performance analysis of measurement methods for interference quantum radar[J]. Journal of Nanjing University (Natural Sciences), 2016, 52(5): 939-945.
[5] Allen EH, Karageorgis M. Radarsystems and methods using entangled quantum particles: US7375802[P/OL].2008-05-20[2019-07-29]. https:∥patents.glgoo.top/patent/US7375802B2/en.
[6] Myers J M. Speaking of sensing in the language of quantum mechanics[J]. Proceedings of SPIE, 2007, 6573: 657302.
[7] LanzagortaM. 量子雷达[M]. 周万幸, 吴鸣亚, 胡明春, 等, 译. 北京: 电子工业出版社2013: 54- 59.
LanzagortaM. Quantum radar[M]. Zhou W X, Wu M Y, Hu M C, et al., Transl. Beijing: Publishing House of Electronics Industry, 2013: 54- 59.
[8] 肖怀铁, 刘康, 范红旗. 量子雷达及其目标探测性能综述[J]. 国防科技大学学报, 2014, 36(6): 140-145.
Xiao H T, Liu K, Fan H Q. Overview of quantum radar and target detection performance[J]. Journal of National University of Defense Technology, 2014, 36(6): 140-145.
[9] Gallagher M W, Connolly P J, Whiteway J, et al. An overview of the microphysical structure of cirrus clouds observed during EMERALD-1[J]. Quarterly Journal of the Royal Meteorological Society, 2006, 131(607): 1143-1169.
[10] 孙贤明. 大气中离散随机介质的波传播和散射特性研究[D]. 西安: 西安电子科技大学, 2007: 21- 30.
Sun XM. Study on wave propagation and scattering characteristics of atmospheric discrete random media[D]. Xi'an:Xidian University, 2007: 21- 30.
[11] 徐世龙, 胡以华, 赵楠翔, 等. 金属目标原子晶格结构对其量子雷达散射截面的影响[J]. 物理学报, 2015, 64(15): 154203.
Xu S L, Hu Y H, Zhao N X, et al. Impact of metal target's atom lattice structure on its quantum radar cross-section[J]. Acta Physica Sinica, 2015, 64(15): 154203.
[12] 李旭, 聂敏, 杨光, 等. 基于纠缠的量子雷达生存性策略及性能仿真[J]. 光子学报, 2015, 44(11): 1127002.
[13] 徐泽华, 李伟, 许强, 等. 基于分光链路模拟的曲面量子雷达散射截面研究[J]. 空军工程大学学报(自然科学版), 2019, 20(1): 90-95.
Xu Z H, Li W, Xu Q, et al. Research on scattering section of surface quantum radar based on simulation of spectral link[J]. Journal of Air Force Engineering University (Natural Science Edition), 2019, 20(1): 90-95.
[14] 王书, 任益充, 饶瑞中, 等. 大气闪烁对相干态量子干涉雷达探测性能的影响[J]. 中国激光, 2018, 45(8): 0810002.
[15] 王书, 任益充, 饶瑞中, 等. 大气损耗对量子干涉雷达的影响机理[J]. 物理学报, 2017, 66(15): 150301.
Wang S, Ren Y C, Rao R Z, et al. Influence of atmosphere attenuation on quantum interferometric radar[J]. Acta Physica Sinica, 2017, 66(15): 150301.
[16] Yang P, Bi L, Baum B A, et al. Spectrally consistent scattering, absorption, and polarization properties of atmospheric ice crystals at wavelengths from 0.2 to 100 μm[J]. Journal of the Atmospheric Sciences, 2013, 70(1): 330-347.
[17] 张琳. 卷云的辐射传输与散射特性研究[D]. 西安: 西安电子科技大学, 2010: 23- 24.
ZhangL. Study on the radiative transmission and scattering properties of cirrus clouds[D]. Xi'an:Xidian University, 2010: 23- 24.
[18] . 107(D13): AAC 7[J]. Key J R. Parameterization of shortwave ice cloud optical properties for various particle habits. Journal of Geophysical Research: Atmospheres, 2002.
[19] 李铁飞, 李伟, 杨峰, 等. 卫星量子通信的光子偏振误差影响与补偿研究[J]. 量子电子学报, 2015, 32(6): 678-685.
[21] Strohbehn J, Clifford S. Polarization and angle-of-arrival fluctuations for a plane wave propagated through a turbulent medium[J]. IEEE Transactions on Antennas and Propagation, 1967, 15(3): 416-421.
[22] 唐守荣, 聂敏, 杨光, 等. 空间尘埃等离子体对量子卫星通信性能的影响[J]. 光子学报, 2017, 46(12): 1206002.
[23] 李凤娇, 邹雪峰, 崔亮, 等. 便携式偏振纠缠量子光源[J]. 激光与光电子学进展, 2019, 56(9): 092701.
[24] 史鹏, 赵士成, 李文东, 等. 水下量子密钥分配的误码率和成码率[J]. 中国海洋大学学报, 2017, 47(4): 114-120.
Shi P, Zhao S C, Li W D, et al. Bit error rate and key generation rate for underwater quantum key distribution[J]. Periodical of Ocean University of China, 2017, 47(4): 114-120.
[26] Hufford G. Amodel for the complex permittivity of ice at frequencies below 1 THz[J]. International Journal of Infrared and Millimeter Waves, 1991, 12(7): 677-682.
Article Outline
聂敏, 张怡心, 杨光, 张美玲, 孙爱晶, 裴昌幸. 冰晶粒子对量子干涉雷达探测性能的影响[J]. 激光与光电子学进展, 2020, 57(13): 130101. Min Nie, Yixin Zhang, Guang Yang, Meiling Zhang, Aijing Sun, Changxing Pei. Influences of Ice Crystal Particles on the Detection Performance of Quantum Interference Radar[J]. Laser & Optoelectronics Progress, 2020, 57(13): 130101.