高功率掺Tm³⁺光纤放大器热效应管理的泵浦方式优化理论研究

高功率掺Tm 3+ 光纤放大器热效应管理的泵浦方式优化理论研究 Theoretical Study of Pumping Method of High Power Tm 3+ -doped Fiber Amplifier for Thermal Effect Management 1 first-author 张轲(1995-), 男, 山西运城人, 硕士研究生, 2017年于中国石油大学(华东)获得学士学位, 主要从事2 μm波段光纤激光器方面的研究。E-mail:17120162@bjtu.edu.cn 张轲 (1995-), 男, 山西运城人, 硕士研究生, 2017年于中国石油大学(华东)获得学士学位, 主要从事2 μm波段光纤激光器方面的研究。E-mail: 17120162@bjtu.edu.cn 17120162@bjtu.edu.cn 1 corresp YAN Feng-ping, E-mail:fpyan@bjtu.edu.cn YAN Feng-ping, E-mail:fpyan@bjtu.edu.cn 延凤平(1966-), 男, 山西兴县人, 博士, 教授, 1996年于北方交通大学获得博士学位, 主要从事光纤激光器、光纤传感器、光纤通信、基于太赫兹超材料等方面的研究。E-mail:fpyan@bjtu.edu.cn 延凤平 (1966-), 男, 山西兴县人, 博士, 教授, 1996年于北方交通大学获得博士学位, 主要从事光纤激光器、光纤传感器、光纤通信、基于太赫兹超材料等方面的研究。E-mail: fpyan@bjtu.edu.cn fpyan@bjtu.edu.cn 1 2 corresp FENG Ting, E-mail:wlxyft@hbu.edu.cn FENG Ting, E-mail:wlxyft@hbu.edu.cn 冯亭(1986-), 男, 河北张家口人, 博士, 副教授, 2014年于北京交通大学获得博士学位, 主要从事光纤激光、光纤传感器及其应用等方面的研究。E-mail:wlxyft@hbu.edu.cn 冯亭 (1986-), 男, 河北张家口人, 博士, 副教授, 2014年于北京交通大学获得博士学位, 主要从事光纤激光、光纤传感器及其应用等方面的研究。E-mail: wlxyft@hbu.edu.cn wlxyft@hbu.edu.cn 北京交通大学光波技术研究所 全光网络与现代通信网教育部重点实验室, 北京 100044 北京交通大学光波技术研究所 全光网络与现代通信网教育部重点实验室, 北京 100044 Key Laboratory of All Optical Network and Advanced Telecommunication Network of Ministry of Education, Institute of Lightwave Technology, Beijing Jiaotong University, Beijing 100044, China Key Laboratory of All Optical Network and Advanced Telecommunication Network of Ministry of Education, Institute of Lightwave Technology, Beijing Jiaotong University, Beijing 100044, China 河北大学物理科学与技术学院 光信息技术创新中心, 河北 保定 071002 河北大学物理科学与技术学院 光信息技术创新中心, 河北 保定 071002 Photonics Information Innovation Center, College of Physics Science & Technology, Hebei University, Baoding 071002, China Photonics Information Innovation Center, College of Physics Science & Technology, Hebei University, Baoding 071002, China 基于主振荡功率放大器结构的高功率掺Tm 3+ 光纤激光器是2 μm波段高功率光纤激光器的主要实现形式，掺Tm 3+ 光纤放大器（Thulium-doped fiber amplifier，TDFA）热效应管理的研究对于其输出激光功率的不断提升具有重要意义。本文主要对TDFA热效应管理的泵浦方式优化方面进行理论研究，利用龙格库塔法以及牛顿迭代法求解不同泵浦方式下TDFA的稳态速率方程，并根据热传导方程，模拟掺Tm 3+ 光纤（Thulium-doped fiber，TDF）温度沿径向和轴向的分布。结合遗传算法理论，研究了分段泵浦方式，经过参数优化，在功率为5 W的2 020 nm输入信号光、总功率为1 000 W的793 nm激光泵浦、TDF吸收系数为3.1 dB/m条件下，将总长度为11 m的TDF分为2.4，2，2，2，2.6 m的5段进行泵浦，得到放大信号激光输出功率为284.5 W、斜率效率为28.45%、光纤外包层边界最高温度为86.28℃且温度总体分布均匀。与传统前向泵浦、双端泵浦方式下的TDFA相比，其热效应有明显改善。 The high-power thulium(Tm 3+ )doped fiber laser(TDFL) based on the structure of master oscillator power-amplifier is the main type of 2 μm band high-power fiber lasers. Study on thermal effect management of the thulium doped fiber amplifier(TDFA) stage means a lot to the continuously power-scaling of the TDFL. This paper mainly focuses on the theoretical research on the optimization of the pumping method of TDFA for thermal effect management. We use the Runge-Kutta method and Newton iteration method to solve the steady-state rate equation of TDFA under different pumping methods. And, on the basis of the heat conduction equation, the temperature distributions of Tm 3+ doped fiber(Thulium-doped fiber, TDF) along the radial and axial direction are simulated. Combining with the genetic algorithm theory, the segmentally pumping method is studied. Through optimizing of the relative parameters, using a 2 020 nm input signal light with a power of 5 W, a total 793 nm laser pump power of 1 000 W and a TDF with an absorption coefficient of 3.1 dB/m and a total length of 11 m which is divided into 5 segments(2.4, 2, 2, 2, 2.6 m), we obtained an amplified signal laser with an output power of 284.5 W and a slope efficiency of 28.45%, and the maximum temperature at the outer cladding boundary of the TDF is only 86.28 ℃ with a uniform temperature distribution along the whole fiber. Compared with a TDFA using traditional forward pumping or dual-end pumping mode, the thermal effect has been significantly improved. 掺Tm 3+ 光纤放大器 热效应管理 分段泵浦 遗传算法 thulium-doped fiber amplifier thermal effect management segmented pump genetic algorithm 引言 1 2 μm波段掺Tm 3+ 光纤激光器(Thulium-doped fiber laser, TDFL)的输出波长处于人眼安全波段，并且同时覆盖了多个重要的吸收带，如水的吸收峰和大气窗口等，因此有望被广泛应用于诸多领域，如激光医疗手术、激光雷达和激光大气监测与传感等 [ 1 1 - 3 3 1-3 ] 。近年来，光纤激光器的输出功率已从几瓦提高至几千瓦量级，大多采用主振荡功率放大器(Master oscillator power-amplifier, MOPA)结构 [ 4 4 - 5 5 4-5 ] 。除种子源激光器外，高功率放大器的性能是决定MOPA系统输出功率量级的关键，且随着诸多理论与实验工作的进一步开展，研究发现强烈的热效应是限制高功率光纤放大器输出功率进一步提升的主要因素之一。然而，目前关于2 μm波段掺Tm 3+ 光纤放大器(Thulium-doped fiber amplifier, TDFA)热效应管理的理论研究报道还比较少。双包层掺Tm 3+ 光纤(Thulium-doped fiber, TDF)增益介质本身的表面积与体积比很大，散热性能良好。但是，对于上千瓦量级的TDFA，其泵浦光转化为信号光时量子亏损引起强烈的热效应，从而引起光纤热透镜效应、应力和折射率变化等问题 [ 6 6 ] 。此外，热效应还会严重影响激光输出特性，如输出波长不稳定、转换效率下降、光束质量变差、噪声高等。因此，研究高功率TDFA的热效应管理问题对其整体性能的提升有着重要意义 [ 5 5 - 6 6 5-6 ] 。 目前，对于高功率光纤放大器热效应的研究主要集中在三个方面 [ 6 6 ] ：增益光纤制作质量提升、泵浦优化及外部散热辅助，研究也多以1 μm波段的掺镱光纤放大器为主 [ 7 7 - 14 14 7-14 ] 。其中，以泵浦方式及泵浦源选择研究为主的泵浦优化方面可以开展大量的实验与理论研究工作。对于TDFA，常用793 nm的半导体激光器(Laser diode, LD)作为TDF的泵浦源，而由于793 nm~2 μm之间的量子亏损比较大，这种泵浦方式可实现的最高斜率效率为60%，泵浦源产生的热沉积成为限制其输出功率进一步提升的主要原因 [ 15 15 - 16 16 15-16 ] 。研究者们尝试使用其他波段激光泵浦TDF，用于实现更高的输出功率。2007年，Meleshkevich等利用中心波长为1 500 nm附近的铒镱共掺光纤激光器作为泵浦源，获得了中心波长为1 940 nm、功率为415 W的连续激光输出 [ 17 17 ] ；2014年，Creeden等利用中心波长为1 908 nm的TDFL泵浦掺铥光纤获得中心波长为1 993 nm、功率为123 W的连续激光输出 [ 18 18 ] 。但是，特殊波长的高功率泵浦源的获取是非常大的难题，对于2 μm波段TDFA的实现，目前仍以高功率793 nm的LD泵浦为主。可见，在泵浦方式方面开展研究工作，是目前高功率TDFA热效应管理中泵浦优化的主要途径。 本文主要针对高功率TDFA的分段泵浦方式进行理论建模和仿真研究，采用龙格库塔法以及牛顿迭代法求解出耦合速率方程，并根据热传导方程得出TDF的径向和轴向温度分布，再结合遗传算法理论对分段泵浦的每段泵浦功率及长度做出优化选择。以5段泵浦为例，在得到TDFA较高转换效率的同时实现TDF最高温度满足实际需要且整体光纤温度均匀分布。 理论模型与方法 1 泵浦方案和理论模型 2 图 1 图 1 为 3 H 6 → 3 H 4 泵浦方案的能级跃迁示意图。该泵浦方案对应于800 nm的吸收峰，泵浦源通常使用中心波长为793 nm的LD激光器。 在 3 H 6 → 3 H 4 的跃迁方式下得到双包层TDFA的稳态速率方程为 [ 19 19 - 20 20 19-20 ] ： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963525&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963525&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963525&type=middle 其中， τ i 为能级 i 的寿命; N 为总的离子数浓度； G 为交叉弛豫因子； W 03 、 W 10 、 W 01 分别为基态泵浦吸收系数、2 μm激光受激发射系数和受激吸收系数，它们可分别表示为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963526&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963526&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963526&type=middle 其中， λ p 、 λ s 分别为泵浦光和信号光的波长， Γ p 和 Γ s 为泵浦光和信号光的重叠因子， h 为普朗克常量， c 为光速， k ijkl 表示能级 i 到能级 j 和能级 k 到能级 l 的能量转移系数。 σ e ( λ p )、 σ α ( λ p )、 σ e ( λ s )、 σ α ( λ s )分别为泵浦光的发射和吸收截面、2 μm激光的发射和吸收截面； P p ± ( z )和 P s ± ( z )分别为前、后向泵浦光功率和前、后向信号光功率，忽略SBS效应的影响，可分别表示为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963532&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963532&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963532&type=middle 其中， δ p 和 δ s 分别为泵浦光和2 μm激光的本征吸收系数，正、负号分别表示沿光纤的正、反方向。 分段泵浦示意图如 图 2 图 2 所示，其中 P s1 + (0)为种子光功率。公式(4)为满足分 n 段的泵浦边界方程： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963538&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963538&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963538&type=middle 其中， P p m + (0)、 P p m － ( L )、 P s m + (0)、 P s m － (L)分别为第 m 段的前、后向泵浦光与信号光功率。 p m + 和 p m － 分别是第 m 段前、后向注入泵浦光功率， μ 为泵浦点泵浦光的泄露比率，取0.123； η 为前、后向激光在经过泵浦点时的损耗比率，取0.01；下标 m 表示相应的光纤段序数。 温度分布和理论模型 2 TDFA中的大部分热量是吸收泵浦光时量子亏损引起的，其在光纤径向的热传导方程为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963544&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963544&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963544&type=middle 其中 r 为光纤径向方向坐标， k 为热传导系数， Q 为热功率密度。假设TDF的纤芯半径为 r 0 ，内包层半径为 r 1 ，外包层半径为 r 2 。将热传导方程与光纤温度分布的边界条件结合，得到TDF中径向温度分布: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963550&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963550&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963550&type=middle 其中， T 0 ( r )、 T 1 ( r )、 T 2 ( r )分别为纤芯、内包层、外包层中坐标 r 处的温度，且有 T 0 = T 0 (0)。单位体积内泵浦功率密度及纤芯轴向温度可分别表示为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963553&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963553&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963553&type=middle http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963555&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963555&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963555&type=middle 其中， β p 为泵浦吸收系数， H 为对流传热系数， T c 为环境温度或制冷温度，且若只考虑量子亏损产生的热量时， η h =( λ s － λ p )/ λ s 。 理论方法 2 龙格库塔法配合牛顿迭代法可求解两点边值问题的微分方程组，在计算分段泵浦的速率方程时可以将多点边值问题转换为两点问题。基本思路是：假设存在3个泵浦点，分别位于 z =0、 z = kL 、 z = L 处，其中 L 为光纤长度；将泵浦点1和2之间的区间通过线性坐标变换 z = z 1 / k 1 由[0, k 1 L ]换成[0, L ]，泵浦点2和3之间光纤通过线性变化同样由[0, k 2 L ]变为[0, L ]， k 1 、 k 2 分别为两段光纤与整段光纤的长度比，则分别有d z 1 = k 1 L d z 和d z 2 = k 2 L d z 。这样就可将三点边值问题转换为两点边值问题 [ 21 21 ] 。使用MATLAB进行数值仿真，计算中所用到的各个参数的选取如 表 1 表 1 所示。 http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963557&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963557&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963557&type=middle 基于 表 1 表 1 中的参数，对于TDFA， 图 3 图 3 给出了前向泵浦方式下，泵浦功率为1 000 W、信号光功率为5 W、纤芯掺杂浓度 N =2.5×10 26 m －3 时光纤径向温度分布以及光纤外包层边界( r =200 μm)处温度沿光纤轴向的分布情况。从 图 3(a) 图 3(a) 中可以看出，TDF中纤芯与外包层外侧处温差为20 ℃，而纤芯与外包层外侧正常工作所允许的最高温度分别为700 ℃与200 ℃，故外包层外侧的温度为主要的限制因素。从 图 3(b) 图 3(b) 中可以看出，泵浦端的外包层外侧温度为693.21 ℃，远高于200 ℃，所以在高功率TDFA中需要进行热管理。 分段泵浦仿真分析与优化 1 泵浦方式对比及分析 2 本文仿真中参数设置为低掺浓度掺杂Tm 3+ 的交叉弛豫系数，所以得到的斜率效率比较低，但是对应的量子亏损更大，产生的热效应也就更为明显，得到的仿真结果更适用于一般情况，更有助于说明分段泵浦方式热效应管理的优势。当采用前向泵浦方式，泵浦光功率为1 000 W、信号光功率为5 W时，转化斜率效率与交叉弛豫系数关系如 图 4 图 4 所示。从 图 4 图 4 中可以看出，斜率效率随着交叉弛豫系数的增大而增大，这是因为增强交叉弛豫过程会增加量子转换效率。 基于 表 1 表 1 中的仿真参数，对于TDFA， 图 5 图 5 给出了在总泵浦功率为1 000 W条件下，2 020 nm种子源信号光功率为5 W、纤芯掺杂浓度 N =2.5×10 26 m －3 时，分别在前向泵浦、双端泵浦、分段泵浦(5段)3种泵浦方式下，TDF中放大信号光功率及光纤外包层边界( r =200 μm)处温度沿光纤轴向分布情况。其中，3种情况下TDF长度均为5 m。 从 图 5(a) 图 5(a) 中可以看出，在光纤长度相同时，3种泵浦方式在光纤末端处的放大信号光功率分别为326.68，322.71，264.73 W，其中前向泵浦的信号光功率最大，其次为双端泵浦，分段泵浦的信号光功率最低。造成这种情况的原因主要有两个：一是泵浦光转化效率不同，因为在不同的泵浦方式下光纤的最佳长度不一样，存在光纤长度过短时泵浦光未被完全吸收现象；二是分段泵浦的泵浦光以及产生激光受每段泵浦点的耦合效率以及传输效率的影响，即泵浦点越多，泵浦光以及信号光损耗的功率就越多。由 图 5(b) 图 5(b) 可以明显看出，在不考虑外部散热条件下，在传统前向泵浦及双端泵浦方案中，当泵浦光总功率达到千瓦级别时，TDF外包层边界最高温度分别为693.21 ℃和360.71 ℃。但一般双包层TDF外包层外起保护作用的涂覆层在温度达到200 ℃时就会被烧毁。分段泵浦的最高温度为173.29 ℃，而且整体分布均匀，可以达到TDF工作时涂覆层对于温度的要求。由此可见，采用分段泵浦的方式可以大大降低光纤外包层边界的温度并且使温度大致均匀地分布在光纤中。另外，前向、双端、分段泵浦的斜率效率分别为32.68%、32.27%、26.47%，可见采用分段泵浦斜率效率下降了6%左右。这有望通过进一步分析得到分段泵浦的最佳光纤长度后得以改善，而且在提高斜率效率的同时，有望进一步减小泵浦光在每个泵浦点的沉积从而降低光纤中的最高温度。 泵浦吸收系数分析 2 除光纤长度外，从公式(7)、(8)可以看出，光纤工作温度也与平均泵浦吸收系数 β p 密切相关，长度为 L 的TDF光纤的 β p 可以定义为: http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963581&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963581&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963581&type=middle 双端泵浦可以看作一种特殊的分段泵浦方式，故以双端泵浦为例，找到兼顾激光器输出性能以及温度的最佳泵浦吸收系数，为接下来分析分段泵浦提供数据依据。在总泵浦功率为1 000 W(前、后向分别为500 W)条件下，泵浦吸收系数为3.1 dB/m时, 泵浦光及放大信号光功率分布及信号光输出功率与光纤总长度的关系分别如 图 6(a) 图 6(a) 、 (b) (b) 所示。 由于双端泵浦方式下增益光纤两端都有泵浦源，所以当光纤长度一定时，信号光功率在光纤上的分布整体呈增长趋势，无法直观判断出最优光纤长度，需要通过比较不同长度的光纤后才能得出结果。如 图 6(b) 图 6(b) 所示，在泵浦吸收系数为3.1 dB/m时，光纤最优长度为5.75 m，其他泵浦吸收系数及其对应的最优光纤长度研究方法与此相同。 光纤最优长度与泵浦吸收系数的关系及光纤最优长度时对应的放大信号激光输出功率、光纤外包层边界最高温度及平均温度与泵浦吸收系数的对应关系如 图 7 图 7 所示。从 图 7(a) 图 7(a) 中可以看出，随着泵浦吸收系数的增大，光纤最优长度减小。 图 7(b) 图 7(b) 、 (c) (c) 显示，在利用 图 7(a) 图 7(a) 中得到的 β p 对应的光纤最优长度时，随着 β p 增大，放大信号激光的输出功率逐渐增大，光纤的最高温度及平均温度也逐渐升高。当 β p ≥3.1 dB/m时，信号输出功率增加逐渐趋于平缓，斜率效率大于32%且不再明显增加，但此时光纤的温度依然处于线性升高状态，所以此时不应过于追求增加斜率效率。由双端泵浦可以得出 β p =3.1 dB/m为泵浦吸收系数的最佳取值，以此为条件来提高分段泵浦的斜率效率。 分段泵浦分析 2 由以上分析可知, 分段泵浦存在一个最佳光纤长度使得光光转换效率最高并且可以减小泵浦光的沉积。 图 8 图 8 表示了在总泵浦功率为1 000 W、泵浦吸收系数为3.1 dB/m条件下，分段泵浦的TDFA中放大信号激光输出功率，以及光纤外包层边界最高温度随光纤长度变化的趋势。图中3条曲线分别表示分段数为4，5，6时的情况，其中每段的光纤长度、注入的泵浦光功率都相等。 从 图 8(a) 图 8(a) 可以看出，随着光纤长度的增加，信号光的输出功率逐渐增大，达到最大功率后逐渐下降，存在一个最佳长度。逐渐增大对应泵浦光转换为信号光的过程，再减小是由于光纤过长时信号光会被吸收而在其他波段产生自发辐射。经过计算可知，3种分段方式下，最佳长度分别为8.5，9.6，11.1 m，而且可以发现段数越多，最佳长度越长，这是因为泵浦点越多就需要更长的光纤来吸收泵浦光；同时可以看出，随着分段数的增加，TDFA信号光的最大功率逐渐减小，分段数为4，5，6时对应的最大信号功率分别为281.64，276.78，272.61 W。这是因为每个泵浦点处泵浦光并不能完全耦合入光纤中，即耦合效率不是100%，存在泵浦光损耗的因素。其次，在分段光纤熔接处，即使熔接质量很好，也会存在一定程度的激光泄露。功率损耗会随着分段数进行叠加，分段越多，损耗越大。所以在光纤温度处于正常时，在满足功率输出要求下，尽可能采取分段少的方式进行热管理。 图 8(b) 图 8(b) 表明随着光纤长度的增加，光纤最高温度逐渐降低，这也验证了之前关于光纤过短时泵浦光在泵浦点累积导致光纤最高温度升高的观点。 接下来以5段泵浦方式为例，分析其在最佳泵浦吸收系数以及最佳光纤长度下功率与温度的分布。 图 9 图 9 为泵浦光与信号光以及光纤温度在 β p =3.1 dB/m对应的最佳长度时沿光纤的轴向分布。 图 9(a) 图 9(a) 表明信号光的输出功率为276.78 W, 大于未在最优长度时的264.73 W，斜率效率为27.6%，相比其他泵浦方式只下降5%。从 图 9(b) 图 9(b) 中可以看出光纤最高温度为100.28 ℃，低于未在最优长度时的173.29 ℃。但是可以发现光纤每段的最高温度分布不均匀，最大为100.28 ℃，最小为76.45 ℃，这主要是因为泵浦光功率分布未优化，泵浦光沉积现象依然存在。 温度分布优化 2 由公式(7)、(8)可以看出在 η h 不变时，解决温度分布不均匀的主要方法有对分段光纤设置非均匀的泵浦吸收系数、不同的泵浦光功率、光纤长度等。本文主要通过对每段光纤设置不同的泵浦光功率以及光纤长度来达到温度分布均匀的目的。具体方法为：使用遗传算法同时对光纤长度和泵浦光功率进行选择，在实现温度均匀分布的同时获得较大的斜率效率。其基本方法为：设置优化参数 L i ， P p, i + ， P p, i － ， i =1, 2, …， n ， P p, i －1 － = P p, i + ，即前后相邻的泵浦光功率设置相等； n 是光纤段数， L 2 = L 3 =…= L n －1 ，除首尾两端外，中间分段光纤的长度相同，这是因为首尾端的光纤要充分吸收前向、后向泵浦功率。遗传算法的目标函数为参数优化后最终要达到的目的，其定义为： http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963597&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963597&type=small http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=2963597&type=middle 其中， T i 为光纤每段中最高温度， T ave 为每段的平均温度， P out 为输出功率。目标函数的第一部分为光纤中最高温度的标准方差，第二部分为输出功率的倒数；遗传算法优化时采用搜索目标函数最大值机制，所以取两项的相反数，同时满足得到最小标准方差和最大输出功率的结果，实现温度的均匀分布及获得较高的斜率效率。对分5段泵浦的情况进行优化计算得到结果为 L 1 =2.4 m， L 5 =2.6 m，其他分段长度为2 m，总长11 m。泵浦功率设置分别为 P p, 1 + =115 W、 P p, 5 － =89 W、 P p, 1 － = P p, 2 + =108 W、 P p, 2 － = P p, 3 + =100 W、 P p, 3 － = P p, 4 + =94 W、 P p, 4 － = P p, 5 + =96 W。 TDF中各功率成分和外包层边界温度沿光纤轴向分布分别如 图 10(a) 图 10(a) 、 (b) (b) 所示。可得TDFA的输出功率为284.5 W，斜率效率28.45%，相比单向及双端泵浦方式仅下降了4%。但光纤中最高温度仅为86.28 ℃，且每段最高温度间最大差值仅为1.97 ℃，实现了温度的均匀分布，并且与优化前相比降低了光纤最高工作温度。 结论 1 本文对基于分段泵浦的高功率TDFA热效应管理进行理论研究，利用龙格库塔法以及牛顿迭代法求解不同泵浦方式下TDFA的稳态速率方程，并根据热传导方程，模拟TDF温度沿径向和轴向的分布，通过对比得出分段泵浦对于热效应管理的优势。结合遗传算法理论，以5段泵浦为例，经过参数优化，当给定功率为5 W的种子光输入时，将吸收系数为3.1 dB/m、总长度为11 m的TDF分为2.4，2，2，2，2.6 m的5段进行泵浦，在总功率为1 000 W的793 nm激光泵浦下，得到激光输出功率为284.5 W、斜率效率为28.45%、光纤外包层边界最高温度为86.28 ℃且温度总体分布均匀，总体热效应较传统前向泵浦、双端泵浦以及未优化前的分段泵浦TDFA均有明显改善。本文的研究工作及研究方法为2 μm波段高功率TDFL的热效应管理及功率的进一步提升提供了重要的研究基础。 参考文献 龙井宇, 白晋涛, 任兆玉, 等. 2μm掺Tm 3+ 石英光纤激光器的泵浦效率分析和研究进展[J].激光杂志, 2009, 30(4):1-3. LONG J Y, BAI J T, REN Z Y, et al .. Pumping efficiency analysis and progress of 2μm Tm 3+ doped silica fiber lasers[J]. Laser J ., 2009, 30(4):1-3. (in Chinese) MA Q L, BO Y, ZONG N, et al .. Light scattering and 2μm laser performance of Tm:YAG ceramic[J]. Opt. Commun ., 2011, 284(6):1645-1647. WANG X, ZHOU P, ZHANG H W, et al .. 100 W-level Tm-doped fiber laser pumped by 1173 nm Raman fiber lasers[J]. Opt. Lett ., 2014, 39(15):4329-4332. 刘江, 刘晨, 师红星, 等. 342 W全光纤结构窄线宽连续掺铥光纤激光器[J].物理学报, 2016, 65(19):194209-1-5. LIU J, LIU C, SHI H X, et al .. 342 W narrow-linewidth continuous-wave thulium-doped all-fiber laser[J]. Acta Phys. Sinica , 2016, 65(19):194209-1-5. (in Chinese) 胡志涛, 何兵, 周军, 等.高功率光纤激光器热效应的研究进展[J].激光与光电子学进展, 2016, 53(8):080002-1-10. HU Z T, HE B, ZHOU J, et al .. Research progress in thermal effect of high power fiber lasers[J]. Laser Optoelectron. Progr ., 2016, 53(8):080002-1-10. (in Chinese) 王郡婕, 冯晓强, 冯选旗, 等.掺铥双包层光纤激光器的热效应分析和实验研究[J].激光杂志, 2013, 34(1):4-6. WANG J J, FENG X Q, FENG X Q, et al .. Thermal effect and experimental research for Tm 3+ -doped double-clad fiber laser[J]. Laser J ., 2013, 34(1):4-6. (in Chinese) WANG Y, XU C Q, PO H. Thermal effects in kilowatt fiber lasers[J]. IEEE Photonics Technol. Lett ., 2004, 16(1):63-65. WANG Y. Heat dissipation in kilowatt fiber power amplifiers[J]. IEEE J. Quantum Electron ., 2004, 40(6):731-740. HE X L, LIAO L, ZHANG F F, et al ..Analysis of laser efficiency and thermal effects in kilowatt fiber lasers based on distributed side-coupled fibers[J]. Opt. Eng ., 2018, 57(9):096108-1-12. 楼祺洪, 周军, 漆云凤, 等.千瓦级双包层光纤激光器冷却方案设计理论和实验研究[J].物理学报, 2008, 57(8):4966-4971. LOU Q H, ZHOU J, QI Y F, et al .. Theoretical and experimental research on cooling scheme design of kilowatt double-clad fiber laser[J]. Acta Phys. Sinica , 2008, 57(8):4966-4971. (in Chinese) 薛冬, 周军, 楼祺洪, 等.高功率双包层光纤激光器热效应及功率极限[J].强激光与粒子束, 2009, 21(7):1013-1018. XUE D, ZHOU J, LOU Q H, et al .. Thermal effect and power limit in high power double-clad fiber laser[J]. High Power Laser Part. Beams , 2009, 21(7):1013-1018. (in Chinese) 陈金宝, 曹涧秋, 潘志勇, 等.全国产分布式侧面抽运光纤激光器实现千瓦输出[J].中国激光, 2015, 42(2):0219002. CHEN J, CAOJ Q, PAN Z Y, et al .. Nationally produced distributed side pumped fiber lasers to achieve kilowatt output[J]. Chin. J. Lasers , 2015, 42(2):0219002. (in Chinese) 陈爽, 冯莹.高功率光子晶体光纤激光器温度分布研究[J].光子学报, 2008, 37(6):1134-1138. CHEN S, FENG Y. Temperature distribution in high power photonic crystal fiber laser[J]. Acta Photon. Sinica , 2008, 37(6):1134-1138. (in Chinese) 代守军, 何兵, 周军, 等.高功率散热技术及高功率光纤激光放大器[J].中国激光, 2013, 40(5):13-18. DAI S J, HE B, ZHOU J, et al .. Cooling technology of high-power and high-power fiber laser amplifier[J]. Chin. J. Lasers , 2013, 40(5):13-18. (in Chinese) TANG Y L, LI F, XU J Q. High peak-power gain-switched Tm 3+ -doped fiber laser[J]. IEEE Photonics Technol. Lett ., 2011, 23(13):893-895. WANG Y, YANG J L, HUANG C Y, et al .. High power tandem-pumped thulium-doped fiber laser[J]. Opt. Express , 2015, 23(3):2991-2998. MELESHKEVICH M, PLATONOV N, GAPONTSEV D, et al .. 415 W single-mode CW thulium fiber laser in all-fiber format[C]. Proceedings of 2007 European Conference on Lasers and Electro-Optics and The International Quantum Electronics Conference , Munich, Germany , 2007: 1. CREEDEN D, JOHNSON B R, RINES G A, et al .. High power resonant pumping of Tm-doped fiber amplifiers in core- and cladding-pumped configurations[J]. Opt. Express , 2014, 22(23):29067-29080. 黎大军, 杜戈果, 闫培光. LD泵浦掺铥(Tm 3+ )光纤激光器的数值分析[J].应用光学, 2007, 28(4):439-444. LI D J, DU G G, YAN P G. Numerical analysis of Tm 3+ -doped silica fiber lasers pumped by LD[J]. J. Appl. Opt ., 2007, 28(4):439-444. (in Chinese) JACKSON S D, KING T A. Theoretical modeling of Tm-doped silica fiber lasers[J]. J. Lightwave Technol ., 1999, 17(5):948-956. 欧攀, 贾豫东, 白明, 等.高等光学仿真:MATLAB版:光波导·激光[M].第2版.北京:北京航空航天大学出版社, 2014:7. OU P, JIA Y D, BAI M, et al .. Advanced Optical Simulation:MATLAB Version:Optical Waveguide·Laser [M]. 2nd ed. Beijing:Beijing University of Aeronautics and Astronautics Press, 2014:7. (in Chinese) 10.37188/fgxb20204107.0839 O437 fgxb-41-7-839 国家自然科学基金（61827818，61620106014，61975049）资助项目 Supported by National Natural Science Foundation of China(61827818, 61620106014, 61975049) 2020-05-03 2020-05-13 2020-09-10 版权所有©《发光学报》编辑部2020 Copyright ©2020 Chinese Journal of Luminescence.All rights reserved. 2020 发光学报 7 41