一种融合总流量与电导探针阵列信号的水平井持水率测量方法
【专利说明】一种融合总流量与电导探针阵列信号的水平井持水率测量 方法 【技术领域】
[0001] 本发明属于多相流检测技术领域，特别是针对中低产液水平井，提出的一 种融合总流量和电导探针阵列信号并结合非线性支持向量回归（Supportvector regression,SVR)模型进行水平井持水率测量的方法。 【【背景技术】】
[0002] 生产测井指在采油井、注水井、观察井等投产至报废的生产过程中进行的地球物 理测井，可分为工程技术测井、产层评价测井和生产动态测井。生产动态测井主要任务是对 产出井的流动剖面进行动态监测，以了解产层出力状况和流体组分含量，为油田挖潜与细 化管理提供依据。由于近年来石油开采的深入，国内大多数陆上注水垂直井已进入高含水 生产阶段，现有测井仪在分辨率与精度方面已难以满足测井要求。相比于垂直井，水平井可 大幅提高单井产油量以及油藏的采收率，在地面及井下受限制的油气现场进行大规模的应 用。目前我国水平井技术相对比较落后，常规测井仪多采用中心采样方法进行持水率的测 量，无法满足高精度的测量需要。因此，水平井中油水两相流持水率测量的研究对推动水平 井测井技术的发展有着非常重要的意义。
[0003]在油气田开发过程中，油水分相流量测量对监测各产层产出、估计油产量、指导油 井注采、提高采收率、节约能耗以及预测油井开发寿命等方面具有重要的作用。在井下难以 通过油水分离分别测量油和水的分相流量，通常是先测得油水两相流总流量和含水率，或 者先测得油和水的分相流速和持水率，然后计算出油和水的分相流量。由于油水分相流速 难以测量，国内一些油田通常采用测量总流量和持水率的方法近似得到油水分相流量。因 此，持水率是产出剖面测井中非常重要的多相流参数。持水率是指在井筒中任一截面内水 相面积所占总面积的百分比，含水率是指单位时间内流过井筒的水相体积流量占多相流总 体积流量的百分比。含水率也是非常重要的多相流参数，但是油水两相流型多变，且两相流 相间存在复杂的界面效应和滑脱现象，致使持水率往往并不等于含水率。因此，水平井持水 率的高精度测量仍然是极具挑战性的工作。
[0004]目前，水平井持水率的测量得到国家的高度重视，国家科技部先后组织开展了水 平井持水率测量的相关课题的研究工作。在"十二五"期间，立项了"油气测井重大技术与 装备"国家科技重大专项，并在子课题"基于电导法的两相流成像测井传感器优化设计及井 下电路设计"中对持水率测井仪样机的研制和测量方法进行了更深入地研究，先后在水平 模拟井和油气现场进行了实验，将持水率的测量精度控制在8%以内。国家知识产权局授 权了一些关于石油生产测井监测领域的发明专利。例如，授权的两项发明专利"一种基于 遗传算法的多环电极阵列传感器结构优化方法"（ZL201210544383.X)和"一种环形持水率 测井传感器阵列的结构优化方法"（ZL201010543247.X)都是对持水率测井仪的结构进行 了研究，并对探针阵列的位置布局和探针结构进行优化，提高了中低产液水平井的持水率/ 含水率测量和流型识别的精度。授权发明专利"一种基于电导探针阵列传感器的水平井多 参数估计方法"（ZL201310193498.3)提出了对电导探针阵列传感器的方位角、水平井的水 相电导率和油水界面高度的估计方法，可用于低产液水平井产出剖面的在线测量和动态监 测，而且测量下限低、精度高、可靠性好。其中，提出的水平井油水界面高度的估计方法用于 水平井的持水率测量，能提高水平井持水率的测量精度。
[0005] 国内外学者对持水率的测量方法积极开展了研究，并提出了大量的方法，如探针 法、射线法、电学层析成像法以及软测量方法等。探针法可分为电导探针法、光纤探针法和 电容探针法三类。电导探针法及电容探针法是通过测量探针处的流体导电性的电压响应值 来确定该点介质分布，进行持水率的测量。光纤探针法则是利用液相的折射率高于气相的 原理来区分液相和气相，通过测量多相流在棱镜上的反射光强度获得探针周围的介质分布 来进行持水率的测量。射线法是比较经过多相流体的射线的衰减程度来获得相持率。射线 层析成像可获取井筒某一横截面相持率的时变信息，实现流动可视化，而且测量精度高，但 是需要一个稳定可靠的射线源，且存在辐射安全问题。电学层析成像法通过对被测对象施 加电激励，检测其边界值的变化，并利用逆问题求解算法逆向推出被测物体内部电学特性 的参数分布，从而实时获得持水率信息。基于信息融合的软测量方法是将传统的传感器测 量技术与人工神经网络和支持向量机结合，该方法也被广泛应用于持水率的检测。
[0006] 传统的测井仪多是通过判断水平井中油水界面高度来测量持水率，这种方法测量 持水率的精度不仅受探针的结构和位置分布的影响，而且受油水两相流的流型及流速的影 响。目前多数的持水率测量方法是对测井仪探针阵列的结构进行优化，并结合人工神经网 络和支持向量机等软测量的方法提高测量精度。由于进行信息融合的主要成分是探针阵列 在油水两相流中的通断信息，限制了水平井持水率的测量精度。 【
【发明内容】
】
[0007] 本发明的目的是针对中低产液水平井提出一种融合总流量与电导探针阵列信号 并结合非线性SVR模型进行水平井持水率测量的方法，提高水平井持水率的测量精度。本 发明提出的持水率测量方法的总体实施方案如图1所示。可收缩式双环电导探针阵列传感 器结构如图2所示，双环电导探针阵列测井仪通过电导测量电路（24)记录各支电导探针的 输出，得到油水两相流的时变电学特性作为测量信号。将油水两相流总流量和电导探针阵 列信号作为非线性SVR模型的输入，在不改变测井仪探针结构的基础上大幅提高水平井持 水率的测量精度。
[0008] 为实现上述目的，本发明提供了一种融合总流量与电导探针阵列信号的水平井持 水率测量方法，采用以下技术方案：
[0009] -种融合总流量与电导探针阵列信号的水平井持水率测量方法，其特征在于，步 骤如下：
[0010] 第一步，建立由水平井油相扇形面弧长到水平井持水率的计算公式。
[0011] 在中低产液水平井中，油水两相流由于重力分异作用而明显出现相对稳定的分层 流动现象，密度较小的油相靠近管道上壁流动，密度较大的水相在管道底部流动，形成相对 稳定的界面高度。水平井油水两相分层流的横截面示意图如图3所示，棕红色代表油（31)， 蓝色代表水（35)，L代表油相扇形面的弧长（32)，r代表水平井内径半径（33)，Θ代表油 相扇形面的圆心角（34)。
[0012] 水平井持水率α是指在多相流管流中任一截面内水相面积占管道截面总面积的 百分比，即
[0014] 式中，圆心角Θ可表示为
[0016] 结合式[1]和式[2]可进一步将水平井持水率的计算公式表示为
[0018] 第二步，获取建立非线性SVR模型的训练样本，获取过程如下：
[0019] (a)在水平井油水两相流管道上，以△F为步长，从最低总流量F_到最高总流量 ?_逐步调节油水两相流的总流量，获得总流量设定值向量{Fi|i= 1，2,…，1}，其中，i为 总流量设定值的序号^为总流量设定值的个数七二？^?，？^?，？^，。、）'?^ 总流量设定值Fji= 1，2,…，I)下，以Δβ为步长，从最低含水率β_到最高含水率β_ 逐步调节油水两相流的含水率，获得含水率设定值向量{βiIj= 1，2,…，J}，其中，j为含 水率设定值的序号，J为含水率设定值的个数，β1=β_，β;=β_，ββ_+(j-l) Αβ，从而得到IXJ种总流量设定值和含水率设定值的组合；
[0020] (b)在每种总流量设定值和含水率设定值的组合下，分别测量24支电导探针 的电压输出，当输出电压V、j大于阈值电压Vth时，判定该探针在油中，将该探针的输出 用逻辑"0"表示，反之用逻辑"1"表示，因此获得电导探针阵列输出的向量为Rt(t= 1，2,…，IXJ)，其中，t为探针阵列输出向量序号，IXJ为探针阵列输出向量的个数，艮为 1行24列的行向量；同时利用高速相机（45)进行摄像，获取设定值组合下的油相扇形面的 弧长Lt(t= 1，2,…，IXJ);
[0021] (c)根据第一步建立的油相扇形面弧长到持水率的计算关系，利用第二步中 步骤（b)得到的油相扇形面弧长值Lt(t= 1，2,…，IXJ)获得水平井持水率at(t =1，2，一，IXJ)，并与总流量设定值Fji= 1，2，···，Ι)和探针阵列输出向量Rt(t= 1，2,…，IXJ)共同构成支持向量回归模型的训练样本；
[0022] 第三步，建立从总流量设定值和电导探针阵列信号到高速摄像获得的持水率参考 值的非线性SVR模型，过程如下：
[0023] (a)将第二步获得的IXJ组训练样本进行50次随机划分，得到各占80%和20% 的训练集和测试集，测试集和训练集的数据格式一样。训练集的一个样本记作
[0024] (爲；，:1^)..·"..1%!.:=尺L()，1jW1
[0025] 式中，xm表示1X25的行向量，即测试集的总流量设定值和24支电导探针阵 列信号，yn表示利用高速相机（45)获得的持水率参考值，1表示训练集的长度，即1 = 0· 8 · (IXJ) 〇
[0026] (b)建立非线性SVR模型。非线性SVR模型的建模过程如图5所示，其非线性超平 面回归函数f(x)的表达式为
[0027]f(χ) =ωTx+b [5]
[0028] 为确定其权重向量ω和偏置b(52)，需要求解正则化风险函数
[0033] 式中，C>0表示惩罚因子，ε>0表示给定精度，f表示松弛变量。为方便求解 公式[6],通过构建Lagrange函数，并引入高斯径向基核函数（54)
[0035]式中，σ表示核半径。将公式[6]写成对偶形式
[0039]式中，λι; 表示非负的Lagrange乘子（55)。解得非线性SVR模型的非线性超 平面回归函数为
[0043]用测试集样本进行交叉验证，从而建立高精度持水率预测的非线性SVR模型。
[0044] 第四步，利用第三步获得的非线性SVR模型用于水平井持水率的在线测量，即
[0045] 将涡轮流量计实际测得的总流量和探针阵列测得的持水率作为非线性SVR模型 的输入，预测更高精度的水平井持水率值。
[0046] 本发明专利提出的一种基于总流量与电导探针阵列信号并结