br>[0039]步骤2,点击确定，程序自动处理。程序内部处理的流程如图1所示的流程图，首先 进行对GR采用FIR滤波，FIR滤波采用31位有限长带通滤波器，设定GR时间采用4us，滤波频 率为 l-8khz。
[0040]步骤3，对滤波的GR离散化曲线，进行二阶求导：
[0041 ] GR" =GR[ i+1 ]+GR[ i-1 ]-2*GR[ i ]
[0042] 式中，i表示当前伽马曲线数据点，i-1表示前一个数据点，i + 1为后一个数据点； GR"为伽马曲线的二阶导数。
[0043] 步骤4,当二阶导数小于0,就是局部最大值的位置，也就是局部GR峰值。
[0044] 步骤5,将局部GR峰值与阈值比较，小于阈值的，则不认为是局部最大值，过滤掉； 大于阈值的，则为所需要的伽马曲线GR局部最大值；
[0045] 其中，程序启动时，读取伽马曲线GR所有值，从中统计出伽马曲线GR最大值的十分 之一作为默认阈值。
[0046] 步骤6,将所有局部最大值对应深度的GR值、伪深度、真深度、时间值写入表格中， 计算流体速度：
[0047] V = 60*(Depth[ i_l ]_Depth[ i] )/(Time[ i_l ]_Time[ i])
[0048] 其中，V表示〇6口1±[;[-1]至深度〇6口1:11[;[]的平均流速，〇6口1:11[;[-1]为第;[-1点所对 应的真深度，Depth [ i ]为第;[点所对用的真深度，Time [ i-1 ]为第i-1点所对应的时间，Time [i]为第i点所对应的时间。
[0049] 步骤7,将计算得到的流体速度点值绘制成曲线，通过下述方式实现：
[0050] Depthl，Depth2的流体速度分别为Vdepthi，Vdepth2，那么深度为Vdepth的流体速度为：
[0051 ] Vdepth = Vdepth2+(Depth-Depth2)*(Vdepthi-VdePth2)/(Depthl-Depth2) 〇
[0052]生成的表格如图4,峰值自动拾取的效果图如图3,波峰对应横线部分是图4表格的 可视化显示，对应的是GR峰值的位置，右侧曲线是追踪深度曲线，对应的是真实的深度值， 计算的是上提方向的峰值，对应的追踪深度曲线是往上深度增加的部分。可以看出识别峰 值位置很准确，没有遗漏。
[0053]根据保存的表格，计算速度曲线及合成曲线；图5为相应的设置对话框，选择生成 的TVAUOl表格，点击"转换曲线"，则生成相应的曲线，方法是对生成的速度表格，进行线性 差值，生成速度曲线。生成的速度效果图如图6所示。
[0054]以上对本发明实施例所提供的实施方案进行了详细介绍，本文中应用了具体个例 对本发明实施例的原理以及实施方式进行了阐述，以上实施例的说明只适用于帮助理解本 发明实施例的原理；同时，对于本领域的一般技术人员，依据本发明实施例，在具体实施方 式以及应用范围上均会有改变之处，综上所述，本说明书内容不应理解为对本发明的限制。
【主权项】
1. 一种生产测井相关流量数据预处理的方法，其特征在于，包括以下步骤： 步骤1，对采用单探头追踪法测量的生产测井相关流量数据，设置深度曲线DEPTH、伽马 曲线GR和时间曲线??ΜΕ; 步骤2，对伽马曲线GR进行FIR滤波，去除毛刺和干扰； 步骤3,对滤波后的伽马曲线GR进行二阶求导： GR"=GR[i+l]+GR[i-l]-2*GR[i] 式中，i表示当前伽马曲线数据点，i-1表示前一个数据点，i + 1为后一个数据点;GR"为 伽马曲线的二阶导数； 步骤4,当二阶导数小于0,即为伽马曲线GR局部最大值的位置，即，伽马曲线GR局部峰 值； 步骤5,将伽马曲线GR局部峰值与系统预设阈值比较，小于阈值的，则不认为是伽马曲 线GR局部最大值，过滤掉;大于阈值的，则为所需要的伽马曲线GR局部最大值； 步骤6,将步骤5获得的所需要的伽马曲线GR局部对应深度的伪深度、真深度和时间值 一并写入表格中，计算流体速度； 步骤7，将计算得到的流体速度点值绘制成曲线。2. 根据权利要求1所述的方法，其特征在于，步骤2中，对伽马曲线GR采用31位有限长带 通滤波器进行FIR滤波，设定GR离散化曲线时间采用4us，滤波频率为l-8khz。3. 根据权利要求1所述的方法，其特征在于，步骤5中，程序启动时，读取伽马曲线GR所 有值，从中统计出伽马曲线GR最大值的十分之一作为默认阈值。4. 根据权利要求1所述的方法，其特征在于，步骤6中，将所有伽马曲线GR局部对应深度 的伪深度、真深度和时间值写入表格中计算流体速度，通过下式得到： V = 60*(Depth[i_l]-Depth[i] )/(Time[i_l]-Time[i]) 其中，V表示〇6口1:11[;[-1]至深度〇6口1:11[;[]的平均流速，〇6口1:11[;[-1]为第;[-1点所对应的 真深度，Depth [ i ]为第;[点所对用的真深度，Time [ i-Ι ]为第i-Ι点所对应的时间，Time [ i ]为 第i点所对应的时间。5. 根据权利要求1或5所述的方法，其特征在于，步骤7中，将计算得到的流体速度点值 绘制成曲线，通过下述方式实现： Depth 1，DePth2的流体速度分别为Vdepthl，VdePth2，那么深度为Vdepth的流体速度为： Vdepth一 Vdepth2+ (Depth_Depth2) * (Vdepthi_Vdepth2) / (Depthl _Depth2) 〇
【专利摘要】本发明提供了一种生产测井相关流量数据预处理的方法，包括对相关流量数据设置深度曲线DEPTH、伽马曲线GR和时间曲线TIME；FIR滤波，对滤波后的伽马曲线GR进行二阶求导；当二阶导数小于0，即为伽马曲线GR局部峰值；局部峰值与系统预设阈值比较，小于阈值的过滤掉，大于阈值的保留；将伽马曲线GR局部对应深度的伪深度、真深度和时间值一并写入表格中，计算流体速度；将计算得到的流体速度点值绘制成曲线。其核心在如何通过算法准确识别GR曲线局部最大值，此处理流程操作简便，自动识别准确性高，处理效果较好，具备一定的推广应用价值。
【IPC分类】E21B47/00
【公开号】CN105626039
【申请号】CN201511032493
【发明人】倪路桥, 韩炜, 杜钦波, 宁卫东, 陈小磊, 王青艳, 左俊林
【申请人】中国石油天然气集团公司, 中国石油天然气集团公司管材研究所
【公开日】2016年6月1日
【申请日】2015年12月31日