获取高分辨率测井曲线的方法及装置.pdf

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定，基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明实施例提供一种获取高分辨率测井曲线的技术，该技术通过对测井曲线的校正使得校正后的测井曲线更加符合目标地层的真实地质状况，从而为准确描述目标地层的地质奠定良好基础。

基于此，本实施例提供一种获取高分辨率测井曲线的方法，如图1所示，图1为本实施例的一种获取高分辨率测井曲线的方法的示意流程图，该流程包括如下步骤：

101.通过岩心实验采集目标岩心的目标地质参数。

本实施例得先采集目标岩心的目标地质参数，该步骤通过岩心实验来采集，其中目标岩心为从目标测井地层采样的地质标本，该标本可以代表该地层的地质特性，比如若要对探测较薄的地层的原始测井曲线进行校正，那就将在该较薄的地层中采集岩心。而目标地质参数是指欲采集的地质参数，该参数是与原始测井曲线地质属性相同的参数，比如当原始测井曲线是描述较薄地层电阻率的，那目标地质参数就是待测电阻率参数；当原始测井曲线是描述较薄地层放射性的(即由自然伽马值组成的曲线)，那目标地质参数就是待测自然伽马值。

而为了较准确的采集地质数据，本步骤在采集目标岩心的目标地质参数时，采样间距不大于0.125米，比如在目标岩心上每隔0.1米采样一次地质数据。

102.根据各个多级反褶积校正公式以及原始测井曲线获取对应各个校正公式的多级校正曲线。

本实施例通过多级反褶积校正公式来获取校正曲线，该多级反褶积公式可以为3+(n-1)×2级反褶积校正公式，其中n为自然数。这样该多级反褶积校正公式至少为3级反褶积，所以本实施例中该各个多级反褶积校正公式至少包括当n为1时的3级的校正公式：

XR(z_i)＝M[-X(z_i-1)+N·X(z_i)-X(z_i+1)]

其中XR(Z_i)表示对应目标岩心上目标采集点的地质参数校正值，X(Z_i)表示原始测井曲线上对应目标采集点的地质参数原始值，该参数可通过原始测井曲线获得，Z_i为目标采集点，Z_i+1和Z_i-1为目标采集点具有一个采样间距的相邻采样点，M＝1/(αΔz)²，α为地质特征参数，Δz为采样间距，N＝2+(αΔz)²。

该公式可由下述方法推导而来，当确定了采集点Z_i后，即可确定距该采集点Z_i一个采样间距Δz的相邻采样点Z_i+1和Z_i-1，由于该地质参数校正值XR(Z_i)可以表示为：

XR(zi)=X(zi)-1α2d2Xdz2|z=zi]]>

所以在离散化二阶导数，通过二阶中心差商可导出公式：

d2Xdz2|z=zi≈X(zi+Δz)-2X(zi)+X(zi-Δz)(Δz)2]]>

将导出公式代入地质参数校正值XR(Z_i)公式后，即可获得3级反褶积校正公式：

XR(z_i)＝M[-X(z_i-1)+N·X(z_i)-X(z_i+1)]

为了更准确的校正原始测井曲线，确定与目标地质参数最接近的校正曲线，本实施例可以通过更多的多级反褶积校正公式来获取校正曲线。比如当n为2时的5级的校正公式：

XR(zi)=14M[(λ-1)X(zi+2)-4λ·X(zi+1)+(4N+6λ+2)X(zi)-4λ·X(zi-1)+(λ-1)X(Zi-2)]]]>

当根据加权平均获得公式：

d2Xdz2|z=zi≈λX(zi+1)-2X(zi)+X(zi-1)(Δz)2+(1-λ)X(zi+2)-2X(zi)+X(zi-2)(2Δz)2]]>

后，将其代入地质参数校正值XR(Z_i)公式即可获得5级的校正公式。其中各个参量的含义与3级的校正公式相同，故不再赘述，而Z_i+2和Z_i-2为目标采集点具有两个采样间距的采样点。

同理，多级反褶积校正公式还包括当n为3时的七级的校正公式：

XR(zi)=136M[-4λ·X(zi+3)-9μ·X(zi+2)-36(1-λ-μ)X(zi+1)+(36N-64λ-64μ)X(zi)]]>

-36(1-λ-μ)X(zi-1)-9μ·X(Zi-2)-4λ·X(Zi-3)]]]>

其中各个参量的含义与3级的校正公式相同，故不再赘述，而Z_i+3和Z_i-3为目标采集点具有三个采样间距的采样点。

这样就可以获取每个反褶积校正公式所对应的离散地质参数校正值，将这些相邻的离散地质参数校正值各个连接即可获得每个反褶积校正公式所对应的多级校正曲线。

203.根据获取到的各个多级校正曲线判断与所述采集到的目标地质参数最接近的一条多级校正曲线。

由于通过岩心实验采集的地质参数是比较准确的，所以越靠近通过岩心实验采集的地质参数，则校正曲线的校正效果越好。本步骤根据步骤202获取到的各个多级校正曲线判断与所述采集到的目标地质参数最接近的一条多级校正曲线。比如当通过岩心实验采集的地质参数几乎都与3级校正公式对应的3级校正曲线吻合(即通过岩心实验采集的地质参数几乎都落在曲线上)，那么该3级校正曲线就是最接近目标地质参数的曲线。对于目标岩心对应的油井，可以通过3级校正曲线来描述该油井的地质状况。而对于该地层上的其他油井，则可以通过3级反褶积校正公式来对测井曲线进行校正。这样可以对不同地层的测井曲线进行有效校正，提高测井曲线的分辨率，从而根据获取到的高分辨率测井曲线准确分析地层的地质状况。

为了本领域技术人员能更清楚的理解本发明实施例所述的方法，下面引用一示例来说明，需要注意的是本示例中的数字及地质参数仅仅是示例性的，不能将其做为限制本发明保护范围的解释，本领域技术人员应清楚了解本发明实施例所述的方法还可以通过其他具体的示例来表达。

如图2所示，图2为岩心实验采集目标岩心自然伽马值数据与校正前后测井曲线的对比图，从图2可以清楚看出目标岩心为2313～2320米的岩心样本，在这段深度内的岩心自然伽马值通过岩心实验采集获得并反映在了图2中，其中采样间距Δz为0.125米。图2中的连续曲线为原始测井曲线，该测井曲线是实地测量的描述该地层放射性的测井曲线。

为了对该原始测井曲线进行校正，本示例计算了三条多级反褶积校正曲线：3级反褶积校正曲线、5级反褶积校正曲线和7级反褶积校正曲线，其中3级反褶积校正公式中取α＝2，Δz＝0.125，所以M＝16，N＝2.0625，将这些参数代入公式，即

XR(z_i)＝16[-X(z_i-1)+2.0625·X(z_i)-X(z_i+1)]

根据原始测井曲线将采集点的地质参数原始值代入后可以算出各采集点的地质参数校正值，将相邻的校正值连接后即为3级反褶积校正曲线。同理，当取λ＝0.6，μ＝0.4时，5级反褶积校正公式为：

XR(zi)=14×16[(0.6λ-1)X(zi+2)-4×0.6X(zi+1)+(4×2.0625+6×0.6+2)X(zi)-4×0.6X(zi-1)+(0.6-1)X(Zi-2)]]]>

从而也可以算出各采集点的地质参数校正值，将相邻的校正值连接后即为5级反褶积校正曲线。当取λ＝0.2，μ＝0.3时，7级反褶积校正公式为：

XR(zi)=136×16[-4×0.2GR(zi+3)-9×0.3GR(zi+2)-36×0.5GR(zi+1)+(36×2.0625-64×0.5)GR(zi)]]>

-36×0.5GR(zi-1)-9×0.3GR(Zi-2)-4×0.2GR(Zi-3)]]]>

从而也可以算出各采集点的地质参数校正值，将相邻的校正值连接后即为7级反褶积校正曲线。将采集到的目标岩心自然伽马值与各个多级校正曲线进行对比，即可判断出与采集到的目标岩心自然伽马值最接近的一条多级校正曲线为7级反褶积校正曲线。从图2中也可以清楚看出GRR7曲线(即7级反褶积校正曲线)与岩心自然伽马能够较好的吻合。

实现上述方法实施例的装置可以是一个独立设备，下面对实现上述功能的装置进行说明。如图3所示，图3为本发明实施例的获取高分辨率测井曲线装置的功能结构框图，从图3中可知，该装置可以包括：

采集单元301，用于通过岩心实验采集目标岩心的目标地质参数；

获取单元302，用于根据各个多级反褶积校正公式以及原始测井曲线获取对应各个校正公式的多级校正曲线；

决策单元303，用于根据获取单元302获取到的各个多级校正曲线判断与采集单元301采集到的目标地质参数最接近的一条多级校正曲线。

存储单元304，用于存储3+(n-1)×2级反褶积校正公式，其中n为自然数。其中所述3+(n-1)×2级反褶积校正公式至少包括n为1时的三级的校正公式：

XR(z_i)＝M[-X(z_i-1)+N·X(z_i)-X(z_i+1)]

其中XR(Z_i)表示对应目标岩心上目标采集点的地质参数校正值，X(Z_i)表示原始测井曲线上对应目标采集点的地质参数原始值，Z_i为目标采集点，Z_i+1和Z_i-1为目标采集点具有一个采样间距的相邻采样点，M＝1/(αΔz)²，α为地质特征参数，Δz为采样间距，N＝2+(αΔz)²。各个多级反褶积校正公式还包括当n为2时的五级的校正公式：

XR(zi)=14M[(λ-1)X(zi+2)-4λ·X(zi+1)+(4N+6λ+2)X(zi)-4λ·X(zi-1)+(λ-1)X(Zi-2)]]]>

其中XR(Z_i)表示对应目标岩心上目标采集点的校正值，X(Z_i)表示原始测井曲线上对应目标采集点的原始值，Z_i为目标采集点，Z_i+1和Z_i-1为目标采集点具有一个采样间距的相邻采样点，M＝1/(αΔz)²，α为地质特征参数，Δz为采样间距，N＝2+(αΔz)²，λ和μ均为权重系数。

各个多级反褶积校正公式还包括当n为3时的七级的校正公式：

XR(zi)=136M[-4λ·X(zi+3)-9μ·X(zi+2)-36(1-λ-μ)X(zi+1)+(36N-64λ-64μ)X(zi)]]>

-36(1-λ-μ)X(zi-1)-9μ·X(Zi-2)-4λ·X(Zi-3)]]]>

其中XR(Z_i)表示对应目标岩心上目标采集点的校正值，X(Z_i)表示原始测井曲线上对应目标采集点的原始值，Z_i为目标采集点，Z_i+1和Z_i-1为目标采集点具有一个采样间距的相邻采样点，M＝1/(αΔz)²，α为地质特征参数，Δz为采样间距，N＝2+(αΔz)²，λ和μ均为权重系数。

为了有效获取多级校正曲线，获取单元302包括：

计算模块321，用于根据存储单元304存储的3+(n-1)×2级反褶积校正公式计算地质参数校正值；

曲线拟合模块322，用于根据计算模块321计算出的地质参数校正值拟合曲线。

本实施例所述的装置可以是计算机等具有信息处理能力的设备，通过对不同地层的测井曲线进行有效校正，提高测井曲线的分辨率，从而根据获取到的高分辨率测井曲线准确分析地层的地质状况。

当然，以上所述的具体实施例，对本发明的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本发明的具体实施例而已，并不用于限定本发明的保护范围，凡在本发明的逻辑和原则之内，所做的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。