CN115165792B - 管线中气液两相流状态检测方法及装置 - Google Patents

Abstract

本发明公开一种管线中气液两相流状态检测方法及装置，该检测方法基于沿待测管线的管径方向布置的至少一对位于所述管线外部的光发射器和光接收器进行检测，包括：实时检测并处理光接收器接收到的光信号，以生成反应光信号强度变化率的功率谱；获取功率谱中当前时刻的峰值；将功率谱的峰值与预设的第一阈值比较，并根据比较结果确定管线中当前气液两相流介质的存在状态；上述检测方法，通过反应光信号强度变化率的功率谱来表征管线中的介质流状态，无论管径大小，也即，即使极小管径的管线，管线中的气液两相流介质的存在状态的变化也会反应在功率谱中，从而有效提高通过光线传感器对小管径中气液两相流介质状态检测的精准度。

Description

管线中气液两相流状态检测方法及装置

技术领域

本发明涉及气液两相流检测技术领域，尤其涉及一种管线中气液两相流状态检测方法及装置。

背景技术

通常情况下，直接将待输送的液体流导入管道进行输送，此时管道中的介质流称之为单相流。另一种情况，为提高液体的输送效率或特殊需要，在将液体导入输送管道的同时通入高压气体流，从而在输送管道中形成气液两相流，因此，气液两相流指的是包括气体和液体的混合介质流，一般用于在管道输送中。

在应用过程中，一般通过红外检测装置来监控输送管道中流动液体的状态，其工作原理是：红外激光发射器沿管径方向向管道投射红外光，红外接收器接收从管道的另一侧透射出的红外光，然后分析接收到的红外光的光强度，以计算红外光通过管道后的衰减率，当红外光的衰减率大于预设值时表示管道内存在有液体介质流。对于这种检测方式，当输送管道的管径特别小(毫米级别) 时，如油气润滑管线(管线内有空气和油液的混合介质)，由于管道内的气液介质流的截面比较单薄，通过上述红外检测处理方式，很难精确检测管道内的介质流的状态，因此，现有技术中，一般通过超声波来检测小管径的液流状态，检测结果不够理想。

相比于单相流，气液两相流在管道内具有特有的流型，如图11，沿流动方向，液体介质膜呈周期性间隔分布，且液体介质膜附着在管壁上。因此，可基于气液两相流的流型，对光传感器测到的数据处理进行改进，以使得可基于光线检测的传感器对气液两相流的状态进行检测，并获得良好效果。

发明内容

本发明的目的是提供一种可用于小管径的并采用光线传感器对管线内的气液两相流介质的状态进行检测的管线中气液两相流状态检测方法及装置。

为了实现上述目的，本发明公开了一种管线中气液两相流状态检测方法，基于沿待测管线的管径方向布置的至少一对位于所述管线外部的光发射器和光接收器进行检测，所述检测方法包括：

实时检测并处理所述光接收器接收到的光信号，以生成反应光信号强度变化率的功率谱；

获取所述功率谱中当前时刻的峰值；

将所述功率谱的峰值与预设的第一阈值比较，并根据比较结果确定所述管线中当前气液两相流介质的存在状态。

较佳地，所述确定所述管线中当前气液两相流介质的存在状态包括：

确定所述功率谱的峰值是否小于第一阈值，如果是，

则，所述管线内不存在气液两相流介质，或所述管线内的气液两相流介质不流动；

如果否，则所述管线内存在流动的气液两相流介质。

较佳地，当所述功率谱的峰值不小于所述第一阈值时，计算前后时刻的功率谱的峰值的频率差，将该频率差与预设的第二阈值比较，并根据比较结果确定所述管线中当前气液两相流介质的流动状态。

较佳地，所述确定所述管线中当前气液两相流介质的流动状态包括：

确定前后时刻的功率谱的峰值的频率差是否大于第二阈值，如果是，

则，所述管线中的气液两相流介质流动不均匀；

如果否，则所述管线中的气液两相流介质流动均匀。

较佳地，对所述光信号进行处理的方法包括：

将所述光信号转换为与其光强度相适配的电信号；

对所述电信号进行噪声去除和分解处理，以得到目标信号；

通过功率谱估计的方法估计所得到的一系列所述目标信号的峰值频率，以得到所述功率谱。

较佳地，所述第一阈值为a，e-5＜a＜e-3。

较佳地，所述第二阈值为b，b≤1Hz。

较佳地，所述管线的外径为c，内径为d，3mm≤c≤5mm，1.5mm≤d≤3mm。

较佳地，所述光发射器包括LED灯。

较佳地，所述气液两相流介质包括空气和油液的混合介质。

较佳地，所述光发射器发出的光线被所述管线的吸收率小于被所述油液的吸收率。

较佳地，所述光发射器发出的光线为波长为923nm的红外光。

本发明还公开一种管线中气液两相流状态检测装置，其包括：

沿待测管线的管径方向布置的至少一对位于所述管线外部的光发射器和光接收器；

信号处理模块，实时检测并处理所述光接收器接收到的光信号，以生成反应光信号强度变化率的功率谱；

第一比较模块，用于比较所述功率谱的峰值与预设的第一阈值的大小；

第一确定模块，用于根据所述第一比较模块的比较结果确定所述管线中当前气液两相流介质的存在状态。

较佳地，还包括计算模块、第二比较模块和第二确定模块；

所述计算模块，用于计算在所述功率谱的峰值不小于第一阈值的情况下，前后时刻的功率谱的峰值的频率差；

所述第二比较模块，用于将所述频率差与预设的第二阈值比较；

所述第二确定模块，用于根据所述第二比较模块的比较结果确定所述管线中当前气液两相流介质的流动状态。

较佳地，所述信号处理模块中设置有光电转换模块、噪声处理模块以及功率谱生成模块；

所述光电转换模块，用于将所述光接收器接收到的光信号转换为与其光强度相适配的电信号；

所述噪声处理模块，用于对所述电信号进行噪声去除和分解处理，以得到目标信号；

所述功率谱生成模块，用于通过功率谱估计的方法估计所得到的一些列所述目标信号的峰值频率，以得到所述功率谱。

较佳地，所述第一阈值为a，e-5＜a＜e-3。

较佳地，所述第二阈值为b，b≤1Hz。

较佳地，所述管线的外径为c，内径为d，3mm≤c≤5mm，1.5mm≤d≤3mm。

较佳地，所述光发射器包括LED灯。

较佳地，所述气液两相流介质包括空气和油液的混合介质。

较佳地，所述光发射器发出的光线为波长为923nm的红外光。

本发明还公开一种管线中气液两相流状态检测系统，其包括：

一个或多个处理器；

存储器；

以及一个或多个程序，其中一个或多个程序被存储在所述存储器中，并且被配置成由所述一个或多个处理器执行，所述程序包括用于执行如上所述的管线中气液两相流状态检测方法的指令。

本发明还公开一种计算机可读存储介质，其包括计算机程序，所述计算机程序可被处理器执行以完成如上所述的管线中气液两相流状态检测方法。

与现有技术相比，本申请通过光线传感器也即沿待测管线的管径方向布置的一对位于管线外部的光发射器和光接收器来检测管线中气液两相流介质的状态，并根据气液两相流在管线中的特有流型对检测到的数据进行处理，以得到反应光信号强度变化率的功率谱，并将该功率谱的峰值与预设的第一阈值比较，根据比较结果即可确定出该管线中当前气液两相流介质的存在状态；由此可知，本申请通过反应光信号强度变化率的功率谱来表征管线中的介质流状态，无论管径大小，也即，即使极小管径的管线，管线中的气液两相流介质的存在状态的变化也会反应在功率谱中，从而有效提高通过光线传感器对小管径中气液两相流介质状态检测的精准度。

附图说明

图1为本发明实施例中对管线进行光学检测的布置结构示意图。

图2为本发明其中一实施例中检测方法流程图。

图3为本发明另一实施例中检测方法流程图。

图4为本发明实施例中其中一功率谱图，其中管内介质不流动。

图5为本发明实施例中另一功率谱图，其中管内不存在介质。

图6为本发明实施例中另一功率谱图，其中管内介质正常流动。

图7为本发明实施例中对油气润滑管线进行检测时采集的功率谱图，其中油液流动均匀。

图8为本发明实施例中对油气润滑管线进行检测时采集的功率谱图，其中油液流动不均匀。

图9为本发明实施例中对油气润滑管线进行检测时采集的功率谱图，其中油液不流动。

图10为本发明实施例中对油气润滑管线进行检测时采集的功率谱图，其中管线内无油液。

图11为本发明实施例中油气混合液的流型图。

图12为本发明其中一实施例中检测装置结构示意图。

图13为本发明另一实施例中检测装置结构示意图。

具体实施方式

为详细说明本发明的技术内容、构造特征、所实现目的及效果，以下结合实施方式并配合附图详予说明。

本实施例公开了一种气液两相流状态检测方法，以对通过管线进行输送的气液两相流介质的状态进行检测，该气液两相流介质为导入管线中的气体和液体的混合介质，如油气润滑管中的油液和空气的混合油液，通过该检测方法，可有效检测油气润滑管中的油气混合油液的状态。本实施例中的检测方法基于光学传感器进行检测，具体地，如图1，沿待测管线的管径方向布置至少一对位于管线外部的光发射器10和光接收器11，光发射器10用于向管线中发射检测光，该检测光沿管线的管径方向透过管线，光接收器11用于接收从管线的管径方向射出的光信号。

为方便理解本实施例中的检测方法，首先基于气液两相流在管线中的特有流型对入射光线的影响予以说明。如图11，对于管线中的气液两相流介质来说，以油液空气混合介质为例，油膜呈周期性断裂，所以，当在管线的某一侧壁处经由管径方向投射检测光时，出射光也应呈周期变化。且，在一个周期内，随着油液向前流动，管线底部油膜的扰动波峰将率先通过检测区域，此时油膜厚度急剧增大，由于出射光强度与油液厚度成负相关，因此，导致出射光强度急剧减小。而后在检测区域液膜厚度逐渐减小的过程中，还会遭遇顶部油膜的扰动波波峰，因此对应的出射光强度将再次急剧减小，之后，就是油膜厚度逐渐减小到0的过程，至此一个周期结束。

基于上述气液两相流在管线中的特有流型对入射光线的影响的原理，如图1 和图2，本实施例中检测方法包括如下步骤：

S1：实时检测并处理光接收器11接收到的光信号，以生成反应光信号强度变化率的功率谱。功率谱，又叫功率谱密度函数，定义了单位频带内的信号功率。在本实施例中，对于任一时刻的功率谱来说，如图4，其横坐标表示信号变化频率，与管线中的介质流的流动速度正相关，纵坐标表示光信号强度变化率的幅值，反应光信号强度的变化情况，其大小也与管线中的介质流的流动速度正相关，也即，管线中的介质流的流动速度越快，光接收器11接收到的光信号强度变化越快，相应地，光信号强度变化率的幅值越大，相反地，管线中的介质流的流动速度越慢，光接收器11接收到的光信号强度变化越慢，相应地，光信号强度变化率的幅值越小。

S2：获取功率谱中当前时刻的峰值。

S3：将功率谱的峰值与预设的第一阈值比较，并根据比较结果确定管线中当前气液两相流介质的存在状态。

具体地，管线中的气液两相流介质的存在状态包括三种，一种状态为管线内不存在气液两相流介质，另一种状态为管线内存在不流动的气液两相流介质，再一种状态为，管线内存在流动的气液两相流介质。对此，确定管线中当前气液两相流介质的存在状态的方法为：

S30：确定功率谱的峰值是否小于第一阈值，如果是，则进入S31，如果否，则进入S32。

S31：管线内不存在气液两相流介质，或管线内的气液两相流介质不流动。

S32：管线内存在流动的气液两相流介质。

上述实施例，通过对光接收器11接收到的光信号的处理，生成反应光信号强度变化率的功率谱，并根据该功率谱中的有效信号的峰值与预设的第一阈值的比较来确定当前管线中的气液两相流的状态。当功率谱中的峰值较小(小于第一阈值)时，表示管线内的介质对检测光强度变化的影响几乎可忽略，因此，可确定管线内不存在气液两相流介质(如图5)，或管线内的气液两相流介质不流动(如图4)。当功率谱中的峰值较大(不小于第一阈值)时，表示管线内的介质对检测光强度变化的影响较大，因此，可确定管线内存在流动的气液两相流介质(如图6)。

在另一实施例中，如图3，当功率谱的峰值不小于第一阈值时，上述检测方法还包括如下步骤：

S4：计算前后时刻的功率谱的峰值的频率差。

S5：将该频率差与预设的第二阈值比较，并根据比较结果确定管线中当前气液两相流介质的流动状态。

在功率谱中，由于功率谱的峰值的频率与气液两相流介质的流动速度相关，因此，当管线中气液两相流介质流速恒定时，理论上前后两时刻的功率谱的峰值的频率差为零，相反地，管线中气液两相流介质流速变化越大，前后两时刻的功率谱的峰值的频率差相应越大。因此，具体地，本实施例中确定管线中当前气液两相流介质的流动状态的方法包括：

S50：确定前后时刻的功率谱的峰值的频率差是否大于第二阈值，如果是，则进入S51，如果否，则进入S52。

S51：管线中的气液两相流介质流动不均匀；

S52：管线中的气液两相流介质流动均匀。

采用上述实施例公开的检测方法对油气润滑管线中的油液和空气混合介质进行检测时，通过光接收器11实时接收光发射器10经由管线的径向射出的光信号，并将该光信号转换为相应的电信号，然后，采用该电信号生成获得反应光信号强度变化率的功率谱。接着，从该功率谱中获取有效信号的峰值，然后，将该峰值与第一预设值比较，如果该峰值小于第一预设值，则给出如下检测结果：管线内不存在油液(如图5、图10)，或管线内的油液不流动(如图4、图 9)。如果该峰值不小于第一预设值，则给出如下检测结果：管线内存在流动的油液(如图6)。另外，当确定管线内存在流动的油液时，从功率谱中取出一组前后时刻的峰值的频率差，并确定该频率差是否大于第二预设值，如果是，则给出如下检测结果：管线内油液流动不均匀(如图8)，如果否，则给出如下检测结果：管线内油液流动均匀(如图7)。由此，工作人员根据检查结果即可确定出油气润滑管线当前的工作状态，当出现故障时，如管线内没有油液，及时做出应对措施。

另一实施例中，对光信号进行处理的方法包括：

首先，将光信号转换为与其光强度相适配的电信号。其次，对电信号进行噪声去除和分解处理，以得到目标信号。然后，通过功率谱估计的方法估计所得到的一系列目标信号的峰值频率，以得到功率谱。

具体地，得到电信号后，首先去基线漂移，去基线漂移主要是为了去除信号中直流分量的干扰。然后进行小波包信号分解，小波分解则是将信号分解为多个频段，通过计算各个频段信号与原信号的互相关系数，来找出实际信号和噪声分布的频段，然后再结合小波阈值降噪进行分组处理，接着，对处理后的各频段进行小波重构。最后，通过功率谱估计的方式估计信号的峰值频率，以得到功率谱。

另一实施例中，第一阈值为a，e-5＜a＜e-3。具体地，a优选为e-4。

再一实施例中，关于第二阈值的大小选择，可在管线形成稳定气液两相流的情况下，分别截取每个气压工况下5个时刻的信号波形，重复数据处理流程，提取峰值频率如下表1-1所示。从中可看出，在管线稳定运行的过程中，信号的峰值频率变化幅度很小，统计各气压下峰值频率的极差，可以看出最大不会超过1Hz。因此，第二阈值为b，b≤1Hz。

表1-1

根据上述实施例公开的检测方法的原理可知，其可应用于极小管径的管线，本实施例中，管线的外径为c，内径为d，3mm≤c≤5mm，1.5mm≤d≤3mm。具体地，c为4mm，d为2.5mm。另外，对于本实施例公开的小管径的管线，由于检测光的路径较短，因此，检测光通过管线的过程中，即使有些许发散，也可被光接收器11有效接收，所以，采用LED灯作为光发射器10，相比激光源，有效降低设备成本。

再一实施例中，为取得较好的检测效果，光发射器10发出的光线被管线的吸收率小于被油液的吸收率，从而可有效区分管线中油液混合介质的状态。具体地，光发射器10发出的光线优选为波长为923nm的红外光。

如图12，本发明还公开一种管线中气液两相流状态检测装置，其包括：

沿待测管线的管径方向布置的至少一对位于管线外部的光发射器10和光接收器11、信号处理模块12、第一比较模块13以及第一确定模块14。

光发射器10用于向管线中发射检测光，光接收器11用于接收从管线的管径方向射出的光信号。

信号处理模块12，用于实时检测并处理光接收器11接收到的光信号，以生成反应光信号强度变化率的功率谱。

第一比较模块13，用于比较功率谱的峰值与预设的第一阈值的大小。

第一确定模块14，用于根据第一比较模块13的比较结果确定管线中当前气液两相流介质的存在状态。

进一步地，如图13，该检测装置还包括计算模块15、第二比较模块16和第二确定模块17。

计算模块15用于计算在功率谱的峰值不小于第一阈值的情况下，前后时刻的功率谱的峰值的频率差。

第二比较模块16用于将频率差与预设的第二阈值比较。

第二确定模块17用于根据第二比较模块16的比较结果确定管线中当前气液两相流介质的流动状态。

更进一步地，信号处理模块12中设置有光电转换模块120、噪声处理模块 121以及功率谱生成模块122。

光电转换模块120，用于将光接收器11接收到的光信号转换为与其光强度相适配的电信号。

噪声处理模块121，用于对电信号进行噪声去除和分解处理，以得到目标信号。

功率谱生成模块122，用于通过功率谱估计的方法估计所得到的一些列目标信号的峰值频率，以得到功率谱。

另外需要说明的是，上述实施例中的检测装置的工作原理和使用方式详见上述检测方法实施例，在此不再赘述。

本发明还公开另一种气液两相流状态检测系统，其包括一个或多个处理器、存储器以及一个或多个程序，其中一个或多个程序被存储在所述存储器中，并且被配置成由所述一个或多个处理器执行，所述程序包括用于执行如上所述的气液两相流状态检测方法的指令。处理模块可以采用通用的中央处理模块 (Central Processing Unit，CPU)，微处理模块，应用专用集成电路(Application Specific Integrated Circuit，ASIC)，或者一个或多个集成电路，用于执行相关程序，以实现本申请实施例的气液两相流状态检测系统中的模块所需执行的功能，或者执行本申请方法实施例的气液两相流状态检测方法。

本发明还公开一种计算机可读存储介质，其包括计算机程序，所述计算机程序可被处理器执行以完成如上所述的气液两相流状态检测方法。该计算机可读存储介质可以是计算机能够存取的任何可用介质或者是包含一个或多个可用介质集成的服务器、数据中心等数据存储设备。该可用介质可以是只读存储器 (read-onlymemory，ROM)，或随机存取存储器(random access memory，RAM)，或磁性介质，例如，软盘、硬盘、磁带、磁碟、或光介质，例如，数字通用光盘(digital versatile disc，DVD)、或者半导体介质，例如，固态硬盘(solid state disk， SSD)等。

本申请实施例还公开了一种计算机程序产品或计算机程序，该计算机程序产品或计算机程序包括计算机指令，该计算机指令存储在计算机可读存储介质中。电子设备的处理器从计算机可读存储介质读取该计算机指令，处理器执行该计算机指令，使得该电子设备执行上述气液两相流状态检测方法。

以上所揭露的仅为本发明的优选实施例而已，当然不能以此来限定本发明之权利范围，因此依本发明申请专利范围所作的等同变化，仍属本发明所涵盖的范围。

Claims (19)

1.一种管线中气液两相流状态检测方法，其特征在于，基于沿待测管线的管径方向布置的至少一对位于所述管线外部的光发射器和光接收器进行检测，所述检测方法包括：

实时检测并处理所述光接收器接收到的光信号，以生成反应光信号强度变化率的功率谱；

获取所述功率谱中当前时刻的峰值；

将所述功率谱的峰值与预设的第一阈值比较，并根据比较结果确定所述管线中当前气液两相流介质的存在状态；确定所述管线中当前气液两相流介质的存在状态包括：

确定所述功率谱的峰值是否小于第一阈值，如果是，

则，所述管线内不存在气液两相流介质，或所述管线内的气液两相流介质不流动；

如果否，则所述管线内存在流动的气液两相流介质；

当所述功率谱的峰值不小于所述第一阈值时，计算前后时刻的功率谱的峰值的频率差；

确定前后时刻的功率谱的峰值的频率差是否大于第二阈值，如果是，

则，所述管线中的气液两相流介质流动不均匀；

如果否，则所述管线中的气液两相流介质流动均匀。

2.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，对所述光信号进行处理的方法包括：

将所述光信号转换为与其光强度相适配的电信号；

对所述电信号进行噪声去除和分解处理，以得到目标信号；

通过功率谱估计的方法估计所得到的一系列所述目标信号的峰值频率，以得到所述功率谱。

3.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，所述第一阈值为a，e-5＜a＜e-3。

4.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，所述第二阈值为b，b≤1Hz。

5.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，所述管线的外径为c，内径为d，3mm≤c≤5mm，1.5mm≤d≤3mm。

6.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，所述光发射器包括LED灯。

7.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，所述气液两相流介质包括空气和油液的混合介质。

8.根据权利要求7所述的管线中气液两相流状态检测方法，其特征在于，所述光发射器发出的光线被所述管线的吸收率小于被所述油液的吸收率。

9.根据权利要求1所述的管线中气液两相流状态检测方法，其特征在于，所述光发射器发出的光线为波长为923nm的红外光。

10.一种管线中气液两相流状态检测装置，其特征在于，包括：

沿待测管线的管径方向布置的至少一对位于所述管线外部的光发射器和光接收器；

信号处理模块，实时检测并处理所述光接收器接收到的光信号，以生成反应光信号强度变化率的功率谱；

第一比较模块，用于比较所述功率谱的峰值与预设的第一阈值的大小；

第一确定模块，用于根据所述第一比较模块的比较结果确定所述管线中当前气液两相流介质的存在状态；

还包括计算模块、第二比较模块和第二确定模块；

所述计算模块，用于计算在所述功率谱的峰值不小于第一阈值的情况下，前后时刻的功率谱的峰值的频率差；

所述第二比较模块，用于将所述频率差与预设的第二阈值比较；

所述第二确定模块，用于在所述功率谱的峰值不小于所述第一阈值时，如果前后时刻的功率谱的峰值的频率差大于所述第二阈值，输出第一确定结果；如果前后时刻的功率谱的峰值的频率差不大于所述第二阈值，输出第二确定结果；所述第一确定结果为：管线中的气液两相流介质流动不均匀；所述第二确定结果为：管线中的气液两相流介质流动均匀。

11.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述信号处理模块中设置有光电转换模块、噪声处理模块以及功率谱生成模块；

所述光电转换模块，用于将所述光接收器接收到的光信号转换为与其光强度相适配的电信号；

所述噪声处理模块，用于对所述电信号进行噪声去除和分解处理，以得到目标信号；

所述功率谱生成模块，用于通过功率谱估计的方法估计所得到的一系列所述目标信号的峰值频率，以得到所述功率谱。

12.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述第一阈值为a，e-5＜a＜e-3。

13.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述第二阈值为b，b≤1Hz。

14.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述管线的外径为c，内径为d，3mm≤c≤5mm，1.5mm≤d≤3mm。

15.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述光发射器包括LED灯。

16.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述气液两相流介质包括空气和油液的混合介质。

17.根据权利要求10所述的管线中气液两相流状态检测装置，其特征在于，所述光发射器发出的光线为波长为923nm的红外光。

18.一种管线中气液两相流状态检测系统，其特征在于，包括：

一个或多个处理器；

存储器；

以及一个或多个程序，其中一个或多个程序被存储在所述存储器中，并且被配置成由所述一个或多个处理器执行，所述程序包括用于执行如权利要求1至9任一项所述的管线中气液两相流状态检测方法的指令。

19.一种计算机可读存储介质，其特征在于，包括计算机程序，所述计算机程序可被处理器执行以完成如权利要求1至9任一项所述的管线中气液两相流状态检测方法。