夜雨聆风
题目：Decoupling the Isotopic Source Code of Complex Mixtures: A Topological Trajectory Strategy for Dynamic Forensic Sourcing
摘 要
在纵火案件调查中，对汽油助燃剂进行溯源代表着一个关键但极为棘手的分析难题，因为其固有的地质特征与炼油工艺产物以及环境降解噪声不可避免地混杂在一起。传统的静态指纹图谱技术缺乏可行的策略来解析这些复杂动态混合物中的混淆因素。为了弥补这一空白，我们提出了一种分层化学计量学框架，将范式从搜索“不变标记物”转变为通过数学方法将地质基线特征与工艺过程引起的分馏效应进行解耦。通过集成嵌套方差分量分析与多变量判别轨迹分析，我们正交分解了混合同位素方差，从而分离出一组稳健的标记物面板，其中包括三个具有“抗瑞利分馏”特性的区域溯源锚点，以及一个记录工艺过程信息的标记物。通过动力学压力测试，我们量化了精确的“法证分析有效时间窗口”，结果表明：区域地质特征能够在超过48小时的风化作用下持续存在；而工艺过程的特定细节则是瞬态的，若要在燃烧条件下捕捉这些信息，需要在燃烧后2分钟内进行即时取样。至关重要的是，多变量判别轨迹分析揭示了一种“拓扑记忆效应”：尽管热力学分馏导致了显著的同位素漂移，但不同来源的演化轨迹仍然保持着严格的拓扑分离性。这项工作为动态溯源提供了一种通用的数学策略。此外，我们引入了一个开源的计算工作流程，确保了分析的透明性和可重复性。这使得法证分析人员能够在备受关注的野火或纵火案件中区分助燃剂来源，而无需进行大量的再培训或具备专业的编程知识，从而弥合了高级化学计量学与一线执法实践之间的差距。
研究背景
确定复杂化学混合物的来源仍是现代分析化学中最棘手的挑战之一，尤其是当目标分析物在开放体系中经历快速的热力学演化时。在火灾调查中，示踪汽油等可燃液体残留物（ILRs，如汽油）正是这一挑战的典型例证。其固有的地质特征与人为炼制产物以及随机环境降解不可避免地混杂在一起。尽管基于气相色谱-质谱联用的模式识别是ILR分类的基石，但这些静态指纹图谱方法本质上依赖于质谱稳定性的假设。在实际场景中，风化效应会极大地扭曲色谱轮廓，导致传统匹配算法失效。这一方法论缺口凸显了对动态化学计量学策略的迫切需求。
追求更高分辨率的溯源方法推动了先进分析仪器与精密计算智能的整合。最先进的方法利用多维色谱最大化复杂混合物的特征回收率，随后通过涉及精细特征选择和集成建模的机器学习流程进行处理。例如，最近有研究表明，递归特征添加与决策树分类器相结合，可以通过捕捉简单模型遗漏的互补特征交互来提高分类准确性。类似地，有研究进行同位素分析，通过优化特征集和融合多个分类器，在区域汽油溯源方面取得了显著成功。
然而，这些先进工作流程主要设计用于对静态或准静态样本进行分类。其成功依赖于在原始化学或同位素特征中最大化判别力，其隐含假设是风化和燃烧仅仅是需要最小化的噪声因素，而非需要建模的动态过程。关于关键燃料组分在真实火灾非平衡条件下的动力学行为，仍然存在深刻的机理空白。这就留下了一个未回答的关键问题：这些经过精心识别的静态样本在法证分析中能保持多长的有效时间？缺乏对这种时间有效性窗口的定量框架，严重限制了现有模型在真实火灾现场中的直接应用，因为在这些现场中，样本总是在经历复杂热力学转化后才被采集。
本研究建立分层化学计量学框架，将溯源从静态匹配转向动态轨迹分析。通过整合嵌套方差分量分析与多变量判别轨迹分析，将混淆的同位素变异性分解为正交分量，分离出双重功能标记物面板：抗瑞利分馏的区域锚点与记录热力学历史的工艺敏感变量。与传统方法不同，本研究利用同位素漂移模拟燃料动力学演化。关键发现是判别空间中拓扑关系的持久性：尽管绝对同位素值发生显著漂移，不同来源间的相对几何方向始终保持不变，为溯源提供稳健机制。此外，首次阐明了烷烃与芳香烃在竞争降解下的化合物特异性同位素分馏机理，通过动力学压力测试建立了分析有效性窗口。该策略将瞬态化学位移转化为矢量轨迹，为非平衡体系复杂样本解混提供了通用方法。
主要发现
1、采样区域
总共采集了60份商业汽油样本，这些样本分布在中国四大主要地理区域（西部、南部、东北部和中部地区）的零售加油站（图1A），在燃烧与风化实验中，从国家样本库中进一步选取了八份汽油样品作为子集（图1B）。
图1 汽油采样点分布图。(A)中国各地采样厂商的地理分布情况，按西部、中部、东北部和南部地区进行分类。(B)北京用于燃烧与风化实验的具体采样点及其对应的汽油等级。
2、层次化特征发现策略
将色谱保真度与嵌套方差解耦技术相结合，从复杂石油烃基质中提取化学计量学上可靠的源追踪标记物需采用系统化方法。首先，通过严格的“同位素保真度纳入标准”筛选出18种具有高色谱纯度的关键诊断化合物。初步的无监督主成分分析表明，仅凭高质量色谱数据不足以有效区分样品，说明地质特征与加工特征存在混淆。为解决此问题，我们应用了嵌套方差成分分析方法（VCA），通过构建线性混合模型（LMM）实现了变异来源的正交分解（图2A）。如图2B所示，方差分解揭示了这18个特征中存在不同的同位素行为特征。通过应用50%贡献度阈值，我们成功识别出13个区域标记物（主要反映地质来源）和3个加工标记物（主要反映炼制分馏过程）。解耦分析表明，部分化合物（如F3）对炼制过程具有抗性，能保留原油特征指纹；而其他化合物则对加工条件高度敏感。
图2 稳定同位素特征的方差解耦与筛选策略。(A)嵌套方差分解的计算框架。 (B)特征化合物的方差成分指纹图谱。
3、示踪剂判别效能的稳健性评估与正交验证
尽管基于VCA的初步筛选将候选区域标记物范围缩小至13个、过程标记物3个，但关键的分析挑战依然存在——即验证地理同位素特征能否经受住不同精炼工艺带来的分馏效应。为解决这一问题，采用了正交偏最小二乘判别分析（OPLS-DA）。与标准主成分分析不同，OPLS -DA提供了一种监督回归机制，能够区分与类别标签（地理来源）相关的变异与正交变异（无关的过程噪声），从而为评估候选示踪剂的判别稳健性提供了严谨的过滤方法。
4、通过VIP过滤技术识别过程不变的区域特征
汽油同位素的“混合信号”特性是源追踪的核心难点，原油固有的地理特征常与精炼装置引入的动力学同位素效应（KIE）相互交织。为检验13个候选标记物在不同制造商工艺扰动下的稳定性，构建了制造商分层的OPLS-DA模型。各分析物的判别能力通过投影变量重要性（VIP）评分进行量化，具体如图3A所示。采用VIP值大于1.0的严格标准，从候选化合物库中筛选出四种“核心区域示踪剂”：3-甲基己烷、2,3-二甲基戊烷、1-甲基-2-丙基苯和2-甲基己烷。较高的VIP值表明，对于这些特定同系物而言，由地质来源引起的区域间变异显著主导了由炼油异质性导致的区域内变异。这说明这些化合物具有“过程无关性”特征，这对可靠的法证溯源至关重要。图3B-E进一步展示了这些优先标记物的同位素分布情况。尽管每个数据框内存在可观测的离散度（反映了不同制造商炼油参数预期存在的波动），但四个区域间的同位素基线在统计学上呈现显著分离（p<0.05）。这种对工艺干扰的高度耐受性证实，所选标记物可作为稳健的区域指纹特征，即使在多来源炼油环境中也能保持其诊断性同位素比值。
图3 通过正交信号校正评估标记稳健性。(A)基于13个候选区域标记的OPLS-DA模型得出的VIP分析结果。红色虚线表示VIP = 1.0。(B-E)箱线图展示了四种已识别核心示踪剂在不同地理区域的δ13C分布情况。
5、高分辨率制造商鉴别：揭示工艺诱导的同位素指纹特征
除了地理来源追溯，可靠的法证方法还需区分同一地区内的不同制造商，即足够的微观分辨率，这构成了重大分析挑战。为此，应用OPLS-DA策略，对VCA标记出的三个潜在工艺标志物进行分析。如图4A所示，VIP分析揭示候选标志物存在明显层级差异：仅特征12的VIP得分>1.0，具有显著区分能力；特征1和14虽在VCA中显示差异，但未对监督分类模型产生显著贡献，其变异可能与其他特征存在共线性或被区域信号稀释。特征12的诊断能力在图4B-E中得到进一步阐明：虽然区域基线决定整体同位素趋势，但特征12在同一区域内不同制造商间表现出显著偏差。例如，中部区域（图4E）中，制造商E和F呈现明显偏离区域均值的独特聚类分布。这证实特征12的碳同位素比值主要受工艺诱导分馏作用（由反应器温度、压力及催化剂类型驱动）而非原料本身主导。因此，特征12可作为“工艺记录器”，实现制造商层面的高分辨率溯源。
图4 通过OPLS-DA技术鉴定过程特异性标志物。(A)三个候选过程标志物的VIP评分显示。插图展示了特征12的分子结构。(B-E)特征12在各区域内的同位素异质性分布。
6、同位素指纹的时空有效性：动力学约束与法证采样窗口
为将实验室鉴定的标记物转化为可操作的法证证据，必须严格界定其在真实犯罪现场场景下的存留能力。同位素溯源面临的主要威胁是由相变（蒸发）或化学转化（燃烧）引起的瑞利分馏效应。为了确定区域示踪剂的热力学稳定性边界，并量化特定过程信息随时间的流失程度，从而构建精确的采样窗口，通过模拟两种极端环境胁迫因素（快速氧化（燃烧）和缓慢挥发（风化）），对候选标记物进行了压力测试。
7、基于瑞利分馏抗性的区域标记物热力学稳健性筛选
理想的区域溯源标记物必须具备“保守性”，即在显著质量损失事件中仍能保持同位素特征不变。为此，我们通过分析四种候选标记物在不同环境胁迫下δ13C演化轨迹的变异系数，评估其对瑞利分馏的抗性。图5A表明，候选标记物稳定性与分子结构密切相关：支链烷烃F2、F3和F4对瑞利分馏表现出卓越抗性，构成稳健的“瑞利抗性组”。这三种化合物在燃烧前8分钟及风化48小时内形成明显“同位素平台期”，期间未检测到显著动力学同位素效应，且变异系数始终低于0.05的严格阈值，证实了其作为核心区域指纹的可靠性。与之相反，芳香族化合物F16表现出截然不同的动力学行为：环境胁迫开始时出现显著正向同位素偏移，变异系数迅速突破安全限值。这种不稳定性可归因于其苄位烷基侧链易发生优先氧化反应，以及高沸点导致的非线性蒸馏效应。基于这些标准，F16被归类为环境敏感变量予以排除。最终确定的核心区域标记物组合包括F2、F3和F4。
图5 动力学稳定性评估及法证溯源窗口的界定。(A)热力学稳健性筛选：区域候选物在燃烧（A部分）与风化（B部分）应力下的δ13C值演变及其变异系数。红色虚线框标示“同位素平台期”，其中F2、F3和F4可作为保守示踪剂。(B)判别力衰减：采用时间分解OPLS-DA模型对过程标志物（乙苯）进行评估。评分图（上图）展示制造商聚类的逐渐消失过程，条形图（下图）量化了有效的预测区间（绿色区域）。
8、基于滑动时间窗口策略的过程区分效能定量评估
与区域标记物要求绝对稳定性不同，工艺标记物（F12，乙苯）的价值在于其在不同制造商间相对分布模式的唯一性，核心挑战是量化该模式在环境噪声掩盖前的稳定时长。为此，我们开发了滑动时间窗口OPLS-DA策略以量化特征12的证据存续期。如图5B所示，基于VIP和Q2指标的模型区分效能演化轨迹表明：不同环境压力下信号衰减模式存在显著差异。燃烧过程中，剧烈的热物理变化导致信号快速衰减，特征12仅在2分钟的极短时间内仍能保持对制造商的显著区分能力（VIP > 1且Q2 > 0.5）。相比之下，在温和的扩散控制性风化条件下，该标记物的工艺特异性同位素分辨率展现出卓越耐久性，有效识别期可达120小时。这些发现确立了一个明确的“法证有效性域”：对于涉及高温暴露的火灾现场残留物，必须在燃烧初期立即进行抢救性采样以捕捉瞬息万变的工艺细节；而对于以蒸发为主的环境泄漏事件，该方法具有更宽的时间容忍度。这一区别为现场证据收集提供了关键科学指导。
9、通过多元判别轨迹分析（MDTA）解析复杂混合物：一种用于动态取证的拓扑学方法
为突破单变量分析在混沌环境中的局限性，我们开发了多元判别轨迹分析（MDTA）框架。该方法将三个区域锚点（F2、F3、F4）与工艺示踪剂（F12）相结合，构建高维“同位素指纹图谱”。MDTA的核心创新在于将分析重点从静态点对点比较转向动态轨迹拓扑结构。图6A展示了计算流程：八家制造商的高维同位素数据通过线性判别分析映射到低维潜在空间，再将同一样本序列的得分关联形成演化轨迹。图6B可视化了极端环境胁迫下的轨迹投影，揭示了显著的拓扑记忆效应。燃烧场景中，所有样本沿LD轴出现显著位移，但八家制造商的演化轨迹仍保持严格的拓扑完整性，表现为平行或发散的路径且无交叉缠绕。理论上，这种拓扑不变性基于热力学分馏速率对源基质不敏感这一假设：尽管绝对同位素值发生变化，制造商质心间的相对距离保持不变。区域稳定性与工艺特异性的协同作用形成了抵御非线性热噪声的缓冲区。在长期风化条件下，轨迹呈现径向演化模式：区域标记物锚定基底同位素变异，工艺标记物驱动制造商间区分能力，正交贡献确保整个时间窗口内实现100%分离。MDTA结果验证了该四标记组合不仅是静态分类器，更是能够通过连贯演化路径保留汽油来源信息的动态记录器，为追踪法证证据提供了数学上稳健的基础。
图6 采用多元判别轨迹分析（MDTA）评估动态溯源性能。(A) MDTA工作流程示意图。该流程旨在评估标记物组合在动态环境胁迫下的综合判别能力。(B)溯源演变轨迹。散点图展示了八家不同制造商（A-H，颜色编码）汽油样本在燃烧（上图）和风化（下图）过程中的LDA投影结果。
10、方法学不确定性评估与适用范围界定
为确保动态追踪模型的法证可接受性，建立了涵盖计量溯源性、统计稳健性及适用范围界定的质量保证框架。首先，采用基于VPDB标度的严格计量学规程验证同位素差异的真实性。分析不确定度（扩展不确定度<0.5‰，k=2）显著小于区域间变异（>1.5‰），提供了足够判别边界以可靠解析源信息。针对样本量较小（n=8）的限制，采用留一法交叉验证和200次迭代置换检验评估过拟合风险；所得指标（Q2 > 0.5）证实分类结果由标记物特有的化学逻辑驱动。最后，尽管概念验证阶段未索引样本可能存在的II类错误，但“拓扑记忆效应”为开放集识别提供了稳健机制：未知样本预期形成独特拓扑聚类而非被误分类，为概率化归因奠定数学基础。然而需特别指出，本文动态轨迹边界由当前试点数据集（n=8）参数化。因此，将该MDTA框架扩展至更多样化的全局精炼工艺或替代原料类型时，可能需要重新校准基线以适应扩大的同位素变异范围，确保统计学稳健性。
11、结论
追溯汽油助燃剂等复杂化学混合物的来源，长期以来一直受到身份识别难题的困扰——天然地质特征与炼油过程产生的伪影及环境降解噪声相互交织。传统的静态指纹识别技术在动态场景中往往失效，导致无法采信的法证证据。本研究通过建立分层化学计量策略突破了这一局限，证明实现可靠的来源溯源并非依赖不可改变的标记物，而是通过数学方法将地质基线与工艺诱导的组分分离。借助VCA和严格的动力学应力测试，我们成功鉴定出包含三个“瑞利抗性”区域锚点（F2、F3、F4）和一个“工艺记录器”（F12）的标记组合。该方法精确界定了“法证有效性窗口”，量化表明：虽然区域特征在风化超过48小时后仍保持稳定，但工艺特异性细节具有瞬时性，因此在燃烧场景中需立即采样（<2分钟）。此外，MDTA的引入揭示了显著的拓扑记忆效应——尽管热力学分馏导致显著的同位素漂移，不同来源的演化轨迹在判别空间中仍保持严格的拓扑分离。这一发现标志着法证化学领域从“静态匹配”向“动态轨迹分析”的范式转变，证明了多标记物组合中的相对几何关系比绝对同位素值更具持久性特征。通过为复杂环境样本的解混提供通用模板，并建立开放集识别的稳健框架，本研究弥合了分析严谨性与法律可采性之间的差距，为下一代概率溯源数据库奠定了基础。
原文链接：
https://doi.org/10.1021/acs.analchem.5c07748
说明：本推送只用于学术交流，如有侵权，请联系删除。
投稿&合作请联系：ecs_pku@163.com
