基于Bloch方程优化的90°脉冲DC极化方法及其在快速高灵敏度Overhauser磁力仪中的应用

《Sensors》：Radio Frequency Interference Suppression for High-Frequency Ocean Remote Sensing Radar with Inter-Pulse Phase Agility Waveform Heng Zhou, Xiongbin Wu, Liang Yu, Fuqi Mo and Xiaoyan Li

本文针对传统质子旋进磁力仪极化时间长、测量效率低的问题，提出一种基于Bloch方程动态分析的渐变式90°脉冲DC极化方法，通过优化脉冲波形指数α(2~30)、强度B0(≥3mT)及脉宽tp(≥3ms)，实现毫秒级质子磁矩高效偏转，使JOM-5SF Overhauser磁力仪在自然与人工磁场环境下均达到0.021nT@1Hz灵敏度，为地球物理勘探与空间磁场探测提供关键技术支撑。

　　

在地球物理勘探、空间环境监测等领域，高精度磁场测量始终是核心需求。传统质子旋进磁力仪虽具备绝对精度优势，却长期受困于“极化效率与速度难以兼得”的矛盾——为获得足够信噪比的Larmor信号，往往需数十秒极化时间，严重制约实时观测能力。更棘手的是，突变式极化脉冲易导致质子磁化矢量绕总磁场作大半径进动，不仅降低有效横向磁化分量，还引发信号幅值随机波动。如何在不牺牲灵敏度的前提下突破极化时间瓶颈，成为仪器研发的关键挑战。

针对这一难题，吉林大学电子科学与工程学院团队在《Sensors》发表研究，创新性地从质子磁矩动力学本质出发，建立基于Bloch方程的磁化演化模型，首次提出“渐变式90°脉冲DC极化”方法，通过精准调控脉冲参数实现质子磁矩高速定向偏转。研究表明，当脉冲波形指数α=3、强度B₀=3mT、脉宽t_p=3ms时，质子磁化矢量可沿总磁场方向完成<90°高效偏转，偏振饱和因子θ降至0.7°，远优于传统方案（θ=81°）。基于此研制的JOM-5SF Overhauser磁力仪，在自然与一级弱磁计量站人工磁场环境中均实现0.021nT@1Hz灵敏度，绝对误差控制在±0.1nT以内，成功将单次测量周期压缩至秒级，为高动态磁场探测提供了全新解决方案。

研究团队首先推导了质子磁化矢量M满足的Bloch方程dM/dt=γ_pB×M（γ_p为质子旋磁比），通过中心差分法将其离散化为迭代计算模型，量化极化脉冲对磁化轨迹的影响；继而定义偏振饱和因子θ（磁化矢量与总磁场夹角）作为评价指标，结合JOM-5SF磁力仪平台（含谐振腔、音频线圈及¹⁵N-D-oxo-TEMPO自由基溶液），采用PWM-DAC混合调制技术生成可编程渐变脉冲；最后通过双机同步比对（自然磁场）、补偿式OPM标准场标定（人工磁场）等多场景验证系统性能，并结合功率谱密度(PSD)分析确认灵敏度极限。

3.1 质子磁化动力学方程建立

基于经典力学与电磁理论，研究团队从单质子磁矩在磁场中的转矩作用出发，推导出宏观磁化矢量M的动力学方程。当忽略弛豫效应时，该方程简化为Bloch方程（式5），揭示磁化矢量运动速率取决于总磁场B（含地磁场B_e与极化场B_p）与旋磁比γ_p的乘积，为后续数值模拟奠定理论基础。

3.2 磁化动力学方程求解与离散化

为实现工程应用，团队引入时间微元Δt近似处理微分项，得到递推关系式（式8）：M(t)≈M(t-2Δt)+2Δtγ_pB(t-Δt)×M(t-Δt)。该模型证明当Δt足够小时，磁化状态仅由前两时刻历史值与瞬时磁场决定，从而将连续动力学问题转化为可编程计算的离散序列。

4.1 90°脉冲波形对极化效率的影响

通过设定极化场B_p(t)=B₀(t/t_p)^α（B₀为峰值强度，t_p为脉宽），研究不同指数α对应的磁化轨迹。仿真显示：α过小（如0.5）时脉冲前沿陡峭，总磁场方向突变过快，磁化矢量无法同步跟踪，最终以44.7°大倾角绕磁场进动（θ=44.7°）；α过大（如>30）时后沿突变同样恶化效率；唯α∈[2,30]时，磁化矢量能以极小半径（θ≈0.005°）紧贴总磁场方向螺旋偏转，实现近完美对齐。

4.2 脉冲强度与持续时间的影响分析

固定α=3、t_p=20ms时，B₀从0.05mT增至5mT过程中，θ变化不显著但最终偏转角β（总磁场与地磁场夹角）持续增大——B₀≥3mT时β≈90°，确保极化结束后磁化矢量垂直地磁场，最大化横向分量。脉宽测试则揭示：t_p<2ms时脉冲突变导致θ急剧上升至27°；t_p≥3ms后θ稳定于<1°，确认最短有效脉宽为3ms。

5.2 不同脉冲参数的实验验证

基于JOM-5SF平台实测显示：当α=0.5时因θ=81°，Larmor信号幅度波动大（1~1.8V）；而α=3时θ=0.7°，信号稳定达2.3V，信噪比提升至44。脉冲强度实验中，B₀=0.1mT时信号幅值较5mT下降0.57V，主因偏转角不足致横向分量减小；脉宽实验中t_p=1ms时信号幅值不稳，与仿真中θ=27°的结果一致，验证理论模型的准确性。

5.3 自然环境下的灵敏度评估

双机同步测量（间距5m）结果显示，夜间低干扰时段磁场差值标准差为0.0217nT，剔除地磁场日变与电磁干扰后的灵敏度达0.021nT，证实该方法在真实环境中的稳定性。

5.4 人工磁场系统的灵敏度标定

在一级弱磁计量站（背景噪声0.005nT）中，JOM-5SF在20~80μT标准场范围内的绝对误差小于±0.1nT，功率谱密度分析进一步确认0.021nT/√Hz的灵敏度量级，与时间域统计结果高度吻合。

本研究通过建立质子磁化动力学的Bloch方程模型，首次系统阐明极化脉冲参数（波形指数α、强度B₀、脉宽t_p）与偏振饱和因子θ的定量关系，揭示“脉冲渐变度匹配磁矩响应速度”是实现高效极化的关键机制。所提出的优化参数组合（α=3、B₀=3mT、t_p=3ms）突破传统方案的时间瓶颈，使Overhauser磁力仪在保持0.021nT高灵敏度的同时，将极化时间压缩至毫秒级，解决了长期存在的“速度-精度”权衡难题。该成果不仅为核旋式磁力仪的工程优化提供普适性设计准则，更为海洋磁测、卫星磁场遥感等对时效性要求极高的应用场景提供核心器件支撑，推动地磁探测技术向高动态、高效率方向跨越发展。