基于粒子群优化算法的宽带混合光纤放大器的设计

巩 译，刘 芳，孟繁轲

（1.北京信息科技大学 信息与通信工程学院，北京 100867；
2.西安邮电大学 通信与信息工程学院，陕西 西安 710121）

波分复用技术（wavelength division multiplexing，WDM）作为有效提升传输系统中继距离与可用带宽的技术手段，随着高速扩张的信息量和日益庞大的网络用户总数逐渐被通信行业广泛投入使用[1]。光放大器为WDM 传输系统中关键性器件[2]，其性能的改善与开发决定着WDM 和现代光波通信的发展速度。WDM 系统中应用最广泛的光放大器是掺铒光纤放大器（erbium doped fiber amplifier，EDFA）和拉曼光纤放大器（Raman fiber amplifier，RFA）。掺铒光纤放大器广泛用于C+L 波段放大，提供高增益系数和成本效益，但工作带宽固定、带宽区间较小[3-6]，增加了自发辐射（ASE）噪声。RFA 可以放大任意波长的信号[7-11]，且增益波段灵活、噪声指数（NF）低，但增加了系统的非线性效应，降低了泵浦光到信号光的功率传输效率[12]。

2 个放大器EDFA 和RFA 各有优缺点，级联2 个放大器可以增加两者的优势。使用混合光放大器（hybrid optical amplifier，HOA）能够降低系统的非线性和噪声系数，提高增益带宽和系统容量[13]。但混合放大器的性能取决于许多参数，包括掺铒光纤放大器和拉曼放大器的光纤长度、泵浦功率及波长等。因此需要优化HOA 的主要参数，以使混合系统的增益最大化。

本文以宽带混合放大器的设计与优化为研究目标，采用能有效改善浓度淬灭、离子聚集效应的铒/镱共掺光纤放大器[14]（erbium-ytterbium doped fiber amplifier，EYDFA）取代EDFA。对EYDFA/RFA混合放大器展开研究，给出了EYDFA/RFA 混合放大器理论模型，并对所设计的混合放大器结构进行合理的参数配置。为了进一步提升混合放大器增益，改善其增益平坦度，建立了二阶多泵浦分布式EYDFA/RFA 混合放大器结构。以增益平坦度和增益值为优化目标设计了二阶多泵浦EYDFA/RFA 混合放大器模型下的粒子群优化算法。

1.1 理论模型

混合放大器的数学模型由RFA 的Raman 耦合波方程与EYDFA 传输方程联合构成，稳态EYDFA的传输方程为

式中：δp为泵浦光散射损耗系数；
δs为信号光散射损耗系数；
NEr*(z)与NYb*(z)分别表示在z 处Er3+、Yb3+的粒子数浓度；
σe(λp)与σe(λs)分别表示泵浦光与信号光辐射过程中波长相关的截面面积；
σα(λp)与σα(λs)分别表示泵浦光与信号光吸收过程中波长相关的截面面积。

结合RFA 功率耦合波方程组，忽略自发Raman散射过程与Rayleigh 散射过程，将研究重心放在信号光与泵浦光以及它自身间的作用过程[15-16]，将完整考虑了RFA 内部各类物理效应的Raman 耦合波微分方程[17-18]简化为如下形式：

式中：Pk、Pi、Pj分别表示第k、i、j路信道中的光信号所对应的功率；
vi、vj分别为信号光和泵浦光的初始输入频率；
Δv为不同信道的频移差；
α(vi)表示光信号频率在vi时所对应的光纤损耗系数；
Keff为偏振相关因子；
Aeff为光纤有效横截面积；
γ(vi)为瑞利散射系数；
gR(vi-vj)为i、j两路信号光之间的拉曼增益系数；
gR(vj-vk)为j、k两路信号光之间的拉曼增益系数；
k和h分别表示波尔兹曼常量和普朗克常量；
[exp(h(vi-vj)/KT)-1]-1为玻色-爱因斯坦因子；
T为光纤绝对温度。

上述理论模型以WDM 系统为基础，EYDFA/RFA 混合放大器的具体求解过程如下：将信号光光功率与EYDFA 的泵浦光功率作为已知的初值代入EYDFA 的传输方程中，构成已知初值的非线性微分方程组，使用四阶龙格-库塔进行求解，得到光信号通过EYDFA 部分后的输出光功率。将所得的信号光输出光功率作为RFA 的输入信号光功率，联合RFA 的泵浦光功率构成第2 个已知初值的非线性微分方程组，继续使用四阶龙格-库塔，所得结果即为混合放大器的输出光功率。

1.2 混合放大器的结构设计

本文采用级联结构的前向EYDFA/RFA 混合放大器模型。

图1 为所设计的二阶EYDFA/RFA 混合放大器的结构示意图。采用前向级联结构的EYDFA/RFA 混合放大器基本结构[19]，以81 路信道间隔为1.125 nm 的信号光组合注入混合放大器，EYDFA部分采用单泵浦输入，RFA 部分使用二阶多泵浦模型。81 路信号光先与EYDFA 泵浦光耦合进入EYDFA 放大，放大完成后经矩形滤波器与RFA 的泵浦光共同耦合进入石英光纤中放大。光谱仪1、2 可用于监测EYDFA 放大前后的信号光谱，光谱仪3、4 则可监测RFA 放大前后的信号光谱。

图1 二阶EYDFA/RFA 混合放大器结构图Fig.1 Structure diagram of second-order EYDFA/RFA hybrid amplifier

粒子群优化算法（particle swarm optimization，PSO）起源于鸟群觅食过程中自发的群体搜索行为。PSO 将个体鸟儿抽象成“粒子”，鸟群寻找的食物则是“目标最优解”，作为一种搜索目标范围内最优解的智能算法[20-23]。使用粒子群优化算法对二阶EYDFA/RFA 混合放大器优化过程如下。

1）粒子群初始化，设置学习因子c1、c2，随机数r1、r2，使种群中的全部粒子得到初始位置xi与初始速度vi。为个体粒子最佳位置Pbest与粒子群最佳位置gbest设置初始值；

2）计算当前模型下各粒子的适应度值fiti；

3）比较此次迭代过程中的粒子适应度值与已获得的个体粒子最佳位置Pbest，若当前适应度值优于已获得的个体极值，以此适应度值替换掉个体极值；
若当前适应度值小于已获得的个体极值，个体极值维持原值；

4）比较此次迭代过程中的粒子适应度值与全局最佳位置gbest，若当前适应度值优于已有的全局极值，以此适应度值替换掉全局极值；
若当前适应度值小于已有的全局极值，全局极值维持原值；

5）将步骤3、4 中所得极值代入（9）式、（10）式，更新当前所得最优解：

6）限制粒子变化时的边界条件，令粒子群在正确的区域内完成迭代，当粒子位置或速度超出此时模型的范围时，将适应度归0；

7）对比此时结果是否满足终止条件，满足条件时，结束循环迭代过程输出所获最优解；
无法满足条件时，返回步骤2 重新进入循环，按照步骤继续依次向下进行迭代。

在进行优化之前，需要对粒子群算法设置一些基本的优化参数，如表1所示。由于采用了5 个泵浦，而每一个泵浦具备泵浦波长和泵浦功率2 个因素，所以总共有5 个泵浦波长和5 个泵浦功率，故一个粒子的维度为10。

表1 粒子群算法参数Table 1 Parameters of particle swarm algorithm

表2 为本次仿真的系统参数，分别确定了泵浦光和信号光的波长、功率的取值范围，以及光纤长度、光纤损耗系数等基本参数。基于表1 和表2所设置的基本参数，按照粒子群优化算法的流程进行优化。在优化过程中，每一次迭代的平均增益和适应度值都可得到，如果增益和增益平坦度满足循环终止条件，则结束算法，获得优化结果；
否则，会对粒子的位置和速度进行更新，而更新后的值将会用于下一次的迭代过程。

表2 二阶EYDFA/RFA 混合放大器仿真参数Table 2 Simulation parameters of second-order EYDFA/RFA hybrid amplifier

完成多组粒子群优化过程后，最终得到了一个拥有90 nm 带宽、平均增益达到38.78 dB 的宽带高增益混合放大器，并且其增益平坦度仅为1.1 dB，所对应的最优泵浦功率与波长的组合方式如表3所示。按照这种组合方式，得到了二阶EYDFA/RFA混合放大器增益谱，如图2所示。

表3 优化后的混合放大器泵浦配置方式Table 3 Optimized pump configuration of hybrid amplifier

图2 EYDFA/RFA 混合放大器增益谱Fig.2 Gain spectrum of EYDFA/RFA hybrid amplifier

利用前向级联结构的EYDFA/RFA 混合放大器基本结构，信号光与EYDFA 泵浦光耦合进入EYDFA 放大，放大完成后信号光与RFA 的泵浦光共同耦合进入石英光纤中进行二次放大，得到了高增益、低平坦度的拉曼增益谱。由图2 可以看出，在90 nm 带宽内，最大增益为39.29 dB，最小增益为38.18 dB，增益曲线在38.78 dB 上下波动，但波动幅度不超过1.1 dB，整体上非常稳定。

为了分析5 个泵浦光在传输过程中的相互作用，研究了泵浦功率随光纤长度的变化情况。从图3 可以看出，RFA 二阶泵浦光刚进入光纤初始端其功率值就开始急速直线下降，当减小到300 mW时开始缓慢变小，随后几乎为0。而其他4 路一阶泵浦光则出现先增大后缓慢减小的趋势。这是因为RFA 二阶泵浦光将能量传递给了4 个一阶泵浦光，这4 路泵浦光在传输过程中从二阶泵浦光处吸收能量，随后因光纤损耗和信号光对能量的吸收又逐渐减小。在整个传输过程中，4 个泵浦光互相传递能量，将二阶泵浦光的作用充分展现，是较优的泵浦光配置。

由图3 可知，RFA 二阶泵浦光波长为1 306.21 nm，初始功率为2 421 mW，为整个RFA 结构提供能量，所以其对EYDFA/RFA 混合放大器的影响也很大。运用控制变量法，即保持系统其他参数不变，不断改变RFA 二阶泵浦光功率值的大小，如图4所示。当RFA 二阶泵浦光功率值从0 开始逐渐增大时，混合放大器的平均增益也在不断增大，当增大到2.4 W 左右时到达峰值，此后随功率值的增大而开始减小，这是因为泵浦功率为2.4 W 时，系统增益达到了饱和，而后当再增加泵浦功率时增益便不再增加，因此应该合理设置泵浦功率值，使混合放大器的性能达到最优。

图3 RFA 泵浦光功率沿光纤长度的变化Fig.3 Variation of RFA pump power along fiber length

图4 RFA 二阶泵浦功率对混合放大器平均增益的影响Fig.4 Effect of RFA second-order pump power on average gain of hybrid amplifier

图5(a)为ASE 噪声和双瑞利散射（DRBS）噪声的信道SNR随光纤长度的变化情况。从图中我们可以看出，光纤输入端系统ASE 噪声和DRBS噪声的信噪比分别为46 dB、38 dB，SNR 随光纤长度的增加不断减小，其中ASE 噪声的信噪比减小比较缓慢，在10 km 之后信噪比大致稳定在32 dB，而DRBS 噪声的信噪比在0～30 km 之间一直在减小。在图5(b)中，输入端系统ASE 噪声和DRBS噪声的噪声系数为1 dB，随着光纤长度的增加，ASE 噪声和DRBS 噪声的噪声系数都成增长趋势。ASE 噪声系数在5 km 之后停止增长，并保持不变，而DRBS 噪声系数随着光纤长度增大而不断增大。结合图5(a)、5(b)我们可知，在输入端时ASE 噪声和DRBS 噪声较小，噪声性能良好，但随着光纤长度的增加，信噪比快速下降、噪声系数快速增加，严重影响了系统的性能。DRBS 噪声与ASE 噪声相比，信噪比下降更快，噪声系数增加更为明显，因此，DRBS 噪声对系统影响更大，在优化设计放大器时应考虑噪声对RFA 性能的影响。

图5 EYDFA/RFA 混合放大器各信道SNR 及增益NF 随光纤长度的变化Fig.5 Variation of SNR and NF of each channel of EYDFA/RFA hybrid amplifier with fiber length

根据以上仿真结果，说明了在1 530 nm～1 620 nm的放大带宽内，当放大器的系统增益达到饱和之后，即使再增加泵浦功率，增益也不再增加。同时，由于噪声的存在使得对系统产生较大影响，因此，在优化设计放大器时应使用优化算法合理设置泵浦功率值，找到达到饱和增益时的最大泵浦功率，使混合放大器的性能达到最优，并考虑噪声对RFA 性能的影响。

建立了二阶多泵浦分布式EYDFA/RFA 模型，为了优化二阶EYDFA/RFA 混合放大器的增益以及平坦度，设计了二阶混合放大器粒子群优化算法流程，并制定了相应的适应度函数，经过粒子群优化过程后，得到了一个带宽90 nm、增益38.78 dB的宽带高增益混合放大器，增益平坦度仅为1.1 dB，工作区间覆盖C+L 波段。同时也分析了二阶泵浦对其增益的影响，以及当光纤长度增加时系统ASE 和DRBS 噪声的变化情况，为混合放大器的设计和优化提供了参考。随着机载、车载以及舰载的各类光电成像系统的广泛应用，该装置可以提高系统检测效率，为光电成像系统设计、装配和性能提升提供计量保障，具有广泛的使用价值。