血球分析仪的阻抗法检测系统和识别阻抗通道堵孔的方法与流程

本申请涉及微粒检测技术领域，具体涉及一种血球分析仪的阻抗法检测系统和识别阻抗通道堵孔的方法。

背景技术：

通过阻抗法检测微粒(细胞)时，将小孔浸入电解质溶液中，在小孔两端通上恒流电源，当微粒(细胞)通过小孔时，小孔两端的电压发生变化，微粒(细胞)的体积越大，电压的变化就越大(电压脉冲越大)。以血液检测领域的应用为例，按照该原理将血液在稀释液(一种试剂具有导电性，且具有生理盐水的特性)中稀释混匀，然后通过两端有电极小孔，就能够得到通过小孔的血细胞的脉冲(体积)。

在使用阻抗法检测微粒(细胞)时，有时会发生堵孔现象，即，小孔堵塞导致待检测的微粒(细胞)无法通过小孔，影响检测的继续进行。血细胞稀释混匀后，在溶液中是均匀分布的，所以单位时间内(比如100ms)所有通过小孔的血细胞的平均脉冲大小应该是一致的。如果某一时刻发生了堵孔，则该时刻前后的细胞平均脉冲会相差很大。因此，有些方法中通过小孔电压曲线的变化来判定是否堵孔，如果小孔的电压在某一时刻后有明显的变化则认为发生了堵孔。但是，当全程堵孔时，小孔的电压可能无明显变化，该方法无法准确报告堵孔。

技术实现要素：

本发明主要解决的技术问题是：现有的方法无法在全程堵孔时，准确报告堵孔。

第一方面，一种实施例中提供一种血球分析仪的阻抗法检测系统，包括：

检测室，包括一小孔；所述检测室用于接收血液样本、溶血剂和稀释液以制备待测定的试样；

压力源部件，用于提供压力以使得检测室中的试样的细胞通过所述小孔；

恒流电源，所述恒流电源的两端分别与所述小孔的两端电连接；

电压脉冲检测装置，所述电压脉冲检测装置与所述恒流电源连接，用于检测细胞通过所述小孔时产生的电压脉冲，并发送所述电压脉冲的数据；

存储器，所述存储器用于接收和存储历史平均脉冲数据，和所述电压电源脉冲检测装置所发送的电压脉冲的数据；以及

处理器，所述处理器用于获取历史平均脉冲数据，对历史平均脉冲数据求均值和标准差；获取本次采样分析的电压脉冲的数据，统计出每个单位时间内通过小孔的细胞的平均脉冲；将每个单位时间前所有单位时间的平均脉冲的均值定义为前均值，依次计算出每个单位时间的前均值，将每个单位时间后所有单位时间的平均脉冲的均值定义为后均值，依次计算每个单位时间的后均值，根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间；计算可能堵孔的单位时间的前均值与可能堵孔的单位时间的后均值间的第一变化量，将所述第一变化量与可能堵孔的单位时间前和/或可能堵孔的单位时间后所有单位时间的平均脉冲的标准差进行比较，第一次判断可能堵孔的单位时间是否发生了堵孔，如果第一次判断发生了堵孔，上报异常，如果第一次判断未发生堵孔，继续执行下一步；第一次判断可能堵孔的单位时间未发生堵孔后，计算本次采样的所有单位时间的平均脉冲的均值，并计算本次采样的所有单位时间的平均脉冲的均值与历史平均脉冲数据的均值间的第二变化量，将所述第二变化量与历史平均脉冲数据的标准差进行比较，第二次判断可能堵孔的单位时间是否发生了堵孔，如果第二次判断发生了堵孔，上报异常，如果第二次判断未发生堵孔，将本次采样的所有单位时间的平均脉冲的均值发送给存储器，并存入历史平均脉冲数据。

一种实施例中，所述处理器获取历史平均脉冲数据，对历史平均脉冲数据求均值和标准差，包括：

获取历史平均脉冲数据的链表premplist，所述历史平均脉冲数据为过去每次采样的所有时间单位的平均脉冲的均值，获取历史平均脉冲数据的均值的计算公式，计算出历史平均脉冲数据的均值mpre，获取历史平均脉冲数据的标准差的计算公式，计算出历史平均脉冲数据的标准差spre；

所述历史平均脉冲数据的均值的计算公式包括：

其中，n为premplist内数据的个数；

所述历史平均脉冲数据的标准差的计算公式包括：

其中，n为premplist内数据的个数。

一种实施例中，所述处理器获取本次采样分析的电压脉冲的数据，统计出每个单位时间内通过小孔的细胞的平均脉冲，包括：

获取总采集时间t，获取单位时间δt，获取单位时间的总个数n的计算公式，计算出单位时间的总个数n，所述单位时间的总个数n的计算公式包括：

n＝t/δt，

获取计算第i个单位时间的脉冲个数的公式，计算出第i个单位时间的脉冲个数，所述计算第i个单位时间的脉冲个数的公式包括：

mi＝count(p(t,v))，

其中，t的取值范围为(i-1)δt≤t＜iδt，i的取值范围为0＜i≤n，i为正整数；

获取计算第i个单位时间的平均脉冲的计算公式，计算出第i个单位时间的平均脉冲，所述计算第i个单位时间的平均脉冲的计算公式包括：

pi＝(∑v)/mi，

其中，v∈p(t,v)，t的取值范围为(i-1)δt≤t＜iδt，i的取值范围为0＜i≤n，i为正整数；

将每个单位时间内的平均脉冲依次存入记数组p中，得到记数组p＝[p1,p2,…,pi,…,pn]。

一种实施例中，所述处理器依次计算出每个单位时间的前均值，依次计算每个单位时间的后均值，包括：

获取记数组p，获取计算每个单位的前均值的计算公式，依次计算出每个单位时间的前均值s1i，所述每个单位的前均值的计算公式包括：

其中，i的取值范围为0＜i≤n-1，i为正整数；

获取计算每个单位的后均值的计算公式，依次计算出每个单位时间的后均值s2i，所述每个单位的后均值的计算公式包括：

其中，i的取值范围为0＜i≤n-1，i为正整数；

获取初始化的链表list1和链表list2，将每个单位时间的前均值s1i存入list1中，得到链表list1[i]，将每个单位时间的后均值s2i存入list2中，得到链表list2[i]。

一种实施例中，所述处理器根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间，包括：

获取链表list1[i]和链表list2[i]，获取差分计算公式，对链表list1[i]内的数据进行差分计算，获取初始化的链表dlist1，将链表list1[i]差分计算的结果存入链表dlist1中，得到链表dlist1[i]，对链表list2[i]内的数据进行差分计算，获取初始化的链表dlist2，将链表list2[i]差分计算的结果存入链表dlist2中，得到链表dlist2[i]，所述差分计算公式包括：

dlist1[i]＝list1[i+1]-list1[i]，

dlist2[i]＝list2[i+1]-list2[i]，

其中，i的取值范围为0＜i≤n-1，i为正整数；

获取可能堵孔的单位时间判定公式，所述可能堵孔的单位时间判定公式包括：

dlist2[i]*dlist1[i+1]<0，

遍历链表dlist2[i]，将链表dlist1[i]和链表dlist2[i]中的数据输入所述可能堵孔的单位时间判定公式中，如果数据符合所述可能堵孔的判定公式，则第i+1个单位时间为可能堵孔的单位时间。

一种实施例中，所述处理器根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间，包括：

获取链表list1[i]和链表list2[i]，获取差分计算公式，对链表list1[i]内的数据进行差分计算，获取初始化的链表dlist1，将链表list1[i]差分计算的结果存入链表dlist1中，得到链表dlist1[i]，对链表list2[i]内的数据进行差分计算，获取初始化的链表dlist2，将链表list2[i]差分计算的结果存入链表dlist2中，得到链表dlist2[i]，所述差分计算公式包括：

dlist1[i]＝list1[i+1]-list1[i]，

dlist2[i]＝list2[i+1]-list2[i]，

其中，i的取值范围为0＜i≤n-1，i为正整数；

获取可能堵孔的单位时间判定公式，所述可能堵孔的单位时间判定公式包括：

dlist1[i]*dlist2[i+1]<0，

遍历链表dlist1[i]，将链表dlist1[i]和链表dlist2[i]中的数据输入所述可能堵孔的单位时间判定公式中，如果数据符合所述可能堵孔的判定公式，则第i+1个单位时间为可能堵孔的单位时间。

一种实施例中，所述处理器第一次判断可能堵孔的单位时间是否发生了堵孔，包括：

获取记数组p，获取第i+1个单位时间为可能堵孔的单位时间，获取前i个单位时间的平均脉冲的均值的计算公式，计算出前i个单位时间的平均脉冲的均值mp1，所述前i个单位时间的平均脉冲的均值的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数；

获取前i个单位时间的平均脉冲的标准差的计算公式，计算出前i个单位时间的平均脉冲的标准差sp1，所述前i个单位时间的平均脉冲的标准差的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数；

获取第i+1到n个单位时间的平均脉冲的均值的计算公式，计算出第i+1到n个单位时间的平均脉冲的均值mp2，所述第i+1到n个单位时间的平均脉冲的均值的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数；

获取第i+1到n个单位时间的平均脉冲的标准差的计算公式，计算出第i+1到n个单位时间的平均脉冲的标准差sp2，所述第i+1到n个单位时间的平均脉冲的标准差的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数；

获取第一次判断的判断公式，所述第一次判断的判断公式包括：

|mp1-mp2|＜sp1，

|mp1-mp2|＜sp2，

将前i个单位时间的平均脉冲的均值mp1、前i个单位时间的平均脉冲的标准差sp1、第i+1到n个单位时间的平均脉冲的均值mp2、第i+1到n个单位时间的平均脉冲的标准差sp2输入第一次判断的判断公式，如果所述第一次判断的判断公式中的至少一个公式成立，则第一次判断发生了堵孔，否则第一次判断未发生堵孔。

一种实施例中，所述处理器第二次判断可能堵孔的单位时间是否发生了堵孔，包括：

获取记数组p，获取本次采样的所有单位时间的平均脉冲的均值的计算公式，计算出本次采样的所有单位时间的平均脉冲的均值mp，所述本次采样的所有单位时间的平均脉冲的均值的计算公式包括：

获取历史平均脉冲数据的均值mpre和历史平均脉冲数据的标准差spre，获取第二次判断的判断公式，将mp、mpre和spre代入第二次判断的判断公式，

如果公式

|mpre-mp|＜3spre

成立，则第二次判断未发生堵孔；

如果公式

|mpre-mp|≥3spre

成立，则第二次判断发生了全程堵孔。

第二方面，一种实施例中提供一种识别阻抗通道堵孔的方法，用于检测血球分析仪的阻抗法检测系统的是否发生了堵孔，包括：

获取历史平均脉冲数据，对历史平均脉冲数据求均值和标准差；

获取本次采样分析的电压脉冲的数据，统计出每个单位时间内通过小孔的细胞的平均脉冲；

将每个单位时间前所有单位时间的平均脉冲的均值定义为前均值，依次计算出每个单位时间的前均值，将每个单位时间后所有单位时间的平均脉冲的均值定义为后均值，依次计算每个单位时间的后均值，根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间；

计算可能堵孔的单位时间的前均值与可能堵孔的单位时间的后均值间的第一变化量，将所述第一变化量与可能堵孔的单位时间前和/或可能堵孔的单位时间后所有单位时间的平均脉冲的标准差进行比较，第一次判断可能堵孔的单位时间是否发生了堵孔，如果第一次判断发生了堵孔，上报异常，如果第一次判断未发生堵孔，继续执行下一步；

第一次判断可能堵孔的单位时间未发生堵孔后，计算本次采样的所有单位时间的平均脉冲的均值，并计算本次采样的所有单位时间的平均脉冲的均值与历史平均脉冲数据的均值间的第二变化量，将所述第二变化量与历史平均脉冲数据的标准差进行比较，第二次判断可能堵孔的单位时间是否发生了堵孔，如果第二次判断发生了堵孔，上报异常，如果第二次判断未发生堵孔，将本次采样的所有单位时间的平均脉冲的均值存入历史平均脉冲数据。

第三方面，一种实施例中提供一种计算机可读存储介质，所述介质上存储有程序，所述程序能够被处理器执行以实现上述的方法。

依据上述实施例的血球分析仪的阻抗法检测系统和识别阻抗通道堵孔的方法，由于增加了将第二变化量与历史平均脉冲数据的标准差进行比较，第二次判断可能堵孔的单位时间是否发生了堵孔，实现了对全程堵孔的情况的识别和报告。

附图说明

图1为库尔特原理的示意图；

图2为阻抗法检测示意图；

图3为一种实施例中血球分析仪的阻抗法检测系统的框架示意图；

图4为一种实施例中识别阻抗通道堵孔的方法的流程示意图；

图5为一种实施例中未堵孔时list2[i]中数据的曲线示意图；

图6为一种实施例中堵孔时list2[i]中数据的曲线示意图。

具体实施方式

下面通过具体实施方式结合附图对本发明作进一步详细说明。其中不同实施方式中类似元件采用了相关联的类似的元件标号。在以下的实施方式中，很多细节描述是为了使得本申请能被更好的理解。然而，本领域技术人员可以毫不费力的认识到，其中部分特征在不同情况下是可以省略的，或者可以由其他元件、材料、方法所替代。在某些情况下，本申请相关的一些操作并没有在说明书中显示或者描述，这是为了避免本申请的核心部分被过多的描述所淹没，而对于本领域技术人员而言，详细描述这些相关操作并不是必要的，他们根据说明书中的描述以及本领域的一般技术知识即可完整了解相关操作。

另外，说明书中所描述的特点、操作或者特征可以以任意适当的方式结合形成各种实施方式。同时，方法描述中的各步骤或者动作也可以按照本领域技术人员所能显而易见的方式进行顺序调换或调整。因此，说明书和附图中的各种顺序只是为了清楚描述某一个实施例，并不意味着是必须的顺序，除非另有说明其中某个顺序是必须遵循的。

本文中为部件所编序号本身，例如“第一”、“第二”等，仅用于区分所描述的对象，不具有任何顺序或技术含义。而本申请所说“连接”、“联接”，如无特别说明，均包括直接和间接连接(联接)。

请参考图1和2，图1为库尔特原理的示意图，将小孔浸入电解质溶液中，在小孔两端通上恒流电源，当微小颗粒通过小孔时，小孔两端的电压发生变化，微小颗粒体积越大，电压的变化就越大(电压脉冲越大)。

图2为阻抗法(库尔特原理)检测示意图，按照库尔特原理将血液在稀释液(一种试剂具有导电性，且具有生理盐水的特性)中稀释混匀，然后通过两端有电极小孔，就能够得到通过小孔细胞的脉冲(体积)。

血细胞稀释混匀后，在溶液中是均匀分布的，所以单位时间内(比如100ms)所有通过小孔的细胞的平均脉冲大小应该是一致的。如果某一时刻发生了堵孔，则该时刻前后的细胞平均脉冲会相差很大。本申请就是基于该原理实现堵孔的判定。

每个通过小孔的细胞，仪器不仅会记录下脉冲的大小，也会记录记录下这个脉冲的时间，即每个脉冲都有两个维度的标记，即分别是脉冲值的大小v，和这个脉冲产生的时间t，记p(t,v)表示一个脉冲。记t为样本采集的总时长(单位s)。

一方面，本实施例提供一种血球分析仪的阻抗法检测系统。

请参考图3，该阻抗法检测系统包括检测室1、压力源部件2、恒流电源3、电压脉冲检测装置4、存储器5和处理器6。

检测室1包括一小孔，检测室1用于接收血液样本、溶血剂和稀释液以制备待测定的试样。压力源部件2用于提供压力以使得检测室1中的试样的细胞通过小孔，恒流电源3的两端分别与小孔的两端电连接。电压脉冲检测装置4与恒流电源3连接，用于检测细胞通过小孔时产生的电压脉冲，并发送电压脉冲的数据。存储器5用于接收和存储历史平均脉冲数据，和电压电源脉冲检测装置4所发送的电压脉冲的数据。

处理器6用于获取历史平均脉冲数据，对历史平均脉冲数据求均值和标准差。获取本次采样分析的电压脉冲的数据，统计出每个单位时间内通过小孔的细胞的平均脉冲。将每个单位时间前所有单位时间的平均脉冲的均值定义为前均值，依次计算出每个单位时间的前均值，将每个单位时间后所有单位时间的平均脉冲的均值定义为后均值，依次计算每个单位时间的后均值，根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间。计算可能堵孔的单位时间的前均值与可能堵孔的单位时间的后均值间的第一变化量，将第一变化量与可能堵孔的单位时间前和/或可能堵孔的单位时间后所有单位时间的平均脉冲的标准差进行比较，第一次判断可能堵孔的单位时间是否发生了堵孔，如果第一次判断发生了堵孔，上报异常，如果第一次判断未发生堵孔，继续执行下一步。第一次判断可能堵孔的单位时间未发生堵孔后，计算本次采样的所有单位时间的平均脉冲的均值，并计算本次采样的所有单位时间的平均脉冲的均值与历史平均脉冲数据的均值间的第二变化量，将第二变化量与历史平均脉冲数据的标准差进行比较，第二次判断可能堵孔的单位时间是否发生了堵孔，如果第二次判断发生了堵孔，上报异常，如果第二次判断未发生堵孔，将本次采样的所有单位时间的平均脉冲的均值发送给存储器，并存入历史平均脉冲数据。

另一方面，本实施例提供一种识别阻抗通道堵孔的方法，该方法用于检测血球分析仪的阻抗法检测系统的是否发生了堵孔。

本实施例的识别阻抗通道堵孔的方法有如下优点：

第一，有的方案通过粒子流的稳定性来判定是否堵孔，如果粒子流在某一时刻后发生跳变则认为发生了堵孔。但是，当全程堵孔时粒子流无明显变化，无法准确报准堵孔。当样本混匀不充分时，粒子流稳定性无法保证，上述方案易误报堵孔。本实施例的方法与上述方案相比，由于增加了将第二变化量与历史平均脉冲数据的标准差进行比较，第二次判断可能堵孔的单位时间是否发生了堵孔，实现了对全程堵孔的情况的识别和报告。也解决了由于混匀不好，导致粒子流跳变误报堵孔的问题。

第二，有的方案通过小孔电压曲线的变化来判定是否堵孔，但是，当全程堵孔时小孔电压可能无明显变化，无法准确报准堵孔。本实施例的方法与上述方案相比，实现了对全程堵孔的情况的识别和报告。

第三，有的方案以开机时的小孔电压为参考，结合小孔电压曲线的变化来判定是否发生堵孔，如果小孔电压无变化但明显高于开机电压，就判定全程堵孔，如果小孔电压在某一时刻后有明显的变化则认为发生了堵孔。但是，当稀释液的温度和开机时的温度有明显差异时，因稀释液在不同温度下导电率不一样使得开机时的参考电压失效，无法准确报堵孔。本实施例的方法与上述方案相比，不易受到温度变化的干扰。

请参考图4，该识别阻抗通道堵孔的方法包括：

请参考图4，步骤100，获取历史平均脉冲数据，对历史平均脉冲数据求均值和标准差。

在一种实施例中，获取历史平均脉冲数据，对历史平均脉冲数据求均值和标准差，包括：

获取历史平均脉冲数据的链表premplist，历史平均脉冲数据为过去每次采样的所有时间单位的平均脉冲的均值，获取历史平均脉冲数据的均值的计算公式，计算出历史平均脉冲数据的均值mpre，获取历史平均脉冲数据的标准差的计算公式，计算出历史平均脉冲数据的标准差spre。

历史平均脉冲数据的均值的计算公式包括：

其中，n为premplist内数据的个数。

历史平均脉冲数据的标准差的计算公式包括：

其中，n为premplist内数据的个数。

请参考图4，步骤200，获取本次采样分析的电压脉冲的数据，统计出每个单位时间内通过小孔的细胞的平均脉冲。

在一种实施例中，获取本次采样分析的电压脉冲的数据，统计出每个单位时间内通过小孔的细胞的平均脉冲，包括：

获取总采集时间t，获取单位时间δt，获取单位时间的总个数n的计算公式，计算出单位时间的总个数n，单位时间的总个数n的计算公式包括：

n＝t/δt，

具体的，本实施例中δt取0.1s，在其他实施例中，δt也可以根据实际需求选择其他合适的值。

获取计算第i个单位时间的脉冲个数的公式，计算出第i个单位时间的脉冲个数，计算第i个单位时间的脉冲个数的公式包括：

mi＝count(p(t,v))，

其中，t的取值范围为(i-1)δt≤t＜iδt，i的取值范围为0＜i≤n，i为正整数。

获取计算第i个单位时间的平均脉冲的计算公式，计算出第i个单位时间的平均脉冲，计算第i个单位时间的平均脉冲的计算公式包括：

pi＝(∑v)/mi，

其中，v∈p(t,v)，t的取值范围为(i-1)δt≤t＜iδt，i的取值范围为0＜i≤n，i为正整数。

将每个单位时间内的平均脉冲依次存入记数组p中，得到记数组(向量)p＝[p1,p2,…,pi,…,pn]。

请参考图4，步骤300，将每个单位时间前所有单位时间的平均脉冲的均值定义为前均值，依次计算出每个单位时间的前均值，将每个单位时间后所有单位时间的平均脉冲的均值定义为后均值，依次计算每个单位时间的后均值。

在一种实施例中，依次计算出每个单位时间的前均值，依次计算每个单位时间的后均值，包括：

获取记数组p，获取计算每个单位的前均值的计算公式，依次计算出每个单位时间的前均值s1i，每个单位的前均值的计算公式包括：

其中，i的取值范围为0＜i≤n-1，i为正整数。

获取计算每个单位的后均值的计算公式，依次计算出每个单位时间的后均值s2i，每个单位的后均值的计算公式包括：

其中，i的取值范围为0＜i≤n-1，i为正整数。

获取初始化的链表list1和链表list2，将每个单位时间的前均值s1i存入list1中，得到链表list1[i]，请参考图5和6，将每个单位时间的后均值s2i存入list2中，得到链表list2[i]。

请参考图4，步骤400，根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间。

在一种实施例中，根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间，包括：

获取链表list1[i]和链表list2[i]，获取差分计算公式，对链表list1[i]内的数据进行差分计算，获取初始化的链表dlist1，将链表list1[i]差分计算的结果存入链表dlist1中，得到链表dlist1[i]，对链表list2[i]内的数据进行差分计算，获取初始化的链表dlist2，将链表list2[i]差分计算的结果存入链表dlist2中，得到链表dlist2[i]，差分计算公式包括：

dlist1[i]＝list1[i+1]-list1[i]，

dlist2[i]＝list2[i+1]-list2[i]，

其中，i的取值范围为0＜i≤n-1，i为正整数。

获取可能堵孔的单位时间判定公式，可能堵孔的单位时间判定公式包括：

dlist2[i]*dlist1[i+1]<0，

遍历链表dlist2[i]，将链表dlist1[i]和链表dlist2[i]中的数据输入可能堵孔的单位时间判定公式中，如果数据符合可能堵孔的判定公式，则第i+1个单位时间为可能堵孔的单位时间。

在另一种实施例中，根据每个单位时间的前均值和后均值的变化情况判断出可能堵孔的单位时间，包括：

获取链表list1[i]和链表list2[i]，获取差分计算公式，对链表list1[i]内的数据进行差分计算，获取初始化的链表dlist1，将链表list1[i]差分计算的结果存入链表dlist1中，得到链表dlist1[i]，对链表list2[i]内的数据进行差分计算，获取初始化的链表dlist2，将链表list2[i]差分计算的结果存入链表dlist2中，得到链表dlist2[i]，差分计算公式包括：

dlist1[i]＝list1[i+1]-list1[i]，

dlist2[i]＝list2[i+1]-list2[i]，

其中，i的取值范围为0＜i≤n-1，i为正整数。

获取可能堵孔的单位时间判定公式，可能堵孔的单位时间判定公式包括：

dlist1[i]*dlist2[i+1]<0，

遍历链表dlist1[i]，将链表dlist1[i]和链表dlist2[i]中的数据输入可能堵孔的单位时间判定公式中，如果数据符合可能堵孔的判定公式，则第i+1个单位时间为可能堵孔的单位时间。

请参考图4，步骤500，计算可能堵孔的单位时间的前均值与可能堵孔的单位时间的后均值间的第一变化量，将第一变化量与可能堵孔的单位时间前和/或可能堵孔的单位时间后所有单位时间的平均脉冲的标准差进行比较，第一次判断可能堵孔的单位时间是否发生了堵孔，如果第一次判断发生了堵孔，上报异常，如果第一次判断未发生堵孔，继续执行下一步。

在一种实施例中，第一次判断可能堵孔的单位时间是否发生了堵孔，包括：

获取记数组p，获取第i+1个单位时间为可能堵孔的单位时间，获取前i个单位时间的平均脉冲的均值的计算公式，计算出前i个单位时间的平均脉冲的均值mp1，前i个单位时间的平均脉冲的均值的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数。

获取前i个单位时间的平均脉冲的标准差的计算公式，计算出前i个单位时间的平均脉冲的标准差sp1，前i个单位时间的平均脉冲的标准差的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数。

获取第i+1到n个单位时间的平均脉冲的均值的计算公式，计算出第i+1到n个单位时间的平均脉冲的均值mp2，第i+1到n个单位时间的平均脉冲的均值的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数。

获取第i+1到n个单位时间的平均脉冲的标准差的计算公式，计算出第i+1到n个单位时间的平均脉冲的标准差sp2，第i+1到n个单位时间的平均脉冲的标准差的计算公式包括：

其中，i的取值范围为0＜i＜n，i为正整数。

获取第一次判断的判断公式，第一次判断的判断公式包括：

|mp1-mp2|＜sp1，

|mp1-mp2|＜sp2，

将前i个单位时间的平均脉冲的均值mp1、前i个单位时间的平均脉冲的标准差sp1、第i+1到n个单位时间的平均脉冲的均值mp2、第i+1到n个单位时间的平均脉冲的标准差sp2输入第一次判断的判断公式，如果第一次判断的判断公式中的至少一个公式成立，则第一次判断发生了堵孔，否则第一次判断未发生堵孔。判断发生了堵孔后，上报白/红细胞通道异常，并下发排堵指令。上述“第一变化量”为公式中的|mp1-mp2|。

请参考图4，步骤600，第一次判断可能堵孔的单位时间未发生堵孔后，计算本次采样的所有单位时间的平均脉冲的均值，并计算本次采样的所有单位时间的平均脉冲的均值与历史平均脉冲数据的均值间的第二变化量，将第二变化量与历史平均脉冲数据的标准差进行比较，第二次判断可能堵孔的单位时间是否发生了堵孔，如果第二次判断发生了堵孔，上报异常，如果第二次判断未发生堵孔，将本次采样的所有单位时间的平均脉冲的均值存入历史平均脉冲数据。

在一种实施例中，第二次判断可能堵孔的单位时间是否发生了堵孔，包括：

获取记数组p，获取本次采样的所有单位时间的平均脉冲的均值的计算公式，计算出本次采样的所有单位时间的平均脉冲的均值mp，本次采样的所有单位时间的平均脉冲的均值的计算公式包括：

获取历史平均脉冲数据的均值mpre和历史平均脉冲数据的标准差spre，获取第二次判断的判断公式，将mp、mpre和spre代入第二次判断的判断公式，

如果公式

|mpre-mp|＜3spre

成立，则第二次判断未发生堵孔。

如果公式

|mpre-mp|≥3spre

成立，则第二次判断发生了全程堵孔。判断发生了堵孔后，上报白/红细胞通道异常，并下发排堵指令。上述“第二变化量”为公式中的|mpre-mp|。

另一方面，本实施例提供一种计算机可读存储介质，该计算机可读存储介质上存储有程序，程序能够被处理器执行以实现上述的方法。

以上应用了具体个例对本发明进行阐述，只是用于帮助理解本发明，并不用以限制本发明。对于本发明所属技术领域的技术人员，依据本发明的思想，还可以做出若干简单推演、变形或替换。