CN106157158A - 建立盐膏岩矿物成分与岩石速度关系模版的方法及装置 - Google Patents

Abstract

本发明涉及一种建立盐膏岩矿物成分与岩石速度关系模版的方法及装置，其中，方法包括：获得单井的盐膏岩矿物成分各自的骨架速度；盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；获得单井的各层段中的盐膏岩矿物成分含量；在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量，纵坐标表示岩石速度。

Description

建立盐膏岩矿物成分与岩石速度关系模版的方法及装置

技术领域

本发明涉及地球物理勘探技术领域，特别涉及一种建立盐膏岩矿物成分与岩石速度关系模版的方法及装置。

背景技术

国内的库车坳陷、江汉盆地等含油气盆地都是含盐盆地。针对含盐盆地地震勘探中，遇到一个重要问题是如何针对含盐盆地构造获得准确成像，这一问题的核心是获得精确的盐膏岩速度问题。由于盐膏岩异常速度层的存在，对地震波能量有强烈的吸收、衰减作用，导致目的层反射能量弱，地震时间剖面出现“假构造”的陷阱，因此，对盐膏岩的速度研究具有重要意义。目前，解决盐膏岩速度问题的主要方法是：解释盐膏岩地震层位，由盐膏岩层位控制，对盐岩和膏岩分别充填一定的速度后，在进行构造成图，消除构造假象，取得一定的效果。但由于盐膏岩地层中一般发育膏岩、盐岩和泥膏岩，包含了石膏、石盐、泥岩和水共4种矿物，不同矿物组分含量岩石的速度差别很大，这需要搞清它们对速度的影响关系。

在岩石速度与矿物组分含量的关系中，一般是针对砂泥岩地层，在实验室测定泥质含量和孔隙度对岩石速度的影响，很少研究盐膏岩的矿物成分含量对盐膏岩速度的影响。

如何确定盐膏岩的矿物成分含量对盐膏岩速度的影响是目前急需解决的问题。

发明内容

为解决现有技术的问题，本发明提出一种建立盐膏岩矿物成分与岩石速度关系模版的方法及装置，通过该关系模块能直观的评估盐膏岩组分变化对速度的影响。

为实现上述目的，本发明提供了一种建立盐膏岩矿物成分与岩石速度关系模版的方法，包括：

获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

优选地，还包括：

把目标区域的盐膏岩矿物组分含量数据投射在所述盐膏岩矿物成分与岩石速度关系模版上，估算目标区域的盐膏岩矿物组分含量对岩石速度的影响程度。

优选地，所述孔隙度含量确定为0、5％、10％、15％、20％、25％。

优选地，所述获得单井的盐膏岩矿物成分含量的步骤包括：

利用多矿物测井成分解释原理，建立多组分测井响应方程；

假设盐膏岩孔隙内部为地层水，且完全充填，根据GR测井曲线和DT测井曲线，对所述多组分测井响应方程进行转化；

根据双水线性模型，对转化后的多组分测井响应方程进行优化；

对优化后的方程组进行求解，获得所述单井的各层段的盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量。

为实现上述目的，本发明还提供了一种建立盐膏岩矿物成分与岩石速度关系模版的装置，包括：

骨架速度获取单元，用于获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

盐膏岩矿物成分含量获取单元，用于获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

岩石速度获取单元，用于在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

盐膏岩矿物成分与岩石速度关系模版获取单元，用于将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

优选地，还包括：

反演单元，用于把目标区域的盐膏岩矿物组分含量数据投射在所述盐膏岩矿物成分与岩石速度关系模版上，估算目标区域的盐膏岩矿物组分含量对岩石速度的影响程度。

优选地，所述岩石速度获取单元将孔隙度的含量确定为0、5％、10％、15％、20％、25％。

优选地，所述盐膏岩矿物成分含量获取单元包括：

多组分测井响应方程建立模块，用于利用多矿物测井成分解释原理，建立多组分测井响应方程；

转化模块，用于假设盐膏岩孔隙内部为地层水，且完全充填，根据GR测井曲线和DT测井曲线，对所述多组分测井响应方程进行转化；

优化模块，用于根据双水线性模型，对转化后的多组分测井响应方程进行优化；

求解模块，用于对优化后的方程组进行求解，获得所述单井的各层段的盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量。

上述技术方案具有如下有益效果：

由于盐膏岩矿物组分与岩石速度无法用一个统一的关系式表达，本技术方案提供不同孔隙度尺度下的矿物组分含量与岩石速度的关系模版，把实际盐膏岩矿物组分含量数据投射在模版上，就可以直观的估算其矿物组分含量对岩石速度的影响程度。

附图说明

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

图1为本发明提出一种建立盐膏岩矿物成分与岩石速度关系模版的方法流程图；

图2为本实施例的盐膏岩矿物成分含量示意图；

图3为石膏、盐、孔隙三元矿物组分含量与岩石速度关系模版示意图；

图4为石膏、盐、泥质三元矿物组分含量与岩石速度关系模版示意图；

图5为盐膏岩矿物成分与岩石速度关系模版示意图之一；

图6为盐膏岩矿物成分与岩石速度关系模版示意图之二；

图7为本发明提出一种建立盐膏岩矿物成分与岩石速度关系模版的装置框图之一；

图8为本发明提出一种建立盐膏岩矿物成分与岩石速度关系模版的装置框图之二；

图9为本发明提出的一种设备框图；

图10为实际运用中编号为Ber-21井石膏层的盐膏岩矿物成分与岩石速度关系模版示意图之一；

图11为实际运用中编号为Ber-21井石膏层的盐膏岩矿物成分与岩石速度关系模版示意图之二。

具体实施方式

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

本技术方案的工作原理为：现有的文献资料涉及的岩石速度与矿物组分含量的关系中，一般是针对砂泥岩地层，在实验室测定泥质含量和孔隙含量对岩石速度的影响，很少研究盐膏岩的矿物成分含量对盐膏岩速度的影响。

本发明适用于不同孔隙度尺度下，石膏、盐、泥质及孔隙四元矿物组分含量对盐膏岩速度的影响程度，形成了多尺度下的四元矿物组分含量与岩石速度关系的模版。并用实际数据校正对应盐膏岩的四种矿物组分的岩石物理参数，得到的适用于实际工区的多尺度下的盐膏岩矿物组分与岩石速度的模版，能直观的评估目标区域内盐膏岩矿物组分含量变化对速度的影响程度。

基于上述工作原理，本发明提出一种建立盐膏岩矿物成分与岩石速度关系模版的方法，如图1所示。包括：

步骤101)：获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

针对石膏矿物骨架速度来说，把盐膏岩的石膏层段中石膏矿物组分含量是99.99％的速度统计出来(由录井取芯或化验分析得出)，一般为48.38ft/us(换算为6299.817m/s)，同理论值相近，可以视其为石膏矿物的骨架速度。

针对盐矿物骨架速度来说，把盐膏岩的盐岩层段中盐矿物组分含量是99.99％的速度统计出来(由录井取芯或化验分析得出)，一般为67ft/us(换算为4549.032m/s)，同理论值相近，可以视其为盐矿物的骨架速度。

针对孔隙矿物骨架速度来说，一般取其理论值为189ft/us(换算为1612m/s)可以视其为孔隙矿物骨架速度。

针对泥质矿物骨架速度来说，由于泥岩是混合物，其骨架速度在不同地区是不同的。通过上下围岩分析、盐膏岩下伏碳酸盐岩目的层骨架分析来获取，本研究区泥岩的骨架速度为85.7ft/us(换算为3556m/s)。

步骤102)：获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

如图2所示，为本实施例的盐膏岩矿物成分含量示意图。在本步骤中，所述获得单井的盐膏岩矿物成分含量的步骤包括：

首先，利用多矿物测井成分解释原理，建立多组分测井响应方程；其中，将岩石组分归类为石膏、盐岩、泥质和孔隙四大部分，利用多矿物测井成分解释原理，建立多组分测井响应方程如下：

$\left\{\begin{array}{c}{V}_1{X}_1+V_2{X}_2+V_3{X}_3+{\mathrm{\ΦX}}_f=X\\ V_1+V_2+V_3+\mathrm{\Φ}=1\end{array}\right.---\left(1\right)$

式中，X为某一测井响应值；V₁、V₂、V₃、Φ分别为石膏、盐岩、泥质和孔隙体积相对含量；X₁、X₂、X₃、X_f分别为石膏、盐岩、泥质和孔隙所含流体的某一测井响应特征值。

然后，假设盐膏岩孔隙内部为地层水，且完全充填，根据GR测井曲线和DT测井曲线，对所述多组分测井响应方程进行转化；其中，能用于矿物成分解释的曲线只有GR测井曲线和DT测井曲线，由于盐膏岩孔隙内没有油气聚集，因此认为盐膏岩孔隙内部为地层水完全充填，包括可动水和束缚水。因此上述方程组可以转化为：

$\left\{\begin{array}{c}{V}_1\·{\mathrm{GR}}_1+V_2\·{\mathrm{GR}}_2+V_3\·{\mathrm{GR}}_3+\mathrm{\Φ}\·{\mathrm{GR}}_w=GR\\ V_1\·{\mathrm{DT}}_1+V_2\·{\mathrm{DT}}_2+V_3\·{\mathrm{DT}}_3+\mathrm{\Φ}\·{\mathrm{DT}}_w=DT\\ V_1+V_2+V_3+\mathrm{\Φ}=1\end{array}\right.---\left(2\right)$

式中，GR_w和DT_w分别为地层水自然伽玛值和地层水声波时差值。GR₁、GR₂、GR₃分别表示石膏、盐、泥质的自然伽玛值；GR表示盐膏岩总的自然伽玛值；DT₁、DT₂、DT₃分别表示石膏、盐、泥质的声波时差值；DT表示盐膏岩总的声波时差值。

再然后，根据双水线性模型，对转化后的多组分测井响应方程进行优化；

在本步骤中，由于研究区盐膏岩有4种主要组分，至少需要4个方程进行联立求解，而以上方程组只有3个方程，不能满足求解需求，因此基于优化矿物组分测井解释方法，引入了双水线性模型概念，结合GR和DT曲线，进行优化矿物成分解释。

双水线性模型方程如下：

$\left\{\begin{array}{c}(V_{U_{bw}}+V_{U_{fw}})(V_{U_{bw}}{C}_{U_{bw}}+V_{U_{fw}}{C}_{U_{fw}})=CT\\ (V_{bw}+V_{fw})(V_{bw}{C}_{bw}+V_{fw}{C}_{fw})={\mathrm{CX}}_O\end{array}\right.---\left(3\right)$

式中，CT为原状地层电导率；Vu_bw+Vu_fw＝V_bw+V_fw；Cu_bw、Cu_fw分别为原状地层束缚水、自由水电导率；Vu_bw、Vu_fw分别为原状地层束缚水和自由水相对体积。

对优化后的方程组进行求解，获得所述单井的各层段的盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量。

优化后的方程组为：

$\left\{\begin{array}{c}(V_{U_{bw}}+V_{U_{fw}})(V_{U_{bw}}{C}_{U_{bw}}+V_{U_{fw}}{C}_{U_{fw}})=CT\\ (V_{bw}+V_{fw})(V_{bw}{C}_{bw}+V_{fw}{C}_{fw})={\mathrm{CX}}_O\\ V_1\·{\mathrm{GR}}_1+V_2\·{\mathrm{GR}}_2+V_3\·{\mathrm{GR}}_3+\mathrm{\Φ}\·{\mathrm{GR}}_w=GR\\ V_1\·{\mathrm{DT}}_1+V_2\·{\mathrm{DT}}_2+V_3\·{\mathrm{DT}}_3+\mathrm{\Φ}\·{\mathrm{DT}}_w=DT\\ V_1+V_2+V_3+\mathrm{\Φ}=1\\ V_1{X}_1+V_2{X}_2+V_3{X}_3+{\mathrm{\ΦX}}_f=X\\ V_{U_{bw}}+V_{U_{fw}}=V_{bw}+V_{fw}\end{array}\right.---\left(4\right)$

对式(4)进行求解，获得V₁、V₂、V₃、Φ，从而获得单井的各层段的盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量。

步骤103)：在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

如图3所示，泥质含量为0，孔隙含量分别为0、5％、10％、15％、20％、25％时，盐膏岩矿物成分含量与岩石速度关系模版。其中，

当孔隙含量为0时，盐的含量和石膏的含量总为100％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；其中，石膏含量也为0时，岩石速度为盐的骨架速度；石膏含量为100％时，岩石速度为石膏的骨架速度。

当孔隙含量为5％时，盐的含量和石膏的含量总为95％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

当孔隙含量为10％时，盐的含量和石膏的含量总为90％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

当孔隙含量为15％时，盐的含量和石膏的含量总为85％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

当孔隙含量为20％时，盐的含量和石膏的含量总为80％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

当孔隙含量为25％时，盐的含量和石膏的含量总为75％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线。

为了获得四元关系模版图，固定孔隙的含量，建立石膏、盐、泥质三种矿物与岩石速度的关系图版。包括：

孔隙的含量为零时，建立石膏、盐、泥质三种矿物与岩石速度的关系图版；如图4所示。其中，

泥质含量为0时，盐的含量和石膏的含量总为100％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；其中，石膏含量也为0时，岩石速度为盐的骨架速度；石膏含量为100％时，岩石速度为石膏的骨架速度；

泥质含量为10％时，盐的含量和石膏的含量总为90％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

泥质含量为20％时，盐的含量和石膏的含量总为80％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

泥质含量为30％时，盐的含量和石膏的含量总为70％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

泥质含量为40％时，盐的含量和石膏的含量总为60％，石膏的含量发生变化，且盐的含量亦随之变化，获得石膏的含量每一次变化时相应的岩石速度；根据石膏的含量多次变化对应的岩石速度，拟合出相应的盐膏岩矿物成分含量与岩石速度关系曲线；

同理，孔隙含量分别为5％、10％、15％时，用上述相同的方式建立石膏、盐、泥质三种矿物与岩石速度的关系图版。

步骤104)：将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

在本实施例中，将孔隙含量分别为0％、5％、10％、15％时所对应的盐膏岩矿物成分含量与岩石速度关系曲线放在同一坐标系中，就构成了盐膏岩矿物成分与岩石速度关系模版。如图5、图6所示。在图5中，横坐标表示石膏的含量，纵坐标表示对应的岩石速度。在图6中，横坐标表示盐的含量，纵坐标表示对应的岩石速度。相同的孔隙含量，坐标点使用相同的标记。

从图5、图6可以看出，石膏含量变化范围是0-100％，盐含量变化范围是0-100％，泥质含量变化范围是0-20％，孔隙变化是0-15％。由于盐膏岩的孔隙度值较低，大多数是低于5％，因此，孔隙度0-5％时，图5、图6所示的模版可以适应于多数的盐膏岩地层，泥质含量的0-20％变化一般也适用多数的盐膏岩地层。

由于盐膏岩矿物组分与岩石速度无法用一个统一的关系式表达，本方法提供不同孔隙度尺度下的矿物组分含量与岩石速度的关系模版，实际盐膏岩地层的孔隙度一般是小于5％，对于阿姆河右岸盐膏岩的孔隙度多数为2％，计算出对应尺度下石膏、盐、泥质、孔隙四组分的含量与岩石速度关系的曲线，把实际盐膏岩矿物组分含量数据投射在模版上，如图10、图11所示。就可以直观的估算其矿物组分含量对岩石速度的影响程度。

应当注意，尽管在附图中以特定顺序描述了本发明方法的操作，但是，这并非要求或者暗示必须按照该特定顺序来执行这些操作，或是必须执行全部所示的操作才能实现期望的结果。附加地或备选地，可以省略某些步骤，将多个步骤合并为一个步骤执行，和/或将一个步骤分解为多个步骤执行。

如图7所示，为本发明提出一种建立盐膏岩矿物成分与岩石速度关系模版的装置框图之一。包括：

骨架速度获取单元701，用于获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

盐膏岩矿物成分含量获取单元702，用于获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

岩石速度获取单元703，用于在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

盐膏岩矿物成分与岩石速度关系模版获取单元704，用于将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

如图8所示，在图7的基础上，本装置还包括：

反演单元705，用于把目标区域的盐膏岩矿物组分含量数据投射在所述盐膏岩矿物成分与岩石速度关系模版上，估算目标区域的盐膏岩矿物组分含量对岩石速度的影响程度。

此外，尽管在上文详细描述中提及装置的若干单元，但是这种划分仅仅并非强制性的。实际上，根据本发明的实施方式，上文描述的两个或更多单元的特征和功能可以在一个单元中具体化。同样，上文描述的一个单元的特征和功能也可以进一步划分为由多个单元来具体化。

基于上述示例性装置和方法，本实施例还提出一种设备，如图9所示。包括：

存储器a，用于存储请求指令；

处理器b，其与所述存储器耦合，该处理器被配置为执行存储在所述存储器中的请求指令，其中，所述处理器被配置的应用程序用于：

获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

本发明实施例还提供一种计算机可读程序，其中当在电子设备中执行所述程序时，所述程序使得计算机在所述电子设备中执行如图1所述的方法。

本发明实施例还提供一种存储有计算机可读程序的存储介质，其中所述计算机可读程序使得计算机在电子设备中执行如图1所述的方法。

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

Claims (8)

1.一种建立盐膏岩矿物成分与岩石速度关系模版的方法，其特征在于，包括：

获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

2.如权利要求1所述的方法，其特征在于，还包括：

把目标区域的盐膏岩矿物组分含量数据投射在所述盐膏岩矿物成分与岩石速度关系模版上，估算目标区域的盐膏岩矿物组分含量对岩石速度的影响程度。

3.如权利要求1所述的方法，其特征在于，所述孔隙度含量确定为0、5％、10％、15％、20％、25％。

4.如权利要求1、2或3所述的方法，其特征在于，所述获得单井的盐膏岩矿物成分含量的步骤包括：

利用多矿物测井成分解释原理，建立多组分测井响应方程；

假设盐膏岩孔隙内部为地层水，且完全充填，根据GR测井曲线和DT测井曲线，对所述多组分测井响应方程进行转化；

根据双水线性模型，对转化后的多组分测井响应方程进行优化；

对优化后的方程组进行求解，获得所述单井的各层段的盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量。

5.一种建立盐膏岩矿物成分与岩石速度关系模版的装置，其特征在于，包括：

骨架速度获取单元，用于获得单井的盐膏岩矿物成分各自的骨架速度；其中，所述盐膏岩矿物成分包括：石膏、盐、泥质和孔隙；

盐膏岩矿物成分含量获取单元，用于获得单井的各层段中的盐膏岩矿物成分含量；其中，所述盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量之和为1；

岩石速度获取单元，用于在孔隙度含量确定的情况下，泥质含量发生变化，相应的石膏和盐的含量也发生变化，根据所述盐膏岩矿物成分各自的骨架速度以及孔隙度含量、泥质含量、石膏含量、盐含量获得相应的岩石速度；

盐膏岩矿物成分与岩石速度关系模版获取单元，用于将盐膏岩中不同的孔隙度含量所对应的岩石速度放在同一坐标系中，构成盐膏岩矿物成分与岩石速度关系模版；其中，所述盐膏岩矿物成分与岩石速度关系模版的横坐标表示盐膏岩矿物成分中石膏的含量或盐的含量；所述盐膏岩矿物成分与岩石速度关系模版的纵坐标表示岩石速度；在模版中，将相同孔隙度含量所对应的岩石速度拟合出盐膏岩矿物成分含量与岩石速度关系曲线。

6.如权利要求5所述的装置，其特征在于，还包括：

反演单元，用于把目标区域的盐膏岩矿物组分含量数据投射在所述盐膏岩矿物成分与岩石速度关系模版上，估算目标区域的盐膏岩矿物组分含量对岩石速度的影响程度。

7.如权利要求5所述的装置，其特征在于，所述岩石速度获取单元将孔隙度的含量确定为0、5％、10％、15％、20％、25％。

8.如权利要求5、6或7所述的装置，其特征在于，所述盐膏岩矿物成分含量获取单元包括：

多组分测井响应方程建立模块，用于利用多矿物测井成分解释原理，建立多组分测井响应方程；

转化模块，用于假设盐膏岩孔隙内部为地层水，且完全充填，根据GR测井曲线和DT测井曲线，对所述多组分测井响应方程进行转化；

优化模块，用于根据双水线性模型，对转化后的多组分测井响应方程进行优化；

求解模块，用于对优化后的方程组进行求解，获得所述单井的各层段的盐膏岩矿物成分中石膏含量、盐含量、泥质含量以及孔隙含量。