4. 稳定性分析的数据采集与实验设计

长期稳定性分析的本质,是用足够长、足够密、足够干净的时间序列去刻画系统输出随时间的慢漂移。它和短时精度测试的最大区别在于:被测对象不再只是仪器本身,还包括环境、供电、时钟、操作流程等一整套“共模条件”。因此实验设计的第一步不是选仪器,而是先明确要回答的问题——是评估某台量子测量系统在恒温实验室内的日漂移,还是评估它在真实机房环境下的月漂移?问题不同,采样频率、时长和环境记录粒度都会完全不同。

一、先确定观测对象与时间尺度

量子测量系统的输出通常可分解为三部分:被测物理量引起的真实变化、系统自身的慢漂移、以及随机噪声。长期监测要捕捉的是第二项,因此必须让第一项尽可能小或可独立记录。以超导量子比特读取链路为例,若关注的是读出谐振频率的长期漂移,就应固定偏置点、固定读取脉冲参数,并在每次采样时记录芯片温度、稀释制冷机各温区读数、磁屏蔽内剩磁和本振参考源状态。若关注的是整机对某标准信号的测量偏差,则需引入一个稳定的参考量(如原子钟秒脉冲或标准电阻)作为“锚点”,否则无法区分是系统漂移还是被测对象变化。

时间尺度上,建议先做一次短周期的预实验(例如连续 72 小时),用 Allan 方差或功率谱密度估计噪声类型和拐点,再据此确定正式实验的采样间隔。经验上,采样间隔应至少比目标漂移特征时间小一个数量级;若要分辨 1 天周期的温度漂移,采样间隔不宜大于 10 分钟。

二、采样频率与时长:用统计量反推,而非拍脑袋

采样频率的选择要同时满足两个约束:一是能覆盖关心的最快变化(抗混叠),二是单次采样本身的测量不确定度要远小于预期漂移量。设单点测量标准差为 σ,预期在时长 T 内产生的漂移为 D,则单点信噪比 D/σ 决定了需要多少次平均。若采用 N 点平均,平均后噪声降为 σ/√N,要求 σ/√N ≤ D/3,即可反推 N。例如某读出链路单次频率估计 σ≈200 Hz,预期 30 天漂移约 5 kHz,则需 N≥(3×200/5000)²≈0.014,说明单点已足够;但若预期漂移只有 300 Hz,则 N≥4,此时可在每个采样时刻做 4 次快速平均再记录。

时长方面,长期稳定性通常要求覆盖至少一个完整的“环境周期”(日周期 + 周周期),并建议不少于 30 天;若关注季节漂移,则需 3~6 个月。时长不足会导致低频分量被截断,Allan 方差在长 τ 处出现虚假下降。一个实用的做法是:正式实验前先设定“最短可接受时长”和“目标时长”两档,前者用于早期判断是否存在异常漂移,后者用于最终统计。

三、环境记录:与主数据同源、同时、同精度

环境记录不是“顺便记一下”,而是漂移归因的核心证据。建议至少记录:实验室空气温度(多点布置,精度 0.1 ℃)、相对湿度、供电电压与频率、冷却水温或制冷机冷头温度、振动(可用加速度计 RMS 值)、磁场(三轴磁通门)。关键要求是环境通道与主测量通道共用同一时间基准,否则后续相关分析会出现时间错位。若条件允许,环境采样频率应不低于主数据采样频率,至少为其 1/2。

建议在数据文件中为每个环境通道保留原始值和时间戳,不要只存分钟平均值。后期做漂移归因时,原始序列可以重采样,而平均值无法还原。

四、数据存储与质量控制

存储设计要保证“可追溯、可复算”。每条记录至少包含:时间戳(UTC,含闰秒标志)、设备编号、通道号、原始读数、单位、量程、增益、状态标志(正常/饱和/超时/人工标记)。建议采用“原始层 + 派生层”两级存储:原始层只追加不修改,派生层存放平均、去趋势、 Allan 方差等结果,并记录所用算法版本和参数。文件命名应包含设备、通道、起止时间,避免依赖人工记忆。

质量控制应贯穿采集全程,而不是事后清洗。可设置三类自动检查:范围检查(读数超出物理合理区间即标记)、跳变检查(相邻点变化超过 k 倍历史标准差即标记)、连续性检查(采样间隔超过标称值 1.5 倍即标记为缺口)。被标记的数据不删除,只打标签,统计时按标签剔除或加权。此外,建议每 7 天插入一次“零输入/参考输入”校验点,用于监测链路自身增益漂移。

长期稳定性监测:采集—记录—质控—存储 数据流 量子测量系统 固定偏置/读取参数 参考锚点 原子钟 / 标准电阻 环境通道 温度 / 湿度 / 供电 振动 / 磁场 / 冷头 与主通道共时钟 采集与时间基准 统一 UTC 时间戳 采样间隔 Δt 单点平均次数 N 量程 / 增益记录 在线质量控制 范围检查 跳变检查(kσ) 连续性 / 缺口标记 7 天参考输入校验 只打标签不删除 状态标志写入记录 两级存储 原始层:只追加 派生层:平均/去趋势 算法版本与参数 设备/通道/起止时间 可追溯、可复算 稳定性分析输出 Allan 方差 / 功率谱密度 / 漂移率 环境相关分析 · 异常时段归因 · 不确定度预算 注:环境通道与主测量通道共用同一时间基准,是漂移归因成立的前提

五、一份可直接套用的实验方案要素

把上述要点落成可执行方案,建议至少写清以下字段(示例参数为合理假设,非实测结果):

要素示例设定设定依据
监测对象读出谐振频率 fr 与读取保真度固定偏置点与读取脉冲
采样间隔10 min(主)、5 min(环境)目标分辨日周期漂移
单点平均N=8 次快速读出σ/√N 小于预期漂移 1/3
时长最短 30 天,目标 90 天覆盖日/周周期,留季节观察余量
环境通道温度 3 点、湿度、供电、振动、磁场覆盖主要漂移诱因
校验点每 7 天一次参考输入监测链路增益漂移
质控规则3σ 跳变、1.5Δt 缺口、量程越界只打标签,不删原始数据
存储原始层 + 派生层,含算法版本保证可追溯与可复算

最后需要强调,长期实验的最大风险往往不是仪器噪声,而是“中途有人动了系统却没人记录”。因此方案中应明确变更日志:任何参数调整、维护、断电、软件升级都必须记录时间、内容和操作人,并与数据标签关联。只有这样,后期在 Allan 方差曲线上看到的拐点,才能被正确解释为系统漂移,而不是一次未被记录的人为干预。