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

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


背景技术：

2.砂岩型铀矿是指铀矿床产于砂岩、砂砾岩等中粗碎屑岩的铀矿。在对砂岩型铀矿进行勘测时，主要采用地震勘探技术。在地震勘探技术中，主要通过观测和分析大地对人工激发地震波的响应，推断出地下岩层性质和形态。
3.地震波可以分为横波和纵波，粒子振动方向与地震波的前进方向平行的波为纵波，粒子振动方向与地震波的前进方向垂直的波为横波，在应用地震勘探技术时，例如在进行叠前地震反演研究时，需要获得横波测井曲线。


技术实现要素：

4.鉴于上述问题，提出了本技术以便提供一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法。
5.本技术的实施例提供一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法，包括如下步骤：s1：测量砂岩型铀矿勘查区，确定测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线；s2：根据测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线，确定计算横波测井曲线；s3：对计算横波测井曲线进行高频滤波处理，确定计算横波的高频测井曲线；s4：对测量纵波测井曲线进行低频滤波处理，确定测量纵波的低频测井曲线；s5：将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线；s6：根据预测横波测井曲线的低频曲线和计算横波的高频测井曲线，确定预测横波测井曲线。
6.本技术实施例中确定砂岩型铀矿勘查区的预测横波测井曲线的方法，可以准确地确定预测横波测井曲线。
附图说明
7.图1为根据本技术的实施例的确定砂岩型铀矿勘查区的预测横波测井曲线的方法的流程图；
8.图2为根据本技术的实施例的确定计算横波测井曲线的方法的流程图；
9.图3为根据本技术的实施例的确定计算横波测井曲线的方法的流程图；
10.图4为根据本技术的实施例的将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线的方法的流程图；
11.图5为根据本技术的实施例的确定横波测井曲线和纵波测井曲线的转换关系式的方法的流程图；
12.图6为根据本技术的实施例确定的松辽盆地海力锦勘查区zkl7-8井中的计算纵波测井曲线和计算横波测井曲线的数据交会图；
13.图7为根据本技术的实施例确定的松辽盆地海力锦勘查区zkl7-8井中的预测横波测井曲线、预测横波测井曲线的低频曲线及计算横波的高频测井曲线的示意图。
14.图中：1、预测横波测井曲线；2、预测横波测井曲线的低频曲线；3、计算横波的高频测井曲线。
15.还应该注意的是，附图只是为了便于描述优选实施例，而不是本技术本身。附图没有示出所描述的实施例的每个方面，并且不限制本技术的范围。
具体实施方式
16.为使本技术的目的、技术方案和优点更加清楚，下面将结合本技术实施例的附图，对本技术的技术方案进行清楚、完整地描述。显然，所描述的实施例是本技术的一个实施例，而不是全部的实施例。基于所描述的本技术的实施例，本领域普通技术人员在无需创造性劳动的前提下所获得的所有其他实施例，都属于本技术保护的范围。
17.需要说明的是，除非另外定义，本技术使用的技术术语或者科学术语应当为本技术所属领域内具有一般技能的人士所理解的通常意义。若全文中涉及“第一”、“第二”等描述，则该“第一”、“第二”等描述仅用于区别类似的对象，而不能理解为指示或暗示其相对重要性、先后次序或者隐含指明所指示的技术特征的数量，应该理解为“第一”、“第二”等描述的数据在适当情况下可以互换。若全文中出现“和/或”，其含义为包括三个并列方案，以“a和/或b”为例，包括a方案，或b方案，或a和b同时满足的方案。
18.本技术的发明人发现，在对砂岩型铀矿进行勘测时，由于横波测井曲线并未被列入常规测井参数中，因此，无法通过勘测而直接获得横波测井曲线。
19.为此，本技术的实施例提供一种准确的确定砂岩型铀矿勘查区的预测横波测井曲线的方法，如图1所示出的根据本技术的实施例的一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法的流程图，其中，包括如下步骤：s1：测量砂岩型铀矿勘查区，确定测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线；s2：根据测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线，确定计算横波测井曲线；s3：对计算横波测井曲线进行高频滤波处理，确定计算横波的高频测井曲线；s4：对测量纵波测井曲线进行低频滤波处理，确定测量纵波的低频测井曲线；s5：将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线；s6：根据预测横波测井曲线的低频曲线和计算横波的高频测井曲线，确定预测横波测井曲线。
20.其中，伽马测井曲线中的伽马值可以通过测量砂岩型铀矿勘查区中的各个位置的伽马射线强度而得到。电阻率测井曲线中的电阻率值可以通过测量砂岩型铀矿勘查区中的各个位置的岩石(包括其中的流体)的电阻率而得到。测量纵波测井曲线可以通过声波测井技术得到。
21.在一些实施例中，在s1步骤中，伽马测井曲线和电阻率测井曲线可以从砂岩型铀矿勘查区的现有资料中收集。
22.在一些实施例中，在s2步骤中，如图2所示出的根据本技术的实施例的确定计算横波测井曲线的方法的流程图，其还包括如下步骤：s21：根据伽马测井曲线和电阻率测井曲线，确定泥质含量测井曲线；s22：根据测量纵波测井曲线及泥质含量测井曲线，确定孔隙度测井曲线；s23：根据孔隙度测井曲线和泥质含量测井曲线，确定计算横波测井曲线。
23.其中，孔隙度测井曲线中的孔隙度，是指岩石中所有孔隙空间体积之和与该岩石体积的比值。
24.砂岩型铀矿勘查区中包括非铀矿化段和铀矿化段，在一些实施例中，在s21步骤中，还包括如下步骤以确定砂岩型铀矿勘查区的泥质含量测井曲线：根据以下表达式，确定非铀矿化段的泥质含量测井曲线：
25.其中，v
sh1i
为非铀矿化段的泥质含量测井曲线中序号为i的值，a值为3至3.7，δi根据以下表达式确定：
26.δi＝(gr
i-gr
min
)/(gr
max-gr
min
)，其中，gri为伽马测井曲线中序号为i的值，gr
min
为伽马测井曲线中的最小值，gr
max
为伽马测井曲线中的最大值；确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式；根据非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式，在铀矿化段将电阻率测井曲线转换为泥质含量测井曲线。通过上述步骤可以方便快捷地计算出整个砂岩型铀矿勘查区中钻孔的泥质含量测井曲线，而不用探测铀矿化段的泥质含量测井曲线。
27.在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以用非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线进行交会分析。在一些实施例中，可以利用excel软件对非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线进行交会，绘制散点图并添加趋势线转换关系式。在一些实施例中，还可以利用origin等其他软件进行交会分析。
28.在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以选择相关性较好的非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率曲线来分析它们的转换关系式，并且可以通过非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的交会图判断相关性，当某非铀矿化段的泥质含量测井曲线和该段的电阻率测井曲线在交会图中存在明显的正相关或负相关时，可以判断出该段的相关性较好。
29.例如，在松辽盆地海力锦勘查区zkl7-8井中，483m～503m段为非铀矿化段，确定该段的泥质含量测井曲线后，将该段的泥质含量测井曲线与该段的电阻率曲线进行交会分析，以确定松辽盆地海力锦勘查区zkl7-8井中的泥质含量测井曲线和电阻率测井曲线的转换关系式，得到的转换关系式为：
[0030]vshi
＝-0.007rti+0.2399，其中，v
shi
为泥质含量测井曲线中第i个值，rti为电阻率测井曲线中第i个值。然后在铀矿化段512-552m，利用转换关系式，将rti转换得到新的v
sh_trani
，最后在铀矿化段用新的转换的v
sh_trani
替换原来该段的v
shi
值，形成新的泥质含量测井曲线v
sh
。
[0031]
据此可以得到松辽盆地海力锦勘查区zkl7-8井中的泥质含量测井曲线。
[0032]
在一些实施例中，在s22步骤中，根据以下表达式确定孔隙度测井曲线：
[0033]
其中，φi为孔隙度测井曲线中第i个值，v
shi
为泥质含量测井曲线中第i个值，δt为纵波时差，并通过以下表达式确定：
[0034]
δt＝1000000/v
p_wave
，其中，v
p_wave
为纵波速度值，δts为砂岩纵波时差，δtc为泥岩纵波时差，δtf为饱和流体纵波时差，c为非压实泥岩和压实泥岩纵波时差之比。
[0035]
例如，在松辽盆地海力锦勘查区zkl7-8井中，砂岩纵波时差为165，泥岩纵波时差为293，饱和流体纵波时差为597，并且非压实泥岩为浅层泥岩，压实泥岩为深层泥岩，该井中的c值为浅层泥岩和深层泥岩的声波时差之比，为1.5625。
[0036]
在一些实施例中，在s23步骤中，如图3所示出的根据本技术的实施例的确定计算横波测井曲线的方法的流程图，其还包括如下步骤：s231：根据孔隙度测井曲线和泥质含量测井曲线，确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线；s232：根据孔隙度测井曲线和泥质含量测井曲线，确定砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线；s233：根据砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线、孔隙度测井曲线、砂岩孔隙度测井曲线和泥岩孔隙度测井曲线，确定计算横波测井曲线。
[0037]
在一些实施例中，在s231步骤中，根据以下表达式确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线：
[0038]
φ
ci
＝v
shi
φi/(1-φi)，其中，φ
ci
为泥岩孔隙度测井曲线中第i个值，φi为孔隙度测井曲线中第i个值，v
shi
为泥质含量测井曲线中第i个值；
[0039]
φi＝φ
si
+φ
ci
，其中，φ
si
为砂岩孔隙度测井曲线中第i个值。
[0040]
在一些实施例中，在s232步骤中，根据以下表达式确定砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线：
[0041]
t
pmi
＝165v
smi
+293v
shmi
，其中，t
pmi
为砂泥岩石混合基质的纵波时差曲线中第i个值，v
smi
根据以下表达式确定：
[0042]vsmi
＝v
si
/(1-φi)，其中，φi为孔隙度测井曲线中第i个值，v
si
根据以下关系式确定：
[0043]vsi
+v
shi
+φi＝1，其中，v
shi
为泥质含量测井曲线中第i个值，v
shmi
为砂泥岩石混合基质的泥质含量测井曲线中第i个值，根据以下表达式确定：
[0044]vshmi
＝v
shi
/(1-φi)；
[0045]
t
smi
＝244v
smi
+609v
shmi
，其中，t
smi
为砂泥岩石混合基质的横波时差曲线中第i个值；
[0046]
ρ
mi
＝2.65v
smi
+2.6v
shmi
，其中，ρ
mi
为砂泥岩石混合基质的密度曲线中第i个值。
[0047]
在一些实施例中，在s233步骤中，根据以下表达式确定计算横波测井曲线：
[0048]vs_sati
＝(μ
si
/ρ
si
)
1/2
，其中，v
s_sati
为计算横波测井曲线中第i个值，ρ
si
为砂岩基质的密度曲线中第i个值，根据以下表达式确定：
[0049]
ρ
si
＝ρ
mi
(1-φi)+1.10958φi，其中，ρ
mi
为砂泥岩石混合基质的密度曲线中第i个值，φi为孔隙度测井曲线中第i个值，并且，μ
si
根据以下表达式确定：
[0050]
μ
si
＝μ
di
，μ
di
根据以下表达式确定：
[0051]
μ
di
＝μ
mi
(1-φ
si
)
qsi
(1-φ
ci
)
qci
，其中，φ
si
为砂岩孔隙度测井曲线中第i个值，φ
ci
为泥岩孔隙度测井曲线中第i个值，并且μ
mi
根据以下表达式确定：
[0052]
μ
mi
＝ρ
mi
(1/t
smi2
)，其中，t
smi
为砂泥岩石混合基质的横波时差曲线中第i个值，ρ
mi
为砂泥岩石混合基质的密度曲线中第i个值，并且，q
si
根据以下表达式确定：
[0053]qsi
＝(2/f
2si
+1/f
4si
+(f
4sif5si
+f
6sif7si-f
8sif9si
)/f
2sif4si
)/5，其中，f
2si
根据以下表达式确定：
[0054]f2si
＝-3(gs+hs)/2+ri(3gs+5hs)/2+bi(3-2ri)-(3b
i-1)(3-4ri)(gs+hs-ri(gs-hs+
2hs2))/2，其中，gs＝-0.0324，hs＝0.197414，bi根据以下表达式确定：
[0055]bi
＝3.17453/3k
mi
，其中，k
mi
为砂泥岩石混合基质的弹性模量曲线中第i个值，根据以下表达式确定：
[0056]kmi
＝ρ
mi
[1/t
pmi2-(4/3)(1/t
smi2
)]，其中，t
pmi
为砂泥岩石混合基质的纵波时差曲线中第i个值，并且，ri根据以下表达式确定：
[0057]ri
＝3μ
mi
/(3k
mi
+4μ
mi
)，并且，f
4si
根据以下表达式确定：
[0058]f4si
＝1-(gs+3hs-ri(gs-hs))/4，并且，f
5si
根据以下表达式确定：
[0059]f5si
＝(gs-ri(gs+hs-4/3))+bihs(3-4ri)，并且，f
6si
根据以下表达式确定：
[0060]f6si
＝ri(gs+hs)-gs+bi(1-h)(3-4ri)，并且，f
7si
根据以下表达式确定：
[0061]f7si
＝2-(3gs+9hs-ri(3gs+5hs))/4+bihs(3-4ri)，并且，f
8si
根据以下表达式确定：
[0062]f8si
＝-(1-2ri+gs(r
i-1)/2+hs(5r
i-3)/2)+b(1-hs)(3-4ri)，并且，f
9si
根据以下表达式确定：
[0063]f9si
＝-(gs(r
i-1)-rihs)+bihs(3-4ri)，并且，q
ci
根据以下表达式确定：
[0064]qci
＝(2/f
2ci
+1/f
4ci
+(f
4cif5ci
+f
6cif7ci-f
8cif9ci
)/f
2cif4ci
)/5，其中，f
2ci
根据以下表达式确定：
[0065]f2ci
＝-3(gc+hc)/2+ri(3gc+5hc)/2+bi(3-2ri)-(3b
i-1)(3-4ri)(gc+hc-ri(gc-hc+2hc2))/2，其中，gc＝-0.00446，hc＝0.073819，并且，f
4ci
根据以下表达式确定：
[0066]f4ci
＝1-(gc+3hc-ri(gc-hc))/4，并且，f
5ci
根据以下表达式确定：
[0067]f5ci
＝(gc-ri(gc+hc-4/3))+bihc(3-4ri)，并且，f
6ci
根据以下表达式确定：
[0068]f6ci
＝ri(gc+hc)-gc+bi(1-h)(3-4ri)，并且，f
7ci
根据以下表达式确定：
[0069]f7ci
＝2-(3gc+9hc-ri(3gc+5hc))/4+bihc(3-4ri)，并且，f
8ci
根据以下表达式确定：
[0070]f8ci
＝-(1-2ri+gc(r
i-1)/2+hc(5r
i-3)/2)+b(1-hc)(3-4ri)，并且，f
9ci
根据以下表达式确定：
[0071]f9ci
＝-(gc(r
i-1)-rihc)+bihc(3-4ri)。
[0072]
在一些实施例中，在s5步骤中，如图4所示出的根据本技术的实施例的将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线的方法的流程图，其还包括如下步骤：s51：根据伽马测井曲线和电阻率测井曲线，确定泥质含量测井曲线；s52：根据测量纵波测井曲线及泥质含量测井曲线，确定孔隙度测井曲线；s53：根据孔隙度测井曲线和泥质含量测井曲线，确定横波测井曲线和纵波测井曲线的转换关系式；s54：根据横波测井曲线和纵波测井曲线的转换关系式，将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线。
[0073]
在一些实施例中，在s51步骤中，还包括如下步骤以确定砂岩型铀矿勘查区的泥质含量测井曲线：根据以下表达式，确定非铀矿化段的泥质含量测井曲线：
[0074]
其中，v
sh1i
为非铀矿化段的泥质含量测井曲线中序号为i的值，a值为3至3.7，δi根据以下表达式确定：
[0075]
δi＝(gr
i-gr
min
)/(gr
max-gr
min
)，其中，gri为伽马测井曲线中序号为i的值，gr
min
为伽马测井曲线中的最小值，gr
max
为伽马测井曲线中的最大值；确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式；根据非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式，在铀矿化段将电阻率测井曲线转换为
泥质含量测井曲线。通过上述步骤可以方便快捷地计算出整个砂岩型铀矿勘查区中的泥质含量测井曲线，而不用探测铀矿化段的泥质含量测井曲线。在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以用非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线进行交会分析。在一些实施例中，可以利用excel软件对非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线进行交会，绘制散点图并添加趋势线转换关系式。在一些实施例中，还可以利用origin等其他软件进行交会分析。
[0076]
在一些实施例中，在确定非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的转换关系式时，可以选择相关性较好的非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率曲线来分析它们的转换关系式，并且可以通过非铀矿化段的泥质含量测井曲线和非铀矿化段的电阻率测井曲线的交会图判断相关性，当某非铀矿化段的泥质含量测井曲线和该段的电阻率测井曲线在交会图中存在明显的正相关或负相关时，可以判断出该段的相关性较好。
[0077]
例如，在松辽盆地海力锦勘查区zkl7-8井中，483m～503m段为非铀矿化段，确定该段的泥质含量测井曲线后，将该段的泥质含量测井曲线与该段的电阻率曲线进行交会分析，以确定松辽盆地海力锦勘查区zkl7-8井中的泥质含量测井曲线和电阻率测井曲线的转换关系式，得到的转换关系式为：
[0078]vshi
＝-0.007rti+0.2399，其中，v
shi
为泥质含量测井曲线中第i个值，rti为电阻率测井曲线中第i个值。然后在铀矿化段512-552m，利用转换关系式，将rti转换得到新的v
sh_trani
，最后在铀矿化段用新的转换的v
sh_trani
替换原来该段的v
shi
值，形成新的泥质含量测井曲线v
sh
。
[0079]
据此可以得到松辽盆地海力锦勘查区zkl7-8井中的泥质含量测井曲线。在一些实施例中，根据以下表达式确定孔隙度测井曲线：
[0080]
其中，φi为孔隙度测井曲线中第i个值，v
shi
为泥质含量测井曲线中第i个值，δt为纵波时差，并通过以下表达式确定：
[0081]
δt＝1000000/v
p_wave
，其中，v
p_wave
为测量纵波测井曲线中数据的平均值，δts为砂岩纵波时差，δtc为泥岩纵波时差，δtf为饱和流体纵波时差，c为非压实泥岩和压实泥岩纵波时差之比。
[0082]
例如，在松辽盆地海力锦勘查区zkl7-8井中，砂岩纵波时差为165，泥岩纵波时差为293，饱和流体纵波时差为597，并且非压实泥岩为浅层泥岩，压实泥岩为深层泥岩，该井中的c值为浅层泥岩和深层泥岩的声波时差之比，为1.5625。
[0083]
在一些实施例中，在s53步骤中，如图5所示出的根据本技术的实施例的确定横波测井曲线和纵波测井曲线的转换关系式的方法的流程图，还包括如下步骤：s531：根据孔隙度测井曲线和泥质含量测井曲线，确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线；s532：根据孔隙度测井曲线和泥质含量测井曲线，确定砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线；s533：根据砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线、孔隙度测井曲线、砂岩孔隙度测井曲线和泥岩孔隙度测井曲线，确定计算横波测井曲线和计算纵波测井曲线；s534：根据计算的横波测井曲线和计算纵波测井曲线，确定横波测
井曲线和纵波测井曲线的转换关系式x。
[0084]
在一些实施例中，在步骤s531中，根据以下表达式确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线：
[0085]
φ
ci
＝v
shi
φi/(1-φi)，其中，φ
ci
为泥岩孔隙度测井曲线中第i个值，φi为孔隙度测井曲线中第i个值，v
shi
为泥质含量测井曲线中第i个值；
[0086]
φi＝φ
si
+φ
ci
，其中，φ
si
为砂岩孔隙度测井曲线中第i个值。
[0087]
在一些实施例中，在步骤s532中，根据以下表达式确定砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线：
[0088]
t
pmi
＝165v
smi
+293v
shmi
，其中，t
pmi
为砂泥岩石混合基质的纵波时差曲线中第i个值，v
smi
根据以下表达式确定：
[0089]vsmi
＝v
si
/(1-φi)，其中，φi为孔隙度测井曲线中第i个值，v
si
根据以下关系式确定：
[0090]vsi
+v
shi
+φi＝1，其中，v
shi
为泥质含量测井曲线中第i个值，v
shmi
为砂泥岩石混合基质的泥质含量测井曲线中第i个值，根据以下表达式确定：
[0091]vshmi
＝v
shi
/(1-φi)；
[0092]
t
smi
＝244v
smi
+609v
shmi
，其中，t
smi
为砂泥岩石混合基质的横波时差曲线中序号为i的值；
[0093]
ρ
mi
＝2.65v
smi
+2.6v
shmi
，其中，ρ
mi
为砂泥岩石混合基质的密度曲线中第i个值。
[0094]
在一些实施例中，在s533步骤中，根据以下表达式确定计算横波测井曲线和计算纵波测井曲线：
[0095]vs_sati
＝(μ
si
/ρ
si
)
1/2
，其中，v
s_sati
为计算横波测井曲线中第i个值，ρ
si
为砂岩基质的密度曲线中第i个值，根据以下表达式确定：
[0096]
ρ
si
＝ρ
mi
(1-φi)+1.10958φi，其中，ρ
mi
为砂泥岩石混合基质的密度曲线中第i个值，φi为孔隙度测井曲线中第i个值；
[0097]vp_sati
＝((k
si
+4μ
si
/3)/ρ
si
)
1/2
，其中，v
p_sati
为计算纵波测井曲线中第i个值，并且，μs根据以下表达式确定：
[0098]
μ
si
＝μ
di
，μ
di
根据以下表达式确定：
[0099]
μ
di
＝μ
mi
(1-φ
si
)
qsi
(1-φ
ci
)
qci
，其中，φ
si
为砂岩孔隙度测井曲线中第i个值，φ
ci
为泥岩孔隙度测井曲线中第i个值，并且μ
mi
根据以下表达式确定：
[0100]
μ
mi
＝ρ
mi
(1/t
smi2
)，其中，t
smi
为砂泥岩石混合基质的横波时差曲线中第i个值，ρ
mi
为砂泥岩石混合基质的密度曲线中第i个值，并且，q
si
根据以下表达式确定：
[0101]qsi
＝(2/f
2si
+1/f
4si
+(f
4sif5si
+f
6sif7si-f
8sif9si
)/f
2sif4si
)/5，其中，f
2si
根据以下表达式确定：
[0102]f2si
＝-3(gs+hs)/2+ri(3gs+5hs)/2+bi(3-2ri)-(3b
i-1)(3-4ri)(gs+hs-ri(gs-hs+2hs2))/2，其中，gs＝-0.0324，hs＝0.197414，bi根据以下表达式确定：
[0103]bi
＝3.17453/3k
mi
，其中，k
mi
为砂泥岩石混合基质的弹性模量中第i个值，根据以下表达式确定：
[0104]kmi
＝ρ
mi
[1/t
pmi2-(4/3)(1/t
smi2
)]，其中，t
pmi
为砂泥岩石混合基质的纵波时差曲线中第i个值，并且，ri根据以下表达式确定：
[0105]ri
＝3μ
mi
/(3k
mi
+4μ
mi
)，并且，f
4si
根据以下表达式确定：
[0106]f4si
＝1-(gs+3hs-ri(gs-hs))/4，并且，f
5si
根据以下表达式确定：
[0107]f5si
＝(gs-ri(gs+hs-4/3))+bihs(3-4ri)，并且，f
6si
根据以下表达式确定：
[0108]f6si
＝ri(gs+hs)-gs+bi(1-h)(3-4ri)，并且，f
7si
根据以下表达式确定：
[0109]f7si
＝2-(3gs+9hs-ri(3gs+5hs))/4+bihs(3-4ri)，并且，f
8si
根据以下表达式确定：
[0110]f8si
＝-(1-2ri+gs(r
i-1)/2+hs(5r
i-3)/2)+b(1-hs)(3-4ri)，并且，f
9si
根据以下表达式确定：
[0111]f9si
＝-(gs(r
i-1)-rihs)+bihs(3-4ri)，并且，q
ci
根据以下表达式确定：
[0112]qci
＝(2/f
2ci
+1/f
4ci
+(f
4cif5ci
+f
6cif7ci-f
8cif9ci
)/f
2cif4ci
)/5，其中，f
2ci
根据以下表达式确定：
[0113]f2ci
＝-3(gc+hc)/2+ri(3gc+5hc)/2+bi(3-2ri)-(3b
i-1)(3-4ri)(gc+hc-ri(gc-hc+2hc2))/2，其中，gc＝-0.00446，hc＝0.073819，并且，f
4ci
根据以下表达式确定：
[0114]f4ci
＝1-(gc+3hc-ri(gc-hc))/4，并且，f
5ci
根据以下表达式确定：
[0115]f5ci
＝(gc-ri(gc+hc-4/3))+bihc(3-4ri)，并且，f
6ci
根据以下表达式确定：
[0116]f6ci
＝ri(gc+hc)-gc+bi(1-h)(3-4ri)，并且，f
7ci
根据以下表达式确定：
[0117]f7ci
＝2-(3gc+9hc-ri(3gc+5hc))/4+bihc(3-4ri)，并且，f
8ci
根据以下表达式确定：
[0118]f8ci
＝-(1-2ri+gc(r
i-1)/2+hc(5r
i-3)/2)+b(1-hc)(3-4ri)，并且，f
9ci
根据以下表达式确定：
[0119]f9ci
＝-(gc(r
i-1)-rihc)+bihc(3-4ri)；
[0120]
并且，ks根据以下表达式确定：
[0121]ksi
＝k
di
+(1-k
di2
/k
mi2
)/(φi/kf+(1-φi)/k
mi-k
di
/k
mi2
)，其中，φi为孔隙度测井曲线中第i个值，k
fi
为孔隙流体的体积模量曲线中第i个值，其为根据岩石物理相关手册查询到的固定值，k
di
为干岩石基质的弹性模量曲线中第i个值，根据以下表达式确定：
[0122]kdi
＝k
mi
(1-φ
si
)
psi
(1-φ
ci
)
pci
，其中，p
si
根据以下表达式确定：
[0123]
p
si
＝f
1si
/f
2si
，并且，p
ci
根据以下表达式确定：
[0124]
p
ci
＝f
1ci
/f
2ci
。
[0125]
在一些实施例中，在s534步骤中，在确定横波测井曲线和纵波测井曲线的转换关系式x时，可以用计算纵波测井曲线和计算横波测井曲线进行交会分析。在一些实施例中，可以利用excel软件对计算纵波测井曲线和计算横波测井曲线进行交会，绘制散点图并添加趋势线转换关系式。在一些实施例中，还可以利用origin等其他软件进行交会分析。
[0126]
例如，在松辽盆地海力锦勘查区zkl7-8井，用计算纵波测井曲线和计算横波测井曲线进行交会分析，得到图6中示出的根据本技术的实施例确定的松辽盆地海力锦勘查区zkl7-8井中的计算纵波测井曲线和计算横波测井曲线的交会图，根据上述交会图，确定横波测井曲线和纵波测井曲线的转换关系式x为：
[0127]vs_sati
＝0.8114v
p_sati-0.8957，其中，v
s_sati
为横波测井曲线中第i个值，v
p_sati
为纵波测井曲线中第i个值。
[0128]
在一些实施例中，在s3步骤中，可以根据以下方法对计算横波测井曲线进行高频率波处理：利用带通滤波器对对计算横波测井曲线进行高频滤波处理，将带通滤波器的截止频率设置为40，40，9999，9999，滤波得到计算横波的高频测井曲线。
[0129]
在一些实施例中，在s4步骤中，可以根据以下方法对测量的纵波测井曲线进行低频滤波处理：利用带通滤波器对测量纵波测井曲线进行低频滤波处理，将带通滤波器的截止频率设置为0，0，35，40，滤波得到测量纵波的低频测井曲线。
[0130]
在一些实施例中，在s6步骤中，根据以下表达式确定预测横波测井曲线：
[0131]vs_predictedi
＝v
s_lowi
+v
s_highi
，其中，v
s_predictedi
为预测横波测井曲线中第i个值，v
s_lowi
为预测横波测井曲线的低频曲线中第i个值，v
s_highi
为计算横波的高频测井曲线中第i个值。
[0132]
在一些实施例中，转换关系式x为：v
s_sati
＝0.8114v
p_sati-0.8957，其中，v
s_sati
为横波测井曲线中第i个值，v
p_sati
为纵波测井曲线中第i个值。利用带通滤波器对测量的纵波测井曲线进行低频滤波处理，将带通滤波器的截止频率设置为0，0，35，40，滤波得到测量纵波的低频测井曲线；利用转换关系式x，将测量纵波的低频测井曲线转换为横波低频测井曲线；利用带通滤波器对计算横波测井曲线进行高频滤波处理，将带通滤波器的截止频率设置为40，40，9999，9999，滤波得到计算横波的高频测井曲线；根据以下表达式确定所述预测横波测井曲线：v
s_predictedi
＝v
s_lowi
+v
s_highi
，其中，v
s_predictedi
为预测横波测井曲线中第i个值，v
s_lowi
为预测横波测井曲线的低频曲线中第i个值，v
s_highi
为计算横波的高频测井曲线中第i个值。
[0133]
图7中示出了根据本技术的实施例确定的松辽盆地海力锦勘查区zkl7-8井中的预测横波测井曲线1、预测横波测井曲线的低频曲线2及计算横波的高频测井曲线3。
[0134]
上面结合附图和实施例对本发明作了详细说明，但是本发明并不限于上述实施例，在本领域普通技术人员所具备的知识范围内，还可以在不脱离本发明宗旨的前提下作出各种变化。本发明中未作详细描述的内容均可以采用现有技术。

技术特征：
1.一种确定砂岩型铀矿勘查区的预测横波测井曲线的方法，其中，包括如下步骤：s1：测量所述砂岩型铀矿勘查区，确定测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线；s2：根据所述测量纵波测井曲线、所述伽马测井曲线以及所述电阻率测井曲线，确定计算横波测井曲线；s3：对所述计算横波测井曲线进行高频滤波处理，确定计算横波的高频测井曲线；s4：对所述测量纵波测井曲线进行低频滤波处理，确定测量纵波的低频测井曲线；s5：将所述测量纵波的低频测井曲线转换成所述预测横波测井曲线的低频曲线；s6：根据所述预测横波测井曲线的低频曲线和所述计算横波的高频测井曲线，确定所述预测横波测井曲线。2.根据权利要求1所述的方法，其中，在s2步骤中，还包括如下步骤：根据所述伽马测井曲线和所述电阻率测井曲线，确定泥质含量测井曲线；根据所述测量纵波测井曲线及所述泥质含量测井曲线，确定孔隙度测井曲线；根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定所述计算横波测井曲线。3.根据权利要求2所述的方法，其中，所述根据所述伽马测井曲线和所述电阻率测井曲线，确定泥质含量测井曲线的步骤中，还包括如下步骤：根据以下表达式，确定非铀矿化段的泥质含量测井曲线：其中，v
sh1i
为所述非铀矿化段的泥质含量测井曲线中第i个值，a值为3至3.7，δ
i
根据以下表达式确定：δ
i
＝(gr
i-gr
min
)/(gr
max-gr
min
)，其中，gr
i
为所述伽马测井曲线中第i个值，gr
min
为所述伽马测井曲线中的最小值，gr
max
为所述伽马测井曲线中的最大值；确定所述非铀矿化段的泥质含量测井曲线和所述非铀矿化段的电阻率测井曲线的转换关系式；根据所述非铀矿化段的泥质含量测井曲线和所述非铀矿化段的电阻率测井曲线的转换关系式，在铀矿化段将所述电阻率测井曲线转换为所述泥质含量测井曲线。4.根据权利要求2所述的方法，其中，根据以下表达式确定所述孔隙度测井曲线：其中，φ
i
为所述孔隙度测井曲线中第i个值，v
shi
为所述泥质含量测井曲线中第i个值，δt为纵波时差，并通过以下表达式确定：δt＝1000000/v
p_wave
，其中，v
p_wave
为纵波速度值，δt
s
为砂岩纵波时差，δt
c
为泥岩纵波时差，δt
f
为饱和流体纵波时差，c为非压实泥岩和压实泥岩纵波时差之比。5.根据权利要求2所述的方法，其中，所述根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定所述计算横波测井曲线的步骤，还包括如下步骤：根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线；根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线；根据所述砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线、所述孔隙度
测井曲线、所述砂岩孔隙度测井曲线和所述泥岩孔隙度测井曲线，确定所述计算横波测井曲线。6.根据权利要求5所述的方法，其中，根据以下表达式确定所述砂岩孔隙度测井曲线和所述泥岩孔隙度测井曲线：φ
ci
＝v
shi
φ
i
/(1-φ
i
)，其中，φ
ci
为所述泥岩孔隙度测井曲线中第i个值，φ
i
为所述孔隙度测井曲线中第i个值，v
shi
为所述泥质含量测井曲线中第i个值；φ
i
＝φ
si
+φ
ci
，其中，φ
si
为所述砂岩孔隙度测井曲线中第i个值。7.根据权利要求5所述的方法，其中，根据以下表达式确定所述砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线：t
pmi
＝165v
smi
+293v
shmi
，其中，t
pmi
为所述砂泥岩石混合基质的纵波时差曲线中第i个值，v
smi
根据以下表达式确定：v
smi
＝v
si
/(1-φ
i
)，其中，φ
i
为所述孔隙度测井曲线中第i个值，v
si
根据以下关系式确定：v
si
+v
shi
+φ
i
＝1，其中，v
shi
为所述泥质含量测井曲线中第i个值，并且，v
shmi
为所述砂泥岩石混合基质的泥质含量测井曲线中第i个值，根据以下表达式确定：v
shmi
＝v
shi
/(1-φ
i
)；t
smi
＝244v
smi
+609v
shmi
，其中，t
smi
为所述砂泥岩石混合基质的横波时差曲线中第i个值；ρ
mi
＝2.65v
smi
+2.6v
shmi
，其中，ρ
mi
为所述砂泥岩石混合基质的密度曲线中第i个值。8.根据权利要求1所述的方法，其中，在步骤s5中，还包括如下步骤：根据所述伽马测井曲线和所述电阻率测井曲线，确定泥质含量测井曲线；根据所述测量纵波测井曲线及所述泥质含量测井曲线，确定孔隙度测井曲线；根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定横波测井曲线和纵波测井曲线的转换关系式；根据所述横波测井曲线和纵波测井曲线的转换关系式，将所述测量纵波的低频测井曲线转换成所述预测横波测井曲线的低频曲线。9.根据权利要求8所述的方法，其中，所述根据所述伽马测井曲线和所述电阻率测井曲线，确定泥质含量测井曲线的步骤中，还包括如下步骤：根据以下表达式，确定非铀矿化段的泥质含量测井曲线：其中，v
sh1i
为所述非铀矿化段的泥质含量测井曲线中第i个值，a值为3至3.7，δ
i
根据以下表达式确定：δ
i
＝(gr
i-gr
min
)/(gr
max-gr
min
)，其中，gr
i
为所述伽马测井曲线中第i个值，gr
min
为所述伽马测井曲线中的最小值，gr
max
为所述伽马测井曲线中的最大值；确定所述非铀矿化段泥质含量测井曲线和所述非铀矿化段的电阻率测井曲线的转换关系式；根据所述非铀矿化段的泥质含量测井曲线和所述非铀矿化段的电阻率测井曲线的转换关系式，在铀矿化段将所述电阻率测井曲线转换为所述泥质含量测井曲线。10.根据权利要求8所述的方法，其中，根据以下表达式确定所述孔隙度测井曲线：
其中，φ
i
为所述孔隙度测井曲线中第i个值，v
shi
为所述泥质含量测井曲线中第i个值，δt为纵波时差，并通过以下表达式确定：δt＝1000000/v
p_wave
，其中，v
p_wave
为纵波速度值，δt
s
为砂岩纵波时差，δt
c
为泥岩纵波时差，δt
f
为饱和流体纵波时差，c为非压实泥岩和压实泥岩纵波时差之比。11.根据权利要求8所述的方法，其中，所述根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定所述横波测井曲线和纵波测井曲线的转换关系式的步骤，还包括如下步骤：根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定砂岩孔隙度测井曲线和泥岩孔隙度测井曲线；根据所述孔隙度测井曲线和所述泥质含量测井曲线，确定砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线；根据所述砂泥岩石混合基质的纵波时差曲线、横波时差曲线和密度曲线、所述孔隙度测井曲线、所述砂岩孔隙度测井曲线和所述泥岩孔隙度测井曲线，确定所述计算横波测井曲线和计算纵波测井曲线；根据所述计算横波测井曲线和所述计算纵波测井曲线，确定所述横波测井曲线和纵波测井曲线的转换关系式x。12.根据权利要求11所述的方法，其中，根据以下表达式确定所述砂岩孔隙度测井曲线和所述泥岩孔隙度测井曲线：φ
ci
＝v
shi
φ
i
/(1-φ
i
)，其中，φ
ci
为所述泥岩孔隙度测井曲线中第i个值，φ
i
为所述孔隙度测井曲线中第i个值，v
shi
为所述泥质含量测井曲线中第i个值；φ
i
＝φ
si
+φ
ci
，其中，φ
si
为所述砂岩孔隙度测井曲线中第i个值。13.根据权利要求11所述的方法，其中，根据以下表达式确定所述砂泥岩石混合基质的纵波时差曲线，横波时差曲线和密度曲线：t
pmi
＝165v
smi
+293v
shmi
，其中，t
pmi
为所述砂泥岩石混合基质的纵波时差曲线中第i个值，v
smi
根据以下表达式确定：v
smi
＝v
si
/(1-φ
i
)，其中，φ
i
为所述孔隙度测井曲线中第i个值，v
si
根据以下关系式确定：v
si
+v
shi
+φ
i
＝1，其中，v
shi
为所述泥质含量测井曲线中第i个值，并且，v
shmi
为所述砂泥岩石混合基质的泥质含量测井曲线中第i个值，根据以下表达式确定：v
shmi
＝v
shi
/(1-φ
i
)；t
smi
＝244v
smi
+609v
shmi
，其中，t
smi
为所述砂泥岩石混合基质的横波时差曲线中第i个值；ρ
mi
＝2.65v
smi
+2.6v
shmi
，其中，ρ
mi
为所述砂泥岩石混合基质的密度曲线中第i个值。14.根据权利要求11所述的方法，其中，所述横波测井曲线和纵波测井曲线的转换关系式x为：实例中为v
s_sati
＝0.8114v
p_sati-0.8957，其中，v
s_sati
为所述横波测井曲线中第i个值，v
p_sati
为所述纵波测井曲线中第i个值。15.根据权利要求14所述的方法，其中，还包括以下步骤：利用带通滤波器对所述测量纵波测井曲线进行低频滤波处理，将所述带通滤波器的截止频率设置为0，0，35，40，滤波得到所述测量纵波的低频测井曲线；
利用所述转换关系式x，将所述测量纵波的低频测井曲线转换为所述预测横波测井曲线的低频曲线；利用所述带通滤波器对所述计算横波测井曲线进行高频滤波处理，将所述带通滤波器的截止频率设置为40，40，9999，9999，滤波得到所述计算横波的高频测井曲线；根据以下表达式确定所述预测横波测井曲线：v
s_predictedi
＝v
s_lowi
+v
s_highi
，其中，v
s_predictedi
为所述预测横波测井曲线中第i个值，v
s_lowi
为所述预测横波测井曲线的低频曲线中第i个值，v
s_highi
为所述计算横波的高频测井曲线中第i个值。

技术总结
本申请的实施例涉及地震勘探领域，具体涉及确定砂岩型铀矿勘查区的预测横波测井曲线的方法，包括如下步骤：测量砂岩型铀矿勘查区，确定测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线；根据测量纵波测井曲线、伽马测井曲线以及电阻率测井曲线，确定计算横波测井曲线；对计算横波测井曲线进行高频滤波处理，确定计算横波的高频测井曲线；对测量纵波测井曲线进行低频滤波处理，确定测量纵波的低频测井曲线；将测量纵波的低频测井曲线转换成预测横波测井曲线的低频曲线；根据预测横波测井曲线的低频曲线和计算横波的高频测井曲线，确定预测横波测井曲线。通过本申请实施例中的方法，可以准确地确定预测横波测井曲线。可以准确地确定预测横波测井曲线。可以准确地确定预测横波测井曲线。


技术研发人员：吴曲波 黄昱丞 曹成寅 李子伟 乔宝平 潘自强
受保护的技术使用者：核工业北京地质研究院
技术研发日：2023.05.31
技术公布日：2023/9/5