氦气成因资源评价方法及装置与流程

1.本技术涉及氦气资源评价技术领域，具体涉及一种氦气成因资源评价方法、一种氦气成因资源评价装置、一种机器可读存储介质及一种处理器。


背景技术：

2.氦气被称为“黄金气体”。发明人发现，由于氦源岩的厚度、运聚系数以及幔源岩的难以确定，目前国内外尚没有成熟的生氦法来开展氦气资源评价方法。发明人发现，氦气百分含量法虽然计算准确，但其依赖于氦气数据点的数量和质量、天然气储量的准确性。同时在成因法计算资源量时候，由于氦源岩的厚度、运聚系数以及幔源岩的难以确定，资源量的计算误差很大。发明人还发现，尽管氦气与天然气聚集在同一圈闭中，由于生氦和生烃机制完全不同，天然气与氦气的运聚系数不存在相关性，氦气藏的运聚有其自身的规律。
[0003][0004]
鉴于上述问题，本技术提出了克服上述问题或者至少部分地解决上述问题的一种氦气成因资源评价方法。


技术实现要素：

[0005]
本技术实施例的目的是提供一种氦气成因资源评价方法、一种氦气成因资源评价装置、一种机器可读存储介质及一种处理器。所述方法能够实现氦气资源的有效评价，对于开展氦气成藏和富集规律研究具有重要指导作用。
[0006]
为了实现上述目的，本技术第一方面提供一种氦气成因资源评价方法，所述方法包括：
[0007]
获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
[0008]
确定目标区域的氦源岩的体积；
[0009]
基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
[0010]
基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
[0011]
基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
[0012]
基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
[0013]
确定目标区域的氦源岩中氦气的运聚系数；
[0014]
基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
[0015]
其中，所述确定目标区域的氦源岩的体积，包括：
[0016]
基于预设温度确定目标区域的氦源岩的底界深度；
[0017]
确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
[0018]
基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
[0019]
基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
[0020]
在本技术实施例中，所述确定目标区域的氦源岩的底界深度，包括：基于预设温度确定目标区域的氦源岩的底界深度。
[0021]
在本技术实施例中，所述基于预设温度确定目标区域的氦源岩的底界深度，包括：
[0022]
将预设温度所对应的地层深度确定为氦源岩的底界深度，所述氦源岩包括：花岗岩和火山岩。
[0023]
在本技术实施例中，所述预设温度的温度范围为：220℃～240℃。
[0024]
在本技术实施例中，所述确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，包括：
[0025]
基于地震资料和钻井资料确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，所述氦源岩包括：煤层、泥页岩。
[0026]
在本技术实施例中，基于地震资料和钻井资料确定氦源岩的分布面积，包括：
[0027]
对于厚度不稳定或者分布不连续的氦源岩，基于地震资料和钻井资料，通过小面元积分法求取氦源岩的分布面积。
[0028]
在本技术实施例中，通过(1)式，将预设温度所对应的地层深度确定为氦源岩的底界深度：
[0029]h底
＝(t
预设温度-t
地表温度
)/k+50
ꢀꢀꢀ
(1)；
[0030]
式中，h
底
表示氦源岩的底界深度，单位为m；
[0031]
t
地表温度
表示地表温度，单位为℃；
[0032]
k表示地温梯度，单位为℃/100m。
[0033]
在本技术实施例中，通过(2)式，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度：
[0034]
h＝h
底-h
顶
ꢀꢀꢀ
(2)；
[0035]
式中，h表示氦源岩的厚度，单位为m；
[0036]h底
表示氦源岩的底界深度，单位为m；
[0037]h顶
表示氦源岩的顶面深度，单位为m。
[0038]
在本技术实施例中，通过(3)式，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积：
[0039]
v＝h
×sꢀꢀꢀ
(3)；
[0040]
式中，v表示氦源岩的体积，单位为m3；
[0041]
h表示氦源岩的厚度，单位为m；
[0042]
s表示氦源岩分布的面积，单位为m2。
[0043]
在本技术实施例中，通过(4)式，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量：
[0044]
p＝v
×
ρ
源岩
×
x
ꢀꢀꢀ
(4)；
[0045]
式中，p表示氦源岩中未衰变的母体同位素的含量，单位为mol；
[0046]
v表示氦源岩的体积，单位为m3；
[0047]
ρ
源岩
表示氦源岩的密度，单位为t/m3；
[0048]
x表示氦源岩中放射性元素的含量，单位为mol/g。
[0049]
在本技术实施例中，通过(5)式，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量：
[0050]
d＝(e
λt-1)
×v×
ρ
源岩
×
x
ꢀꢀꢀ
(5)；
[0051]
式中，d表示氦源岩中已衰变的子体同位素的含量，单位为mol；
[0052]
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
[0053]
t表示氦源岩的绝对地质年龄，单位为百万年；
[0054]
v表示氦源岩的体积，单位为m3；
[0055]
ρ
源岩
表示氦源岩的密度，单位为t/m3；
[0056]
x表示氦源岩中放射性元素的含量，单位为mol/g。
[0057]
在本技术实施例中，氦源岩中的放射性元素包括：238u和232th；
[0058]
针对238u衰变，通过(6)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
[0059]
n＝d
×
8/238
ꢀꢀꢀ
(6)；
[0060]
针对232th衰变，通过(7)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
[0061]
n＝d
×
6/232
ꢀꢀꢀ
(7)；
[0062]
式中，n表示氦源岩中氦气的物质的量，d表示氦源岩中已衰变的子体同位素的含量。
[0063]
在本技术实施例中，通过(8)式，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量：
[0064]vhe生
＝n
×
m/ρ
he
ꢀꢀꢀ
(8)；
[0065]
式中，v
he生
表示生氦总量，即生氦资源量，单位为m3；
[0066]
n表示氦源岩中氦气的物质的量；
[0067]
m表示生氦总量的摩尔质量；
[0068]
ρ
he
表示氦气的密度，取值0.1786g/l。
[0069]
在本技术实施例中，通过(9)式，确定目标区域的氦源岩中氦气的运聚系数：
[0070]
φ＝φ1
×
φ2
ꢀꢀꢀ
(9)；
[0071]
式中，φ1为氦源岩排氦率，单位为％；φ2为水溶氦脱气率，单位为％。在本技术实施例中，通过(10)式，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量：
[0072]vhe聚
＝φ
×vhe生
ꢀꢀꢀ
(10)；
[0073]
式中，v
he聚
表示氦气聚集量，即氦气聚集资源量，单位为m3；
[0074]
φ表示运聚系数，单位为％；
[0075]vhe生
表示生氦总量，单位为m3。
[0076]
在本技术实施例中，所述氦气包括3he和4he，所述方法还包括：
[0077]
根据二元复合公式，通过3he/4he的比值，确定幔源氦气和壳源氦气的资源量。
[0078]
本技术第二方面提供一种氦气成因资源评价装置，所述装置包括：
[0079]
获取模块，用于获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
[0080]
第一确定模块，用于确定目标区域的氦源岩的体积；
[0081]
第二确定模块，用于基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
[0082]
第三确定模块，用于基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
[0083]
第四确定模块，用于基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
[0084]
第五确定模块，用于基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
[0085]
第六确定模块，用于确定目标区域的氦源岩中氦气的运聚系数；
[0086]
第七确定模块，用于基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
[0087]
其中，所述第一确定模块，包括：
[0088]
底界深度确定模块，用于基于预设温度确定目标区域的氦源岩的底界深度；
[0089]
顶面深度以及分布面积确定模块，用于确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
[0090]
厚度确定模块，用于基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
[0091]
体积确定模块，用于基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
[0092]
本技术第三方面提供一种处理器，被配置成执行上述的氦气成因资源评价方法。
[0093]
本技术第四方面提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令在被处理器执行时使得所述处理器被配置成执行上述的氦气成因资源评价方法。
[0094]
与现有技术相比，本发明的上述技术方案具有如下有益效果：
[0095]
本技术提出的一种氦气成因资源评价方法、装置、存储介质及处理器，所述方法通过氦源岩中已衰变的子体同位素的含量确定出氦源岩中氦气的物质的量，然后基于氦源岩中氦气的物质的量确定出氦源岩的生氦总量，最后基于氦源岩中氦气的运聚系数，确定氦气聚集资源量。本技术通过上述方法能够实现氦气资源的有效评价，对于开展氦气成藏和富集规律研究具有重要指导作用。
[0096]
本技术实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
[0097]
附图是用来提供对本技术实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本技术实施例，但并不构成对本技术实施例的限制。在附图中：
[0098]
图1示意性示出了根据本技术实施例的氦气成因资源评价方法的应用环境示意图；
[0099]
图2示意性示出了根据本技术实施例的氦气成因资源评价方法的整体流程示意
图；
[0100]
图3示意性示出了根据本技术实施例的氦气成因资源评价方法中，计算目标区域的氦源岩的体积的流程示意图；
[0101]
图4示意性示出了根据本技术实施例的he、n2和ch4亨利常数随温度的变化曲线图；
[0102]
图5示意性示出了根据本技术实施例的氦源岩有效厚度示意图；
[0103]
图6示意性示出了根据本技术实施例的氦气成藏模式及过程示意图；
[0104]
图7示意性示出了根据本技术实施例的氦气成因资源评价装置的结构框图；
[0105]
图8示意性示出了根据本技术实施例的计算机设备的内部结构图。
具体实施方式
[0106]
为使本技术实施例的目的、技术方案和优点更加清楚，下面将结合本技术实施例中的附图，对本技术实施例中的技术方案进行清楚、完整地描述，应当理解的是，此处所描述的具体实施方式仅用于说明和解释本技术实施例，并不用于限制本技术实施例。基于本技术中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本技术保护的范围。
[0107]
需要说明，若本技术实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本技术要求的保护范围之内。
[0108]
本技术提供的氦气成因资源评价方法，可以应用于如图1所示的应用环境中。其中，终端102通过网络与服务器104通过网络进行通信。其中，终端102可以但不限于是各种个人计算机、笔记本电脑、智能手机、平板电脑和便携式可穿戴设备，服务器104可以用独立的服务器或者是多个服务器组成的服务器集群来实现。
[0109]
图2示意性示出了根据本技术实施例的氦气成因资源评价方法的流程示意图。如图2所示，在本技术一实施例中，提供了一种氦气成因资源评价方法，本实施例主要以该方法应用于上述图1中的终端102(或服务器104)
[0110]
来举例说明，包括以下步骤：
[0111]
步骤110，获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数。
[0112]
在本实施例中，目标区域指的是氦气资源研究区域。在本实施例中，氦源岩包括花岗岩、火山岩和沉积岩等。花岗岩或者火山岩的年龄，可以根据测定岩石中锆石u-pb的同位素年龄确定。沉积岩的年龄主要通过测定岩石中古生物年代确定。氦源岩的密度一般是常数，如花岗岩的密度为2.3g/cm3。氦源岩中放射性元素包括：铀(
235
u、
238
u)和钍(
232
th)。可以通过采集不同层次、不同构造部位具有代表性的氦源岩样品，来测定氦源岩中铀或者钍等放射性元素的含量，可以多次测量铀或者钍的含量，同一样品误差小于5％，求取平均值。
235
u的衰变常数为9.8485
×
10-10
/a，
238
u的衰变常数为1.55125
ꢀ×
10-10
/a，
232
th的衰变常数为4.9745
×
10-11
/a。
[0113]
步骤120，确定目标区域的氦源岩的体积。
[0114]
图3示意性示出了根据本技术实施例的氦气成因资源评价方法中，计算目标区域的氦源岩的体积的流程示意图。如图3所示，在本实施例中，所述目标区域的氦源岩的体积通过步骤121～步骤124确定：
[0115]
步骤121，基于预设温度确定目标区域的氦源岩的底界深度。
[0116]
示例性地，计算氦源岩的体积，首先，需要计算氦源岩的底界深度。
[0117]
若氦源岩为花岗岩或者火山岩，则所述氦源岩的底界深度可以根据220℃～240℃所对应的地层深度确定。
[0118]
图4示意性示出了根据本技术实施例的he、n2和ch4亨利常数随温度的变化曲线。如图4所示，依据亨利曲线，氦气(he)在180℃时与甲烷(ch4)、氮气(n2)溶解度、分压基本趋于一致，再进一步释放氦气的潜力基本消失， 240℃时氦源岩释放潜力趋于零，因此，本实施例将研究区220℃～240℃所对应的地层深度确定为花岗岩和火山岩等氦源岩的下限深度，即氦源岩的底界深度。
[0119]
图5示意性示出了根据本技术实施例的氦源岩有效厚度示意图。
[0120]
示例性地，若氦源岩为花岗岩或者火山岩，则氦源岩的底界深度可以通过(1)式计算：
[0121]h底
＝(t
预设温度-t
地表温度
)/k+50
ꢀꢀꢀ
(1)；
[0122]
式中，h
底
表示氦源岩的底界深度，单位为m；
[0123]
t
预设温度
表示预设温度，t
预设温度
的取值范围为220℃～240℃。
[0124]
t
地表温度
表示地表温度，单位为℃；
[0125]
k表示地温梯度，单位为℃/100m。
[0126]
若氦源岩为沉积岩，例如，煤层或者泥页岩，则所述氦源岩的底界深度可以根据地震资料和钻井资料确定。
[0127]
步骤122，确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积。
[0128]
在本实施例中，可以通过地震资料和钻井资料确定氦源岩的顶面深度以及氦源岩的分布面积。对于厚度不稳定或者分布不连续的氦源岩，可以通过小面元积分法求取氦源岩的分布面积。
[0129]
步骤123，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度。
[0130]
示例性地，氦源岩的厚度可以通过(2)式计算：
[0131]
h＝h
底-h
顶
ꢀꢀꢀ
(2)；
[0132]
式中，h表示氦源岩的厚度，单位为m；
[0133]h底
表示氦源岩的底界深度，单位为m；
[0134]h顶
表示氦源岩的顶面深度，单位为m。
[0135]
步骤124，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
[0136]
示例性地，氦源岩的体积可以通过(3)式计算：
[0137]
v＝h
×sꢀꢀꢀ
(3)；
[0138]
式中，v表示氦源岩的体积，单位为m3；
[0139]
h表示氦源岩的厚度，单位为m；
[0140]
s表示氦源岩的分布面积，单位为m2。
[0141]
本实施例提出基于亨利定律的生氦法，解决了氦源岩厚度参数确定的关键难点。常规方法氦源岩的下限采用主观推断法，根据氦源岩的露头厚度来推断地下氦源岩厚度，但地下氦源岩厚度变化大，导致氦源岩的体积计算不准确。本技术以氦气成藏理论为指导，将氦气“死亡线”220℃-240℃作为氦原岩深度的下限温度，将该温度对应的深度作为氦源岩的深度下限，由此确定氦源岩厚度，有效提高了氦源岩体积和生氦资源量准确性，误差下降10％。本实施例有效解决了现有技术中存在的氦源岩的厚度难以确定的问题。
[0142]
步骤130，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量。
[0143]
在本实施例中，可以通过(4)计算氦源岩中未衰变的母体同位素的含量：
[0144]
p＝v
×
ρ
源岩
×
x
ꢀꢀꢀ
(4)；
[0145]
式中，p表示氦源岩中未衰变的母体同位素的含量，单位为mol；
[0146]
v表示氦源岩的体积，单位为m3；
[0147]
ρ
源岩
表示氦源岩的密度，单位为t/m3；
[0148]
x表示氦源岩中放射性元素的含量，单位为mol/g。
[0149]
步骤140，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量。
[0150]
在本实施例中，可以通过(5)计算氦源岩中已衰变的子体同位素的含量：
[0151]
d＝(e
λt-1)
×v×
ρ
源岩
×
x
ꢀꢀꢀ
(5)；
[0152]
式中，d表示氦源岩中已衰变的子体同位素的含量，单位为mol；
[0153]
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
[0154]
t表示氦源岩的绝对地质年龄，单位为百万年；
[0155]
v表示氦源岩的体积，单位为m3；
[0156]
ρ
源岩
表示氦源岩的密度，单位为t/m3；
[0157]
x表示氦源岩中放射性元素的含量，单位为mol/g。
[0158]
步骤150，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量。
[0159]
在本实施例中，氦源岩中的放射性元素包括：238u和232th；
[0160]
针对238u衰变，通过(6)式，计算氦源岩中氦气的物质的量：
[0161]
n＝d
×
8/238
ꢀꢀꢀ
(6)；
[0162]
针对232th衰变，通过(7)式，计算氦源岩中氦气的物质的量：
[0163]
n＝d
×
6/232
ꢀꢀꢀ
(7)；
[0164]
式中，n表示氦源岩中氦气的物质的量，d表示氦源岩中已衰变的子体同位素的含量。
[0165]
步骤160，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量。
[0166]
在本实施例中，通过(8)式，计算氦源岩的生氦总量：
[0167]vhe生
＝n
×
m/ρ
he
ꢀꢀꢀ
(8)；
[0168]
式中，v
he生
表示生氦总量，即生氦资源量，单位为m3；
[0169]
n表示氦源岩中氦气的物质的量；
[0170]
m表示生氦总量的摩尔质量；
[0171]
ρ
he
表示氦气的密度，取值0.1786g/l。
[0172]
示例性地，针对238u衰变，通过(8.1)式，计算氦源岩的生氦总量：
[0173]vhe生
＝n
×
m/ρ
he
[0174]
＝d
×m×
8/238/ρ
he
[0175]
＝(e
λt-1)
×v×
ρ
源岩
×
x
×m×
8/238/ρ
he
ꢀꢀꢀ
(8.1)；
[0176]
示例性地，针对232th衰变，通过(8.2)式，计算氦源岩的生氦总量：
[0177]vhe生
＝n
×
m/ρ
he
[0178]
＝d
×m×
6/232/ρ
he
[0179]
＝(e
λt-1)
×v×
ρ
源岩
×
x
×m×
6/232/ρ
he
ꢀꢀꢀ
(8.2)；
[0180]
式(9)-(10)中，v
he生
表示研究区生氦资源量，即氦源岩的生氦总量，单位为m3；
[0181]
λ表示研究区铀或者钍同位素的衰变常数，无量纲；
[0182]
t表示研究区氦源岩的绝对地质年龄，单位为百万年；
[0183]
ρ
源岩
表示研究区氦源岩的密度，单位为t/m3；例如，花岗岩的密度取值2.5t/m3；
[0184]
ρ
he
表示研究区氦气的密度，取值0.1786g/l；
[0185]
v表示研究区氦源岩的体积，单位为m3；
[0186]
x表示研究区氦源岩中铀或者钍的含量，单位为g/mol。
[0187]
步骤170，确定目标区域的氦源岩中氦气的运聚系数。
[0188]
图6示意性示出了根据本技术实施例的氦气成藏模式及过程示意图。如图6所示，本实施例根据研究区域构造活动期次和排氦特征，确定不同成藏模式下的氦源岩排氦率和水溶氦脱气率，两者相乘即可得到运聚系数。其中，不同成藏模式是指，不同来源的氦气在不同的圈闭中的聚集形式：按照气体成分，分为含氦氮气藏、含氦co2气藏、含氦天然气藏等；按圈闭类型分为背斜、断层以及岩性；按照来源分为幔源、壳源或混源等等。
[0189]
已经排氦的氦源岩与全部氦源岩体积之比即为氦源岩排氦率。通过测定刻度区已知水溶氦含量和气藏含氦量，计算已经聚集的氦气量与水溶氦的体积，两者之比即为水溶氦脱气率。
[0190]
在本实施例中，通过(9)式，计算目标区域的氦源岩中氦气的运聚系数：
[0191]
φ＝φ1
×
φ2
ꢀꢀꢀ
(9)；
[0192]
式中，φ1为氦源岩排氦率，单位为％；φ2为水溶氦脱气率，单位为％。
[0193]
其中，氦源岩排氦率φ1通过已经排氦的氦源岩与全部氦源岩体积之比来计算。测定刻度区已知水溶氦含量和气藏含氦量，计算已经聚集的氦气量与水溶氦的体积，两者之比即为水溶氦脱气率φ2。
[0194]
在不同的成藏模式下，氦气的初次运移和二次运移系数千差万别。已有的方法没有考虑不同成藏模式的差异，也没有给出具体的计算公式。本实施例将运聚系数与烃类成藏时间相匹配，解决运聚系数难以确定的问题。
[0195]
步骤180，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
[0196]
本实施例根据研究区花岗岩、火成岩、沉积岩等所有生氦总量，与运聚系数相乘，即可得到研究区氦气地质资源量。
[0197]
在本实施例中，通过(10)式，计算氦气聚集资源量：
[0198]vhe聚
＝φ
×vhe生
ꢀꢀꢀ
(10)；
[0199]
式中，v
he聚
表示氦气聚集量，即氦气聚集资源量，单位为m3；
[0200]
φ表示运聚系数，单位为％；
[0201]vhe生
表示生氦总量，单位为m3。
[0202]
在本实施例中，所述氦气包括3he和4he，所述方法还包括：
[0203]
确定壳源氦和幔源氦中3he/4he的标准值，标准值如下：
[0204]
壳源氦取值标准：3he/4he＝2.0
×
10-8
；
[0205]
幔源氦取值标准：3he/4he＝1.1
×
10-8
；
[0206]
壳幔混源氦取值标准：1.1
×
10-8
＜3he/4he＜2.0
×
10-8
。
[0207]
在本实施例中，所述氦气包括3he和4he，所述方法还包括：根据二元复合公式，通过3he/4he比值计算幔源氦气和壳源氦气的资源量，壳源氦的比例计算公式如下：
[0208][0209]
幔源氦资源量比例的计算公式如下：
[0210]i4he
＝1-i
3he
ꢀꢀꢀ
(12)；
[0211]
式中，3he表示氦气中幔源氦气所占的比例，单位为％l；
[0212]4he表示氦气中壳源氦气所占的比例，单位为％l；
[0213]
r表示实测氦气值体积，单位为m3；
[0214]
rc表示壳源氦气体积，单位为m3；
[0215]rm
表示幔源氦气体积，单位为m3。
[0216]
图2为一个实施例中氦气成因资源评价方法的流程示意图。应该理解的是，虽然图2的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
[0217]
图7示意性示出了根据本技术实施例的氦气成因资源评价装置的结构框图。
[0218]
在一个实施例中，如图7所示，提供了一种氦气成因资源评价装置200，包括获取模块210、第一确定模块220、第二确定模块230、第三确定模块 240、第四确定模块250、第五确定模块260、第六确定模块270以及第七确定模块280，其中：
[0219]
获取模块210，用于获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
[0220]
第一确定模块220，用于确定目标区域的氦源岩的体积；
[0221]
第二确定模块230，用于基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
[0222]
第三确定模块240，用于基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
[0223]
第四确定模块250，用于基于所述氦源岩中已衰变的子体同位素的含量确定氦源
岩中氦气的物质的量；
[0224]
第五确定模块260，用于基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
[0225]
第六确定模块270，用于确定目标区域的氦源岩中氦气的运聚系数；
[0226]
第七确定模块280，用于基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
[0227]
其中，所述第一确定模块220，包括：
[0228]
底界深度确定模块，用于基于预设温度确定目标区域的氦源岩的底界深度；
[0229]
顶面深度以及分布面积确定模块，用于确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
[0230]
厚度确定模块，用于基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
[0231]
体积确定模块，用于基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
[0232]
所述氦气成因资源评价装置包括处理器和存储器，上述获取模块210、第一确定模块220、第二确定模块230、第三确定模块240、第四确定模块 250、第五确定模块260、第六确定模块270以及第七确定模块280等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序模块中实现相应的功能。
[0233]
处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来实现氦气成因资源评价方法。
[0234]
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器 (ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flash ram)，存储器包括至少一个存储芯片。
[0235]
本技术实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现上述氦气成因资源评价方法。
[0236]
图8示意性示出了根据本技术实施例的计算机设备的内部结构图。
[0237]
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器a01、网络接口a02、显示屏a04、输入装置a05和存储器(图中未示出)。其中，该计算机设备的处理器a01用于提供计算和控制能力。该计算机设备的存储器包括内存储器a03和非易失性存储介质a06。该非易失性存储介质 a06存储有操作系统b01和计算机程序b02。该内存储器a03为非易失性存储介质a06中的操作系统b01和计算机程序b02的运行提供环境。该计算机设备的网络接口a02用于与外部的终端通过网络连接通信。该计算机程序被处理器a01执行时以实现一种氦气成因资源评价方法。该计算机设备的显示屏a04可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置 a05可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
[0238]
本领域技术人员可以理解，图8中示出的结构，仅仅是与本技术方案相关的部分结构的框图，并不构成对本技术方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
[0239]
在一个实施例中，本技术提供的氦气成因资源评价装置可以实现为一种计算机程序的形式，计算机程序可在如图8所示的计算机设备上运行。计算机设备的存储器中可存储组成该氦气成因资源评价装置的各个程序模块，比如，图7所示的获取模块210、第一确定模块220、第二确定模块230、第三确定模块240、第四确定模块250、第五确定模块260、第六确定模块270 以及第七确定模块280。各个程序模块构成的计算机程序使得处理器执行本说明书中描述的本技术各个实施例的氦气成因资源评价方法中的步骤。
[0240]
图8所示的计算机设备可以通过如图7所示的氦气成因资源评价装置中的获取模块210执行步骤110，计算机设备可通过第一确定模块220执行步骤120，计算机设备可通过第二确定模块230执行步骤130，计算机设备可通过第三确定模块240执行步骤140，计算机设备可通过第四确定模块250 执行步骤150，计算机设备可通过第五确定模块260执行步骤160，计算机设备可通过第六确定模块270执行步骤170，计算机设备可通过第七确定模块280执行步骤180。
[0241]
本技术实施例提供了一种设备，设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序，处理器执行程序时实现以下步骤：
[0242]
步骤110，获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
[0243]
步骤120，确定目标区域的氦源岩的体积；
[0244]
步骤130，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
[0245]
步骤140，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
[0246]
步骤150，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
[0247]
步骤160，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
[0248]
步骤170，确定目标区域的氦源岩中氦气的运聚系数；
[0249]
步骤180，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
[0250]
在一个实施例中，所述步骤120，包括：
[0251]
步骤121，基于预设温度确定目标区域的氦源岩的底界深度；
[0252]
步骤122，确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
[0253]
步骤123，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
[0254]
步骤124，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
[0255]
在一个实施例中，所述确定目标区域的氦源岩的底界深度，包括：
[0256]
基于预设温度确定目标区域的氦源岩的底界深度。
[0257]
在一个实施例中，所述基于预设温度确定目标区域的氦源岩的底界深度，包括：
[0258]
将预设温度所对应的地层深度确定为氦源岩的底界深度，所述氦源岩包括：花岗岩和火山岩。
[0259]
在一个实施例中，所述预设温度的温度范围为：220℃～240℃。
[0260]
在一个实施例中，所述确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，包括：
[0261]
基于地震资料和钻井资料确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，所述氦源岩包括：煤层、泥页岩。
[0262]
在一个实施例中，所述基于地震资料和钻井资料确定氦源岩的分布面积，包括：
[0263]
对于厚度不稳定或者分布不连续的氦源岩，基于地震资料和钻井资料，通过小面元积分法求取氦源岩的分布面积。
[0264]
在一个实施例中，通过(1)式，将预设温度所对应的地层深度确定为氦源岩的底界深度：
[0265]h底
＝(t
预设温度-t
地表温度
)/k+50
ꢀꢀꢀ
(1)；
[0266]
式中，h
底
表示氦源岩的底界深度，单位为m；
[0267]
t
地表温度
表示地表温度，单位为℃；
[0268]
k表示地温梯度，单位为℃/100m。
[0269]
在一个实施例中，通过(2)式，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度：
[0270]
h＝h
底-h
顶
ꢀꢀꢀ
(2)；
[0271]
式中，h表示氦源岩的厚度，单位为m；
[0272]h底
表示氦源岩的底界深度，单位为m；
[0273]h顶
表示氦源岩的顶面深度，单位为m。
[0274]
在一个实施例中，通过(3)式，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积：
[0275]
v＝h
×sꢀꢀꢀ
(3)；
[0276]
式中，v表示氦源岩的体积，单位为m3；
[0277]
h表示氦源岩的厚度，单位为m；
[0278]
s表示氦源岩的分布面积，单位为m2。
[0279]
在一个实施例中，通过(4)式，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量：
[0280]
p＝v
×
ρ
源岩
×
x
ꢀꢀꢀ
(4)；
[0281]
式中，p表示氦源岩中未衰变的母体同位素的含量，单位为mol；
[0282]
v表示氦源岩的体积，单位为m3；
[0283]
ρ
源岩
表示氦源岩的密度，单位为t/m3；
[0284]
x表示氦源岩中放射性元素的含量，单位为mol/g。
[0285]
在一个实施例中，通过(5)式，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量：
[0286]
d＝(e
λt-1)
×v×
ρ
源岩
×
x
ꢀꢀꢀ
(5)；
[0287]
式中，d表示氦源岩中已衰变的子体同位素的含量，单位为mol；
[0288]
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
[0289]
t表示氦源岩的绝对地质年龄，单位为百万年；
[0290]
v表示氦源岩的体积，单位为m3；
[0291]
ρ
源岩
表示氦源岩的密度，单位为t/m3；
[0292]
x表示氦源岩中放射性元素的含量，单位为mol/g。
[0293]
在一个实施例中，氦源岩中的放射性元素包括：238u和232th；
[0294]
针对238u衰变，通过(6)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
[0295]
n＝d
×
8/238
ꢀꢀꢀ
(6)；
[0296]
针对232th衰变，通过(7)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
[0297]
n＝d
×
6/232
ꢀꢀꢀ
(7)；
[0298]
式中，n表示氦源岩中氦气的物质的量，d表示氦源岩中已衰变的子体同位素的含量。
[0299]
在一个实施例中，通过(8)式，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量：
[0300]vhe生
＝n
×
m/ρ
he
ꢀꢀꢀ
(8)；
[0301]
式中，v
he生
表示生氦总量，即生氦资源量，单位为m3；
[0302]
n表示氦源岩中氦气的物质的量；
[0303]
m表示生氦总量的摩尔质量；
[0304]
ρ
he
表示氦气的密度，取值0.1786g/l。
[0305]
在一个实施例中，通过(9)式，确定目标区域的氦源岩中氦气的运聚系数：
[0306]
φ＝φ1
×
φ2
ꢀꢀꢀ
(9)；
[0307]
式中，φ1为氦源岩排氦率，单位为％；φ2为水溶氦脱气率，单位为％。
[0308]
在一个实施例中，通过(10)式，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量：
[0309]vhe聚
＝φ
×vhe生
ꢀꢀꢀ
(10)；
[0310]
式中，v
he聚
表示氦气聚集量，即氦气聚集资源量，单位为m3；
[0311]
φ表示运聚系数，单位为％；
[0312]vhe生
表示生氦总量，单位为m3。
[0313]
在一个实施例中，所述氦气包括3he和4he，所述方法还包括：根据二元复合公式，通过3he/4he的比值，确定幔源氦气和壳源氦气的资源量。
[0314]
本领域内的技术人员应明白，本技术的实施例可提供为方法、系统、或计算机程序产品。因此，本技术可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本技术可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、cd-rom、光学存储器等)上实施的计算机程序产品的形式。
[0315]
本技术是参照根据本技术实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/ 或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装
置。
[0316]
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
[0317]
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
[0318]
在一个典型的配置中，计算设备包括一个或多个处理器(cpu)、输入/ 输出接口、网络接口和内存。
[0319]
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器 (ram)和/或非易失性内存等形式，如只读存储器(rom)或闪存(flash ram)。存储器是计算机可读介质的示例。
[0320]
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体，可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(pram)、静态随机存取存储器(sram)、动态随机存取存储器(dram)、其他类型的随机存取存储器(ram)、只读存储器(rom)、电可擦除可编程只读存储器(eeprom)、快闪记忆体或其他内存技术、只读光盘只读存储器 (cd-rom)、数字多功能光盘(dvd)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
[0321]
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个
……”
限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
[0322]
以上仅为本技术的实施例而已，并不用于限制本技术。对于本领域技术人员来说，本技术可以有各种更改和变化。凡在本技术的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本技术的权利要求范围之内。