第4章 · 4.8

4.8 噪声、伪影与污染

Noise, Artefacts, and Contamination
本页范围与学习目标
弱小谱线可能来自随机噪声、有限采样与处理失误,也可能对应残留溶剂、水或真实耦合。本节从扫描平均、时间截断和复数采样解释异常的触发条件,给出相位与基线检查顺序,并以有明确介质和样品条件的杂质数据说明多证据指认。排查采用保存原始数据、提出预测、单独改变参数和复核结果的流程,使可疑区域进入结构解释前有可追溯的判断依据。
先读图解和核心解释;需要细节时再展开。
先看图 · 阅读入口
打开本节交互图

异常信号先找机制,再找修正方法

伪影和外来物信号的排查路径不同,应按机制逐项检查。

异常信号先找机制,再找修正方法 步骤 1 判断重现性 随机噪声还是固定异常 比较重复记录 步骤 2 回到采集处理 检查折叠、截断、相位 定位数据形成的步骤 步骤 3 再查物质来源 空白、溶剂和另制样 检验污染解释 异常信号先找机制,再找修正方法 步骤 1 判断重现性 随机噪声还是固定异常 比较重复记录 步骤 2 回到采集处理 检查折叠、截断、相位 定位数据形成的步骤 步骤 3 再查物质来源 空白、溶剂和另制样 检验污染解释
图 4.8-1 本页概念关系图。箭头表示本页采用的检查顺序。

读图时保留这项条件:真实小峰可能是污染或卫星;固定异常也可能来自处理。

动手看一看 · 教学模型

区分峰本身与相位、基线造成的外观

先观察参考峰,再加入相位或基线变化,比较对称窗口中的面积。

区分峰本身与相位、基线造成的外观吸收色散
接收投影与斜基线
接收投影与斜基线响应-0.5-0.2250.050.3250.6-8-4048相对频率 / Hz
中心响应0.4547
对称窗口面积0.9444
相位偏差0°
0 °
0
区分峰本身与相位、基线造成的外观吸收色散

图 4.8-2 异常由明确的教学参数产生;实际污染、噪声和截断还需按各自机制排查。

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}\):

\[ \mathrm{SNR}=\frac{h}{\sigma},\qquad \sigma_{\mathrm{avg}}=\frac{\sigma_1}{\sqrt{N_{\mathrm S}}},\qquad \frac{\mathrm{SNR}_2}{\mathrm{SNR}_1}=\sqrt{\frac{N_{\mathrm S,2}}{N_{\mathrm S,1}}}. \]
(4.8-1)

这里 \(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\),并采用负号傅里叶核,有限采集给出:

\[ S_T(f)=A_0\int_0^{T_{\mathrm{acq}}} e^{-t/T_2^*}e^{-i2\pi(f-f_0)t}\,\mathrm dt =A_0\frac{1-e^{-[1/T_2^*+i2\pi(f-f_0)]T_{\mathrm{acq}}}} {1/T_2^*+i2\pi(f-f_0)}. \]
(4.8-2)

\(T_{\mathrm{acq}}\) 和 \(T_2^*\) 以秒表示,\(f\) 与 \(f_0\) 以Hz表示。分子中的边界项随频率振荡,其相位重复的频率间隔为 \(1/T_{\mathrm{acq}}\)。当FID末端幅度接近零时,边界项很小。相同采集时间下,窄峰衰减较慢,可能明显振铃,而宽峰已经衰减,不一定出现同样异常(Keeler,2004,第4章)。

图 4.8-3:时间截断与频率振铃,手绘示意而非实验或数值模拟 时间域:只画衰减包络 频率域:强峰周围出现振荡 采集结束 时间 / s 末端尚有幅度,边界突然归零 交替起伏 频率 / Hz 振荡尺度随采集时间改变 延长采集:末端幅度更低 适度衰减窗:旁瓣减弱,谱线展宽 示意;曲线仅表达边界与旁瓣,不用于读取峰位或振铃幅度。
图 4.8-3 截断振铃的机制示意。虚线表示未被记录的衰减延续。振荡与强峰关联,实际形状还取决于衰减、相位及窗函数。
展开说明与推导

检查先回到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}\)。这里使用滤波、抽取后的有效复数采样间隔,不能把厂商记录的单实数点间隔不加核对地代入。

展开说明与推导
\[ e^{i2\pi(f+kf_{\mathrm s})n\Delta t}=e^{i2\pi fn\Delta t},\qquad f_{\mathrm{alias}}=f-kf_{\mathrm s}\in[-f_{\mathrm s}/2,f_{\mathrm s}/2). \]
(4.8-3)

\(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。该峰容易被误认成脂肪氢,尽管其原频率位于低场。

图 4.8-4:复数采样的整谱宽回卷及扩大谱宽检查 教学设定:400 MHz,中心 4 ppm,频率坐标向右增加 减去 4000 Hz,即 10 ppm −1 1 4 9 11 ppm 窄谱:记录区间 −1 至 9 ppm 折叠峰 真实频率 扩大至 8000 Hz:11 ppm 被纳入记录区间 示意;下图只对照所需频率,不按新谱宽比例重新标尺。
图 4.8-4 折叠位置的可核验预测。上图按标称400 MHz近似换算Hz;下图表示扩大谱宽后的纳入关系,并非实测谱或完整新频率轴。
展开说明与推导

排查应核对谱宽、中心和预期位移范围,再扩大谱宽重新采集并正确定标。上述例子扩大至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)}\) 后取实部,得到:

展开说明与推导
\[ R(f)=A(f)\cos\phi(f)-D(f)\sin\phi(f),\qquad \phi(f)=\phi_0+\phi_1\frac{f-f_{\mathrm p}}{\mathrm{SW}}. \]
(4.8-4)

\(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及实验部分)。

表 4.8-1 CD₂Cl₂中的少量参考信号,条件见正文
候选来源参考位移 / 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)
样品中非目标少量组分,其指认需要介质、峰组和来源证据。

参考文献与延伸阅读

  1. Hornak JP. The Basics of NMR. Rochester Institute of Technology,在线教材,查阅于2026-10-11。第8章:Practical Considerations;第7章:NMR Hardware。用于SNR定义、扫描平均与接收滤波的基本说明。
  2. Traficante DD. 1991. Time averaging. Does the noise really average toward zero? Concepts in Magnetic Resonance, 3: 83–87. 原始论文与DOI。用于扫描平均的统计解释。
  3. Keeler J. 2004. Understanding NMR Spectroscopy,公开课程讲义,第4章:Fourier transformation and data processing. University of Cambridge. 作者公开讲义;课程版本说明。用于相位、窗函数与截断机制。
  4. Claridge T. 2002. Processing Techniques. University of Oxford NMR Facility,署名版权年份为1998、2002。官方处理讲义。用于加窗与线宽、噪声的取舍;当地软件命令不作为跨平台默认值。
  5. 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;大学托管原文;作者机构托管补充材料。仅引用本节必要的数据及条件,未转载整张原表。