[发明专利]一种提取示踪测井资料形成流动成像图的方法

权利要求书

1.一种提取示踪测井资料形成流动成像图的方法，其特征在于，包括下述步骤：

步骤101：全要素提取基于时间驱动的测井数据，将测井数据区分为上测和下测数据序列，分别解析上测和下测测量段内的每一深度点的示踪剂测量值，提取每个深度点不同测量时间测量的示踪幅度；

所述全要素提取基于时间驱动的测井数据包括以下步骤：

S01：加载数据驱动测井数据；

S02：区分上下测示踪追踪曲线；

S03：剔除本底值的影响；

所述步骤S03中，剔除本底值采用如下式(2)所示：

ΔGR＝IT-GR (2)

式中，GR为伽玛基线，API；IT为示踪伽玛，API；

S04：归一化示踪幅度；

所述步骤S04中归一化示踪幅度处理，通过公式(3)计算：

NGR＝ΔGR/max(ΔGR) (3)

式中，NGR为归一化的示踪幅度；ΔGR为示踪幅度，API；

S05：提取各深度点的测试序列；

所述提取每个深度点不同测量时间测量的示踪幅度，包括以下步骤：

S06：提取时间驱动数据中的深度D_i、时间T_i数据和伽玛GR_i，1≤i≤n；

S07：创建上、下测数据缓冲区，并设置当前处理索引i＝1；

S08：判断i是否小于等于n-1；如果不是，就退出循环，否则执行步骤S 09；

S09：判断(D_i+1-D_i)ⅹ(T_i+1-T_i)是否大于0；如果是，下测数据缓冲区增加一行，把D_i、T_i和GR_i数据增加到下测数据缓冲区中；如果不是，判断(D_i+1-D_i)ⅹ(T_i+1-T_i)是否小于0，如果是，将上测数据缓冲区增加一行，把D_i、T_i和GR_i数据增加到上测数据缓冲区中；

S10：处理索引i增加1，返回步骤S08；

所述提取同一深度示踪测量时间谱包括以下步骤：

S11：获得上测没有重复的深度序列D_i，1≤i≤n；

S12：初始化测井矩阵数据缓冲区并设当前处理索引i＝1；

S13：判断i是否小于等于n；如果不是，处理完毕并退出循环；如果是，在上测测井数据缓冲中，提取深度为D_i的不同时间点的示踪剂伽玛测量值并将当前深度Di和不同时间点以及对应的伽玛测量值保存到测井矩阵数据缓冲区中；

S14：当前处理索引增加1，返回步骤S13；

步骤102：反演井筒内示踪剂载体随时间变化的运移成像，对测量深度范围内的每一个深度点的示踪测试序列，采用线性内插的方式对测试时间范围内的示踪幅度进行内插，实现示踪成像重建，直观展示将示踪剂流动形态；

所述步骤102，内插公式见公式(4)，形成测量范围内不同时间的示踪剂运移数据阵列：

NGR＝NGR_i+(T-T_i)×(NGR_j-NGR_i)/(T_j-T_i) (4)

式中，NGR_i为i时刻的归一化示踪幅度；NGR_j为j时刻的归一化示踪幅度；T_i为i时刻，s；T_j为j时刻，s；T为T时刻，s；

步骤103：提取井筒流体的流动速度，通过选取相同示踪形态和示踪峰值计算井筒流体流速；

所述步骤103，采用公式(5)计算井筒流体流速：

式中，V_a为流体流速，m/s；

Q_f＝V_a×S×86400 (6)

式中，Q_f为注入井或产出井特定深度流体日注入量或产出量，m³/d；V_a为流体流速，m/s；S为流体流动横截面积，m²。