异常信号先找机制,再找修正方法
伪影和外来物信号的排查路径不同,应按机制逐项检查。
读图时保留这项条件:真实小峰可能是污染或卫星;固定异常也可能来自处理。
区分峰本身与相位、基线造成的外观
先观察参考峰,再加入相位或基线变化,比较对称窗口中的面积。
4.8.1 随机噪声与可重复异常
本节以第3章的裂分知识及本章的采样、相位和积分为先修。学习任务是在一张问题谱上判断应先重处理、重采集还是检查样品。将峰表用于结构推断前,先保留原始FID和处理参数;仅有一幅压缩图片,通常无法检验异常的来源。
展开说明与推导
随机噪声(random noise)表现为不能逐点重复的强度起伏;伪影(artefact)是采集或处理产生、不能按通常样品共振解释的谱图特征。污染物产生真实核磁共振信号。真实多重峰、同位素卫星和强耦合线形也有物理来源,不能因其弱小或难以指认就归入伪影。异常是否重复、是否依附于强峰,以及怎样响应参数改变,比外观是否规则更有诊断价值。
本节采用峰高型信噪比(signal-to-noise ratio, SNR):正确相位谱中,扣除局部基线后的指定峰高 \(h\),除以无信号区噪声强度的标准差 \(\sigma\)。二者须取自同一谱、同一强度单位。峰峰噪声、均方根噪声与软件特定的SNR输出不能未经换算直接比较。
展开说明与推导
选择噪声区时要避开峰尾、振铃和未识别的小峰,并记录区间与点数。平坦且零均值时,噪声均方根与标准差接近;若基线仍有斜坡,直接求均方根会把基线偏差算成噪声。不同窗函数改变峰高、线宽和噪声相关性,因此不能以两张任意处理谱的峰高型SNR评价采集灵敏度(Hornak,《The Basics of NMR》第8章)。
扫描平均的条件可以用一个简单模型说明。设第 \(j\) 次扫描在固定位置的值为 \(x_j=s+\epsilon_j\),信号 \(s\) 稳定,噪声 \(\epsilon_j\) 均值为零、方差相同且扫描间独立。平均后信号仍为 \(s\),噪声方差由方差可加性得到 \(\sigma_1^2/N_{\mathrm S}\):
这里 \(N_{\mathrm S}\) 是扫描数,\(\sigma_1\) 是单次扫描噪声标准差。若保存的是累加信号,信号幅度随扫描数增长,噪声标准差随其平方根增长,所得SNR规律相同。噪声并未在每个数据点精确抵消;有限次数平均仍留下随机起伏(Traficante,1991)。比较时还要固定脉冲、重复时间、接收设置及处理方法,并满足相位循环要求。样品降解、温度漂移或相关干扰会破坏模型。例如,两次独立采集中的弱小突起若不能在同一频率复现,可以暂按噪声处理;一个稳定的窄线则需继续检查。增加扫描数可降低随机不确定度,却不会按上述规律消除与信号同步的振铃,更不能消除真实杂质。
习题 4.8-1 扫描数与测量时间
教学设定:同条件下16次扫描的指定峰高为0.200强度单位,无信号区标准差为0.025。若希望SNR达到24,应使用多少次扫描?若一次扫描周期不变,总采集时间相对原实验增加多少?为什么换用很强的展宽窗后SNR上升不能直接验证预测?
参考解答原SNR为 \(0.200/0.025=8\)。目标提高到3倍,所需扫描数为 \(16(24/8)^2=144\),符合常见相位循环要求时可选144次;否则应向上取允许的扫描数。扫描周期相同且忽略固定准备时间,时间为原来的 \(144/16=9\) 倍。平均表示下预期峰高仍为0.200,噪声标准差降至 \(0.025/3=0.00833\)。强展宽窗同时改变峰高、线宽和噪声,违背固定处理条件,必须用同一窗函数重处理两组FID后比较。该计算还假定样品和信号稳定、噪声独立且同方差。
4.8.2 截断振铃:时间边界留下的谱线
截断振铃(truncation ringing)常在强而窄的峰附近形成交替起伏。触发机制是采集结束时FID仍有可观幅度,记录在时间边界突然变为零。频域结果相当于原谱与矩形时间窗的傅里叶变换卷积;后者具有振荡旁瓣。它可以遮蔽邻近弱峰,也会被寻峰程序列成一组小峰。
展开说明与推导
对一个衰减时间常数为 \(T_2^*\) 的孤立共振,取初始幅度 \(A_0\)、频率偏置 \(f_0\),并采用负号傅里叶核,有限采集给出:
\(T_{\mathrm{acq}}\) 和 \(T_2^*\) 以秒表示,\(f\) 与 \(f_0\) 以Hz表示。分子中的边界项随频率振荡,其相位重复的频率间隔为 \(1/T_{\mathrm{acq}}\)。当FID末端幅度接近零时,边界项很小。相同采集时间下,窄峰衰减较慢,可能明显振铃,而宽峰已经衰减,不一定出现同样异常(Keeler,2004,第4章)。
展开说明与推导
检查先回到FID,放大末端并与噪声水平比较;再记录强峰两侧起伏的Hz间隔。若有条件重采,延长采集时间并保持其他条件可比,预测起伏间距缩小、边界影响减弱。同一FID上增加零填充只使曲线显示更密,不补回被截掉的信号。增加扫描数甚至可能让振铃更清楚。
复制处理结果,在副本上施加适度衰减窗,比较邻近弱峰及积分。指数窗 \(w(t)=e^{-\pi\mathrm{LB}t}\) 中 \(\mathrm{LB}\) 以Hz表示;对理想指数FID,它增加相同数值的洛伦兹线宽。窗函数可减轻末端不连续,但会改变分辨率及SNR,应保留未加窗版本(Claridge,2002)。
展开说明与推导
不能只凭窗函数使一条弱线消失就判其为伪影:真实窄裂分也会被展宽合并。较有力的证据是起伏间距随采集时间改变,并同时得到FID末端、强峰关联及重处理响应的支持。
习题 4.8-2 核验截断程度与处置代价
教学设定:孤立峰的 \(T_2^*=2.0\,\mathrm s\),采集时间为0.50 s。求末端包络相对初始值及振荡相位重复间隔。若采集时间增至1.00 s,二者怎样变化?若仅对原FID采用 \(\mathrm{LB}=2.0\,\mathrm{Hz}\) 的指数窗,末端相对幅度是多少?
参考解答原末端比例为 \(e^{-0.50/2.0}=0.7788\),振荡间隔为 \(1/0.50=2.0\,\mathrm{Hz}\)。延长后比例为 \(e^{-1.00/2.0}=0.6065\),间隔为1.0 Hz,仍有明显截断。原FID加窗后末端比例为 \(e^{-0.25}e^{-\pi\times2.0\times0.50}=0.03366\)。它明显下降,但未必低于实际噪声,须看测量数据。理想洛伦兹模型下原线宽为 \(1/(\pi T_2^*)=0.1592\,\mathrm{Hz}\),加窗后约2.1592 Hz,因此与主峰相距2 Hz的真实弱裂分可能难以分辨。以上是设定参数的解析计算,没有把示意曲线当作实测谱。
4.8.3 频率折叠:采样不能区分的频率
频率折叠(aliasing)发生在超出记录频率区间的分量仍进入采样数据时。对有效采样间隔 \(\Delta t\) 的复数正交数据,定义 \(f_{\mathrm s}=1/\Delta t\),可选记录区间为 \([-f_{\mathrm s}/2,f_{\mathrm s}/2)\),其全宽为 \(f_{\mathrm s}\)。这里使用滤波、抽取后的有效复数采样间隔,不能把厂商记录的单实数点间隔不加核对地代入。
展开说明与推导
\(n\) 为整数采样序号,\(k\) 为使频率落入所选区间的整数。相差整倍采样频率的分量具有相同离散序列,因此处理程序无法仅靠这份数据决定真实频率。复数正交采样保留正负偏置;本例采用整谱宽回卷,不沿用单通道实数采样的镜像图像。实际接收滤波会衰减区间外信号,所以也可能表现为缺峰(Hornak,第7章)。
教学设定:标称400 MHz仪器,中心为4 ppm,全谱宽4000 Hz,约为10 ppm,记录范围为−1至9 ppm。一个11 ppm共振相对中心的偏置约为 \(7\times400=2800\,\mathrm{Hz}\)。它若未被滤除,会回卷至 \(2800-4000=-1200\,\mathrm{Hz}\),显示在1 ppm。该峰容易被误认成脂肪氢,尽管其原频率位于低场。
展开说明与推导
排查应核对谱宽、中心和预期位移范围,再扩大谱宽重新采集并正确定标。上述例子扩大至8000 Hz、仍以4 ppm为中心时,范围约为−6至14 ppm,11 ppm共振应回到其真实位置;真实1 ppm共振则仍在1 ppm。改变谱宽会同时改变采样间隔,需调整点数以保持足够采集时间。仅在软件里拉宽显示坐标无效。
没有看到折叠峰也不能证明全谱已被覆盖。带外共振可能被接收滤波或有限激发带宽抑制;样品包含醛、酸性氢或异常位移环境时,要回查采集范围。峰在扩大谱宽后仍不出现,还须检验激发、交换与灵敏度。
4.8.4 相位与基线:先恢复线形,再检查面积
相位误差(phase error)把吸收与色散分量混入显示谱。设正确复数谱为 \(A(f)+iD(f)\),乘以残余相位因子 \(e^{i\phi(f)}\) 后取实部,得到:
展开说明与推导
\(A\) 和 \(D\) 是吸收及色散强度;\(\phi_0\)、\(\phi_1\) 为本式定义的弧度参数,\(f_{\mathrm p}\) 是枢轴频率,\(\mathrm{SW}\) 是全谱宽。软件的相位符号和一阶参数归一化可能不同。常数误差使全谱出现相近的偏斜;频率相关误差则可能使枢轴附近峰正常、远端峰出现相反方向的拖尾。接收延迟、数字滤波处理或脉冲带宽可参与形成频率相关误差(Keeler,2004)。
检查应选两个相隔较远、较孤立的已知吸收峰。在副本上先调零阶相位,再以固定枢轴调一阶相位,来回核验峰两侧。若同一调整同时恢复多个峰的对称性,支持处理相位误差;某条线单独保持负强度则不能强行拉正。带编辑或相干选择的实验可以有设计规定的正负号,常规一维氢谱的检查标准不能无条件移用。
基线偏差(baseline distortion)是无共振区域偏离零强度的缓慢变化或波动。其来源可以涉及FID早期点、接收恢复、强溶剂信号或残余相位;它与随机噪声不同。先修正相位,再检查多个无峰区,避免用弯曲的拟合基线吸收色散尾部。选择基线锚点时应排除宽峰,采用能解释无峰区的低复杂度模型,并对照校正前后谱。
恒定残余基线 \(b\) 在积分频宽 \(W\) 内造成面积偏差 \(bW\)。因此,相同高度的基线偏差在宽积分区影响更大。若更换积分边界时面积持续增长,而新增区间没有可辨识信号,应检查基线和宽峰尾部。对真实宽OH、NH或交换信号,自动校正可能把整个宽峰削掉;此时应保留原谱和不同校正方案,把积分敏感性记入结论。
4.8.5 溶剂、水与污染:真实信号的来源证据
残余溶剂信号(residual solvent signal)来自氘代介质中含氢的同位素体,例如CD₂Cl₂中的CDHCl₂。后处理可以标注或屏蔽其区域,却没有因此去除溶液里的核。痕量杂质(trace impurity)还可来自反应残余、萃取溶剂、油脂、管壁及移液器具。其可重复性与真实样品峰相似,须增加来源证据。
展开说明与推导
下表只摘取与本节案例相关的少数数据。原文采用商业氘代溶剂,298 K、300、500或600 MHz仪器;典型杂质测量向约0.4 mL介质加入3 μL标准混合液,以TMS参考。残余溶剂和水在加入杂质前另行定标。这里明确选用CD₂Cl₂列,不把其他氘代介质中的数值混入(Fulmer et al.,2010,表1及实验部分)。
| 候选来源 | 参考位移 / ppm | 需要一并检查的证据 |
|---|---|---|
| 介质含氢同位素体 | 5.32 | 溶剂空白、核种耦合及定标记录 |
| 水 | 1.52 | 含水量、温度、溶质及交换条件 |
| 丙酮 | 2.12,单峰 | 使用或清洗历史、另制样的响应 |
| 乙酸乙酯 | 1.23,三重峰;2.00,单峰;4.08,四重峰 | 对应积分3∶3∶2、乙基裂分约7 Hz |
展开说明与推导
水峰尤其依赖氢键、温度、含水量及溶质,1.52 ppm仅是上述介质与条件下的参考。观察值接近它只能产生候选。D₂O检验使信号改变或消失,说明其参与交换;样品的OH和NH也会如此,不能据此唯一认定水。加入D₂O还会改变组成、溶解度和交换速率,应在另取样品上验证,并保留原样。
例如,CD₂Cl₂样品在1.5 ppm附近出现宽峰,先用同批溶剂和洁净管采集空白。空白出现相近信号支持水由溶剂或器具带入;空白没有该峰,仍需考虑样品带水。再取小份样品,用适合该样品的干燥介质另制样,并比较目标分子的非交换峰。宽峰减少支持水候选,但样品OH、NH也可能因交换条件改变而移动或展宽,应保留这项歧义。水峰的宽窄本身也不足以确定来源。
教学设定:CD₂Cl₂样品中同时出现1.23 ppm三重峰、2.00 ppm单峰和4.08 ppm四重峰。先确认各区域没有重叠,再检验乙基两组谱线的Hz间距及3∶3∶2积分。若制备曾使用乙酸乙酯,候选得到来源支持;若另取一份样品重新除去挥发性残留后三组信号同步减小,则证据进一步增强。单独匹配2.00 ppm单峰不足以确认乙酸乙酯。
展开说明与推导
可逆检查优先包括溶剂空白、洁净管另制样和重处理。需要除溶剂、加水或加候选杂质时,使用独立小份样品,保持原样用于对照。加标必须观察预期的整组信号响应,还要排除重叠及反应;加标使一个已有峰增高,仍不证明原峰全部来自该杂质。
4.8.6 把排查结果写成可复核的证据
逐项排查应围绕一个可被检验的预测展开。强峰周围的规则起伏先检查FID末端;可疑峰位先检查记录区间;全谱偏斜先检查相位;面积随积分区变动先检查基线。这样可避免在未经修复的峰表上直接数氢。一次只改变一个主要因素,另存结果,并记录在该改变下哪些谱线移动、变宽、消失或保持原状。
展开说明与推导
一张谱可以同时含有数种问题。设一个强溶剂峰伴随振铃,邻近还叠加真实杂质;加窗后振铃减弱,杂质也可能变得不易分辨。此时需要联合FID末端、另一采集时间和杂质其他信号,而不能把整个区域一起删除。
小峰重复出现也仅排除了部分随机噪声解释,同步干扰与真实共振都可能重复。
已在3.8节讨论的同位素卫星,应检查对称位置、与母峰的耦合关系及整体强度,而不按噪声剔除。真实多重峰在改变显示点数后保留其频率关系;强耦合导致的不对称线形还需4.9节的自旋分析。若资料不足,记录“来源未定”及尚需检查的条件,比给每条弱线强行安排一个氢环境更合适。
提交结构证据时,列明采用哪一处理版本、排除区域和排除理由;受溶剂覆盖的信号应记作不可独立判断。保留原始FID、噪声区间、谱宽与中心、采集时间、窗函数、相位及基线参数,以及另制样的组成。后续重处理时才能逐项复核判断。
习题 4.8-3 一张问题谱的逐项排查
教学设定:400 MHz、中心4 ppm、全谱宽4000 Hz的复数采样谱中,1 ppm有一条异常峰;扩大至8000 Hz重采后,该线转至11 ppm,已知3 ppm线保持原位。另一峰更正相位后两侧拖尾减弱;其40 Hz积分区仍残留恒定基线0.020强度单位,真实面积为8.00强度单位·Hz。CD₂Cl₂中还见表4.8-1乙酸乙酯候选的三组信号,积分为6.0∶6.0∶4.0。请分别给出计算、分类、下一项检查及结论限度。
参考解答①11 ppm对中心的偏置约为2800 Hz,窄谱回卷至−1200 Hz,对应1 ppm;扩大谱宽后真实频率被纳入,且3 ppm不动,支持频率折叠。应再核对两次参考与采集参数。仅有最初1 ppm峰本身不能确定真实频率,相差整倍谱宽的其他来源也可回卷。
②多峰同时恢复吸收线形支持相位误差。残余基线面积为 \(0.020\times40=0.800\) 强度单位·Hz,原积分为8.800,相对真实面积偏高10%。相位修正并未自动消除基线,下一步用不含峰的锚点检查低复杂度基线模型,同时确认40 Hz区域没有遗漏真实宽信号。
③候选积分除以2得到3∶3∶2,与表中乙酸乙酯的氢数匹配。应进一步核对三重峰与四重峰的约7 Hz耦合、样品制备历史、溶剂空白和另制样时三组线的同步变化。位移、裂分和积分构成候选证据,尚不足以排除重叠或具有相同片段的其他组分。此处三组线属于可能的真实污染物共振,不能统一归为处理伪影。
关键术语
- 随机噪声 (random noise)
- 重复测量间不能逐点复现的随机强度起伏。
- 伪影 (artefact)
- 由采集或处理造成的异常特征,需与真实样品共振区分。
- 信噪比 (signal-to-noise ratio)
- 本节定义为扣除局部基线的指定峰高与无信号区噪声标准差之比。
- 截断振铃 (truncation ringing)
- FID在尚未充分衰减时终止所引起的频域振荡。
- 频率折叠 (aliasing)
- 不同连续频率产生同一离散采样序列而被记录为区间内频率。
- 相位误差 (phase error)
- 残余相位使显示实部混入吸收与色散分量。
- 基线偏差 (baseline distortion)
- 无共振区域偏离零强度的系统性变化,能够改变积分面积。
- 残余溶剂信号 (residual solvent signal)
- 氘代介质中含氢同位素体产生的真实共振。
- 痕量杂质 (trace impurity)
- 样品中非目标少量组分,其指认需要介质、峰组和来源证据。
参考文献与延伸阅读
- Hornak JP. The Basics of NMR. Rochester Institute of Technology,在线教材,查阅于2026-10-11。第8章:Practical Considerations;第7章:NMR Hardware。用于SNR定义、扫描平均与接收滤波的基本说明。
- Traficante DD. 1991. Time averaging. Does the noise really average toward zero? Concepts in Magnetic Resonance, 3: 83–87. 原始论文与DOI。用于扫描平均的统计解释。
- Keeler J. 2004. Understanding NMR Spectroscopy,公开课程讲义,第4章:Fourier transformation and data processing. University of Cambridge. 作者公开讲义;课程版本说明。用于相位、窗函数与截断机制。
- Claridge T. 2002. Processing Techniques. University of Oxford NMR Facility,署名版权年份为1998、2002。官方处理讲义。用于加窗与线宽、噪声的取舍;当地软件命令不作为跨平台默认值。
- Fulmer GR, Miller AJM, Sherden NH, Gottlieb HE, Nudelman A, Stoltz BM, Bercaw JE, Goldberg KI. 2010. NMR Chemical Shifts of Trace Impurities: Common Laboratory Solvents, Organics, and Gases in Deuterated Solvents Relevant to the Organometallic Chemist. Organometallics, 29: 2176–2179. 原始论文DOI;大学托管原文;作者机构托管补充材料。仅引用本节必要的数据及条件,未转载整张原表。