测井曲线解释方法及装置与流程

本发明涉及石油地球物理勘探技术领域，尤其涉及测井曲线解释方法及装置。

背景技术：

本部分旨在为权利要求书中陈述的本发明实施例提供背景或上下文。此处的描述不因为包括在本部分中就承认是现有技术。

在石油勘探开发过程中，测井曲线揭示了井孔周围地层的各种物理特征，是一种直接反映储层特性的地球物理测量方法。通常情况下，常规测井曲线包括声波测井曲线、密度测井曲线、自然伽玛测井曲线(即gr曲线)、自然电位测井曲线(即sr曲线)、各种电阻率曲线、井径曲线及孔隙度曲线等，这些曲线分别反映了储层的岩性、物性以及含流体性质等信息。随着石油勘探开发的不断深入，地震数据与测井曲线结合的更加紧密，一些原本基于地震数据的特色技术开始被应用到测井曲线解释领域，地震数据反映了地下介质的物性特征，而测井曲线中的声波和密度测井曲线与地震数据具有物理一致性，测井曲线与地震数据二者的区别在于研究尺度的差异。

岩石物理理论研究认为纵横波速比属性能够有效判别储层的含流体性质，这使得叠前地震反演受到越来越多的重视，通过叠前地震反演获得纵横波速度比与纵波阻抗，联合应用可以判断储层孔隙流体的分布。然而这项技术在测井曲线解释领域应用却不能得到充分利用，根本原因在于横波测井曲线的缺失或者不可靠。现阶段横波预测主要通过理论岩石的物理模型计算得到，而预测结果的可靠程度严重依赖物理模型的合理性以及基础参数的准确性，由于地下介质复杂多变、非均质性强等原因导致依靠理论模型预测的横波信息往往不能满足要求，导致预测储层孔隙流体的分布的准确度和精度较低。

在地球岩石物理领域，地震波速度具有频散的特征，即随着频率变化地震波速度也发生变化，而导致这种变化的原因则是储层中存在流体的流动。因此，地震波速度频散属性也被认为是一种有效的储层含流体性质的检测因素。当地震波在含流体的岩石中传播时，高频信息会迅速衰减，能量减弱，而低频信息能量则基本保持。同样的，对于声波测井得到的速度测井曲线也存在这种物理现象，因此基于声波测井利用速度测井曲线的速度频散属性能够帮助预测井孔周围储层含流体性质。但是，速度测井曲线在测量的过程中记录了储层的压实趋势信息，而这种压实趋势信息的幅值范围一般较广，使得对细微的速度频散特征起到了淹没的效果，同样导致预测储层孔隙流体的分布的准确度和精度较低。

因此，现有的测井曲线解释在预测储层孔隙流体的分布时，存在准确度和精度较低的问题。

技术实现要素：

本发明实施例提供一种测井曲线解释方法，用以提高预测储层孔隙流体分布的准确度和精度，该方法包括：

获取目标储层的纵波速度测井曲线和密度测井曲线；

分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线；

对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息；

利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线；

根据高频速度测井曲线确定目标储层的流体特征曲线；

根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。

本发明实施例还提供一种测井曲线解释装置，用以提高预测储层孔隙流体分布的准确度和精度，该装置包括：

获取模块，用于获取目标储层的纵波速度测井曲线和密度测井曲线；

值域重构模块，用于分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线；

滤波模块，用于对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息；

高频速度测井曲线获取模块，用于利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线；

流体特征曲线确定模块，用于根据高频速度测井曲线确定目标储层的流体特征曲线；

孔隙流体分布确定模块，用于根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。

本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述测井曲线解释方法。

本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述测井曲线解释方法的计算机程序。

本发明实施例中，获取目标储层的纵波速度测井曲线和密度测井曲线，分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线，对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，根据高频速度测井曲线确定目标储层的流体特征曲线，根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。本发明实施例，在横波测井曲线缺失的情况下，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，进一步根据高频速度测井曲线确定目标储层的流体特征曲线，密度测井曲线的直流分量表明了目标储层的压实趋势信息，流体特征曲线表明了目标储层的速度频散特性，去除目标储层的压实趋势信息后的流体特征曲线，可以提高预测目标储层孔隙流体分布的准确度和精度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。在附图中：

图1为本发明实施例提供的测井曲线解释方法的实现流程图；

图2为本发明实施例提供的测井曲线解释装置的功能模块图；

图3为本发明实施例提供的实际储层重构值域后的纵波速度测井曲线示意图；

图4为本发明实施例提供的实际储层重构值域后的密度测井曲线示意图；

图5为本发明实施例提供的实际储层重构值域后的密度测井曲线直流分量的示意图；

图6为本发明实施例提供的实际储层高频速度测井曲线直流分量的示意图；

图7为本发明实施例提供的实际储层流体特征曲线的示意图。

具体实施方式

为使本发明实施例的目的、技术方案和优点更加清楚明白，下面结合附图对本发明实施例做进一步详细说明。在此，本发明的示意性实施例及其说明用于解释本发明，但并不作为对本发明的限定。

虽然本发明提供了如下述实施例或附图所示的方法操作步骤或装置结构，但基于常规或者无需创造性的劳动在所述方法或装置中可以包括更多或者更少的操作步骤或模块单元。在逻辑性上不存在必要因果关系的步骤或结构中，这些步骤的执行顺序或装置的模块结构不限于本发明实施例或附图所示的执行顺序或模块结构。所述的方法或模块结构的在实际中的装置或终端产品应用时，可以按照实施例或者附图所示的方法或模块结构进行顺序执行或者并行执行。

针对现有技术中测井曲线解释存在的预测储层孔隙流体分布的准确度和精度较低的缺陷，本发明的申请人提出了一种测井曲线解释方法及装置，其通过在横波测井曲线缺失的情况下，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，进一步根据高频速度测井曲线确定目标储层的流体特征曲线。鉴于密度测井曲线的直流分量表明了目标储层的压实趋势信息，流体特征曲线表明了目标储层的速度频散特性，去除目标储层的压实趋势信息后的流体特征曲线，可以达到提高预测目标储层孔隙流体分布的准确度和精度的目的。

图1示出了本发明实施例提供的测井曲线解释方法的实现流程，为便于描述，仅示出了与本发明实施例相关的部分，详述如下：

如图1所示，测井曲线解释方法，包括：

步骤101，获取目标储层的纵波速度测井曲线和密度测井曲线；

步骤102，分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线；

步骤103，对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息；

步骤104，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线；

步骤105，根据高频速度测井曲线确定目标储层的流体特征曲线；

步骤106，根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。

目标储层是指需要确定储层中孔隙分布规律的储层。为了实现对目标储层孔隙流体的分布进行预测的目的，针对目标储层，首先获取目标储层的纵波速度测井曲线和密度测井曲线。

声波速度测井是指在钻孔中通过研究目标储层中声波传播速度来确定目标储层性质的一种测井方法。声波速度包括横波速度和纵波速度两个方面。纵波速度反映的是目标储层的纵波速度属性特征，对应的测井曲线称为目标储层的纵波速度测井曲线。

密度测井曲线，是确定目标储层岩性和岩石密度的重要测井方法，与声波测井、中子测井组合构成岩性孔隙度测井系列。密度反映的是目标储层的密度属性特征，对应的测井曲线称为目标储层的密度测井曲线。

在储层正常压实段内，目标储层的岩层孔隙度随地层埋深变化的关系反映了地层压实的规律，即地层的压实趋势。在发明实施例中，密度测井曲线包含有目标储层的压实趋势信息。将压实趋势信息从纵波速度测井曲线中去除，有助于更加准确的确定目标储层孔隙流体的分布。

鉴于纵波速度测井曲线和密度测井曲线表示的纵波速度和密度的度量存在差异性，因此需要分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，以使目标储层的纵波速度测井曲线和密度测井曲线表示的纵波速度和密度的度量统一，进而获得重构值域后的纵波速度测井曲线和密度测井曲线。

在本发明的一实施例中，重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线的值域范围在0至1之间。该实施例是将纵波速度测井曲线和密度测井曲线分别进行归一化，以使归一化后的纵波速度测井曲线和密度测井曲线表示的纵波速度和密度的度量相一致。

在获得重构值域后的密度测井曲线后，对重构值域后的密度测井曲线进行长波长滤波，进而获得目标储层的压实趋势信息，而该目标储层的压实趋势信息，为密度测井曲线的直流分量，表征的正是目标储层的压实趋势信息。因此，利用重构值域后的纵波速度测井曲线减去重构值域后的密度测井曲线的直流分量，将压实趋势信息从纵波速度测井曲线中去除，获得目标储层的高频速度测井曲线，进而根据高频速度测井曲线对目标储层孔隙流体的分布进行预测，最终确定目标储层孔隙流体的分布。

在根据高频速度测井曲线对目标储层孔隙流体的分布时，首先根据高频速度测井曲线确定目标储层的流体特征曲线，及将目标储层的高频速度属性特征转化为目标储层的流体属性特征，进而根据目标储层的流体特征曲线确定目标储层孔隙流体的分布，实现对目标储层孔隙流体分布规律的预测。

流体特征曲线反映了目标储层含水饱和度的变化趋势，在根据目标储层的流体特征曲线确定目标储层孔隙流体的分布时，流体特征曲线值越大，表明目标储层中流体特征曲线值对应的该处位置孔隙中含水饱和度越小(一般认为含水饱和度不大于20％)、含烃越多，表明了流体特征曲线值对应的该处位置很可能是油气层，即指示了油气层的存在。

流体特征曲线值越小，表明目标储层中流体特征曲线值对应的该处位置孔隙中含水饱和度越大(一般认为含水饱和度不小于80％)、含烃越少，表明了流体特征曲线值对应的该处位置很可能是水层，即指示了水层的存在。

在本发明实施例中，获取目标储层的纵波速度测井曲线和密度测井曲线，分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线，对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，根据高频速度测井曲线确定目标储层的流体特征曲线，根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。本发明实施例，在横波测井曲线缺失的情况下，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，进一步根据高频速度测井曲线确定目标储层的流体特征曲线，密度测井曲线的直流分量表明了目标储层的压实趋势信息，流体特征曲线表明了目标储层的速度频散特性，去除目标储层的压实趋势信息后的流体特征曲线，可以提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，测井曲线解释方法中，步骤102，分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线，包括：

步骤，确定纵波速度测井曲线的最大值和最小值；

步骤，根据纵波速度测井曲线及其最大值和最小值，对目标储层的纵波速度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线。

在对纵波速度测井曲线进行值域的重构时，首先确定纵波速度测井曲线的最大值vmax和最小值vmin，进而根据纵波速度测井曲线及其最大值vmax和最小值vmin，对目标储层的纵波速度测井曲线进行值域的重构，进而获得重构值域后的纵波速度测井曲线。

具体可以通过下述公式对纵波速度测井曲线进行值域的重构：

v2(j)＝[v1(j)-vmin]/[vmax-vmin]；

其中，v2(j)表示重构值域后的纵波速度测井曲线，v1(j)表示值域重构前的纵波速度测井曲线，j表示纵波速度测井曲线的采样点，vmax表示值域重构前的纵波速度测井曲线的最大值，vmin表示值域重构前的纵波速度测井曲线的最小值。

在本发明实施例中，确定纵波速度测井曲线的最大值和最小值，根据纵波速度测井曲线及其最大值和最小值，对目标储层的纵波速度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，测井曲线解释方法中，步骤102，分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线，包括：

步骤，确定密度测井曲线的最大值和最小值；

步骤，根据密度测井曲线及其最大值和最小值，对目标储层的密度测井曲线进行值域的重构，获得重构值域后的密度测井曲线。

在对密度测井曲线进行值域的重构时，首先确定密度测井曲线的最大值ρmax和最小值ρmin，进而根据密度测井曲线及其最大值和最小值，对目标储层的密度测井曲线进行值域的重构，进而获得重构值域后的密度测井曲线。

具体可以通过下述公式对密度测井曲线进行值域的重构：

ρ2(j)＝[ρ1(j)-ρmin]/[ρmax-ρmin]；

其中，ρ2(j)表示重构值域后的密度测井曲线，ρ1(j)表示值域重构前的密度测井曲线，j表示密度测井曲线的采样点，ρmax表示值域重构前的密度测井曲线的最大值，ρmin表示值域重构前的密度测井曲线的最小值。

在本发明实施例中，确定密度测井曲线的最大值和最小值，根据密度测井曲线及其最大值和最小值，对目标储层的密度测井曲线进行值域的重构，获得重构值域后的密度测井曲线，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，测井曲线解释方法中步骤103，对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，包括：

步骤，构建滤波器，确定滤波器的滤波算子；

步骤，利用滤波器的滤波算子对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息。

为了确定目标储层的压实趋势信息，需要利用滤波器对重构值域后的密度测井曲线进行长波长滤波。其中长波长的范围包括1310纳米至1550纳米波长。因此，构建滤波器，确定滤波器的滤波算子。其中，该滤波器为长波长滤波器。最终确定的滤波器的滤波算子为：

其中，t(j)为滤波算子，g(j)为采样点j到频谱中心的距离，g0为滤波器的通带半径，n表示滤波器的幂系数，其决定了滤波器的形状，频谱中心的位置为

在本发明的一实施例中，滤波器的通带半径g0＝10，滤波器的幂系数n＝2，频谱中心的位置为(0,150)，即m＝300。

在构建滤波器，并确定滤波器的滤波算子后，利用滤波器的滤波算子对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息。具体的，可以通过下述公式对重构值域后的密度测井曲线进行长波长滤波：

d(j)＝ρ2(j)×t(j)；

其中，d(j)表示重构值域后的密度测井曲线的直流分量，其表征了目标储层的压实趋势信息，ρ2(j)为重构值域后的密度测井曲线，t(j)为滤波器的滤波算子。

在本发明实施例中，构建滤波器，确定滤波器的滤波算子，利用滤波器的滤波算子对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

在确定目标储层的压实趋势信息后，即重构值域后的密度测井曲线的直流分量d(j)，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线。具体的，可以通过下述公式确定目标储层的高频速度测井曲线：

v3(j)＝v2(j)-d(j)；

其中，v3(j)为去掉压实趋势信息后的目标储层的高频速度测井曲线，v2(j)表示重构值域后的纵波速度测井曲线，d(j)表示重构值域后的密度测井曲线的直流分量。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，测井曲线解释方法中步骤105，根据高频速度测井曲线确定目标储层的流体特征曲线，包括：

步骤，对高频速度测井曲线进行拉普拉斯变换，获得不同频率对应的速度测井曲线；

步骤，根据不同频率对应的速度测井曲线确定目标储层的流体特征曲线。

拉普拉斯变换，又称为拉式变换(s变换)，是一个线性的积分变换。在本发明实施例中，对高频速度测井曲线运用拉普拉斯变换，得到不同频率对应的速度测井曲线。具体可以通过如下公式对高频速度测井曲线进行拉普拉斯变换：

其中，v3(j,f)表示不同频率对应的速度测井曲线，x表示目标储层中某位置的深度v3(x)表示深度x对应的速度曲线值，λa表示拉普拉斯变换的调节因子，p表示拉普拉斯变换的调节幂系数，f表示频率。

在本发明的一实施例中，调节因子λa＝1.2，调节幂系数p＝1。

在获得不同频率对应的速度测井曲线后，具体可以通过下述公式确定目标储层的流体特征曲线：

令则根据对高频速度测井曲线进行拉普拉斯变换的公式，对其进行一阶泰勒公式展开，可以得到：

其中，a表示目标储层的流体特征曲线，v3(j,f)表示不同频率对应的速度测井曲线，v3(j,f0)表示参考频率对应的速度测井曲线，f表示频率，f0表示参考频率。

在本发明的一实施例中，参考频率f0＝5hz。

进而可以根据获得的目标储层的流体特征曲线，确定目标储层孔隙流体的分布。

在本发明实施例中，对高频速度测井曲线进行拉普拉斯变换，获得不同频率对应的速度测井曲线，根据不同频率对应的速度测井曲线确定目标储层的流体特征曲线，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

本发明实施例中还提供了一种测井曲线解释装置，如下面的实施例所述。由于这些装置解决问题的原理与测井曲线解释方法相似，因此这些装置的实施可以参见方法的实施，重复之处不再赘述。

图2示出了本发明实施例提供的测井曲线解释装置的功能模块，为便于说明，仅示出了与本发明实施例相关的部分，详述如下：

参考图2，所述测井曲线解释装置所包含的各个模块用于执行图1对应实施例中的各个步骤，具体请参阅图1以及图1对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述测井曲线解释装置包括获取模块201、值域重构模块202、滤波模块203、高频速度测井曲线获取模块204、流体特征曲线确定模块205及孔隙流体分布确定模块206。

获取模块201，用于获取目标储层的纵波速度测井曲线和密度测井曲线。

值域重构模块202，用于分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线。

滤波模块203，用于对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息。

高频速度测井曲线获取模块204，用于利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线。

流体特征曲线确定模块205，用于根据高频速度测井曲线确定目标储层的流体特征曲线。

孔隙流体分布确定模块206，用于根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。

在本发明实施例中，获取模块201获取目标储层的纵波速度测井曲线和密度测井曲线，值域重构模块202分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线，滤波模块203对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，高频速度测井曲线获取模块204利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，流体特征曲线确定模块205根据高频速度测井曲线确定目标储层的流体特征曲线，孔隙流体分布确定模块206根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。本发明实施例，在横波测井曲线缺失的情况下，高频速度测井曲线获取模块204利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，进一步流体特征曲线确定模块205根据高频速度测井曲线确定目标储层的流体特征曲线，密度测井曲线的直流分量表明了目标储层的压实趋势信息，流体特征曲线表明了目标储层的速度频散特性，去除目标储层的压实趋势信息后的流体特征曲线，可以提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线的值域范围在0至1之间。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，所述测井曲线解释装置中值域重构模块202所包含的各个单元用于执行测井曲线解释方法中对应实施例中的各个步骤，具体请参阅测井曲线解释方法对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述测井曲线解释装置中值域重构模块202包括第一确定单元和第一获取单元。

第一确定单元，用于确定纵波速度测井曲线的最大值和最小值。

第一获取单元，用于根据纵波速度测井曲线及其最大值和最小值，对目标储层的纵波速度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线。

在本发明实施例中，第一确定单元确定纵波速度测井曲线的最大值和最小值，第一获取单元根据纵波速度测井曲线及其最大值和最小值，对目标储层的纵波速度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，所述测井曲线解释装置中值域重构模块202所包含的各个单元用于执行测井曲线解释方法中对应实施例中的各个步骤，具体请参阅测井曲线解释方法对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述测井曲线解释装置中值域重构模块202包括第二确定单元和第二获取单元。

第二确定单元，用于确定密度测井曲线的最大值和最小值。

第二获取单元，用于根据密度测井曲线及其最大值和最小值，对目标储层的密度测井曲线进行值域的重构，获得重构值域后的密度测井曲线。

在本发明实施例中，第二确定单元确定密度测井曲线的最大值和最小值，第二获取单元根据密度测井曲线及其最大值和最小值，对目标储层的密度测井曲线进行值域的重构，获得重构值域后的密度测井曲线，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，所述测井曲线解释装置中滤波模块203所包含的各个单元用于执行测井曲线解释方法对应实施例中的各个步骤，具体请参阅测井曲线解释方法对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述测井曲线解释装置中滤波模块203包括滤波算子确定单元和滤波单元。

滤波算子确定单元，用于构建滤波器，确定滤波器的滤波算子；

滤波单元，用于利用滤波器的滤波算子对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息。

在本发明实施例中，滤波算子确定单元构建滤波器，确定滤波器的滤波算子，滤波单元利用滤波器的滤波算子对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

在本发明的一实施例中，为了进一步提高预测目标储层孔隙流体分布的准确度和精度，所述测井曲线解释装置中流体特征曲线确定模块205所包含的各个单元用于执行测井曲线解释方法对应实施例中的各个步骤，具体请参阅测井曲线解释方法对应实施例中的相关描述，此处不再赘述。本发明实施例中，所述测井曲线解释装置中流体特征曲线确定模块205包括拉普拉斯变换单元和流体特征曲线确定单元。

拉普拉斯变换单元，用于对高频速度测井曲线进行拉普拉斯变换，获得不同频率对应的速度测井曲线；

流体特征曲线确定单元，用于根据不同频率对应的速度测井曲线确定目标储层的流体特征曲线。

在本发明实施例中，拉普拉斯变换单元对高频速度测井曲线进行拉普拉斯变换，获得不同频率对应的速度测井曲线，流体特征曲线确定单元根据不同频率对应的速度测井曲线确定目标储层的流体特征曲线，可以进一步提高预测目标储层孔隙流体分布的准确度和精度。

以下结合某地区的一实际储层，对本发明的实施和功能原理进行阐述：

首先获取该实际储层的纵波速度测井曲线和密度测井曲线，进而分别对实际储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线。在本发明实施例中，该实际储层的深度范围为2000米至2300米的深度范围。

其中，图3为该实际储层重构值域后的纵波速度曲线，图4为该实际储层重构值域后的密度测井曲线。从图3和图4中可以看出，重构值域后实际储层纵波速度曲线的值域范围和重构值域后密度测井曲线的值域范围均处于(0，1)之间。

进而对图4中实际储层重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，即重构值域后的密度测井曲线的直流分量。

图5为该实际储层的压实趋势信息，即重构值域后的密度测井曲线的直流分量的形态，从图5中可以看出来，经过长波长滤波处理后的密度测井曲线的高频细节被完全压制，只保留反映实际储层压实趋势信息的低频部分，即重构值域后的密度测井曲线的直流分量。由于地层沉积持续且缓慢，造成密度测井曲线的直流分量变化平缓且不受该实际储层中岩石微观性质差异的影响。

然后利用图3中的该实际储层重构值域后的纵波速度曲线，减去目标储层的压实趋势信息，即图5所示的重构值域后的密度测井曲线的直流分量，得到图6所示的目标储层的高频速度测井曲线。

从图6中可以看出，高频速度测井曲线中没有反应该实际储层压实趋势信息的直流分量，只包含微观孔隙以及孔隙流体等物性差异带来的高频速度振幅异常响应，通过图6中这些高频速度振幅异常响应，可以推断该实际储层的孔隙流体的分布规律以及孔隙流体的类型等。

在获得图6所示的高频速度测井曲线后，对高频速度测井曲线进行拉普拉斯变换，得到实际储层的流体特征曲线。

图7为获得的实际储层的流体特征曲线。图7中横坐标为流体特征曲线的幅值，纵坐标为深度。从图7中可以看出来，流体特征曲线在为2150米至2250米的深度范围内存在振幅异常，该深度范围内的流体特征曲线的幅值明显大于上层及下层，表明该深度段可能是油气层，而上层和/或下层则为水层或者致密层。

最后在根据目标储层的流体特征曲线确定目标储层孔隙流体的分布时，利用流体特征曲线对目标储层进行解释，解释结论认为在2150米至2250米的深度范围内存在油气层，并且由于2150米至2200米的深度范围内，流体特征曲线的幅值明显要比2200米至2250米的深度范围内的幅值大，根据岩石理论认为：储层含气相对于储层含油能够对地层产生更加强烈的吸收衰减效应，据此推测2150米至2200米的深度范围应该是气层，而2150米至2200米的深度范围应该是油层，并根据测井解释结果提交了两处射孔建议，分别为2223米(#1)和2188米(#2)，结果显示2223米(#1)和2188米(#2)分别为油层和气层，与预测结论完全一致。后续又提交了一处浅层的试油2050米(#3)，结果显示为水层。三次试油结果与基于本流体特征曲线预测的储层孔隙流体分布规律的情况完全吻合和一致，验证了本发明的有效性，即本发明实施例提出的测井解释方法和装置，可以准确预测储层孔隙流体的分布规律，提高储层预测及孔隙流体检测的准确性，降低油气勘探成本和风险。

本发明实施例还提供一种计算机设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现上述测井曲线解释方法。

本发明实施例还提供一种计算机可读存储介质，所述计算机可读存储介质存储有执行上述测井曲线解释方法的计算机程序。

综上所述，本发明实施例获取目标储层的纵波速度测井曲线和密度测井曲线，分别对目标储层的纵波速度测井曲线和密度测井曲线进行值域的重构，获得重构值域后的纵波速度测井曲线和重构值域后的密度测井曲线，对重构值域后的密度测井曲线进行长波长滤波，获得目标储层的压实趋势信息，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，根据高频速度测井曲线确定目标储层的流体特征曲线，根据目标储层的流体特征曲线确定目标储层孔隙流体的分布。本发明实施例，在横波测井曲线缺失的情况下，利用重构值域后的纵波速度测井曲线减去目标储层的压实趋势信息，获得目标储层的高频速度测井曲线，进一步根据高频速度测井曲线确定目标储层的流体特征曲线，密度测井曲线的直流分量表明了目标储层的压实趋势信息，流体特征曲线表明了目标储层的速度频散特性，去除目标储层的压实趋势信息后的流体特征曲线，可以提高预测目标储层孔隙流体分布的准确度和精度。

本领域内的技术人员应明白，本发明的实施例可提供为方法、系统、或计算机程序产品。因此，本发明可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本发明可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。

本发明是参照根据本发明实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。

这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。

这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。

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