C2测井曲线的自动分层问题 - 图文

图3 层内差异法粗分层结果

图4 聚类分析并层后结果

5.2.3 聚类分析

我们通过上一步，确定了几条分层性质良好的曲线。我们选取GR、AC、POR、SW、PORW、RILD、YW进行聚类分析。经我们多次调试，聚类分析中马氏距离判定常数Dmin取0.18。得到1号井的分层结果，如图4所示。图中横坐标代表分层点的个数，纵坐标代表分层界点的深度。

5.2.4 对模型一分析讨论

从模型一的自动分层结果来看，人工分层参考点都在结果中有所体现，但是分层点过于密集。如果调整参数让分层点变少，则有些必要的层界点会消失。也就是说，模型一很难做到两全其美。

5.3 模型二求解

根据前文所述的模型，我们将模型二的分层过程总结为如下图图5所示流程图。

12

开始导入测井数据前期处理测井曲线归一化平滑滤波算出活度函数，做出活度曲线对曲线上每一点求差分活度函数峰值检测di?xi?1?xidi?0且di+1?0是这个点是峰值，即层界面检测这些峰值是否大于设定的门限值T否不是层界面xi?T假界面剔除否并层是检测相邻两层之间的层均值是否大于设定的均值检验标准Sxk?1?xk?S是分层有效输出所有的层界面点结束否并层图5 模型二流程图

根据流程图，我们先对1号井GR曲线归一化处理、平滑滤波。处理之后的曲线见下图图6。图中蓝线表示滤波前的图像，红线表示平滑滤波之后的结果（图6为放大的

13

局部图像）。

图6 平滑滤波结果（局部）

我们选用1号井由5.1的主成分分析得到的6条测井曲线及它们各自的权重。门限值取0.0005，层间差判定值D取0.1，厚度门限T取10。得到的分层曲线如下图图7所示。图中横坐标代表分层点的个数，纵坐标代表深度。蓝色的横线是人工分层结果，红点是通过上述活度函数法得到的分层点。

可见，与人工层界点相比较，误差已经非常小，但是还存在这样的问题：自动分层的界点会多于人工分层界点，此问题可以通过修改使用的测井曲线、调节他们的权重、适当修改参数等方法完善。

图7 模型二自动分层结果

14

5.4 模型一、二求解结果对比

用CarbonsoftIntellExplore测井绘图软件将1号井用两种模型得到的自动分层结果连同人工分层结果一起绘于同一张图上以便对比；同时将两种模型分别用到的DEN、SP、GR、AC、RML、RMN等测井曲线的原始数据也绘制上。绘制结果如下图图8所示，图中第一列为深度，第二列为使用模型一的自动分层结果，第三列为使用模型二的自动分层结果，第四列为人工分层结果，第五到十列为两个模型中用到的DEN、GR、AC、SP、CNL、RT等测井曲线的原始数据。

图8 1号井两种模型自动分层结果

由图可见，模型一、模型二均能与人工分层结果大致吻合。然而模型一自动分层的层界面过于密集，这将给后续进一步分层带来无可避免的困难；而模型二得到的分层界面数目适中，并且绝大多数点与人工分层界面吻合得很好。

另一方面，因为模型一采用的是数理统计方法，因此计算量大，运行效率低。而模型二的效率大大提高。

总之，在对模型一、模型二进行比较后，我们得出结论：模型二的结果明显优于模型一。所以对2~13号测井的自动分层我们采用模型二。

5.5 2~7号井自动分层结果

我们用模型二（即活度函数法）将2~7号井进行自动分层。将分层结果整理成表格形式以方便阅读和对比。

2号井 3号井 4号井 井层 人工分层 自动分层 人工分层 自动分层 人工分层 自动分层 长31 257 252 213 218 276 280.3 长32 290 291.3 251 247 311 315.9 15