C11
手机比色法测定天然水体正磷酸盐：色度校正策略在真实基质下的失效边界与偏差预算
1 · 研究问题
在含浊度与腐殖质底色的真实水样中，手机 RGB 比色法配合色卡内标校正测定正磷酸盐（orthophosphate），相对国标钼锑抗分光光度法的偏差能否稳定控制在 ±10% 以内？若不能，该偏差能否由手机自身可测的两个量（样品空白区的明度 L* 与色相角 h）提前预警，从而给出「何时必须离心或稀释」的可操作判据？
2 · 研究背景与空白
正磷酸盐是水体富营养化（eutrophication）的关键控制指标。国标《水质 总磷的测定 钼酸铵分光光度法》（GB 11893-89）的化学基础是：磷酸根与钼酸铵在酸性条件下生成磷钼杂多酸（phosphomolybdic acid），经抗坏血酸还原为钼蓝（molybdenum blue），于 700 nm 测吸光度，线性范围约 0.01–0.6 mg/L（以 P 计）。手机比色的物理本质，是用相机三通道的宽带积分响应替代单波长吸光度：红通道（约 600–700 nm）恰好覆盖钼蓝的吸收带，因此 −log(R/R₀) 可作为 A₇₀₀ 的替代量。但这个替代关系成立有前提——样品中其他吸收与散射在该通道内可忽略。真实水样恰恰不满足这个前提。
已有工作沿"色卡内标 + 色空间变换 + 机器学习"这条主线推进得相当远。2026 年 2 月 Analytica Chimica Acta 的一篇工作用嵌入式色卡参考加集成学习做亚硝酸盐/铵离子的手机比色，把跨设备 R² 从 0.76–0.89 提升到 0.89–0.95、MSE 下降 25–75%；Biosensors and Bioelectronics（2024）系统比较色空间，指出把 gamma 编码的 sRGB 线性化到 CIEXYZ 再转 CIELAB 可显著改善跨平台一致性；PLOS ONE（2020）则论证了受控光箱下"设备无关"（device-independent）校正的可行性。这些工作的共同形态是：自变量是设备与光源，因变量是拟合优度，基质基本是纯水标准溶液。
空白在——没有工作把评估维度换成真实基质下的失效边界：给定基质底色与浊度，各校正策略在哪个点上开始失效，以及能否用手机自身的读数事先判断该失效。这个空白适合一年期中学课题：仪器只需可见分光光度计与手机；参考方法有国标可核；且结论正负都成立——若校正策略在真实基质下依然稳健，是对现有方法适用范围的一次独立扩展检验；若失效，产出的失效判据本身就是可用成果。
3 · 可检验假设
H1 在色卡归一化的 CIELAB 色差（ΔEab）特征下，当样品浊度替代量 A₆₆₀ ≤ 0.05 且基质底色 b ≤ 5 时，手机法与国标法在 0.05–0.60 mg/L(P) 区间的相对偏差 ≤ 10%（Bland–Altman 95% 一致性界限宽度 ≤ 0.12 mg/L）；越过任一阈值后偏差随该量单调上升，在 A₆₆₀ = 0.20 处 > 25%。
H2 抗光源鲁棒性与抗基质鲁棒性不能由同一特征同时提供：HSV 的色相 H 对光源色温变化最不敏感（三种光源间偏差 ≤ 5%）但对基质底色最敏感，色卡归一化 ΔE*ab 则相反。若 H2 被否证（存在同时最优的特征），该特征即为本课题推荐的默认特征，结论同样成立。
4 · 量化验收标准
- 方法学校验（硬门槛）：按 GB 11893-89，用磷酸二氢钾基准试剂配制 0.00 / 0.05 / 0.10 / 0.20 / 0.40 / 0.60 mg/L(P) 系列，分光光度法标准曲线 R² ≥ 0.995；对自来水与河水各做 3 个浓度水平的加标回收，回收率 95–105%、RSD ≤ 5%。此门槛不过，后续全部手机数据无效。
- 重复口径写死：每个浓度点 ≥ 3 次独立配制并独立显色（不是同一比色皿连拍三次），报均值与 RSD；日内 RSD ≤ 5%、日间（≥ 3 个不同日期）RSD ≤ 8%。
- 检出限按 3σ/斜率法计算，σ 取 11 份独立空白的标准差；定量限按 10σ/斜率。分光光度法与手机法分别报告，并明确写「手机法 LOD 高于国标法不是失败，而是本研究要量化的量」。
- 实验设计：3 台不同品牌手机 × 3 种光照（D65 自制灯箱 / 教室荧光灯 / 窗边散射日光）× 4 种基质（去离子水、自来水、河水、腐殖酸+高岭土合成干扰基质，浊度 5 级梯度）全因子，单元内 n = 3。
- 统计口径：方法一致性一律用 Bland–Altman 偏差与 95% 一致性界限，明确写「不用相关系数 r 代替一致性，因为 r 对系统偏差不敏感」；因子效应用三因素方差分析（ANOVA），前提检查预先写死（Shapiro–Wilk 检验残差正态性、Levene 检验方差齐性；任一不满足改用 Kruskal–Wallis 并做 Bonferroni 校正）。失效阈值用分段线性回归（breakpoint regression）估计，断点置信区间由 1000 次自助抽样（bootstrap）给出。
- 可复现性：全部原始照片（含色卡）、ROI 提取与色空间变换脚本（Python + Pillow/NumPy，纯 CPU）、拟合与作图脚本开源，第三方可一键重跑。
5 · 数据与工具
| 用途 | 来源 / 工具 |
|---|---|
| 参考方法 | GB 11893-89《水质 总磷的测定 钼酸铵分光光度法》全文（国家标准全文公开系统 openstd.samr.gov.cn；生态环境部标准库）。仅用于校验与对比，不计入本项目的数据贡献 |
| 显色试剂 | 钼酸铵（AR 25 g，≈60 元）、酒石酸锑钾（AR 5 g，≈40 元）、抗坏血酸（AR 25 g，≈25 元）、浓硫酸（教学装，≈20 元） |
| 基准物 | 磷酸二氢钾（基准试剂 25 g，≈45 元） |
| 干扰基质 | 腐殖酸钠（100 g，≈35 元）、高岭土（500 g，≈20 元），用于配制可控浊度/底色的合成基质 |
| 真实样品 | 校园周边河湖水、自来水、含磷洗衣粉稀释液（自采，记录采样点坐标与时间） |
| 定值仪器 | 可见分光光度计（721/722 型，学校常规配置），700 nm；浊度替代量用同机 660 nm 消光 |
| 成像 | 3 台学生自有手机（记录型号与系统版本，关闭 HDR 与画质增强，锁定 ISO 与曝光，RAW 优先；无 RAW 则记录 JPEG 参数）；24 色标准色卡（国产仿 ColorChecker，≈50 元）；纸箱自制灯箱 + D65 LED（≈60 元） |
| 图像处理 | Python 3 + Pillow + NumPy + scikit-image（sRGB → 线性 RGB → CIEXYZ → CIELAB）。纯 CPU，单张 < 0.1 s |
| 统计 | Python + SciPy + statsmodels（ANOVA、Bland–Altman、分段回归）；自助法自写并开源 |
| 方法对照文献 | Anal. Chim. Acta（2026）参考校正 + 集成学习；Biosens. Bioelectron.（2024）最优色空间比较；PLOS ONE（2020）设备无关手机比色。仅作方法对照与新颖性界定，不计入本项目贡献 |
6 · 方法路径
- 装配试剂并跑通参考方法：按 GB 11893-89 完成完整标准曲线与加标回收，通过验收标准第 1 条后方可进入下一步。
- 建立成像标准作业程序（SOP）：固定拍摄距离、锁定曝光与白平衡、色卡与样品同框、逐张记录光源与时间。ROI（感兴趣区）选取须由脚本按色卡定位自动裁切，不允许人工目视框选，判据须可复现。
- 核实工具能力边界（第 8–10 周完成，不要边做边发现）：确认每台手机能否导出 RAW、能否锁定 ISO/曝光、厂商色彩增强能否关闭。不可关闭的项写入不确定度预算，不剔除该机型。
- 在纯水基质上建立三套校正特征（原始 R 通道消光、HSV 的 H、色卡归一化 ΔE*ab），各自建标准曲线并与分光光度法做 Bland–Altman 对照。
- 加入合成干扰基质做浊度/底色梯度扫描，用分段线性回归定位每套特征的失效断点，检验断点能否由 L*、h 两个手机自测量预测。
- 在真实河水 / 自来水 / 洗衣粉废水盲样上做外部验证：先用判据预测「是否需前处理」，再实测验证判据命中率（报灵敏度与特异度）。
- 独立交叉校验：对同一批盲样用标准加入法（standard addition）在真实基质上独立核定回收率，作为不依赖外部仪器的第二条验证通道。
7 · 新颖性边界
- 本课题不声称提出新的显色化学、新的手机比色装置或新的机器学习模型。钼锑抗法、色卡内标、色空间变换均为已发表方法。
- 已有工作完成到：内标色卡 + 集成学习把亚硝酸盐/铵的跨设备手机比色 R² 提到 0.89–0.95（Anal. Chim. Acta, 2026）；系统比较色空间并推荐线性化后的 CIELAB（Biosens. Bioelectron., 2024）；受控光箱下设备无关校正的可行性（PLOS ONE, 2020）。它们的自变量是设备与光源，因变量是拟合优度，基质基本是纯水标准溶液。
- 本项目的贡献是评价维度的转换，且构成主结论：自变量换成真实基质的浊度与底色，因变量换成相对国标法的系统偏差与一致性界限，交付物是一条可操作的失效判据（用手机自身读数预警何时必须前处理）与一份完整的不确定度预算。
- 价值：手机比色的实际使用场景（野外、课堂、家庭）恰恰是有基质干扰的场景，而现有文献的适用范围声明基本停在受控条件。给出失效边界，比再提高一次 R² 更能决定这类方法能否被真的用起来。
- 若结果显示「各策略在真实基质下均稳健、差异不显著」，这同样是有效结论，但必须用误差棒证明本研究有能力分辨 10% 量级的偏差（由验收标准第 2、3 条保证）。
8 · 决策门槛（go / no-go）
- 第 4 周末：验收标准第 1 条必须通过（R² ≥ 0.995 且加标回收 95–105%）。未通过 → 降级路径 A：把分析物换成亚甲基蓝或市售合成食用色素的水溶液（显色稳定、无显色反应变量），研究问题、假设结构与统计口径全部保留，只是从体系中移除"化学显色不确定度"；主结论（色度校正策略的基质失效边界）不受影响。
- 第 10 周末：确认至少 2 台手机可锁定曝光与白平衡。3 台全部不可控 → 降级路径 B：改为「固定灯箱 + 单机多日重复」设计，把"跨设备"因子换成"跨日期 / 跨批次试剂"因子，全因子设计与断点分析框架不变。
- 第 20 周末：合成干扰基质上必须已观察到至少一套特征出现可定位的失效断点（断点自助置信区间宽度 ≤ 梯度跨度的 40%）。若所有特征都不失效 → 主动加大干扰强度（腐殖酸至 20 mg/L、浊度至 50 NTU），并把「在中学可及浓度范围内未观察到失效」写成定量结论（附最小可分辨偏差）。
- 第 34 周末：真实盲样验证完成，判据灵敏度应 ≥ 0.80。低于 0.80 则不修补判据，如实报告适用范围并做失败案例分析——这个困难本身就是研究内容。
- 第 44 周数据冻结，转入英文报告与英文 PPT；第 48 周完成实验视频剪辑。
- 实验视频要求：固定机位记录三段——(a) 国标法标准曲线制备与比色皿读数；(b) 灯箱内一次完整的手机拍摄 SOP；(c) 真实水样前处理与盲样测定。每段含操作者出镜、试剂标签特写与仪器读数特写，累计 ≥ 8 分钟，原始素材保留至答辩结束。
- 安全与废液：酒石酸锑钾属有毒品，须由教师配制母液并锁柜管理，学生只接触稀释后工作液，全程手套 + 护目镜；浓硫酸稀释遵循「酸入水」，在通风橱或开窗处操作。含钼、锑的显色废液严禁排入下水道，按「含重金属显色废液」单独收集、贴签注明成分与日期，交学校危废渠道。不涉及剧毒品、易制毒/易制爆化学品、氢氟酸、高压反应釜与明火有机溶剂操作。
- 预算裁剪顺序：先砍手机台数（3 → 2，跨设备因子降为两水平，H1 仍可检验）；再砍标准色卡（改自印色块并用分光光度计标定反射率，精度下降但框架不变）；灯箱与显色试剂不可砍。
- 适配学生：适合愿意做大量重复性定量操作、能写基础 Python 图像处理的学生。不适合追求"造一个新器件"的学生。