摘要:量子通信具有覆盖面广、安全保密的优势, 是当前通信领域国内外的研究热点. 在自由空间量子通信过程中, 光量子信号需要在地表上空一定高度进行传输, 因此各种环境因素, 例如降雪、沙尘暴、降雨、雾霾、浮尘等, 不可避免地会影响量子通信性能. 然而, 迄今为止, 降雪对地表附近自由空间量子信道影响的研究尚未展开. 为此, 根据降雪的强度, 将降雪分为小雪(${S_{\rm{1}}}$)、中雪(${S_{\rm{2}}}$)、大雪(${S_{\rm{3}}}$) 和暴雪(${S_{\rm{4}}}$) 四个等级. 由于空中正在飘落雪花对光量子信号具有能量吸收作用, 称为消光效应, 不同强度的降雪, 其消光效应对自由空间光量子信号的影响不同. 本文首先建立了不同等级降雪对光量子信号消光效应的数学模型; 然后建立了因降雪导致的自由空间消光衰减定量关系, 信道极限生存函数、不同降雪强度下的信道容量和量子误码率等性能参数受降雪影响的变化情况; 最后建立了降雪强度、传输距离与链路衰减、幅值阻尼信道容量、信道生存函数以及信道误码率的数学模型. 仿真结果表明, 当降雪强度为2.1 mm/d (${S_{\rm{1}}}$), 传输距离为2.2 km时, 通信链路衰减为0.0362, 信道容量为0.7745, 信道生存函数为0.2329, 信道误码率为0.0105. 当降雪强度为3.8 mm/d (${S_{\rm{2}}}$), 传输距离为3.5 km时, 通信链路衰减为0.1326, 信道容量为0.4922, 信道生存函数为0.2099, 信道误码率为0.019. 由此可见, 降雪对量子通信性能有不同程度的影响. 所以在实际应用中应根据降雪强度大小, 自适应调节量子通信相关参数, 提高量子通信的可靠性.
关键词: 量子通信/
降雪强度/
消光系数/
生存函数

English Abstract

全文HTML

--> --> -->

量子通信目前是国内外的研究前沿和热点, 在国防、金融等领域有着重要的作用, 并具有显著的优点, 例如量子密钥分发 (quantum key distribution, QKD) 协议具有较强的安全性^[1]. 2016年, Smania等^[2]提出了三层量子通信系统的通信协议, 并且实现了三方量子通信的信息加密共享. 2020年, B?uml等^[3]通过线性程序实现了量子互联网上的量子纠缠与密钥分发, 使得量子互联网的实现更进一步. 2020年Bhaskar等^[4]使用17, 18, 19纳米光子的金刚石谐振器, 实现了异步光子钟态测量, 有效地使量子通信以超过理想的等效损耗直接传输. 2010年, 潘建伟团队^[5]实现了16 km的量子隐形传态实验. 2016年, “墨子号”成功发射, 全球瞩目^[6]. 2017年, 潘建伟团队^[7]完成了53 km的自由空间密钥分发. 2020年, 潘建伟团队等^[8]利用“墨子号”量子科学实验卫星, 通过基于纠缠的QKD协议, 在国际上首次实现无中继千公里级的量子保密通信, 大幅度提升了量子通信距离.
然而, 量子通信在进行传输时, 不可避免地会受到自由空间各种大气环境因素的影响, 如图1所示. 文献[9]研究了三种不同形状的气溶胶粒子、不同大气相对湿度背景下的量子通信性能. 文献[10]研究了地球中纬度地区背景下的偶发E层对量子通信的影响. 文献[11]研究了中尺度沙尘暴对量子卫星通信的影响. 文献[12]根据降雨的分布模型, 提出了变色龙算法, 通过调整每脉冲光子数, 提高量子通信质量. 文献[13]根据自然环境下对量子通信的干扰因素, 提出了基于袋鼠纠缠跳跃模型的自适应通信策略, 提高了量子通信在雾霾沙尘暴等干扰下的性能. 这些都是不同环境下, 量子通信受到自然环境因素的影响以及所采取的抗干扰策略, 对未来的研究有着重要的参考意义.
图 1 自由空间量子通信简图
Figure1. Schematic diagram of free space quantum channel.

由于降雪会产生消光作用^[14], 对量子通信质量有一定的影响. 因此, 本文根据降雪的谱分布函数以及雪对光量子的吸收和散射所引起的消光特性, 研究了降雪消光系数和降雪强度之间的关系, 进一步研究了光量子能量与降雪强度之间的关系, 建立了降雪强度与幅值阻尼信道容量、信道生存函数和信道误码率之间的关系, 通过仿真研究, 为降雪条件下地表附近自由空间量子通信提供参考.

2

-->

2.1.理论模型

根据24 h内的降雪量划分标准^[15], 可将降雪强度按等级划分, 划分结果如表1所列.

降雪量/mm	降雪强度等级
< 2.5	小雪(${S_{\rm{1}}}$)
2.5—5	中雪(${S_{\rm{2}}}$)
5—10	大雪(${S_{\rm{3}}}$)
> 10	暴雪(${S_{\rm{4}}}$)

表1降雪强度划分标准
Table1.Criteria for classification of snowfall intensity.

由于雪花一般都是不规则形状, 因此常用Γ分布描述降雪谱分布函数^[16], 可表示为

$N(D) = {N_0}{D^m}{{\rm{e}}^{ - {\varLambda} D}},$

式中, $D$表示雪花的等效直径; ${N_{\rm{0}}}$表示降雪浓度参数; $m$表示雪的形状因子, 降雪一般服从指数分布, 所以$m$可取0; ${\varLambda} $表示降雪谱分布斜率, 与降雪强度$S$的关系为^[17]

${\varLambda} = 25.5{S^{{\rm{ - }}0.48}}.$

已知降雪谱分布函数为$N\left( D \right)$, 由Mie散射可得雪的消光系数为^[18]

$\sigma = \frac{\pi }{4}\int_{\rm{0}}^\infty {{D^{\rm{2}}} \cdot {Q_{{\rm{ext}}}}(D,\lambda ,m) \cdot N(D)\operatorname{d} D} ,$

式中, ${Q_{{\rm{ext}}}}$表示雪的消光率因子, 由于雪的最小直径大于10 μm, 远大于可见光的波长, 所以可认为${Q_{\rm {ext}}} \approx 2$^[19]; $\sigma $的单位为dB/km.
2 -->

2.2.仿真分析

根据 (1)—(3) 式, 对降雪强度$S$和消光系数$\sigma $进行仿真分析, 结果如图2所示. 图中横坐标表示降雪强度$S$, 纵坐标表示消光系数$\sigma $. 根据图2可知, 随着降雪量增加, 消光系数逐渐上升, 当降雪强度为2 mm/d (${S_{\rm{1}}}$) 时, 消光系数为4.11 dB/km. 据此可知, 降雪会影响光量子信号衰减.
图 2 降雪强度和消光系数之间的关系
Figure2. Relationship between extinction coefficient and snowfall intensity.

2

-->

3.1.降雪对通信链路的影响分析

由Bougure定律^[20]可知, 消光系数$\sigma $和光量子能量$E$的关系式可表示为

$E = {E_0}\exp ( - \sigma \cdot L),$

式中, ${E_0}$表示光量子初始能量, $L$表示光量子传输距离. 取对数可得链路衰减因子为

${L_{\rm{att}}} = 10 \cdot \frac{{\pi {N_0}}}{{{{\left( {25.5{S^{ - 0.48}}} \right)}^3}}} \cdot \lg {\rm{e}} \cdot L . $

对链路衰减${L_{{\rm{att}}}}$、降雪强度S和光量子传输距离$L$之间的关系进行仿真, 结果如图3所示. 由图3可以看出: 当降雪强度S为0时, 此时为理想状态, 通信链路几乎无损耗; 随着降雪强度增加, 通信链路急剧衰减; 当降雪强度为4.1 mm/d (${S_2}$), 传输距离3.7 km时, 通信链路为0.1494. 据此可知, 降雪会对通信链路传输产生较大影响.
图 3 链路衰减与降雪强度、传输距离的关系
Figure3. Photon energy versus snowfall intensity and transmission distance.

2 -->

3.2.降雪对信道容量的影响分析

在量子通信传输时, 降雪下的环境量子态和光量子态的复合系统经受联合幺正演化, 从而导致了消相干现象. 本文研究幅值阻尼信道容量变化.
幅值阻尼信道运算算子可以表示为^[21]

${T_k} = \sum\limits_n {\sqrt {\left( \!{\begin{array}{*{20}{c}} n \\ k \end{array}}\! \right)} \sqrt {{{\left( {1 - \eta } \right)}^{n - k}}{\eta ^k}} \left| {n - k} \right\rangle \left\langle n \right|} ,$

式中, $\left| k \right\rangle , \left\langle n \right|$分别表示降雪环境与量子湮灭算子的本征态; $\eta $为受到降雪丢失一个光量子的概率, 可表示为

$\eta = \frac{{{E_0} - E}}{{{E_0}}} = 1 - \exp \left( { - \sigma \cdot L} \right). $

单光子比特状态下, 有

$ \mathit{\boldsymbol{\rho}} = \left(\! {\begin{array}{*{20}{c}} a&b \\ {{b^*}}&c \end{array}}\! \right),$

通过幅值阻尼信道后变为

$\mathit{{\varepsilon}} \left( \mathit{{\rho}} \right) = \left[ {\begin{array}{*{20}{c}} {1 - \left( {1 - \eta } \right)\left( {1 - a} \right)}&{b\sqrt {1 - \eta } } \\ {b*\sqrt {1 - \eta } }&{c\sqrt {1 - \eta } } \end{array}} \right]. $

若输入字符${\rho _1} = \left| 0 \right\rangle \left\langle 0 \right|, {\rho _2} = \left| 1 \right\rangle \left\langle 1 \right|$, 有

$\begin{split} &{ {\varepsilon}} \left( {\sum\limits_i {{p_i}{\rho _i}} } \right) = { {\varepsilon}} \left[ {{p_1}{\rho _1} + \left( {1 - {p_1}{\rho _2}} \right)} \right] \\ =\;& {p_1}\left[ {\begin{array}{*{20}{c}} 1&0 \\ 0&0 \end{array}} \right] + \left( {1 - {p_1}} \right)\left[ {\begin{array}{*{20}{c}} \eta &0 \\ 0&{1 - \eta } \end{array}} \right] \\ =\; & \left[\! {\begin{array}{*{20}{c}} {{p_1} + \left( {1 - {p_1}} \right)\eta }&0 \\ 0&{\left( {1 - {p_1}} \right)\left( {1 - \eta } \right)} \end{array}}\! \right], \end{split} $

其中${P_i}$表示信源量子符为${\rho _i}$时的概率. (10) 式对应的冯诺依曼熵为^[22]

$\begin{split} & s\left[ {{ \varepsilon} \left( {\sum\limits_i {{p_i}{\rho _i}} } \right)} \right] \\ =& - \big\{ \left[ {{p_1} + \left( {1 - {p_1}} \right)\eta } \right]\log \left[ {{p_1} + \left( {1 - {p_1}} \right)\eta } \right] \\ &{\rm{ }} + \left( {1 - {p_1}} \right)\left( {1 - \eta } \right)\log \left( {1 - {p_1}} \right)\left( {1 - \eta } \right)\big\} . \end{split} $

可得幅值阻尼信道的信道容量为

$\begin{split} C &= - \Big\{ {p_1} + \left( {1 - {p_1}} \right)\eta \log \left[ {{p_1} + \left( {1 - {p_1}} \right)\eta } \right] \\ &{\rm{\;\;\;\;}} + \left( {1 - {p_1}} \right)\left( {1 - \eta } \right)\log \left( {1 - {p_1}} \right)\left( {1 - \eta } \right) \Big\} \\ &{\rm{\;\;\;\;}} - \left( {1 - {p_1}} \right){H_2}\left( \eta \right), \\ \end{split} $

式中, ${H_2}\left( \eta \right)$为二元香农熵, 且${p_1} = \dfrac{{t\left( {1 - \eta } \right) - \eta }}{{\left( {1 + t} \right)\left( {1 - \eta } \right)}}, {\rm{ }}t = {2^{\frac{{{H_2}\left( \eta \right)}}{{1 - \eta }}}}$.
对信道容量、传输距离以及降雪强度进行仿真, 如图4所示. 由图4可知: 当降雪量$S$以及传输距离$L$增加时, 信道容量$C$逐渐减小; 信道容量为1时为理想信道, 且降雪量和传输距离均为0; 当降雪强度为2.2 mm/d, 传输距离为7.2 km时, 信道容量衰减到0.51. 由此可知, 降雪对量子通信的信道容量影响较大, 因此需要根据降雪强度的大小, 自适应调节量子通信的各项参数, 以满足人们对量子通信的需求.
图 4 信道容量与降雪强度、传输距离的关系
Figure4. Channel capacity versus snowfall intensity and transmission distance.

2 -->

3.3.降雪对量子通信系统生存函数的影响分析

生存函数^[23]可表示为

$B = \xi \cdot F,$

式中, $\xi $表示降雪背景下的信道生存系数, $F$表示降雪背景下信道的传输保真度. 生存系数可定义为

$\xi = \frac{{5809}}{{5810}} - \frac{{9\sigma }}{{2905}},$

保真度$F$可表示为^[24]

$F = F\left( {{ \rho} ,{ \rho} '} \right) = {\left( {{\rm tr}\sqrt {{{ \rho} ^{1/2}}{ \rho} '{{ \rho} ^{1/2}}} } \right)^2},$

那么降雪情况下的信道生存函数可定义为

$\begin{split} B =\;& \xi \cdot F = \xi \cdot F\left[ {\sum\limits_i {{p_i}{\rho _i}} ,\varepsilon \left( {\sum\limits_i {{p_i}{\rho _i}} } \right)} \right] \\ & = \xi \cdot {\rm{tr}}\Bigg\{ \Bigg[ {{\left( {\sum\limits_i {{p_i}{\rho _i}} } \right)}^{{\rm{1/2}}}} \cdot \varepsilon \left( {\sum\limits_i {{p_i}{\rho _i}} } \right) \\ & \qquad \times {{\left( {\sum\limits_i {{p_i}{\rho _i}} } \right)}^{{\rm{1/2}}}} \Bigg]^{{\rm{1/2}}} \Bigg\}. \\[-15pt] \end{split} $

根据幅值阻尼信道的性质, 可得

${ \varepsilon} \left( {\sum\limits_i {{p_i}{\rho _i}} } \right) = \left[ {\begin{array}{*{20}{c}} {{p_1} + \left( {1 - {p_1}} \right)}&0 \\ 0&{\left( {1 - {p_1}} \right)\left( {1 - \eta } \right)} \end{array}} \right],$

所以信道生存函数可表示为

$\begin{split}B &= \xi \cdot F = \xi \cdot {\left\{ {{p_1}\left[ {{p_1} + \left( {1 - {p_1}} \right)\eta } \right]} \right\}^{1/2}} \\&+ \left( {1 + {p_1}} \right)\sqrt {1 - \eta } ,\end{split}$

其中p₁表示信源取${\rho _1}$的概率, 可进一步写为

$\begin{split} &B = \xi \cdot \bigg[ {p_1}\bigg( {p_1} + \left( {1 - {p_1}} \right) \\ & \times \bigg\{ {1 - \exp \bigg[ { - \frac{{{\rm{π}} {N_0}}}{{{{( {25.5{S^{ - 0.48}}} )}^3}}} \cdot L} \bigg]} \bigg\} \bigg) \bigg]^{1/2} \\ & + (1 \!+\! p_1)\sqrt{1 - \bigg\{\! {1 - \exp \bigg[ { - \frac{{{\rm{π}} {N_0}}}{{{{( {25.5{S^{ - 0.48}}})}^3}}} \cdot L} \bigg]}\! \bigg\}} . \\ \end{split} $

对降雪强度、传输距离及信道生存函数进行仿真, 如图5所示. 由图5可以看出: 当降雪量以及传输距离增加时, 生存函数逐渐下降; 当降雪强度为4.4 mm/d (${S_2}$), 传输距离为7.7 km时, 生存函数为0.2004. 由此可知, 降雪对量子通信的信道生存函数有着较为明显的影响, 因此需要根据降雪强度的大小, 调整量子通信相关参数, 保证量子通信质量.
图 5 生存函数与降雪强度、传输距离的关系
Figure5. Channel survival ability versus snowfall intensity and transmission distance.

2 -->

3.4.降雪对量子通信系统信道误码率的影响分析

基于BB84协议的量子密钥分发系统, 定义量子误码率${P_{\rm{B}}}$为接收到的误码比特数${P_{\rm{e}}}$与总比特率${P_{\rm{s}}}$的比值^[25], 表示为

${P_{\rm{B}}} = {{{P_{\rm{e}}}}}/{{{P_{\rm{s}}}}}. $

令${P_{{\rm{snow}}}}$为由降雪导致的误码量子比特数, (20)式变成

${P_{\rm{B}}} = {{{P_{{\rm{snow}}}}}}/{{{P_{\rm{s}}}}}.$

由文献[26]得

${P'_{\rm{B}}} = \frac{1}{2} \frac{{S\tau n}}{{1 - \exp \left( { - \mu {T_{\rm{c}}}{P_{\rm{a}}}{T_{\rm{a}}}{\eta _{\rm{d}}}{F_{\rm{m}}}} \right)}},$

式中${T_{\rm c}}$为信道传输率, 且有

${T_{\rm{c}}} = \exp \left( { - \sigma \cdot L} \right);$

其他参数所表示的意义和取值如表2所列.

参数	含义	取值
$\tau $	量子探测器时间窗口	1
$n$	探测器数目	1
$\mu $	平均光子数	1
${P_{\rm a}}$	单光子捕获率	0.5
${T_{\rm a}}$	系统装置传输率	1
${\eta _{\rm d}}$	单光子探测器效率	0.65
${F_{\rm m}}$	测量因子	1

表2信道误码率各参数含义和取值
Table2.Meaning and values of the parameters of the channel bit error rate.

根据表2取值, 可得

$P_{\rm B}^\prime \!=\! \frac{1}{2} \frac{S}{{1 \!-\! \exp \bigg\{\!- 0.325 \exp \bigg[ - \dfrac{\text{π} {N_0}}{( 25.5 S^{-0.48})^3}\cdot L\bigg]\bigg\}}}.$

对误码率${P_{\rm{B}}}$, 传输距离$L$和降雪强度$S$进行仿真分析, 结果见图6. 当降雪量和传输距离增加, 信道误码率急剧上升. 当降雪强度为5.5 mm/d, 传输距离为3.7 km时, 信道误码率达到了0.0275. 据此可得, 降雪会影响量子通信误码率.
图 6 误码率与降雪强度、传输距离的关系
Figure6. Quantum bit error rate versus snowfall intensity and transmission distance.

本文对降雪背景下地表附近自由空间量子信道性能进行了分析. 根据降雪的谱分布函数以及降雪的消光特性, 研究并分析了降雪与链路衰减、幅值阻尼信道容量、生存函数以及信道误码率之间的关系. 由仿真分析可得, 降雪会影响自由空间量子通信性能. 因此在进行量子通信时, 降雪的环境因素必须考虑, 需要根据降雪强度, 适当调整相关参数, 以减小降雪对量子通信的影响.