我最早接触光纤布拉格光栅FBG的时候第一反应是这不就是光纤上刻了一排“印痕”嘛为什么能滤光为什么传感器都爱用它后来真正上手搭系统、做解调、写标定程序才意识到如果没把“波导结构”和“布拉格波长”这两个底子抠透后续看谱线漂移、算温度灵敏度、设计光栅阵列的时候全是靠猜。所以这个系列我打算从源头梳理第一篇就聚焦两件事FBG 的波导结构到底长什么样布拉格波长公式又是怎么一步步推出来的。这篇笔记适合三类人刚接触光纤光学、想搞清楚 FBG 原理的学生在实验室里看过光谱仪上那个反射峰、但没系统推过公式的工程师以及做传感系统集成、需要根据公式估算光栅工作波段和灵敏度的开发者。我会尽量把公式推导过程写完整把每一步的物理含义和工程含义都点出来不跳步。1. 先从“一根光纤”讲起FBG 的波导结构与关键参数1.1 光纤的结构是 FBG 的物理舞台标准单模光纤比如最常见的 SMF-28横截面从里到外分三层最中心的纤芯core、包着纤芯的包层cladding以及最外层的涂覆层coating。纤芯直径通常在 8~10 μm 左右包层直径统一是 125 μm涂覆层是 250 μm 左右的丙烯酸酯或聚酰亚胺材料。真正让光被“关”在纤芯里传播的是纤芯和包层之间的折射率差。石英玻璃本身折射率大约在 1.44~1.46 之间纤芯里掺了锗Ge之后折射率会略高一点典型差值 Δn ≈ 0.003~0.005。这个差值看着很小但就是它决定了光纤能够支持哪些模式、光的能量在横截面上怎么分布。FBG 的本质是在这段纤芯内部引入一个“周期性”的折射率调制。不是把光纤切断也不是在表面开槽而是利用紫外光通过相位掩模板照射在掺锗光纤上引起纤芯材料局部折射率发生永久性改变。于是沿光纤轴向看过去纤芯折射率分布变成了n(z) n_core δn · sin(2πz / Λ)这里的 Λ 就是光栅周期通常用纳米nm来度量典型值在几百纳米到几微米之间。δn 是折射率调制的幅度通常在 10⁻⁵ ~ 10⁻³ 这个量级。这个周期性的折射率扰动就是后续一切“反射”和“透射”现象的来源。1.2 光栅“写”在波导里和写在平面上完全不同如果只看到“周期性结构”很多人会想到平面衍射光栅。区别在于平面光栅是在自由空间里衍射FBG 是在波导里衍射。光在光纤里传播时并不是一束无限平面波而是被波导边界约束的导模。这个约束条件直接影响相位匹配条件。更准确地说光纤的导模有自己的传播常数 β (2π / λ) · n_eff其中 n_eff 是模式的有效折射率。有效折射率不是一个单纯的“材料折射率”它是模式在纤芯、包层中能量分布的综合体现。因为导模的能量不是 100% 都在纤芯里有相当一部分尤其在长波长段会渗透到包层中所以 n_eff 的数值会介于包层折射率和纤芯折射率之间并且随波长变化。这种色散特性在推导布拉格条件时不能忽视。FBG 的波导结构还有一层特殊性因为折射率调制区域只在纤芯里所以它主要影响纤芯导模的传播状态对包层模的影响相对弱。这也解释了为什么我们在实验中看到的 FBG 反射谱比理论上的“单峰”要复杂一些因为还有包层模耦合、辐射模耦合等次要过程参与进来。后面在推导完基本公式后我会专门提一下哪些项在工程上可以忽略、什么条件下不能忽略。1.3 从波导结构的视角看“光栅强度”和“光栅长度”除了周期 Λ 和调制深度 δn波导结构里还有两个参数决定了光栅的反射能力光栅长度 L 和有效折射率调制幅度 δn_eff。把光栅理解成无数个弱反射面的叠加每个折射率突变界面都会反射一小部分光虽然单个界面的反射率极低典型在 10⁻⁶ 以下但在某个特定波长下所有界面的反射光相位一致就能相干叠加成很强的反射峰。光栅越长参与叠加的界面就越多反射率越高同时反射谱带宽越窄。这个关系很像天线阵列阵元越多定向性越强。在弱光栅δn 很小近似下反射率峰值可以用公式表达R_max tanh²(π · δn_eff · L / λ_B)当光栅比较弱时tanh(x) ≈ xR_max 近似正比于 (δn_eff · L)²。所以做光纤光栅的时候想提高反射率加长光栅比单纯增大 δn 更容易实现而且对谱线宽度的影响是两回事。这个“长度换反射率”的思路在传感阵列设计中非常重要。2. 模式耦合到底是谁和谁在“耦合”从波导光学看光栅2.1 为什么不能用“光的反射”一句话糊弄过去很多科普文章说 FBG 就是“特定波长的光被反射回来”这句话没错但不够。既然标题里带了“波导结构”就必须回到模式耦合的框架里去理解入射的纤芯导模在周期性折射率扰动下会与哪些模式发生能量交换。耦合模理论Coupled Mode Theory是分析周期性波导的标准工具。核心思路是把折射率调制看作一个微扰导致原本正交的模式之间发生能量转移。在 FBG 中最主要的过程是正向传播的纤芯基模LP01与反向传播的同一模式之间的耦合。也就是说光栅像一个“镜子”把一部分正向能量转化为同波长的反向能量。这个正向和反向模式之间的耦合专业性上叫“Bragg 反射”或“导模-导模反向耦合”。相位匹配条件决定了哪个波长能被反射而这个条件又由波导结构的传播常数决定。用传播常数写出来就是β_forward β_backward 2π / Λ · m其中 m 是光栅的衍射级次。由于 β_backward -β_forward同一模式反向传播时传播常数方向反了所以上式变成 2β 2πm / Λ。2.2 传播常数和有效折射率的关系把传播常数 β 表示成有效折射率的形式是最常用的套路β (2π / λ) · n_eff这个式子看着简单但要理解n_eff 不是材料的固有属性它是模式在某一个具体波长下的“等效”折射率取决于纤芯尺寸、折射率差、波长等因素。同样的光纤不同波长下的 n_eff 不一样。这正是 FBG 反射峰波长不能简单用“光纤材料折射率”来算的原因必须用有效折射率。举个数值例子SMF-28 在 1550 nm 附近纤芯折射率约 1.4504包层折射率约 1.4440但基模的有效折射率大约是 1.4468 到 1.4475 之间具体数值要看波长和仿真精度。如果我们错用纤芯折射率 1.4504 去算布拉格波长会差出好几个纳米这在滤波和传感场景里都是不可接受的误差。2.3 耦合系数和光栅的“强度特征”耦合模理论里有个关键参数叫耦合系数 κ它衡量的是折射率调制对模式耦合的贡献强度。对于均匀正弦调制光栅κ 可以近似表达为κ π · δn_eff / λ这里的 δn_eff 是有效折射率调制的幅度不是峰值。κ 的单位是 1/m它直观的意义是“单位长度内模式能量转换的速率”。κ 越大光在光栅里走一小段距离就能完成明显的能量交换反射峰自然更强、谱线也相应更宽。这个 κ 还与“光栅强度”的分类直接相关。工程上定义了一个无量纲参数 κL耦合系数乘以光栅长度用来区分弱光栅和强光栅κL 1弱光栅反射率低适用于需要小扰动、窄带宽、甚至需要倾斜光栅的情况。κL ≈ 1~3中等强度光栅传感系统里最常用反射率适中谱线质量好。κL 3强光栅反射率接近饱和适合做滤波器、激光器腔镜但传感解调时谱线偏宽。之所以要讲这部分是因为后续推导布拉格波长时很多教材默认是“弱光栅”条件直接忽略高阶修正。实际用强光栅时反射峰中心会有轻微移动这种系统误差在做高精度传感标定时必须考虑进去。3. 布拉格条件是怎么一步步导出来的从相位匹配到波长公式3.1 用“弹簧队列”理解什么叫相位匹配我想先用一个生活化类比把相位匹配说明白然后再给数学推导。想象一排人站在一条跑道上每个人都在原地匀速拍手。你从跑道一端走过去每遇到一个拍手的人就听到一次掌声。如果所有人的拍手节奏相同你可能在某些位置听到的掌声特别密集、特别响因为多个人拍手的声波在某处叠加了而在另一些位置掌声互相抵消。周期性结构对光波的作用就是让微弱的散射波在某一个方向上同相叠加。对 FBG 来说入射光在每一个折射率周期边界上都会被微弱地反射。如果这些反射波在某一个波长下恰好相位相同它们就会形成一束很强的反向光。这个“恰好”就是布拉格条件。3.2 严格一点的数学推导不跳步版推导的目标是求出反射峰的中心波长 λ_B。我尽量从两个角度推这样你可以挑自己习惯的方式理解。方法一一维相位匹配条件设正向传播模式的电场的相位因子为 exp(iβz)其中 β 是传播常数。纤芯里折射率周期调制产生的“极化强度扰动”中包含一个相位因子与光栅波矢相关。考虑正向光通过每一个周期边界时产生的反向散射。为了让每个周期边界的反向散射同相叠加相邻周期边界之间的往返相位差必须是 2π 的整数倍。相邻两个周期界面的间距是 Λ。正向波从界面 A 传到界面 B相位变化是 βΛ在界面 B 被反向散射后反向波从 B 回到 A相位变化又是 βΛ。所以总的往返相位差是 2βΛ。为了让散射波在 A 处同相叠加2βΛ 2πm这里 m 是整数衍射级次。把 β 2πn_eff / λ 代入2 · (2πn_eff / λ) · Λ 2πm化简得到mλ_B 2n_effΛ这就是光纤布拉格光栅的基本方程。对于最常见的 m 1即一阶光栅λ_B 2n_effΛ这个式子干净漂亮但背后的物理过程其实包含了好几层假设只考虑纤芯基模的正向/反向耦合、忽略包层模、假设弱调制下有效折射率不随光栅本身改变。工程里大部分 FBG 设计都以这个公式为首位近似。方法二从光的动量守恒推导物理学里还有一种更直观的推导思路就是把光子和光栅都看作具有动量的对象。光子动量是 ħk ħ(2π / λ)光栅的周期调制相当于一个动量源动量为 ħ(2π / Λ)针对一阶过程。正向光子的动量是 ħ(2πn_eff / λ)反向散射后的光子动量是 -ħ(2πn_eff / λ)。动量守恒要求正向光子动量 光栅提供的动量 反向光子动量用向量写就是ħ(2πn_eff / λ) ħ(2π / Λ) -ħ(2πn_eff / λ)略去 ħ 然后整理2 · (2πn_eff / λ) 2π / Λ两边约掉 2π得到2n_eff / λ 1 / Λ再整理成 λ 2n_effΛ。这跟方法一的结果完全一致。动量守恒和相位匹配本质上是同一个物理的不同表述但动量视角更容易理解为什么周期越小、反射波长越短光栅周期相当于给光子“借出”的动量越大需要匹配的波长就越短。3.3 高阶衍射级次和模式类型的影响前面公式里的 m 如果大于 1比如 m 2对应波长条件变为λ_B(m) 2n_effΛ / m这意味着同样的光栅周期下还存在更短波长的高阶布拉格共振。但实际上一阶共振m1最容易观察到因为它的反射强度最高。高阶共振要么波长落在材料吸收区要么与包层模耦合纠缠在一起在实验中往往表现为短波长侧的宽包络或杂散峰。另外正向传播的模式如果耦合到反向传播的“包层模”相位匹配条件会变成λ_coupling (n_eff_core n_eff_cladding_mode) · Λ这时反射峰或者损耗峰的位置比布拉格波长更短。在光纤光栅的透射谱上你可能看到主反射峰左侧有一串小损耗峰那就是包层模耦合造成的不是光栅缺陷。设计时要分清。3.4 公式里的每一项和工程单位换算做实际设计的时候公式 λ_B 2n_effΛ 里各物理量的单位必须统一。n_eff 无量纲Λ 如果用纳米nmλ_B 自然也是纳米。举个例子假设要做一个 C 波段 1550 nm 的 FBGn_eff 在室温下是 1.447那么Λ λ_B / (2n_eff) 1550 / (2 × 1.447) ≈ 535.6 nm这个周期低于可见光波长但比 X 射线晶体衍射的晶格常数大得多。用相位掩模法制备 FBG 时掩模板的周期就是 535.6 nm 附近紫外光写入时产生半周期调制对应一阶光栅。如果用 1064 nm 波段的光纤激光器做 FBGn_eff 大约在 1.46 附近因为波长更短模式更集中在纤芯有效折射率略高那么Λ 1064 / (2 × 1.46) ≈ 364.4 nm这个尺度在工程上依然可以加工但意味着相位掩模板需要重新定制不能拿 1550 nm 的模板通用。光栅设计和“目标波长”强绑定不是所有 FBG 都能在一根光纤上随意刻出来这个问题我在实际项目中踩过不少坑。4. 这个公式里面藏着哪些工程“命门”中心波长、有效折射率与谱线形态4.1 为什么布拉格波长的“温度灵敏度”也能从公式里导出来公式 λ_B 2n_effΛ 最有价值的地方不是它能算波长而是它能告诉我们波长受哪些物理量的影响。对温度 T 求导(1/λ_B) · (dλ_B/dT) (1/n_eff) · (dn_eff/dT) (1/Λ) · (dΛ/dT)右边第一项是热光效应材料的折射率随温度变化第二项是热膨胀效应光纤长度和光栅周期随温度变化。对于石英光纤在 1550 nm 附近热光系数约 6~8 × 10⁻⁶ /°C热膨胀系数约 0.5~0.6 × 10⁻⁶ /°C。两项相加FBG 的温度灵敏度大约为 10~11 pm/°C。也就是说温度升高 100°C中心波长漂移约 1 nm。这个数值范围我在实验里验证过很多次。使用普通单模光纤刻写的 FBG在 1550 nm 波段温度标定结果通常在 10.3~11.2 pm/°C 之间。具体数值取决于光纤掺杂浓度和有效折射率色散所以工程上做温度传感器时必须逐支标定不能只套理论值。4.2 轴向应变灵敏度推导材料力学和光学公式联手对轴向应变 ε 求导得到(1/λ_B) · (dλ_B/dε) (1/n_eff) · (dn_eff/dε) (1/Λ) · (dΛ/dε)第二项其实就是 1因为 dΛ/dε Λ代表几何伸缩的贡献第一项被称为弹光效应通常是负值大约 -0.22~-0.25 之间对于石英光纤。所以总灵敏度系数约为 0.75~0.78对应单位应变 (με) 的波长变化K_ε ≈ λ_B × 0.78 × 10⁻⁶ / με在 1550 nm 波段K_ε ≈ 1.2 pm/με。这意味着 1 με百万分之一应变相当于 1 米光纤伸长 1 微米会引起约 1.2 pm 的波长偏移。这个量级很微妙如果解调仪的光谱分辨率不够高就很难区分微小应变和噪声。4.3 有效折射率的色散效应算出来的波长为什么和实测差一点前面所有推导默认 n_eff 在感兴趣的波长范围内是一个常数但严格来说 n_eff 是波长 λ 的函数。光纤的波导色散会导致如果先初定一个 λ_B再用这个 n_eff 去算 Λ然后反推 λ_B会跟原始设定值有一点点偏差。这个偏差通常在几十皮米量级对于普通传感应用无所谓但对于高精度波长锁定或者窄线宽滤波场景就不能忽略。解决办法是迭代计算先假设一个中心波长 λ₀用仿真或经验公式求出对应 n_eff(λ₀)再用 λ 2n_effΛ 算出新的 λ₁把 λ₁ 带回去重新求 n_eff如此迭代两三次就能收敛。这个流程在写 FBG 设计程序时一定要内置否则批量设计不同波长光栅时会碰到系统偏差。4.4 折射率调制形状对谱线形态的影响布拉格公式只告诉我们“中心波长在哪”但反射峰的旁瓣、带宽、滚降特性取决于折射率调制沿光纤轴向的包络形状。如果光栅区域内的折射率调制幅度是均匀的即矩形包络那么反射谱会有明显的高斯旁瓣结构中心主峰两侧有一系列次级峰。这是因为矩形函数傅里叶变换对应的是 sinc 函数。这些旁瓣在工程上很讨厌会造成串扰和波长解算误判。为了压制旁瓣工程上常用“切趾”apodization技术让折射率调制幅度从光栅两端向中间逐渐增大形成类似高斯或升余弦的包络。切趾后旁瓣被大幅压低但主峰也会略宽、峰值反射率略降。这是一个典型的“灵敏度”和“选择性”之间的权衡。我在实际制作传感器阵列时通常要求光栅反射谱旁瓣低于主峰 15~20 dB否则同一个波长附近两个反射峰很容易被误判成同一个。切趾不是可选项而是必需的工艺闭环。4.5 光栅反射率和带宽的经验公式速查以下这套近似公式在项目预研阶段非常有用虽然精度不如数值仿真但足够判断方案可行性反射率 R_maxR_max tanh²(κL)其中 κ πδn_eff / λ。反射带宽 Δλ半高全宽弱光栅近似Δλ ≈ λ_B² / (π · n_eff · L) · √( (κL)² π² )当光栅很长且弱时Δλ 随 L 增长而减小。这个关系直接决定了传感器的最小可分辨波长变化。我把 1550 nm 波段几组典型参数算给大家看如果 L10 mm、δn_eff1×10⁻⁴λ_B1550 nmn_eff1.447则 κ≈π×1e-4/1550e-9 ≈ 203 m⁻¹κL≈2.03R_max≈tanh²(2.03)≈0.823 dB 带宽大约在几百 pm 量级。这个配置做温度传感配合高分辨率解调仪温度分辨率能做到 0.1°C 以下。如果 L 缩短到 1 mmδn_eff 不变那么 κL≈0.2R_max≈0.04反射率太低传感信号弱解调就难。所以做短光栅必须配合大 δn这受限于光纤的光敏性和写入工艺。5. 推导完了回到实际温度应变传感和滤波设计里怎么用它5.1 温度传感器设计波长-温度标定必须做且不能只看斜率基于布拉格公式的温度传感思路很简单先用温控箱做逐点标定得到波长-温度曲线再在实际测量中反推温度。但大多数初学者会犯一个错误只用两个点比如 20°C 和 80°C算出一条直线斜率然后线性外推。实际上 FBG 的温度响应存在轻微非线性尤其在宽温区-40°C 到 150°C使用时会明显。石英光纤的热光系数本身随温度变化所以标定至少需要 5~7 个温度点用二次多项式拟合。温度范围越大非线性修正越重要。我在做海洋温度链的时候用水浴锅做了 -2°C 到 40°C 的逐点标定最终拟合残差能控制在 ±0.05°C 以内。如果不做二次项修正残差会扩大到 ±0.2°C 左右对某些海洋观测场景就不可接受了。5.2 应变传感器设计温度补偿是一个非常麻烦的交叉问题布拉格公式同样表明应变和温度同时影响 λ_B。如果我们想在应变测量中剔除温度干扰一般有两条路线双光栅补偿法在同一个位置并排放两个光栅一个感受应变温度另一个只感受温度不受力。从两个波长漂移量中相减可以分离应变和温度。这个方法简单但“不受力”的光栅封装工艺很考验人。参考光栅法在传感器附近放一个温度参考 FBG认为它和应变光栅处在相同温度场中。这个方法要求布点时温度梯度尽可能小否则补偿不干净。如果你要做高精度分布式应变测量还需要考虑光栅的应变传递效率。光纤粘贴在结构件表面时胶层、基体材料都会分担或衰减应变实测灵敏度往往低于理论裸光纤值。这个“传递系数”需要在实验室里用等强度梁标定出来不能直接套理论 1.2 pm/με。5.3 在光纤激光器里当波长选择反射镜除了传感FBG 最常见的应用是光纤激光器的腔镜。比如线形腔光纤激光器两端各放一个 FBG它们的反射波长决定了激光输出波长。用布拉格公式可以预先算出所需周期和有效折射率。但激光器应用对光栅谱线要求更高不仅要中心波长准还要求反射谱带宽窄、旁瓣低、反射率高。否则激光器可能多纵模起振或者由于旁瓣反射造成额外损耗。此时在设计阶段就要用切趾手段压制旁瓣并结合光栅长度分析纵模间隔与带宽是否匹配。5.4 解调精度和波长分辨率的工程推算FBG 传感不像普通应变片可以直接读电阻值它需要把波长漂移量测准。市面上的商用解调仪分辨率通常在 1~5 pm 左右高端产品能做到 0.1 pm。这个指标直接决定了传感灵敏度下限。以温度传感为例如果解调分辨率是 1 pm温度灵敏度是 10.5 pm/°C那么理论上最小可分辨温度约为 0.1°C。如果你需要分辨 0.01°C解调分辨率必须到 0.1 pm 量级这时候已经不是光栅本身的问题了而是整套光学系统的信噪比、光谱仪热漂移的问题。做系统设计时这些都要在方案阶段算清楚别等光栅刻完了再发现分辨率不够。5.5 关于封装布拉格公式里没有告诉你的那些事理论上光栅裸光纤的灵敏度就是前面推的那些数值但实际工程里没人会拿裸光纤直接扔在现场用。封装会改变应变传递、温度传递甚至改变有效折射率。举个例子把 FBG 贴在不锈钢基片上做温度传感器基片的热膨胀会和光纤互相拉扯。如果基片热膨胀系数比石英大那传感器测到的“温度”其实还叠加了热应力造成的应变响应。这会让温度灵敏度偏离裸光纤的理论值。所以在封装设计时最好选择热膨胀系数接近石英的材料或者用悬空结构让光纤能够自由伸缩。如果你要做高温 FBG 传感器300°C 以上普通掺锗光纤的光栅会被高温擦除需要改用 regenerated FBG 或使用耐高温涂覆层和特殊光纤。这些限制不在布拉格公式里体现但决定了公式在实际项目中能不能落地。6. 我踩过的几个坑写给正在复现推导和做实验的你光看公式推导永远感觉什么都会上手一做全是问题。我把这些年做 FBG 实验时最常见的几个坑列出来希望能帮你避开。第一写光栅前一定要确认光纤的截止波长和模场直径。很多 d 同学拿一根多模光纤或者特种光纤就当单模用刻完光栅以后反射峰又宽又乱因为高阶模也参与了布拉格反射。只有单模光纤在目标波段工作在单模状态下才能用简化的 2n_effΛ 公式。第二光谱仪测得的反射峰中心不等于布拉格波长。如果你用的是宽带光源光谱仪方案反射峰中心会受到光源光谱形状、系统损耗波长特性的影响产生几百皮米的偏移。真要精确标定布拉格波长需要用可调谐激光器扫描或者做参考光栅校准。第三光栅在高温下会“褪色”或中心波长漂移。普通 FBG 在 200°C 以下长期稳定但超过这个温度光栅的折射率调制会逐渐衰减导致反射率下降和中心波长轻微移动。如果你的应用场景是发动机监测、油井高温环境一定提前确认光栅类型和退火工艺。第四不要忽视光纤弯曲带来的有效折射率变化。光纤在铺装过程中一旦弯得太厉害会影响模式传播状态导致 n_eff 略微改变进而影响布拉格波长。虽然一般弯曲半径大于几厘米时影响可以忽略但紧密缠绕时可能产生几十 pm 的波长偏移对高精度测量很致命。最后我再分享一个我常用的验证方法拿到一支新刻好的 FBG不要先急着接入系统而是先用可调谐激光器光功率计直接扫出透射谱或反射谱记录中心波长和带宽。然后用手捏住光栅体温加热观察波长是否向长波方向漂移。这一步能快速验证光栅是否真的写入了、是否正常工作。如果捏住光栅波长纹丝不动多半是光栅区域根本没刻在有效位置或者头尾搞反了。布拉格波长推导看起来就一个公式但它连接着波导光学的模式理论、材料的热光弹光特性、还有实际系统设计的精度预算。把这个公式彻底吃透后续读 FBG 相关的论文、做传感器标定、设计光栅阵列你会发现自己有了一个非常坚固的起点。