确定砂岩型铀矿勘查区的预测横波测井曲线的方法

技术领域

本申请的实施例涉及地震勘探领域，具体涉及一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法。

背景技术

砂岩型铀矿是指铀矿床产于砂岩、砂砾岩等中粗碎屑岩的铀矿。在对砂岩型铀矿进行勘测时，主要采用地震勘探技术。在地震勘探技术中，主要通过观测和分析大地对人工激发地震波的响应，推断出地下岩层性质和形态。

地震波可以分为横波和纵波，粒子振动方向与地震波的前进方向平行的波为纵波，粒子振动方向与地震波的前进方向垂直的波为横波，在应用地震勘探技术时，例如在进行叠前地震反演研究时，需要获得横波测井曲线。

发明内容

鉴于上述问题，提出了本申请以便提供一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法。

本申请的实施例提供一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法，包括如下步骤：S1：测量砂岩型铀矿勘查区，确定测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线；S2：根据测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线，确定计算横波测井曲线；S3：对计算横波测井曲线进行高频滤波处理，确定计算横波的高频测井曲线；S4：对测量纵波测井曲线进行低频滤波处理，确定测量纵波的低频测井曲线；S5：将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线；S6：根据预测横波测井曲线的低频曲线和计算横波的高频测井曲线，确定预测横波测井曲线。

本申请实施例中确定砂岩型铀矿勘查区的预测横波测井曲线的方法，可以准确地确定预测横波测井曲线。

附图说明

图1为根据本申请的实施例的确定砂岩型铀矿勘查区的预测横波测井曲线的方法的流程图；

图2为根据本申请的实施例的确定计算横波测井曲线的方法的流程图；

图3为根据本申请的实施例的确定计算横波测井曲线的方法的流程图；

图4为根据本申请的实施例的将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线的方法的流程图；

图5为根据本申请的实施例的确定横波测井曲线和纵波测井曲线的转换关系式的方法的流程图；

图6为根据本申请的实施例确定的松辽盆地海力锦勘查区ZKL7-8井中的计算纵波测井曲线和计算横波测井曲线的数据交会图；

图7为根据本申请的实施例确定的松辽盆地海力锦勘查区ZKL7-8井中的预测横波测井曲线、预测横波测井曲线的低频曲线及计算横波的高频测井曲线的示意图。

图中：1、预测横波测井曲线；2、预测横波测井曲线的低频曲线；3、计算横波的高频测井曲线。

还应该注意的是，附图只是为了便于描述优选实施例，而不是本申请本身。附图没有示出所描述的实施例的每个方面，并且不限制本申请的范围。

具体实施方式

为使本申请的目的、技术方案和优点更加清楚，下面将结合本申请实施例的附图，对本申请的技术方案进行清楚、完整地描述。显然，所描述的实施例是本申请的一个实施例，而不是全部的实施例。基于所描述的本申请的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。

需要说明的是，除非另外定义，本申请使用的技术术语或者科学术语应当为本申请所属领域内具有一般技能的人士所理解的通常意义。若全文中涉及“第一”、“第二”等描述，则该“第一”、“第二”等描述仅用于区别类似的对象，而不能理解为指示或暗示其相对重要性、先后次序或者隐含指明所指示的技术特征的数量，应该理解为“第一”、“第二”等描述的数据在适当情况下可以互换。若全文中出现“和/或”，其含义为包括三个并列方案，以“A和/或B”为例，包括A方案，或B方案，或A和B同时满足的方案。

本申请的发明人发现，在对砂岩型铀矿进行勘测时，由于横波测井曲线并未被列入常规测井参数中，因此，无法通过勘测而直接获得横波测井曲线。

为此，本申请的实施例提供一种准确的确定砂岩型铀矿勘查区的预测横波测井曲线的方法，如图1所示出的根据本申请的实施例的一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法的流程图，其中，包括如下步骤：S1：测量砂岩型铀矿勘查区，确定测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线；S2：根据测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线，确定计算横波测井曲线；S3：对计算横波测井曲线进行高频滤波处理，确定计算横波的高频测井曲线；S4：对测量纵波测井曲线进行低频滤波处理，确定测量纵波的低频测井曲线；S5：将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线；S6：根据预测横波测井曲线的低频曲线和计算横波的高频测井曲线，确定预测横波测井曲线。

其中，伽马测井曲线中的伽马值可以通过测量砂岩型铀矿勘查区中的各个位置的伽马射线强度而得到。电阻率测井曲线中的电阻率值可以通过测量砂岩型铀矿勘查区中的各个位置的岩石(包括其中的流体)的电阻率而得到。测量纵波测井曲线可以通过声波测井技术得到。

在一些实施例中，在S1步骤中，伽马测井曲线和电阻率测井曲线可以从砂岩型铀矿勘查区的现有资料中收集。

在一些实施例中，在S2步骤中，如图2所示出的根据本申请的实施例的确定计算横波测井曲线的方法的流程图，其还包括如下步骤：S21：根据伽马测井曲线和电阻率测井曲线，确定泥质含量测井曲线；S22：根据测量纵波测井曲线及泥质含量测井曲线，确定孔隙度测井曲线；S23：根据孔隙度测井曲线和泥质含量测井曲线，确定计算横波测井曲线。

其中，孔隙度测井曲线中的孔隙度，是指岩石中所有孔隙空间体积之和与该岩石体积的比值。

砂岩型铀矿勘查区中包括非铀矿化段和铀矿化段，在一些实施例中，在S21步骤中，还包括如下步骤以确定砂岩型铀矿勘查区的泥质含量测井曲线：根据以下表达式，确定非铀矿化段的泥质含量测井曲线：

δ

在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以用非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线进行交会分析。在一些实施例中，可以利用excel软件对非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线进行交会，绘制散点图并添加趋势线转换关系式。在一些实施例中，还可以利用origin等其他软件进行交会分析。

在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以选择相关性较好的非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率曲线来分析它们的转换关系式，并且可以通过非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的交会图判断相关性，当某非铀矿化段的泥质含量测井曲线和该段的电阻率测井曲线在交会图中存在明显的正相关或负相关时，可以判断出该段的相关性较好。

例如，在松辽盆地海力锦勘查区ZKL7-8井中，483m～503m段为非铀矿化段，确定该段的泥质含量测井曲线后，将该段的泥质含量测井曲线与该段的电阻率曲线进行交会分析，以确定松辽盆地海力锦勘查区ZKL7-8井中的泥质含量测井曲线和电阻率测井曲线的转换关系式，得到的转换关系式为：

V

据此可以得到松辽盆地海力锦勘查区ZKL7-8井中的泥质含量测井曲线。

在一些实施例中，在S22步骤中，根据以下表达式确定孔隙度测井曲线：

Δt＝1000000/V

例如，在松辽盆地海力锦勘查区ZKL7-8井中，砂岩纵波时差为165，泥岩纵波时差为293，饱和流体纵波时差为597，并且非压实泥岩为浅层泥岩，压实泥岩为深层泥岩，该井中的C值为浅层泥岩和深层泥岩的声波时差之比，为1.5625。

在一些实施例中，在S23步骤中，如图3所示出的根据本申请的实施例的确定计算横波测井曲线的方法的流程图，其还包括如下步骤：S231：根据孔隙度测井曲线和泥质含量测井曲线，确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线；S232：根据孔隙度测井曲线和泥质含量测井曲线，确定砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线；S233：根据砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线、孔隙度测井曲线、砂岩孔隙度测井曲线和泥岩孔隙度测井曲线，确定计算横波测井曲线。

在一些实施例中，在S231步骤中，根据以下表达式确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线：

Φ

Φ

在一些实施例中，在S232步骤中，根据以下表达式确定砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线：

T

V

V

V

T

ρ

在一些实施例中，在S233步骤中，根据以下表达式确定计算横波测井曲线：

V

ρ

μ

μ

μ

Q

F

B

K

R

F

F

F

F

F

F

Q

F

F

F

F

F

F

F

在一些实施例中，在S5步骤中，如图4所示出的根据本申请的实施例的将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线的方法的流程图，其还包括如下步骤：S51：根据伽马测井曲线和电阻率测井曲线，确定泥质含量测井曲线；S52：根据测量纵波测井曲线及泥质含量测井曲线，确定孔隙度测井曲线；S53：根据孔隙度测井曲线和泥质含量测井曲线，确定横波测井曲线和纵波测井曲线的转换关系式；S54：根据横波测井曲线和纵波测井曲线的转换关系式，将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线。

在一些实施例中，在S51步骤中，还包括如下步骤以确定砂岩型铀矿勘查区的泥质含量测井曲线：根据以下表达式，确定非铀矿化段的泥质含量测井曲线：

δ

在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以选择相关性较好的非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率曲线来分析它们的转换关系式，并且可以通过非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的交会图判断相关性，当某非铀矿化段的泥质含量测井曲线和该段的电阻率测井曲线在交会图中存在明显的正相关或负相关时，可以判断出该段的相关性较好。

例如，在松辽盆地海力锦勘查区ZKL7-8井中，483m～503m段为非铀矿化段，确定该段的泥质含量测井曲线后，将该段的泥质含量测井曲线与该段的电阻率曲线进行交会分析，以确定松辽盆地海力锦勘查区ZKL7-8井中的泥质含量测井曲线和电阻率测井曲线的转换关系式，得到的转换关系式为：

V

据此可以得到松辽盆地海力锦勘查区ZKL7-8井中的泥质含量测井曲线。在一些实施例中，根据以下表达式确定孔隙度测井曲线：

Δt＝1000000/V

例如，在松辽盆地海力锦勘查区ZKL7-8井中，砂岩纵波时差为165，泥岩纵波时差为293，饱和流体纵波时差为597，并且非压实泥岩为浅层泥岩，压实泥岩为深层泥岩，该井中的C值为浅层泥岩和深层泥岩的声波时差之比，为1.5625。

在一些实施例中，在S53步骤中，如图5所示出的根据本申请的实施例的确定横波测井曲线和纵波测井曲线的转换关系式的方法的流程图，还包括如下步骤：S531：根据孔隙度测井曲线和泥质含量测井曲线，确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线；S532：根据孔隙度测井曲线和泥质含量测井曲线，确定砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线；S533：根据砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线、孔隙度测井曲线、砂岩孔隙度测井曲线和泥岩孔隙度测井曲线，确定计算横波测井曲线和计算纵波测井曲线；S534：根据计算的横波测井曲线和计算纵波测井曲线，确定横波测井曲线和纵波测井曲线的转换关系式X。

在一些实施例中，在步骤S531中，根据以下表达式确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线：

Φ

Φ

在一些实施例中，在步骤S532中，根据以下表达式确定砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线：

T

V

V

V

T

ρ

在一些实施例中，在S533步骤中，根据以下表达式确定计算横波测井曲线和计算纵波测井曲线：

V

ρ

V

μ

μ

μ

Q

F

B

K

R

F

F

F

F

F

F

Q

F

F

F

F

F

F

F

并且，K

K

K

P

P

在一些实施例中，在S534步骤中，在确定横波测井曲线和纵波测井曲线的转换关系式X时，可以用计算纵波测井曲线和计算横波测井曲线进行交会分析。在一些实施例中，可以利用excel软件对计算纵波测井曲线和计算横波测井曲线进行交会，绘制散点图并添加趋势线转换关系式。在一些实施例中，还可以利用origin等其他软件进行交会分析。

例如，在松辽盆地海力锦勘查区ZKL7-8井，用计算纵波测井曲线和计算横波测井曲线进行交会分析，得到图6中示出的根据本申请的实施例确定的松辽盆地海力锦勘查区ZKL7-8井中的计算纵波测井曲线和计算横波测井曲线的交会图，根据上述交会图，确定横波测井曲线和纵波测井曲线的转换关系式X为：

V

在一些实施例中，在S3步骤中，可以根据以下方法对计算横波测井曲线进行高频率波处理：利用带通滤波器对对计算横波测井曲线进行高频滤波处理，将带通滤波器的截止频率设置为40，40，9999，9999，滤波得到计算横波的高频测井曲线。

在一些实施例中，在S4步骤中，可以根据以下方法对测量的纵波测井曲线进行低频滤波处理：利用带通滤波器对测量纵波测井曲线进行低频滤波处理，将带通滤波器的截止频率设置为0，0，35，40，滤波得到测量纵波的低频测井曲线。

在一些实施例中，在S6步骤中，根据以下表达式确定预测横波测井曲线：

V

在一些实施例中，转换关系式X为：V

图7中示出了根据本申请的实施例确定的松辽盆地海力锦勘查区ZKL7-8井中的预测横波测井曲线1、预测横波测井曲线的低频曲线2及计算横波的高频测井曲线3。

上面结合附图和实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。