1.一种利用阵列流量测井获取水平井持水率的方法,其特征在于:它包括如下步骤;
①、根据套管内径、仪器外径、仪器在井筒中的位置和各个阵列流量探头的方位角,确定流体过流截面和哥哥阵列流量探头在过流截面的位置;
②、根据过流截面的几何结构,将过流截面以轴心为原点建立直角坐标系,沿横向x轴和纵向y轴,将井筒半径等分为n份,将过流截面划分为若干网格,对描述井筒内粘性流体的纳维-斯托克斯方程(即n-s方程)在过流截面的网格内进行有限差分;将井筒截面沿x轴和y轴上内径等分为n等份,井筒截面划分为若干网格,垂向和横向步长为
h=δx=δy=2r/n(1)
对n-s方程进行有限差分,差分形式为
③、首先,预设流体压力pz和过流截面油水分界面高度h’参数,过流截面油水分布示意图;针对过流截面划分的网格,采用高斯-塞德尔迭代算法计算每个网格点的流速v(i,j);
然后,根据步骤1)中每个阵列流量探头的位置,读取其位置点的流速vical(i=1,2,3,4,5,6);
最后,对比每个阵列流量探头位置处迭代计算的流体速度值vcali(i=1,2,3,4,5,6)与阵列流量探头的测量值vfi(i=1,2,3,4,5,6),建立最优化目标函数
式中v=(pz,h)t为一个列向量;f(v)、ci(v)均为v的函数。ci(v)为约束条件;
通过在约束条件内流体压力pz和过流截面油水分界面高度h’参数的调优变化,采用最小二乘优化算法对目标函数进行求解,只到模型计算的每个阵列流量探头位置出的流体速度值与探头的测量值接近,得到油水分界面的高度h;
④、分界面高度和持水率计算方法;
具体计算公式如下;由步骤③计算得到油水分界面高度h,由油水分界面高度h减去套管半径r得到油水分界面到井筒中心o的距离ho,则
式中aδ为油水分界面相对于井筒中心o的三角形面积;则油相所占截面积为
井筒过流截面的面积为
a=π(r2-r2)(7)
如果油水分界面位于仪器之上,则油相所占截面积为ao为,仪器全部在水中,则持水率为
如果油水分界面位于仪器中间,则油相所占截面积为ao为,仪器部分在油中,部分在水中,则
仪器在油相所占截面积为
则持水率为
如果油水分界面位于仪器之下,则油相所占截面积为ao为,仪器全部在油中,仪器在油中的面积为
ato=πr2(12)
则持水率为
其中:r为井筒半径;仪器半径为r;ho为油水分界面相对于井筒中心o的高度(在井筒中心o点之上为正值,之下为负值);α为油水分界面与井筒交点相对于井筒中心o的夹角;θ为油水分界面与仪器外边界交点相对于井筒中心o的夹角。
2.根据权利要求1所述的一种利用阵列流量测井获取水平井持水率的方法,其特征在于:油水分界面高度h即到井筒底部的高度。