★ FEATURED ARTICLE 简介光束质量是衡量激光器性能的核心指标之一常通过M²因子、光束参数积BPP和数值孔径等参数描述。M²因子反映实际光束偏离理想高斯衍射极限的程度BPP则直接关联光纤耦合的容限判断。与单模光纤激光器不同半导体激光器因发光区结构导致快慢轴束腰位置分离形成像散使单一M²数值难以完整表征其聚焦能力。因此工程中需分轴测量M²并记录快慢轴束腰位置差同时借助CCD成像、刀口法或狭缝扫描法沿光轴采集光斑序列通过双曲线拟合重建束腰。这些评价方法广泛应用于泵浦源验收、光纤耦合方案验证及直接半导体加工等场景。针对半导体激光器光束质量评价中的常见陷阱如CCD饱和、杂散光干扰、像散遗漏与热漂移文章结合实践给出修正思路帮助工程人员建立可靠的光束质量验收流程。1. 半导体激光器光束质量评价功率表之外的另一个硬指标一台 50W 的半导体激光器模块加同一只聚焦镜一个能在焦点打出锐利的加工点另一个在焦面附近就摊成一条亮带。当时第一反应是重测功率结果两台功率读数几乎持平差异不在功率而在光束质量。半导体激光器发光区狭长、发散角大快慢轴严重不对称单看功率没法解释聚焦能力和光纤耦合效率必须把光束质量的评价维度补上。这篇笔记就围绕“半导体激光器光束质量评价方法”展开M² 因子怎么定义、像散为什么让单数值失效、实际测量中 CCD 和刀口法怎么选、从光斑序列到 M² 的数据处理怎么做以及五个高频翻车点。无论做泵浦源验收、光纤耦合方案验证还是直接半导体加工应用这套评价思路都直接对应产线和研发的验收环节。2. 把半导体激光器的束腰拆开快慢轴、M² 因子和像散的三角关系2.1 从发光区到远场先看清“束腰”天生不止一个半导体激光器的发光区是矩形波导结构快轴方向垂直于 pn 结的发光区宽度通常只有 12 μm慢轴方向平行于 pn 结可能是几十到几百微米。这个几何差异直接决定了远场光斑形态快轴发散角半角往往在 25°35°慢轴半角通常在 3°8°。用数值孔径来体会快轴 NA≈0.5这已经是相当强聚焦系统才有的指标裸芯片出光后几毫米内就能扩散开所以绝大多数应用都先做快轴准直。准直之后快慢轴的光斑尺寸和发散角仍然不是同一个量级。以常见的百微米级发光区单管为例准直后慢轴方向的光斑宽度可能是快轴方向的 510 倍但慢轴发散角小得多。于是评价光束质量时就必须分轴讨论快轴接近衍射极限M² 值可以做到 1.21.5慢轴受发光区宽度限制且伴随模式较多M² 值落在 1025 属于常态。把一个标量 M² 套在快慢轴上等于把两个不同物理过程的数字平均既掩盖了快轴的优越性也低估了慢轴对聚焦的拖累。从测量角度讲分轴意味着束腰位置和束腰半径要分别拟合两条曲线。慢轴束腰不在快轴束腰的同一个 z 坐标上这就是后面要说的像散。先把这条根因记下来半导体激光器的“束腰”不是单数而是快轴和慢轴两套参数强行合并只会让报告失去工程意义。2.2 M² 因子、BPP 和数值孔径一个数字说不清全部信息M² 因子的定义是实际光束的束腰半径与远场发散角的乘积除以同波长理想高斯光束的对应乘积。公式通常写成M² π · ω₀ · θ / λ其中 ω₀ 是束腰半径半宽θ 是远场发散角半角λ 是波长。对理想高斯光束 M²1对实际光束 M²1。单模光纤激光器 M² 可做到 1.051.1而半导体激光器的慢轴动辄十几甚至二十几两者相差一个数量级所以不能把半导体激光器当成“接近高斯”的光源来处理。工程上更常用的是 BPP光束参数积定义为束腰半径乘以发散角半角BPP ω₀ · θBPP 的单位是 mm·mrad。光纤耦合方案里光纤的芯径和数值孔径决定了可接受的最大 BPP比如 50 μm 芯径、NA0.15 的光纤可通过的 BPP 约等于 25 mm·mrad。如果半导体激光器慢轴 BPP 超过了这个值耦合效率就会断崖式下跌。这就是为什么做半导体激光器应用的人要同时看 M² 和 BPPM² 回答“光束离衍射极限多远”BPP 回答“能不能塞进这根光纤”。两者都可以由束腰半径和发散角推导但工程语境里 BPP 更接近系统容限判断。还有一个容易混的指标是数值孔径 NA。实际测量时发散角一般不直接给 NA因为半导体激光器远场并不是理想的朗伯分布快轴接近余弦分布慢轴通常有旁瓣。直接用半角正弦值估算 NA 只在近轴假设下成立遇到大发散角快轴时误差明显。我一般会记录一组完整的远场强度分布而不是只记一个 “30° 发散角” 的笼统数字。2.3 像散不处理聚焦判断就会失真像散可以理解为快轴和慢轴的等效发光面不在同一个轴向位置。物理原因是波导结构的不对称性和材料折射率分布的差异常见的等效失焦量为 550 μm。这个数值对裸芯片来说看似微小但在聚焦光学系统里它等效于在焦点处将两个轴的光斑错开一段距离结果是焦面上快轴聚焦良好时慢轴还在扩散中或者反过来。测量中最典型的表现是用一面观察镜从 0 到峰值找最小光斑快慢轴最小值出现在不同位置如果用一个方向的最小光斑位置当整个光束的焦点另一个方向的发散就会掩盖实际聚焦能力。更隐蔽的是准直系统自带像差或安装倾斜会使测量像散值进一步放大。所以评价方案里必须有像散这个维度并且要在报告里给出“快轴束腰位置 — 慢轴束腰位置”的差值而不是只给一个 M²。另一个需要提醒的是半导体激光器在低功率和高功率下热透镜效应不同束腰位置和像散都会轻微漂移。测量功率要和实际使用功率一致否则测出来的 M² 只代表工作点 A不代表工作点 B。这一条后面在避坑章节还会展开。3. 搭一套 Z 扫描光束质量测量台方法选型、硬件参数和最小步骤3.1 CCD、刀口法和狭缝扫描三种测量原理怎么选光束质量测量常见的思路是沿传播方向z 向采集一系列光斑宽度再拟合成双曲线。区别只在“宽度怎么测”。CCD 法直接用相机拍光斑的强度分布通过灰度图像计算二阶矩来定义宽度。优点是能看到完整的功率分布形状旁瓣、干涉条纹、像散都能一眼看出来适合研发阶段做诊断缺点是动态范围有限对微弱旁瓣不敏感像素分辨不足时测量结果偏大。刀口法用一片刀刃横向扫过光束记录透过的功率随位置的变化用功率-位置曲线求微分得到强度剖面。优点是动态范围可以做到很大因为用的是功率计没有像素饱和问题适合高功率半导体激光器缺点是只能得到一维投影信息像散分布要分别在两个轴上扫描而且刀刃位置的对准精度直接影响结果。狭缝扫描法类似刀口法但用窄缝代替刀刃直接测窄缝透过的功率积分。它对杂散光抑制比刀口法更好但狭缝宽度要对光斑尺寸做匹配通常要求缝宽小于光斑直径的 1/10 甚至 1/20测量窄光束时容易遇见信号太弱的问题。三种方法选型时可以按这张表取舍方法动态范围二维形状信息窄光束适配关键不稳因素CCD 成像一般8/12/16bit完整像素分辨率受限饱和、杂散光、像元响应不均刀口法高功率计无需逐轴扫好但定位要求高刀刃平行度、机械回程差狭缝法中高无需逐轴扫好对缝宽敏感缝宽与光斑匹配、衍射干扰我个人的习惯是研发阶段用 CCD 法加上单独一维的刀口复核产线快速验收则用刀口法或狭缝法因为不需要为每台设备调试相机曝光。3.2 Z 扫描平台的硬件参数行程、取样密度和衰减片怎么定Z 扫描的核心是让光束在被测光学系统后形成收敛束腰然后在束腰前后不同 z 位置测量光斑宽度。先选聚焦透镜。透镜的焦距决定瑞利长度的实测尺度瑞利长度 zR π·ω₀²/(M²·λ)数值越小测量窗口越窄。对半导体激光器快轴束腰半径可能只有几十微米zR 可能只有几毫米慢轴束腰半径数百微米zR 在几十毫米量级。一个 100 mm 焦距的聚焦镜在这个场景下是常见起点焦深适中气动行程需求不至于过大。行程范围通常要求覆盖束腰前后至少 35 倍瑞利长度。如果快轴 zR 只有 5 mm慢轴 zR 有 40 mm行程取 ±60 mm 才能把慢轴的双曲线尾巴采集完整。建议用行程 150 mm 左右的电动位移台留出余量。取样密度也分轴在束腰附近 2 倍 zR 以内每 12 mm 取一个点在两侧远区每 510 mm 取一个点。总采样点建议在 1525 个之间太少拟合不稳太多数据冗余影响测量时长。功率控制方面半导体激光器尤其是高功率单管裸出光就能熔化刀刃或打伤 CCD衰减是第一优先级。建议先用反射式衰减器或衰减片堆把功率降到 1 mW 以下再用可调衰减片微调到 CCD 灰度的 60%80% 动态范围。避免把衰减片放在聚焦透镜和探测器之间因为衰减片表面的散射光会被误判为光束外围的功率分布我见过有人因此把 M² 从 12 测到 25。3.3 一段可复用的最小采集流程位移台、相机和功率计的协同控制这套流程的核心是让位移台走到指定位置等系统稳定再触发相机或功率计读数。下面是一段精简的 Python 控制思路不是某个厂商的 SDK 绑定代码仅体现流程骨架import numpy as np import time # 采集参数 z_start, z_end, n_points -60.0, 60.0, 21 z_list np.linspace(z_start, z_end, n_points) # 预扫一次找到光斑所在横向位置避免漂出相机视野 stage_x.move_to(0.0) stage_z.move_to(0.0) frame camera.grab(exp_time2.0) max_ymax np.unravel_index(np.argmax(frame), frame.shape) stage_x, stage_y 的基础定位先对着光斑中心 collected [] for z in z_list: stage_z.move_to(z) # 移动位移台 time.sleep(0.15) # 等待机械定位稳定 frame camera.grab(exp_time1.0) # 根据衰减情况调整曝光 collected.append((z, frame))这里有几个参数值得解释。z_list 按等差数列分布但实际操作时我会在焦点附近加密先用大步长粗扫定位最小光斑位置再把 z 换成加权分布。时间等待 0.15 秒是机械位移台稳定到微米级所需的最小值如果用的是步进电机加丝杠时间要更长。曝光时间的选取遵循“单帧不过曝”原则宁可多采几帧求平均也不要让中心区域饱和。采集完成后立即存原始图像不要只存计算后的宽度值后面如果发现可疑结果还能重新处理。采集流程里还有一个隐蔽环节——横向对准。半导体激光器经准直后快轴方向光斑可能在移动过程中有微小平移如果位移台的直线度不够好光斑会在相机上左右移动导致测得的宽度虚高。建议在每个 z 位置都重新计算光斑质心并在后续计算中把质心偏移记录下来。质心动辄偏移几个像素宽度数据就会失真这一条放在避坑篇里继续讲。4. 从光斑序列到 M² 值Python 拟合与双轴像散重建4.1 背景扣除与二阶矩为什么不能直接数亮斑像素拿到一组 z 位置的光斑图像第一步是背景扣除。相机暗电流和光学系统杂散光叠加在图像上如果不扣二阶矩会把微小但分布范围很广的背景当作光束外围功率导致宽度计算严重偏大。一个实用的做法是取图像四个角落的灰度中值作为背景估计再做逐像素相减并强制负值为零。def preprocess(frame, border15): # frame: 原始灰度图 # border: 取图像边角区域的宽度 corner np.concatenate([ frame[:border, :border].ravel(), frame[:border, -border:].ravel(), frame[-border:, :border].ravel(), frame[-border:, -border:].ravel() ]) bg np.median(corner) img frame.astype(np.float64) - bg img[img 0] 0 return img背景中值法的前提是四角确实不存在光束外围分布。对窄光束和大视场相机这个假设通常成立如果光束发散到占满全画面需要先用刀口法做一个总功率剖面确认背景区域。做完背景扣除后计算质心和二阶矩就有意义了。下面是基于灰度二阶矩的宽度计算def beam_width(img, pixel_size5.5e-6): ys, xs np.indices(img.shape) total img.sum() if total 0: return 0.0, 0.0, 0.0 cx (img * xs).sum() / total cy (img * ys).sum() / total sigma_x2 (img * (xs - cx)**2).sum() / total sigma_y2 (img * (ys - cy)**2).sum() / total # d4sigma 宽度像素尺寸换算成米 w_x 4.0 * np.sqrt(sigma_x2) * pixel_size w_y 4.0 * np.sqrt(sigma_y2) * pixel_size return cx, w_x, w_y这里要说明为什么用 4 倍根号二阶矩。对高斯分布二倍二阶矩半径点和光强衰减到 1/e² 处的半径相等直径就是 4σ这是 ISO 11146 系列标准约定俗成的定义方式。实际光束往往不是理想高斯二阶矩定义的宽度更严格但也会把远场旁瓣计入因此对半导体激光器慢轴算出来的宽度经常比“肉眼可见光斑直径”大 20%30%这属于正常现象不是程序 bug。4.2 双曲线拟合束腰半径、焦点位置和 M² 值的相互关系采集到每个 z 位置的宽度 w(z) 后用双曲线模型拟合w(z) w₀ · sqrt(1 ((z − z₀) / z_R)²)其中 w₀ 是束腰半径z₀ 是束腰位置z_R 是瑞利长度。拟合得到这三个值后M² 由下式计算M² π · w₀² / (z_R · λ)用 Python 的 scipy 可以这样写from scipy.optimize import curve_fit def caustic(z, w0, z0, zR): return w0 * np.sqrt(1.0 ((z - z0) / zR)**2) # 假设 z_list 是位置(m)w_list 是对应宽度(m) p0 [1e-4, 0.0, 0.02] popt, _ curve_fit(caustic, z_list, w_list, p0p0, maxfev5000) w0_fit, z0_fit, zR_fit popt M2 np.pi * w0_fit**2 / (zR_fit * lambda_m)参数初值的选择会影响拟合稳定性。p0 的 w0 可以用最小宽度估算z0 用最小宽度对应的 z 位置zR 用半宽扩展一倍的区间跨度估算。如果把 zR 初值给成 1 mm 而实际是 50 mm拟合很容易陷入局部极小。拟合时的权重也有讲究光斑越大宽度测量的绝对误差通常越大但如果直接对 w(z) 做最小二乘大光斑数据点权重更大。更合理的做法是对 w²(z) 拟合因为大光斑处的 w² 动态范围大而焦点附近的窄光斑在平方空间里更有辨识度。实际中我会先做不加权拟合再看残差分布若有明显趋势再改成对平方数据的加权拟合。另一个来自实践的建议半导体激光器的 CCD 图像在焦点附近快轴方向可能接近衍射极限光斑直径只有几微米如果像素尺寸 5.5 μm一个光斑只占约一个像素这时二阶矩计算会有巨大量化误差。这种情况要优先采用刀口法复核或把相机光路加显微放大系统。不要盲目相信像素级光斑算出的 M²那往往是系统误差主导的结果。4.3 快慢轴分离拟合再把像散写进报告半导体激光器的图像是椭圆光斑宽度计算要将图像旋转到快慢轴方向上分别对两个方向做投影和拟合。旋转角度怎么定先在离束腰较远的位置采集一帧光斑对二维灰度做主轴分析用质心协方差矩阵的特征向量确定快慢轴方向。之后所有帧都沿这个固定方向计算宽度避免在每帧上独立做主成分分析导致角度抖动。def fit_m2_per_axis(z_list, w_list, lam): # 对单个轴做双曲线拟合 popt, _ curve_fit(caustic, z_list, w_list, p0[3e-5, 0.0, 0.01]) w0, z0, zR popt m2 np.pi * w0**2 / (zR * lam) return m2, w0, z0, zR # 假设 widths_fast 和 widths_slow 是两个轴各自的宽度序列 m2_fast, w0_fast, z0_fast, _ fit_m2_per_axis(z_list, widths_fast, lam) m2_slow, w0_slow, z0_slow, _ fit_m2_per_axis(z_list, widths_slow, lam) astigmatism z0_fast - z0_slow # 快慢轴束腰位置差单位米 print(f快轴 M²{m2_fast:.2f}, 慢轴 M²{m2_slow:.2f}, 像散{astigmatism*1e6:.1f} μm)这段代码输出像散差单位是微米。报告里除了给 M² 和 BPP还应该给束腰位置差因为光纤耦合对准时这个差值直接影响耦合透镜的设计。部分半导体激光器模块在出厂时会内置快轴准直微透镜这样快轴束腰被移到透镜后某个位置慢轴束腰还留在芯片附近像散差值可能从数微米放大到数百微米。测量结果和大大小小的像散差可能先让人怀疑是测量错误但实际上这就是系统的真实状态后续光路设计必须考虑它。5. 半导体激光器光束测量中的五个高频翻车点现象、原因与修正5.1 CCD 饱和导致束腰假性变窄中心出现平顶现象是图像中心灰度一片平顶宽度计算偏小尤其是光斑中心的“质量”被均匀化后二阶矩分母变大而分子被削减M² 拟合结果异常偏低。原因很简单衰减没做够或者曝光时间过长相机在光斑峰值区域达到满量程。解决方法是衰减片先粗调到中心灰度约为 60% 满量程再对着目标工作点微调。如果已经拍了饱和帧不要寄希望于软件削峰恢复直接重新采集。我的习惯是在自动采集流程里加入一个“保护动作”每帧判断最大灰度是否超过 95% 满量程是则自动降低曝光并打标记。5.2 杂散光抬高 M²误判为光束质量变差现象是 M² 结果显著大于同类模块的参考值且重复性良好但各次差异明显。原因是聚焦透镜、衰减片、刀刃表面的散射光进入了探测器视场这些杂散光在光斑外围形成微弱但宽广的背景。CCD 法对这种情况尤其敏感因为背景扣不干净二阶矩会把边角区域的微弱噪声当成真实功率分布。解决方案是在聚焦透镜前加小口径光阑限制有效光束范围并在测量光路上用消光材料遮挡不必要的反射路径。另一道保险是给 z 轴平台加“异议复核”用刀口法在同一位置重新测一遍如果刀口法 M² 明显低于 CCD 法基本可以实锤杂散光问题。5.3 焦点两侧采样不对称拟合出的 z0 偏移现象是双曲线拟合残差在一侧系统性偏高z0 拟合值偏离肉眼观察到的最小光斑位置。原因是位移台行程没覆盖到束腰两侧各 23 倍瑞利长度等效于拟合曲线缺了一侧的“肩膀”。比如行程范围是 −40 mm 到 20 mm焦点在 0 附近那么一侧数据量足够另一侧只有一两个点拟合的自由度被单向拉偏。解决方法是先用粗扫找到 z0 的大致位置再把行程范围调整成围绕 z0 对称。这一点在快轴瑞利长度小、慢轴瑞利长度大的系统里尤其容易犯为了省时只做快轴的扫描慢轴数据远没收敛。5.4 像散被当作系统误差报告里不给位置差现象是两个模块标称 M² 相同实际光纤耦合效果差异巨大回查数据发现一个像散接近零另一个像散 300 μm。原因是测量团队只处理了快慢轴各自的 M²把束腰位置差放在表格角落没有当作关键指标。如果有自动化的拟合脚本哪怕不删像散也应该生成一张包含 z0 差的记录。验收入口加一条硬性规则对同一个被测模块像散差超过设计容差就判别不合格不要只看 M² 单值。5.5 热漂移让束腰位置在测量半小时内慢慢移动现象是同一台模块连续测三次M² 一致但 z0 一次比一次偏向某个方向偏移量逐渐增大。原因是半导体激光器芯片在通电后温度升高热透镜效应改变了等效发光面位置导致束腰位置随时间漂移。更隐蔽的是这种漂移在低占空比测试时几乎不存在在连续大功率工作时才明显。解决方法是测量前让模块在额定电流下预稳定至少 10 分钟并把热沉温度控制在允许范围内高功率测试时考虑用短脉冲模式来抑制热透镜同时记录环境温度和热沉温度。如果发现 z0 漂移超过 0.5 倍瑞利长度本次数据作废重新进入热稳定流程。6. 用一把光学刀口做现场快速复核十分钟建立起对测量台的信心最后分享一个我一直在用的现场复核技巧不需要重新搭一套完整测量设备只需要一个刀口和一个功率计就可以在十分钟内验证 CCD 数据的可信度。具体做法是在 CCD 测出的束腰附近选三个位置z0、z0−zR、z0zR。用刀口分别沿快轴和慢轴方向扫描求扫描方向上的光斑直径和 CCD 结果对不上就说明 CCD 数据的某一步处理有问题。这个习惯帮我在现场抓住过不少问题。有一次某测试台上 CCD 结果正常但同一台激光器换到另一个衰减片后慢轴 M² 从 14 跳到 22刀口复核显示真实值接近 15——问题出在新衰减片的楔角把我的光束横向错位导致后续的质心跟踪失效。从那以后我在每批数据开始时都花十分钟做这个对照确认整套光路的状态没有在更换部件后改变。在做刀口复核时要注意刀刃切入方向必须垂直于被测轴否则扫描宽度会被余弦误差放大。刀口法宽度定义通常取 10%90% 透射功率区间换算和 CCD 的二阶矩定义不同直接比较数值没有意义。比较的正确方式是双方都换算成等效高斯半径或者按同一套双曲线模型拟合出的 M²。我的做法是把刀口数据单独拟合成一组宽度序列和 CCD 数据的拟合结果比较 M² 值本身而不是比较某个单一位置的宽度。对一套测量台的最终检查通常有两条源验证路径。第一条用一根单模光纤激光器做标定光源M² 应落在 1.01.1 区间第二条用已知 BPP 的多模光纤跳线做模拟光源测量结果与光纤参数推导值偏差应在 10% 以内。如果这两条都不满足优先排查衰减片、聚焦镜脏污和位移台直线度再去怀疑光源本身的光束质量。落到验收环节我习惯把报告拆成三层第一层给 M² 快轴/慢轴双值第二层给 BPP 和束腰位置第三层给像散差和测量时的热状态。买方能直接拿来跟光纤耦合容限对比研发也能拿回去做光学设计的输入参量。坚持做这套验收方法之后我唯一后悔的事是没有早点把像散差往上提一个级别——希望这套思路对你也有同样的启发希望帮到你。本文还有配套的精品资源点击获取
阅读完成 · 觉得有帮助？
