一种生产测井相关流量数据预处理的方法
【技术领域】
[0001 ]本发明属于石油开发和工程领域，是一种在生产测井注入剖面解释中，相关流量 数据预处理的方法。
【背景技术】
[0002] 随着生产测井技术的发展，相关流量技术越来越多的应用在注水剖面的解释中。 相关流量适用于中、低流量井，一般在注水井使用，在生产井使用由于流型变化会使分辨率 下降。在不能使用涡轮流量计测量的井中，一般采用相关流量计。相关流量计在井筒中确定 流量的方法是速度法。测量根据井况安装仪器，仪器停在两个射孔层之间向井筒中喷射示 踪剂，然后测量示踪剂在两个点传递所需要的时间，一般指两个探测器间或喷射器至探测 器间的时间，由此确定每一个解释层的视流速。对于生产井，喷射点应靠近射孔层夹层的底 部，对于注入井，则选择夹层的顶部。根据测量方式，速度法包括两种方法，一种是静止测量 法;另一种是追踪法。
[0003] 对于伽马射线探测器，由于喷射示踪剂的时间是变化的，因此精确确定喷射示踪 剂到达探头的时间较为困难，因此一般采用连续追踪方法。在这一情况下，在一个夹层内， 要进行至少三次测量。由于测量的同时示踪桥塞也在流动，因此必须保证有足够高的流速 测量完整的示踪段塞。如果发现第一次的位移较大，则应加快测速，反之则降低测速。流速 的计算方法为
[0005] 式中，△ H两次测量示踪剂段塞位移的距离(峰值的深度差）；
[0006] Δ t为段塞位移所需的时间；
[0007] 相关流量预处理所做的就是找到GR曲线的局部最大值，也就是波峰，然后根据波 峰的伪深度，找到Depth曲线对应的真深度及Time曲线对应的时间，然后就可以计算两点的 流速。
【发明内容】
[0008] 本发明的目的是针对目前没有较好的通过程序自动识别GR曲线局部最大值、手工 交互拾取又存在费时费力的问题，寻找一种能够通过程序和设定相应的阈值自动识别局部 最大值进而计算速度曲线的相关流量预处理方法。本发明通过程序自动识别GR曲线局部最 大值，根据局部最大值对应的伪深度，获取相应的真深度、时间，从而计算流速，本发明操作 简便，自动识别准确性高，具有一定的工程应用价值。
[0009] 本发明的目的是通过下述技术方案来实现的。
[0010] -种生产测井相关流量数据预处理的方法，包括以下步骤：
[0011] 步骤1，对采用单探头追踪法测量的生产测井相关流量数据，设置深度曲线DEPTH、 伽马曲线GR和时间曲线??ΜΕ;
[0012]步骤2，对伽马曲线GR进行FIR滤波，去除毛刺和干扰；
[0013]步骤3，对滤波后的伽马曲线GR进行二阶求导：
[0014] GR" =GR[ i+1 ]+GR[ i-1 ]-2*GR[ i ]
[0015]式中，i表示当前伽马曲线数据点，i-1表示前一个数据点，i + 1为后一个数据点； GR"为伽马曲线的二阶导数；
[0016] 步骤4,当二阶导数小于0,即为伽马曲线GR局部最大值的位置，即，伽马曲线GR局 部峰值；
[0017] 步骤5,将伽马曲线GR局部峰值与系统预设阈值比较，小于阈值的，则不认为是伽 马曲线GR局部最大值，过滤掉;大于阈值的，则为所需要的伽马曲线GR局部最大值；
[0018] 步骤6,将步骤5获得的所需要的伽马曲线GR局部对应深度的伪深度、真深度和时 间值一并写入表格中，计算流体速度；
[0019] 步骤7，将计算得到的流体速度点值绘制成曲线。
[0020] 进一步，步骤2中，对伽马曲线GR采用31位有限长带通滤波器进行FIR滤波，设定GR 离散化曲线时间采用4us，滤波频率为l-8khz。
[0021 ]进一步，步骤5中，程序启动时，读取伽马曲线GR所有值，从中统计出伽马曲线GR最 大值的十分之一作为默认阈值。
[0022]进一步，步骤6中，将所有伽马曲线GR局部对应深度的伪深度、真深度和时间值写 入表格中计算流体速度，通过下式得到：
[0023] V = 60*(Depth[ i_l ]_Depth[ i] )/(Time[ i_l ]_Time[ i])
[0024] 其中，V表示〇6口1±[;[-1]至深度〇6。1:11[;[]的平均流速，〇6。1:11[;[-1]为第;[-1点所对 应的真深度，Depth [ i ]为第;[点所对用的真深度，Time [ i-1 ]为第i-1点所对应的时间，Time [i]为第i点所对应的时间。
[0025]进一步，步骤7中，将计算得到的流体速度点值绘制成曲线，通过下述方式实现： [0026] Depthl，Depth2的流体速度分别为VdepthI，Vdepth2，那么深度为V depth的流体速度为：
[0027] Vdepth = Vdepth2+(Depth-Depth2)*(Vdepthi-VdePth2)/(Depthl-Depth2) 〇
[0028] 相对于现有技术，本发明的有益效果在于：
[0029] 本发明提供了一种生产测井相关流量数据预处理的方法。其核心在如何通过算法 准确识别GR曲线局部最大值，本发明采用先对GR曲线进行FIR滤波，然后对曲线二次求导， 通过二次导数小于〇即为局部最大值，将局部最大值和设定阈值相比，小于的滤除掉，然后 根据局部最大值位置获取相应的真深度、时间，从而计算速度曲线。此处理流程具备创新 性。现场操作简便，处理效果较好，具备一定的推广应用价值。
【附图说明】
[0030] 图1为处理流程示意图。
[0031] 图2为参数设置对话框。
[0032] 图3为识别效果图，红色表示识别的局部最大值。
[0033]图4为以表格形式存储的局部最大值对应的曲线数据。
[0034]图5为计算速度曲线对话框。
[0035]图6为计算的速度曲线效果图。
【具体实施方式】
[0036]下面结合具体相关流量测井数据，对本发明的具体的实施方式进行说明。
[0037]如图1所示，本发明一种生产测井相关流量数据预处理的方法，包括下述步骤：
[0038]步骤1，设置相关参数，如图2所示，对采用单探头追踪法测量的生产测井相关流量 数据，设置输入的深度曲线DEPTH、伽马曲线GR、时间曲线TIME，设置幅度系数1为取的波峰， 程序自动计算一个输入阈值大小为603。设置方向为上提测井，则程序只统计计算上提方向 的最大值。深度取向为深，指的是幅度取的不是波峰或波谷时，而是波峰乘以幅度系数时， 这时的取值取的是上半部分还是下半部分。<