递归特性分析在油气水三相流中的应用

50石油仪器PETROLEUMINSTRUMENTS2011年12月方法研究递归特性分析在油气水三相流中的应用王冬梅黄春辉谭1大庆油田测试技术服务分公司黑龙江大庆平王延军李雷2燕山大学信息科学与工程学院河北秦皇岛摘要先用CC算法确定非线性时间序列相空间重构的最佳嵌入维数和最优时间延迟，再利用递归图对油气水三相流各流型进行表征、利用递归定量指标对递归图进行量化分析。试验结果表明水基分散泡状流向水基分散段塞流转变中，递归图的细纹理特征可简单直观快速的识别流型，递归定量指标对分相流量变化能较好的定量分析。因此，递归特性分析技术具有较强的识别能力，为三相流流型辨识提供了一种有效方法。关键词油气水三相流；CC算法；递归图；递归定量指标中图法分类号TE33文献标识码B文章编号100491342011060050030引言油气水三相流动广泛存在于化工、核能、冶金及油气输运过程中，其相问相互作用和相对运动特性复杂，致使流型转化非线性动力学机制尚未完全研究明白。测井传感器响应及资料解释严重依赖对流动特性的了解，所以，开展油气水三相流流型转化的非线性动力学特性研究具有重要意义。自然界中状态的递归性很早就被发现并被研究，但是人们对递归行为的动力学模型在高维相空间M3上缺乏合适的分析手段。ECKMANJ提出的递归图RECURRENCEPLOT和ZBILUTJ提出递归定量分析法RECURRENCEQUANTIFICATIONANALYSIS可以对非平稳时间序列进行处理，近年来在多相流检测3J和生物物理J等领域得到了应用，已成为一种新兴的混沌时间序列分析工具。本文对垂直上升管中油气水三相流不同流型下的递归图结构及递归定量指标进行考察，研究其随流量变化的规律。试验结果表明，递归特性分析技术用于复杂动力学过程的模式识别，具有快速、直观的特点，为油气水三相流流型识别提供了一种有效方法。L递归图和递归定量指标11递归图递归图定义是一种通过高维相空间轨迹在二维平面投影递归行为，来研究高维空间动力学特性的方法。算法如下1在混沌时间序列分析中，相空间重构是必不可少的一步。本文依据CC算法_5同时计算出时间序列相空间重构的最佳嵌入维数M和最优时间延迟。2基于TAKENS嵌入定理，重构原始时间序列X，2，的相空间，相空间向量为XIXI，3CIR，一1RII，N一MIT13计算时间序列的递归值RJ0EL_I，J1，凡一一120为HEAVISIDE函数，NSTD为阈值，STDX为标准差，口为阈值系数，一般取A01025。4画递归图定义递归矩阵R1的位置为递归点，在MM坐标平面相应位置，J处描点，这样就可得到时间序列经过相空间重构后的递归图，因LL置墨L10，故在重构后的递归图中总可以看到一条主对角线。12递归定量指标递归定量分析法是通过定量的表征递归图中基本图形和线段的分布，来达到对时间序列进行定量分析的目的。不同的递归量化参数描述了系统不同方面的动力学递归行为，本文研究如下4种递归定量指标。1递归率RR递归图中所有递归点占总覆盖点数的比率，表征相空间中相邻轨道彼此靠近的相空间点矢量所占的比例。基金项目秦皇岛市科学技术研究与发展计划200901A029201001A062第一作者简介王冬梅，女，1974年生，助理工程师，在大庆油田测试技术服务分公司，主要从事仪器仪表检测方面的研究。邮编1631002011年第25卷第6期王冬梅等递归特性分析在油气水三相流中的应用RRRI,JM3I，12确定性DET构成45。方向线段递归点数占总递归点数的比率，表征递归点中确定性信号所占的比重。M一1MIDETZP1R4LIMF一式4中，PZ是长度为Z的线段数量，只有Z大于预先给定的下限Z时才开始计数，且主对角线递归点不统计。3平均对角线长度L对角线方向线段长度的加权平均值，表征了系统的平均周期。M一1MLZP1PZ5ZLL1M4熵ENTR递归图主对角线方向线段长度的SHANNON熵，表征了系统确定性结构的复杂程度。MLM一1ENTR一【PZP】FF盯一JIN【PZPZ6F2油气水三相流型递归特性分析实验是在大庆油田测试技术服务分公司三相流动模拟试验装置上进行的，透明垂直模拟井筒长8ITI，内径为125MM。实验参数范围油水总流量20IND601T1D，气流量3M0D50MD，含水率KW75。对垂直上升管中油气水三相流流型划分方法见文献E6。依据CC算法得到各信号的时间延迟和嵌入维数，结果表明，最小嵌入维数为24，最优时延为114。试验发现嵌入维数的大小只影响递归点的稀疏，不影响递归图的拓扑结构，所以为了在同一维度空间中观察动力系统的递归特性，本文嵌入维数均取M3，而最优时延R则分别选取，考虑到工况数据中噪声的存在，选择阈值系数A025，各流型递归图如图1所示。水基分散泡状流如图1A所示，由于流道中存在大量小气泡，并伴有气泡产生和破碎，且小气泡能量很小，泡群运动轨迹随机可变，递归图表现为分散的孤立无序点状细纹理结构。水基分散块状流如图1B所示，随着气相流量比例增加，气相与液相间相互搅动加剧，使得流型发生转变，递归图表现为主对角线结构大多较短，长度分布集中，且能隐约看见黑色块状结构。水基分散段塞流如图1C所示，由于存在明显的气塞与液塞周期性交替变化，一些小气泡积聚成塞状大气泡，其直径近似于管内径，长度相当于数倍管径，大气泡之间由液相隔开，油滴以泡状存在于连续相水中，递归图的主对角线结构长短变化较为丰富，长度分布较广，且伴有断续的黑色块，验证了此流型具有一定周期信号特征。因此，递归图对各流型显示了很好的纹理特征，可简单直观的识别流型。囊。，一；彰享。BQ_FL2QG27CHMAW图1不同流型的递归图提取各信号递归图的递归定量指标，结果如图2低，流型的动力学特性变得相对简单。但当液体流量所示。从图2中得出，低液体流量下的指标值总体比不变时，随着气流量的增加，各指标值趋于增大，系统高液体流量的大，且各流型下值能较好区分开。例如相空间的靠近程度、确定性成份和平均周期都趋于增对于递归率水基分散泡状流处于O017100312，大，信息损失速率增大，混乱程度增加，流型的动力学水基分散块状流处于0038300492，水基分散段特性变得复杂，表明分散泡状流向分散段塞流转变时塞流处于0074601281。水基分散泡状流向水基递归特性逐渐增强。因此，各递归定量指标为三相流分散段塞流的转变过程中，递归特征量值整体趋于增流型辨识提供了有效方法。大。当气流量不变时，随着液体流量的增加，各指标值，趋于减小，表明系统相空间的靠近程度、确定性成份和一平均周期都趋于减小，信息损失速率降低，混乱程度降1油气水三相流水基分散泡状流递归图，呈现类石油仪器PETR0LEUMINSTRUMENTS2011年12月图2不同流型下的递归定量指标似于随机信号的分散孤立无序点状细纹理结构；水基2ZBILUTJP，WEBBERCLJREMBEDDINGSANDDELAYSAS分散块状流递归图，主对角线结构大多较短，长度分布DERIVEDFROMQUANTIFICATIONOFRECURRENCEPLOTSJPHYS集中，且能隐约看见黑色块状结构；水基分散段塞流递LETA，1992,7【1归图，主对角线结构长短变化较为丰富，长度分布较3金宁德，郑桂波，陈万鹏气液两相流电导波动信号的广，且伴有断续的黑色块。递归图对各流型显示了很混沌递归特性分析化工学报，2007，5好的纹理特征，可简单直观、快速的识别流型。4TAN，TZ，ZHE“一“GPE“D。一SH。NGWANG2递归率、确定性、平均对角线长度及熵对分相流QAN曲。ATI。ALYI。B。D。HA0G锄二EP。NTATI。量变化能较好的定量分析，有助于理解流型转化非线JJOURNALOFTH。TIALBI。1。，2009，257性动力学机制。递归特性分析技术为油气水三相流流ES吕金虎，陆君安，陈士华混沌时间序列分析及其应用型辨识提供了一种有效方法。EM2武汉大学出版社，2001垒垃士6于莉娜，李英伟，郭学涛垂直上升管内油气水三相流丐X瓢数值模拟方法研究J油气田地面工程，2010，2931ECKMANJP，KAMPHORSTS0，RUELLEDRECURREN