4. 控制图原理与生物过程适用性

Shewhart控制图的基本原理

Shewhart控制图的核心思想是把过程变异分解为两类：由普通原因（common cause）造成的固有随机波动，以及由特殊原因（special cause）造成的异常偏移。控制图并不直接测量“产品是否合格”，而是通过统计量在时间轴上的轨迹，判断过程是否仍处于受控状态。其数学基础是：当过程稳定时，样本统计量围绕中心线CL随机波动，且绝大多数点落在上下控制限UCL、LCL之间。

以均值-极差图（X̄-R图）为例，若子组大小为n，则：

CL = X̄
UCL = X̄ + A2·R̄
LCL = X̄ − A2·R̄
R图：CL = R̄，UCL = D4·R̄，LCL = D3·R̄

其中A2、D3、D4为与子组大小n相关的控制图系数。例如n=5时，A2=0.577，D3=0，D4=2.114。控制限不是规格限，它由过程自身数据计算得到，反映“过程当前实际能做什么”，而规格限反映“客户或工艺要求必须做到什么”。

变差原因：普通原因与特殊原因

普通原因是过程固有的、由大量微小随机因素叠加形成的变异，例如培养基批次内微小差异、电极噪声、环境温度缓慢波动。普通原因造成的分布通常近似稳定，若只存在普通原因，过程处于统计受控状态。特殊原因是可识别、可追溯的异常因素，例如补料泵校准漂移、pH电极污染、灭菌温度异常、操作人员误操作。特殊原因往往造成均值偏移、方差增大或趋势性变化。

控制图判异规则正是围绕这两类原因设计的。常用Western Electric规则包括：1点超出3σ控制限；连续9点在中心线同侧；连续6点递增或递减；连续14点上下交替；连续3点中有2点落在2σ以外；连续5点中有4点落在1σ以外等。这些规则提高对特殊原因的灵敏度，但也带来虚发警报风险。

虚发警报与漏发警报

虚发警报（α风险，第一类错误）指过程实际受控，但控制图却发出异常信号。若仅使用“1点超出3σ”规则，单点虚发概率约为0.27%，即平均每370点出现一次误报。若同时使用多条判异规则，整体虚发概率会显著上升。例如同时使用8条Western Electric规则时，单点误报概率可累积到约1%～2%。

漏发警报（β风险，第二类错误）指过程已经失控，但控制图未能及时报警。漏发概率与偏移大小、子组大小、采样频率有关。偏移越小、子组越小，越难检出。例如均值偏移0.5σ、n=5时，单点检出概率约30%～40%；偏移达到2σ时，检出概率可超过95%。

因此控制图设计需要在虚发与漏发之间权衡：放宽控制限可减少虚发但增加漏发；收紧控制限或增加规则可提高灵敏度但增加误报。生物过程通常不希望频繁误停批，因此常采用“主规则+辅助规则”的分层策略。

生物过程的自相关影响

生物过程数据常具有显著自相关性：当前时刻的pH、溶氧、温度、产物浓度与前一时刻高度相关。这是因为生物反应器具有热惯性、传质滞后、细胞代谢记忆和反馈控制回路。自相关会破坏Shewhart控制图“观测独立”的基本假设，导致两个后果：一是控制限被低估，虚发警报增多；二是对真实偏移的检出延迟，漏发风险上升。

判断自相关可用自相关函数ACF和偏自相关函数PACF。若滞后1阶自相关系数ρ₁显著大于0，例如ρ₁>0.5，则不宜直接使用标准Shewhart图。处理方式包括：

• 增大采样间隔，降低时间相关性；
• 对数据差分或拟合ARIMA模型，用残差作图；
• 采用EWMA或CUSUM控制图，它们对小幅持续偏移更敏感；
• 使用批次间而非批次内数据构建控制图，减少自相关影响。

例如某批次CHO细胞培养中，溶氧DO每5分钟采样一次，ρ₁=0.82。若直接做X̄-R图，控制限过窄，频繁报警。改为每30分钟取一个子组均值，ρ₁降至0.35，再结合EWMA图，误报明显减少，且能及时捕捉到补料泵故障导致的DO持续下降。

生物过程控制图的选型建议

对于生物过程，建议按以下逻辑选择控制图：若数据近似独立且子组合理，优先使用X̄-R图或X̄-s图；若存在明显自相关，先尝试增大采样间隔或差分，再考虑EWMA或CUSUM；若关注批次间一致性，可使用批次均值图或T²多元控制图。无论选择哪种图，都应先验证数据独立性，再计算控制限，最后结合判异规则和工艺知识综合判断，避免把统计信号直接等同于工艺故障。