本发明涉及立管集输系统气液两相流流型判别系统和方法，属于石油工程多相流检测
技术领域：
：。
背景技术：
：：随着陆地油气资源日益枯竭，海洋油气资源的开发力度日渐加大。海洋俨然变成了我国油气资源开发的重要潜在区域。在海洋油气开采中，考虑到自然环境的恶劣，处于经济和环保的考虑，往往不宜就地安装气液分离器或同时铺设两条管线分别输送油气，所以常采用气液混输方式，将原油、天然气和水等产物自海底举升到海面浮式装置。集输立管系统是浮式生产系统用于向船或平台传送流体的基本装置，同时也是薄弱、易损的结构之一。集输立管系统的流型与立管系统油气安全输送密切相关。当流型为严重段塞流等不稳定流型时时，立管出口处气液的间断出流，产量剧烈波动，对下游设备产生影响，使得在下游的分离器中产生断流或溢流现象。液气喷发阶段产生强烈的振动会加剧腐蚀管道系统，特别是立管部分，管道的强烈振动会加剧管道接头、支柱、立管系统的破坏，对其安全和稳定造成影响，也会加剧增压设备的汽蚀。气体从立管出口处喷发时会产生节流效应，使管内温度降低，产生结蜡和水合物，堵塞管道。以严重段塞流为代表的不稳定流型危害影响生产效率和生产安全，导致产量、压力剧烈波动，管线剧烈振动。因此，必须对流型进行监控，保证生产安全，对危害流型进行控制。现有的流型识别专利包括如下：中国发明专利(公布号cn102878431b)、中国发明专利(公布号cn104807589a)和中国发明专利(公布号cn105823617a)，分别采用管线入口的单个压力传感器、管线顶部的压力传感器和集输立管系统气液两相流流型三个不同位置的压差传感器来进行流型判别，上述专利虽然能对立管的流型进行大类的识别，但不能对立管的流型进行更细类的识别，无法提高海上勘探的安全性和高效性。技术实现要素：为了解决现有技术的不足，本发明提供了一种一种集输立管系统气液两相流流型判别方法，采用以下步骤：(1)采集立管顶部与底部的压差信号；(2)将所述压差信号做无量纲化处理；(3)采用统计学和时频分析所述压差无量纲化信号的统计学和时频特征，所述统计学特征包括压差信号的概率密度函数、累积分布函数，所述时频域特征包括压差信号的时域值和功率谱密度；(4)根据统计学和时频特征进行流型的判别。进一步的，步骤(1)结束后，首先采用小波算法过滤所述压差信号中的噪声信号，然后执行步骤(2)。进一步的，当所述压差信号的时域值在设定时间段内的波动幅值大于立管内充满液体工质时静压差的设定百分比时，判别此时流型为不稳定流型，否则流型为稳定流型。进一步的，当流型为不稳定流型时，且所述压差信号的概率密度函数存在波峰，则判别此时流型为严重段塞流型i或严重段塞流型iii，否则流型为严重段塞流型ii或段塞流型。进一步的，当流型为严重段塞流型i或严重段塞流型iii时，且所述压差信号的累积分布函数中横坐标小于第一设定值且大于第二设定值，则判别此时流型为严重段塞流型i，否则流型为严重段塞流型iii。进一步的，当流型为严重段塞流型ii或段塞流型时，且所述压差信号的功率谱密度存在一阶峰值，则判别此时流型为严重段塞流型ii，无波峰为段塞流型。进一步的，当所述压差信号的概率密度函数无波峰时，则判别此时流型为波动气泡流型，否则流型为稳定气泡流型或者乳沫状流型。进一步的，当流型为稳定气泡流型或者乳沫状流型时，且所述压差信号的功率谱密度无波峰时，则判别此时流型为乳沫状流型，否则流型为稳定气泡流型。本发明还提出了一种基于上述方法的集输立管系统气液两相流流型判别系统，包括设置于立管顶部与底部的压力测量装置，压力测量装置依次连接有数据采集模块、数据预处理模块、数据分析模块和流型判别模块：所述数据采集模块，用于采集立管顶部与底部的压差信号；所述数据预处理模块，用于将所述压差信号做无量纲化处理；所述数据分析模块，用于对无量纲的压差信号进行统计学和时频分析，得到统计学和时频特征；所述流型判别模块，基于所述统计学和时频特征对立管进行流型判别；进一步的，本系统还包括流型输出模块，输出并展现集输立管系统气液两相流流型的流型判别结果。进一步的，本系统还包括数据过滤模块，用于采用小波算法过滤所述压差信号中的噪声信号。与现有技术相比，本发明的有益效果是：本发明基于立管顶部与底部的压差信号进行流型判别，通过对采集到的压差信号进行预处理，然后进行分析得到压差信号的pdf((probabilitydistributionfunction概率密度函数)、cdf(cumulativedistributionfunction累积分布函数)和psd(powerspectraldensity功率谱密度)特征。通过压差信号曲线和相应的统计学特征可以对集输立管系统的流型进行划分，且可以得到较细致的流型，提供了一种集输立管流型判别的新思路，对海洋立管系统生产安全有重大意义。附图说明构成本申请的一部分的说明书附图用来提供对本申请的进一步理解，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。图1是集输立管系统气液两相流型判别系统示意图；图2是集输立管系统的各种流型类别图；图3是典型集输流型立管中气液分布示例图；图4是不同流型无量纲压差信号曲线；图5是不同流型无量纲压差信号的pdf曲线；图6是不同流型无量纲压差信号的cdf曲线；图7是不同流型无量纲压差信号的psd曲线；图8是本发明的流型判别模型。图中，1-压差变送器，2-垂直立管，3-压差变送器，4-悬链式立管。具体实施方式应该指出，以下详细说明都是例示性的，旨在对本申请提供进一步的说明。除非另有指明，本文使用的所有技术和科学术语具有与本申请所属
技术领域：
：的普通技术人员通常理解的相同含义。需要注意的是，这里所使用的术语仅是为了描述具体实施方式，而非意图限制根据本申请的示例性实施方式。如在这里所使用的，除非上下文另外明确指出，否则单数形式也意图包括复数形式，此外，还应当理解的是，当在本说明书中使用术语“包含”和/或“包括”时，其指明存在特征、步骤、操作、器件、组件和/或它们的组合。本发明的一种典型的实施例如图1所示：一种集输立管系统气液两相流流型判别系统，包括设置于立管顶部与底部的压力测量装置，压力测量装置依次连接有数据采集模块、数据过滤模块、数据预处理模块、数据统计学分析模块和流型判别模块。一阶峰值的定义：高于某个设定值的极大值只有一个叫一阶波峰，极大值有两个就叫二阶波峰，一阶峰值即只有一个波峰。本实施例中，压力测量装置采用压差变送器，左图为一种悬链式立管，右图为一种垂直立管，采用压差变送器可获得立管顶部和底部之间的压差信号，数据采集模块具体采用数据采集卡获取上述压差信号，然后通过数据过滤模块对压差信号进行小波3层滤波，滤波后的压差信号送入数据预处理模块进行无量纲处理，具体的无量纲处理方法为：对上述压差信号除以测压装置充满液体时的静压差。然后对上述无量纲的压差信号进行分析，得到统计学和时频特征，即压差信号随时间变化曲线和相应的pdf((probabilitydistributionfunction概率密度函数)、cdf(cumulativedistributionfunction累积分布函数)和psd(powerspectraldensity功率谱密度)特征，通过对不同流型无量纲压差信号及统计学特征进行流型判别，输出流型结果。此种流型判别系统提供了一种立管系统流型判别的新思路，可以进行流型现场判别。最后将判别结果发送至流型输出模块，本实施例中的流型输出模块可采用彩色具有3d功能的显示器。本发明的再一实施例是一种集输立管系统气液两相流流型判别方法，在上述实施例的基础上，首先介绍一下压差信号的pdf((probabilitydistributionfunction概率密度函数)、cdf(cumulativedistributionfunction累积分布函数)和psd(powerspectraldensity功率谱密度)特征。其中概率密度函数表征瞬时幅值落在某指定范围内的概率，因此是幅值的函数，它随所取范围的幅值而变化；累积分布函数则是概率密度函数的积分，表示小于某个值的概率；功率谱密度则是表征不同频率信号的能量分布。概率密度函数(pdf)随机测得的数据x(t)落在x-x+δx范围内的概率为：上式为x(t)的概率分布函数。x(t)落在x-x+δx之间的概率为：称为x(t)的概率密度函数，其中tx为t时间内x(t)落在x-x+δx区间的总时间，对应的幅值在x-x+δx范围内的概率也可以用tx/t来表示。累积分布函数(cdf)累积分布函数是概率密度函数的积分，表示随机数值x(t)，小于某个值x的概率：c(x)＝p[x(t)≤x](3)功率谱密度(psd)采用周期图法计算压力信号的功率谱密度：它是把随机序列x(n)的n个观测数据视为一能量有限的序列，直接计算x(n)的离散傅立叶变换得出它在每一个频率上的谱值：然后再取其幅值的平方，并除以n，作为序列x(n)真实功率谱的估计：在利用上述统计学特征对集输立管系统的流行做出判别的基础上，我们先把集输立管系统的流型大致分为以下几个大类，如图2所示，通过实验室进行立管系统气液两相流实验，得到不同工况下立管系统的7种流型，分别为严重段塞流i(severesluggingi,ssi)、严重段塞流ii(severesluggingii,ssii)、严重段塞流iii(severesluggingiii,ssiii)、波动气泡流(unstablebubbleflow,ubf)、稳定气泡流(stablebubbleflow,sbf)、乳沫流(churnflow,cf)和段塞流(slugflow,sf)。其中ssi、ssii、ssiii和段塞流(sf)称为不稳定流(unstableflow,usf),波动气泡流(ubf)、稳定气泡流(sf)和乳沫流(cf)称为稳定流(stableflow,sf)。各流型时，立管系统中的典型的气液分布如图3所示。各流型的描述如下：(1)严重段塞流i(ssi)：发生时的气液流量都处于较低的水平，此时，下倾管内一直保持分层流，液塞间歇性堵管，使气体在下倾管内累积，当气体压力超过立管内液体的静压力时气体进入立管与液体一起喷发，形成的气泡在上升过程中不断膨胀，推动液体加速喷出。喷发结束后，少量液体回落，造成堵管，流动进入下一周期。(2)严重段塞流ii(ssii)：提高入口的气体流量可以观察到，当液塞头部未达到上倾管顶端时，底部气体就已经进入到上倾管中，有气体喷发和回落现象。(3)严重段塞流iii(ssiii)：提高入口液体的流量可以观察到，液体对气泡的携带作用变强，气体间歇进入上倾管有一定的周期性。(4)波动气泡流(ubf))：在严重段塞流iii的基础上继续增加折算液速，会使得气泡进入上倾管段变得剧烈，并且周期性消失。(5)稳定气泡流(sbf：在气相流量较小时，继续增加液体流量，气泡逐渐稳定连续的进入到上倾管段，此时水平管段持液率很大，流动为分层流或气团流。(6)乳沫状流(cf)：上倾管中的流动为乳沫状，水平管中的流型为气团流或段塞流。(7)段塞流(sf)：当气液流量很高时，下倾段为段塞流，液塞速度很快，对立管的冲击强烈。根据无量纲压差信号和对应的pdf、cdf、psd特征进行流型判别。流型判别模型如图5所示。其中的6条判别标准如下：cr1：将流型划分为稳定流(sf)和不稳定流(usf)，立管压差波动幅值大于立管内充满液体工质时静压差的20％，此时流型为不稳定流(usf)，否则流型为稳定流(sf)。cr2：根据pdf曲线是否存在波峰(此处pdf波峰定义：大于0.15的峰值)将不稳定流(usf)分为两类：(1)ssi、ssiii；(2)ssii、段塞流(sf)。若存在波峰则流型为ssi或ssiii，否者为ssii或段塞流(sf)。cr3：此标准用于区分ssi和ssiii,根据cdf曲线中渐变段横坐标来判别，若起始点横坐标(压差信号无量纲数)小于0.5且终止点横坐标大于0.9则为ssi，否者为ssiii。cr4：此标准用于区分ssii和段塞流(sf),通过psd特征来进行区分判别，psd特征曲线有一阶波峰(此处波峰值定义：峰值为极大值且大于6x10-7)为ssii,无波峰为sf。cr5：pdf无波峰(此处pdf波峰定义：大于0.15的峰值)则为波动气泡流(ubf),否则为sbf或cf。cr6：psd无峰值(此处波峰值定义：峰值为极大值且大于1.5x10-7)则为cf,否则为sbf。本发明基于立管顶部与底部的压差信号进行流型判别，通过对采集到的压差信号进行预处理，然后进行统计学和时频分析得到压差信号的pdf((probabilitydistributionfunction概率密度函数)、cdf(cumulativedistributionfunction累积分布函数)和psd(powerspectraldensity功率谱密度)特征。通过压差信号曲线和相应的统计学特征可以对集输立管系统的流型进行划分，且可以得到较细致的流型，提供了一种集输立管流型判别的新思路，对海洋立管系统生产安全有重大意义。以上所述仅为本申请的优选实施例而已，并不用于限制本申请，对于本领域的技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本申请的保护范围之内。当前第1页12当前第1页12