4. 取样频率与样本量的确定

取样频率和样本量不是凭经验“拍脑袋”定的两个数字,而是由三个变量共同约束的结果:过程本身的波动水平、该质量特性失控后造成的后果严重度、以及我们希望在多长时间内发现异常。SPC 提供了把这三者联系起来的定量工具,风险等级决定了我们对漏检的容忍度,工艺稳定性则决定了取样间隔能否拉长。本节围绕“怎么算、怎么定、怎么调”展开。

一、先明确两个不同的量:样本量 n 与取样频率 f

样本量 n 指一次取样中抽取的独立样品数(或一次检测中重复测定的次数),它主要影响对“当前时刻过程水平”估计的精度,以及控制图的检出力。取样频率 f 指两次取样之间的时间间隔或批次间隔,它影响的是发现异常的速度,即平均运行长度 ARL 意义上的响应延迟。二者不能互相替代:n 再大,若两次取样间隔是 8 小时,中间 8 小时的漂移仍然看不见;f 再密,若每次只取 1 个样品且测量噪声大,单点判异会频繁误报。

因此工程上通常先定 n(保证单点判异可靠),再定 f(保证异常暴露及时),最后用风险等级做双向修正。

二、用 SPC 原理确定最小样本量

对于计量型数据(如含量、pH、水分、粒径 D50),控制图基于样本均值分布。若单次测量标准差为 σ,则 n 个样品的均值标准差为 σ/√n。控制限宽度为 ±3σ/√n,因此 n 越大,控制限越窄,对微小偏移越敏感。常用的定量依据是检出力(Power):希望在过程均值偏移 Δ 时,控制图能以 1−β 的概率报警。

以常用的 X̄ 图为例,偏移量用标准差单位表示为 δ = Δ/σ。单点超出 3σ 控制限的检出力近似为:

Power ≈ Φ(δ·√n − 3) + Φ(−δ·√n − 3)

示例(合理假设):某中间体含量 σ = 0.8%,希望检出 Δ = 1.0% 的偏移,
要求检出力 ≥ 0.80。
δ = 1.0 / 0.8 = 1.25
试 n = 4:δ√n = 2.50,Φ(2.50−3)=Φ(−0.5)≈0.309,Power≈0.31
试 n = 9:δ√n = 3.75,Φ(0.75)≈0.773,Power≈0.77
试 n = 10:δ√n = 3.95,Φ(0.95)≈0.829,Power≈0.83  → 满足
结论:最小样本量 n = 10(单次取样合并测定或分批测定后取均值)

若采用计数型数据(如外观不合格品率),则依据二项分布,样本量需满足在给定 p 下能稳定检出变化,常用规则是 np ≥ 5 且 n(1−p) ≥ 5 才能用正态近似建立 p 图;若预期不合格率 p = 2%,则 n 至少 250 件,否则控制限会出现负值而下限截断为 0,图形失真。中试阶段批量小,往往改用累积计数或单位缺陷数 u 图,以“每千克物料缺陷数”为统计量,避免样本量不足。

三、用风险等级修正取样频率

风险等级来自“该特性失控对后续工序、成品质量、安全或注册申报的影响”。可用一个简单的分级矩阵把风险映射为频率系数:

风险等级典型特性示例失控后果频率修正
高(关键质量属性 CQA)有关物质、残留溶剂、无菌/微生物影响安全性或注册放行每批必取,必要时批内多点
中(重要工艺参数关联项)含量、pH、水分、粒度分布影响收率与下游加工每批或每 2 批,批内 1 次
低(一般监控项)外观、颜色、堆密度可返工或影响观感每 3~5 批或每周 1 次

注意:高风险的“每批必取”是底线,不能因为工艺稳定就降低;能调整的是批内取样点数和是否合并测定。

四、用工艺稳定性决定间隔能否拉长

工艺稳定性用控制图的判稳结果和过程能力指数 Cpk 共同判断。中试阶段建议采用“先密后疏”的阶梯策略:

  1. 导入期(前 10~15 批):按最小样本量 n 执行,频率取最密档(如每批、每 2 小时),目的是尽快积累数据建立控制限,并暴露批内漂移。
  2. 稳定期判定:连续 20 点无判异规则触发(无点超出 3σ、无连续 7 点同侧、无连续 6 点递增递减),且 Cpk ≥ 1.33,方可进入放宽评估。
  3. 放宽期:频率可降一档(如每 2 批→每 5 批),但 n 不降低,且必须保留“首批/末批/异常批”加密取样的触发条件。
  4. 回退条件:一旦出现判异、Cpk 跌破 1.0、更换原料供应商或放大设备,立即回到导入期频率,重新积累 10 批数据。
中试阶段批量小、批次数少,切勿直接套用商业化生产的放宽频率。中试的核心目的之一就是验证工艺稳定性,频率过疏会导致“看起来稳定”的假象——因为根本没有足够数据去发现不稳定。

五、把三者合成一张可执行的取样计划

下面这张图把“SPC 计算 → 风险修正 → 稳定性调整”的决策关系画成一条可追溯的判定链,便于在取样计划表中逐项填写依据。

取样频率与样本量确定逻辑链 输入 1:过程波动 σ 历史批数据 / 预试验 目标偏移 Δ、检出力 1−β 输入 2:风险等级 CQA / 重要项 / 一般项 失控后果严重度 输入 3:工艺稳定性 判稳规则 / Cpk 批次数积累情况 计算最小样本量 n Power = Φ(δ√n−3)+Φ(−δ√n−3) 计数型:np≥5 且 n(1−p)≥5 取满足检出力要求的最小 n 风险映射频率档 高:每批必取(底线) 中:每批 / 每 2 批 低:每 3~5 批 / 每周 稳定性阶梯判定 导入期:最密档 稳定期:20 点无判异且 Cpk≥1.33 放宽期:降一档,n 不降 合成取样计划:n(样本量)× f(频率)× 触发条件 n 由检出力决定,f 由风险档与稳定性档共同决定,触发条件兜底 首批 / 末批 / 换料 / 放大 / 判异 → 自动加密 回退触发:判异 / Cpk<1.0 / 换供应商 / 设备放大 立即回到导入期频率,重新积累 10 批数据

六、一个可填写的取样计划示例

以下为某中试车间“结晶工序中间体”的取样计划片段,参数为教学示例(合理假设),实际执行需以本品种的方法验证与历史数据替换。

质量特性风险等级σ(示例)目标 Δ最小 n导入期频率放宽条件
含量(HPLC)中0.8%1.0%10(合并测定)每批 1 次20 点判稳且 Cpk≥1.33 → 每 2 批
有关物质高0.05%0.05%每批必取,不合并每批 1 次不降频,仅可评估批内点数
水分(KF)中0.15%0.20%6每批 1 次每 3 批
外观/颜色低——1每 3 批每周 1 次
填写要点:n 一栏必须写明“单点测定”还是“合并后测定”,因为合并会掩盖批内不均匀性;频率一栏必须同时写明放宽的量化门槛(判稳点数、Cpk 阈值),否则“工艺稳定后可放宽”就是一句无法审计的空话。

最后提醒一个中试常见误区:把“取样频率高”当成质量保证的手段,却忽略了每次取样本身引入的干扰——频繁开盖取样可能影响结晶体系温度与溶剂挥发,反而制造波动。合理的做法是在保证检出力与响应速度的前提下,用最少的取样次数获得足够的判定信息,这正是 SPC 与风险分级结合的价值所在。