抗漂移人工嗅觉：只用一种材料，靠"结构"区分气体

背景与原理

挥发性有机化合物（VOC）是医学诊断、环境安全与工业监测中的重要"潜在指纹"。传统电子鼻模仿生物嗅觉，依赖化学多样性范式——用多种敏感材料组成阵列来获得选择性，但代价是各传感器基线漂移不同步、校准负担沉重。

本 Protocol 提供另一条路线：时空物理编码。只用一种敏感材料（多级多孔 CuP@G 墨水），靠激光加工的程序化微结构"设计"气体扩散路径，把 VOC 吸附转化为可预测、可重复的时空延迟信号，再用深度学习解码。原文标准化了四个逐步演化的平台：

创新点：定制微腔把层流流场（>0.1 m/s）与边界层扩散动力学调制成可预测的时空延迟信号（空间描述符 Sc、Sp、Sc/Sp）；最后用 ResNet-34 或 MTL-CNN 结合 SHAP 可解释性分析，从原始电阻信号直接输出分类结果，破解电子嗅觉的"黑箱"难题。

实验前准备（Before you begin）

准备一：硬件与计算环境

• CO₂ 激光平台：开机前确认排风系统完全运转（即时抽走汽化产物）；针对目标基底平面校准 Z 轴焦距，保证空间分辨率、防止电极断线。
• 动态配气与采集系统：高纯 N₂ 彻底吹扫所有气路；开启 Keithley 3706A 数字多用表并配置多路电阻采集。
• 计算环境：Python 3.11.5 + PyTorch 2.3.0，配 CUDA 兼容 GPU（CUDA Toolkit 12.5）；安装 SHAP（可解释 AI）与 PCA 相关库。

准备二：前驱体试剂与标准气体

• GO 储备液（5 mg/mL）：125 mg GO 粉末分散于 25 mL 去离子水。
• 0.5 M NaOH：2.0 g NaOH 颗粒溶于 100 mL 去离子水，搅拌至完全溶解。
• PDMS 预聚物：基胶与固化剂按 10:1 质量比混合，剧烈搅拌 5 min。
• VOC 标准气体：鼓泡瓶装入目标 VOC 液体（如 1-己醇），恒温 25°C 产生饱和蒸气；按 Antoine 方程计算饱和蒸气压；经两级动态稀释（MFC，Horiba SEC-Z500）定量至目标浓度（如 10 ppm）。

材料与设备（关键配方）

配方表 1：P/G 墨水水热反应混合液

配方表 2：配体功能化 P/G 墨水（可选，用于 OHSA 扩展）

配方表 3：CuP@G 墨水水热反应混合液

分步方案（Step 1–15，按四大模块归组）

路线选择：做 OHSA 用 P/G 基准墨水（Step 1，可选 Step 2 功能化）；做 CISS/SCM/RHC 必须用多级多孔 CuP@G 墨水（Step 3），以保证优异的气体扩散与抗漂移性能。

模块 A：导电墨水合成

Step 1｜P/G 墨水（PDA@rGO）

1. 溶解：25 mg DA·HCl 溶于 20 mL 去离子水，调 pH 至 8.5；与 5 mL GO 悬浮液（5 mg/mL）混合。
2. 水热：转入聚四氟乙烯内衬反应釜，密封，160°C 加热 12 h；冷却至 25°C（pH 通常降至 ~7.5–8.0）。
3. 洗涤：去离子水与无水乙醇交替离心洗涤数次（10,000 × g，10 min，25°C）；末次上清 pH 应达 ~7.0；收集沉淀。
4. 重分散：沉淀分散于无水乙醇，即得 P/G 墨水。可暂停
5. 表征：SEM/TEM 确认 PDA 纳米颗粒均匀锚定于褶皱 rGO 片上；Raman 检查 I_D/I_G；四探针测各基底（纸、Kapton、PDMS 等）方阻。

Step 2｜配体功能化 P/G 墨水（可选，扩展 OHSA 传感库）

1. P/G 墨水分散于无水乙醇，用碱性乙醇溶液调 pH ~8.5。
2. 加入巯基或 NH₂ 类功能化配体，25°C 轻柔搅拌 12 h（Michael 加成或 Schiff 碱/Michael 加成反应）。
3. 无水乙醇离心洗涤数次，重分散于无水乙醇备用。

Step 3｜多级多孔 CuP@G 墨水（CISS/SCM/RHC 必需）

1. 称量：5 mg DA·HCl + 0.4818 mmol（~82.1 mg）CuCl₂·2H₂O。
2. 引发：加 0.2 mL 0.5 M NaOH，摇匀至初始 pH 4.25；再加 1 mL 5 mg/mL GO 溶液混匀。
3. CRITICAL冰浴超声 15 min（保证 Cu²⁺ 与多巴胺均匀配位）。
4. 水热：转入密封瓶（维持自生蒸气压），65°C 加热 24 h，形成多孔 Cu²⁺-聚多巴胺聚集结构。
5. 洗涤：冷至 25°C（预期 pH ~3.5–4.0）；离心洗涤沉淀多次。CRITICAL见 Troubleshooting 1。
6. 重分散：收集沉淀，加 2.5 mL 无水乙醇，振荡+超声 15 min，得多级多孔 CuP@G 墨水。可暂停
7. 表征：SEM 确认 rGO 上蓬松 PDA 层；EDS 验证 C、N、O、Cu 均匀分布；BET 确认比表面积 ~50 m²/g；XRD 验证非晶 rGO 与混合铜相（CuO/Cu₂O）；XPS 证实 Cu²⁺/Cu⁺ 共存。

模块 B：传感器阵列与微腔加工

Step 4｜OHSA（纸基对称折纸结构）

1. 设计模板：11 条独立列，每列 8 个 2 × 2 cm 方块，中部条带区 0.5 × 1 cm；普通打印机打印并裁成 11 列。
2. 沿方块间折线预折成对称蛇形折纸，再展开铺平。
3. 用改装瓶式 DIY 笔（配 10 μL 移液器头，点样精度 ~0.5–1.0 mm）把 P/G 墨水交替写在展开层的中部条带区，形成六通道传感阵列。CRITICAL保持墨道连续；见 Troubleshooting 2。
4. 25°C 完全干燥。可暂停墨道与铝箔条用导电银胶连接成电极，建立良好欧姆接触。

Step 5｜GLIGIE 阵列与模块化微腔（CISS 用，PI 膜）

1. 裁一块 250 μm 厚 PI 膜；CO₂ 激光刻蚀 2 × 5 矩形凹槽（参数见原文表 1）。CRITICAL控制功率防止过烧。
2. 机械去除凹槽内残留 LIG，去离子水+无水乙醇洗净。
3. 刻蚀叉指电极（3 对指，指宽 200 μm、间距 300 μm）。CRITICAL见 Troubleshooting 3。
4. CO₂ 激光加工 PMMA 微腔；用 0.5 mm 硅橡胶圈将微腔与 GLIGIE 阵列耦合，保证优异气密性。

激光参数摘录（原文表 1，额定 25 W）：PI 膜——绿色线刻蚀两遍（首遍 50%、次遍 30% 功率，150 mm/s）；红线刻蚀 20%、200 mm/s；蓝线切割 70%、150 mm/s。PMMA——红线刻蚀 50%、200 mm/s；蓝线切割 70%、50 mm/s。硅橡胶切割 30%、200 mm/s。

Step 6–7｜LIGPIE 阵列（SCM 用）与 PLIGIE 径向阵列（RHC 用，均为 PF 纸）

1. CO₂ 激光在 6 × 6 cm PF 纸（0.09 mm 厚）上刻蚀电极（参数见原文表 2/3）。CRITICAL见 Troubleshooting 4/5。
2. 6.25 W、100 mm/s 切割电极外缘（LIGPIE 成 2 × 2 阵列；PLIGIE 为径向对称）。
3. PMMA 底板刻圆形凹槽（LIGPIE：3.5 × 3.5 cm 板、半径 9.5 mm；PLIGIE：5 × 5 cm 板、半径 14 mm，均为 6.25 W、50 mm/s）；顶板刻同样凹槽并加硅胶密封圈槽位。
4. 0.5 mm 硅胶片上切割密封圈（2.5 W、100 mm/s）。可暂停

模块 C：组装与封装

1. Step 8｜OHSA：沿折线把已贴电极的纸条折成紧凑 2 × 2 cm 方块；胶带密封侧边，双侧入口型留 2 个未封口、单侧型留 1 个。
2. Step 9｜CISS：每个 GLIGIE 单元滴 2 μL CuP@G 墨水，CRITICAL30 s 内转入 80°C 真空干燥箱（Troubleshooting 6）；再滴 2 μL 不同浓度 PDMS 甲苯溶液构建厚度递增的五个传感器（S1–S5）；80°C 烘烤 3 h。可暂停固化阵列与 PMMA 微腔集成，硅橡胶圈准确入槽，螺钉夹紧成最终芯片。
3. Step 10｜SCM 与 RHC：叉指区滴 2 μL CuP@G 墨水，25°C 真空干燥 30 min；制备 4 片传感器逐层堆叠，每层相对前一层顺时针旋转 90°；堆叠体对准定制 PCB，与四正方向引脚电极连接；密封圈装入底腔凹槽，顶盖（"Gas in"）、PCB 与底座（"Gas out"）螺孔垂直对齐后拧紧。CRITICAL层间同心偏差 <50 μm（Troubleshooting 9）。

模块 D：气敏测试与数据分析

Step 11｜OHSA 测试（含可选多刺激）

1. OHSA 固定于电路板，放入程序控温不锈钢测试腔，接多用表/采集系统。
2. 基线：抽真空至 0.06 Torr，去除吸附物，稳定基线电阻 R₀。
3. 施加刺激：注入指定浓度（如 10 ppm）VOC（如庚醛、二甲苯异构体）或预设质量比混合气；温度 25–51°C（±0.5°C）；光强 6.5–310 mW/cm²；湿度 10%–80%、气压 793–868 Torr（湿度精度 ±2% RH）。CRITICAL多参数须逐一解耦。
4. 记录各层（L1–L6）随时间变化的电阻。

Step 12｜CISS/SCM/RHC 标准化动态测试

1. 连接配气系统（MFC 组）。CRITICAL全气路先用高纯 N₂ 吹扫。
2. 纯 N₂ 吹扫稳定基线 R₀；切换至目标 VOC（如 10 ppm）暴露，记录动态响应；随后切回纯 N₂ 恢复。
3. CRITICAL每个 VOC 样品做 6 次独立"暴露/恢复"循环，保证统计显著性。

Step 13｜OHSA 数据分析

1. 电阻计算：不要用单个数据点代表 R₀ 或 R_end；按 ΔR/R₀(%) 计算每层相对电阻变化；每样品 6 循环，弃去初始平衡循环，以均值 ± SD 报告。
2. 速率常数：归一化后空气暴露段用 one-phase association、真空清洗段用 one-phase decay 拟合提取 K。CRITICAL拟合优度 R² 须达标。
3. PCA 区分响应模式（2D/3D 得分图）；用第一主成分（PC1）对预配混合物做成分溯源验证。

Step 14｜CISS 数据分析

1. 特征提取：ΔR/R₀(%) 全程曲线；t_res = 响应达到平衡值 90% 所需时间；t_rec = N₂ 清洗段信号恢复到最大平衡响应 10% 所需时间。
2. 动力学：one-phase association/decay 拟合提取 K 与理论最大响应"Span"；响应热图（OriginPro/MATLAB/Matplotlib）呈现醇、醛、酮类指纹；PCA 取 PC1/PC2 画散点并加椭圆圈簇。
3. （可选）ResNet-34 深度学习：原始时序信号统一标准化为 ΔR/R₀(%)；随机打乱后按 60:20:20 严格划分训练/验证/测试集；仅对 60% 训练集用自编码器做数据增广，验证与测试集绝不增广；验证集调参后再用训练+验证集重训（Troubleshooting 8）；测试集上同时评估质量（A/B/C）与产地（P1–P5）：混淆矩阵、accuracy、macro-P/R/F1、ROC-AUC、softmax 箱线图。CRITICAL测试集完全隔离。

Step 15｜SCM 与 RHC 数据分析

1. 预处理：导出全部传感层（L1–L4）电阻-时间数据；按 3 min 暴露 + 3 min N₂ 清洗循环取 R₀、R_end，计算 ΔR/R₀(%)；6 次重复取均值 ± SD。
2. RHC 专属：在 t = 180 s（暴露结束）同时提取空间位置信号，量化层内中心（Sc）与外周（Sp）差异及层间等级差异（Sc/Sp 比值）。CRITICAL
3. 动力学与时间分析：提取 t_res/t_rec；浓度梯度曲线 one-phase association 拟合；按 3σ 法从低浓度线性拟合计算 LOD。
4. PCA：行为样品、列为四层响应特征；GraphPad Prism 9.5 出热图（时空指纹），Origin 2024 降维，输出 Sc、Sp、Sc/Sp 组合 2D PCA 图。
5. MTL-CNN 多任务分类（PyTorch 2.3.0 + CUDA 12.5，GTX 1650 即可）：两个共享卷积块（Share Layer I/II）生成紧凑表征；分叉为并行任务分支（SCM：Origin 5 类 + Quality 3 类；RHC：Origin 5 类 + Level 3 类），分支内用额外 Conv2d 作空间滤波器重加权，Softmax 全连接输出。
6. 性能评估：5 折交叉验证；混淆矩阵与 PR 曲线、各类 AP 分数；SHAP 可解释性——Sankey 图量化信息流与各层贡献、弦图刻画传感器关联与桥接作用（Troubleshooting 10）。

关键注意事项速查

Troubleshooting（10 个常见问题）

预期结果

1. 多级多孔传感墨水的表征（验证关卡）

CuP@G 墨水在无水乙醇中呈均匀分散；SEM 显示 rGO 片上附着 PDA 层；EDS 显示 C、N、O、Cu 均匀分布；N₂ 吸附确认多级多孔结构与高比表面积（~50.09 m²/g）；XRD/XPS 证实混合铜相与 Cu²⁺/Cu⁺ 共存。

2. 物理编码架构的验证

OHSA：目检确认蛇形折纸边缘密封严实，建立严格的顺序扩散屏障。CISS：CFD 模拟与瞬态输运仿真验证微腔内层流流场（局部流速 >0.1 m/s）；目标气在传感界面 ~4.5 s 内快速达到平衡、~9.5 s 完成基线吹扫。

3. 气敏性能与时空指纹

标准化测试条件下基线漂移严格 <±0.5%。SCM 传感器对 25 ppm 1-己醇连续 12 循环响应稳定、8 周长程稳定；不同弯曲角度与反复弯折下信号稳定；10 ppm 下 180 天长程稳定（不用时存放于真空干燥器）。OHSA 可将温度、湿度与 VOC 多重刺激解耦为可目视识别的空间指纹。

4. 高分辨率区分与真实世界应用

常规桌面验证的硬指标：提取多维时空参数后，PCA/LDA 必须在 11–12 类主要 VOC 上产生互不重叠的清晰聚类——达标即解锁文献级分辨能力：

• OHSA：区分二甲苯三种异构体及其混合物（10 ppm）；识别 2-辛醇、2-丁醇对映体。
• CISS：22 种 VOC（10 ppm）区分；1-己醇/2-己酮/2-庚酮 5–25 ppm 浓度梯度；每次测试后基线完全恢复至初始 R₀；烟草质量（A/B/C）与产地（P1–P5）ResNet-34 分类。
• SCM：12 种 VOC、动态浓度与二元混合物分辨；MTL-CNN 多任务（Origin/Quality）分类 + SHAP 可解释（Sankey/弦图）。
• RHC：Sc、Sp、Sc/Sp 多描述符浓度依赖 PCA 聚类（5–25 ppm）；Origin/Level 多任务分类，5–25 ppm 内可重复。

局限性

1. 依赖可逆气固相互作用：时空编码完全建立在可逆吸附上；不可逆化学吸附会破坏编码的输运动力学，最优分析物限于物理吸附或弱可逆配位物种。
2. 对反应性无机气体交叉敏感：墨水活性位点对强氧化还原活性物种（如 NH₃）本质上有响应，会干扰 ΔR/R₀ 与 R₀。
3. 动力学延迟：完整时空"指纹"依赖气体逐级穿透结构屏障，不适合毫秒级瞬时过程分析。
4. 制造公差严苛：高精度微结构需严格一致性控制才能保证空间描述符可重复。
5. 算法依赖重：特征提取离不开机器学习，需辅以 SHAP 等可解释 AI 框架做严格的物理化学解释。

资源与数据

通讯作者：Min Zhang（华东师范大学，mzhang@chem.ecnu.edu.cn）；所有独特/稳定试剂可向通讯作者索取。数据与代码：原始数据与标准化基准数据集公开存放于 Zenodo（10.5281/zenodo.21259172）。资助：国家自然科学基金 22274053；上海市科委 24140711700。