CN115685377A - 氦气成因资源评价方法及装置 - Google Patents
氦气成因资源评价方法及装置 Download PDFInfo
- Publication number
- CN115685377A CN115685377A CN202211271154.5A CN202211271154A CN115685377A CN 115685377 A CN115685377 A CN 115685377A CN 202211271154 A CN202211271154 A CN 202211271154A CN 115685377 A CN115685377 A CN 115685377A
- Authority
- CN
- China
- Prior art keywords
- helium
- source rock
- helium source
- rock
- determining
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 229910052734 helium Inorganic materials 0.000 title claims abstract description 679
- 239000001307 helium Substances 0.000 title claims abstract description 677
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 title claims abstract description 676
- 239000007789 gas Substances 0.000 title claims abstract description 68
- 238000011156 evaluation Methods 0.000 title claims abstract description 38
- 239000011435 rock Substances 0.000 claims abstract description 495
- 238000000034 method Methods 0.000 claims abstract description 49
- 238000003860 storage Methods 0.000 claims abstract description 17
- 238000004220 aggregation Methods 0.000 claims abstract description 13
- 238000004519 manufacturing process Methods 0.000 claims abstract description 13
- 230000005012 migration Effects 0.000 claims abstract description 13
- 238000013508 migration Methods 0.000 claims abstract description 13
- 230000002776 aggregation Effects 0.000 claims abstract description 12
- 230000002285 radioactive effect Effects 0.000 claims description 50
- 238000009826 distribution Methods 0.000 claims description 42
- 239000000126 substance Substances 0.000 claims description 20
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 claims description 14
- 239000010438 granite Substances 0.000 claims description 11
- 238000005553 drilling Methods 0.000 claims description 9
- 238000007872 degassing Methods 0.000 claims description 7
- 239000003245 coal Substances 0.000 claims description 4
- 150000001875 compounds Chemical class 0.000 claims description 4
- 230000010354 integration Effects 0.000 claims description 4
- NAWXUBYGYWOOIX-SFHVURJKSA-N (2s)-2-[[4-[2-(2,4-diaminoquinazolin-6-yl)ethyl]benzoyl]amino]-4-methylidenepentanedioic acid Chemical compound C1=CC2=NC(N)=NC(N)=C2C=C1CCC1=CC=C(C(=O)N[C@@H](CC(=C)C(O)=O)C(O)=O)C=C1 NAWXUBYGYWOOIX-SFHVURJKSA-N 0.000 claims description 2
- 238000009825 accumulation Methods 0.000 abstract description 10
- 238000011160 research Methods 0.000 abstract description 10
- 238000010586 diagram Methods 0.000 description 17
- 238000004590 computer program Methods 0.000 description 12
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 description 10
- 241000894007 species Species 0.000 description 10
- 230000006870 function Effects 0.000 description 7
- 238000004364 calculation method Methods 0.000 description 6
- 230000008569 process Effects 0.000 description 6
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 6
- 229910052776 Thorium Inorganic materials 0.000 description 5
- 229910052770 Uranium Inorganic materials 0.000 description 5
- 239000003345 natural gas Substances 0.000 description 4
- 238000012545 processing Methods 0.000 description 4
- IJGRMHOSHXDMSA-UHFFFAOYSA-N Atomic nitrogen Chemical compound N#N IJGRMHOSHXDMSA-UHFFFAOYSA-N 0.000 description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 2
- 230000015572 biosynthetic process Effects 0.000 description 2
- 230000008859 change Effects 0.000 description 2
- 238000005516 engineering process Methods 0.000 description 2
- 229930195733 hydrocarbon Natural products 0.000 description 2
- 150000002430 hydrocarbons Chemical class 0.000 description 2
- 238000011835 investigation Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 229910052757 nitrogen Inorganic materials 0.000 description 2
- 230000003287 optical effect Effects 0.000 description 2
- -1 poly (phenylene vinylene) Polymers 0.000 description 2
- 229920000553 poly(phenylenevinylene) Polymers 0.000 description 2
- 238000006116 polymerization reaction Methods 0.000 description 2
- 241000203069 Archaea Species 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000007796 conventional method Methods 0.000 description 1
- 229910001873 dinitrogen Inorganic materials 0.000 description 1
- 238000007599 discharging Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 239000010931 gold Substances 0.000 description 1
- 229910052737 gold Inorganic materials 0.000 description 1
- 150000002371 helium Chemical class 0.000 description 1
- 230000006872 improvement Effects 0.000 description 1
- 230000000155 isotopic effect Effects 0.000 description 1
- 230000000670 limiting effect Effects 0.000 description 1
- 239000004973 liquid crystal related substance Substances 0.000 description 1
- 230000007246 mechanism Effects 0.000 description 1
- 238000010943 off-gassing Methods 0.000 description 1
- 230000036961 partial effect Effects 0.000 description 1
- 230000002829 reductive effect Effects 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 description 1
- 229910052845 zircon Inorganic materials 0.000 description 1
- GFQYVLUOOAAOGM-UHFFFAOYSA-N zirconium(iv) silicate Chemical compound [Zr+4].[O-][Si]([O-])([O-])[O-] GFQYVLUOOAAOGM-UHFFFAOYSA-N 0.000 description 1
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01V—GEOPHYSICS; GRAVITATIONAL MEASUREMENTS; DETECTING MASSES OR OBJECTS; TAGS
- G01V11/00—Prospecting or detecting by methods combining techniques covered by two or more of main groups G01V1/00 - G01V9/00
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02E—REDUCTION OF GREENHOUSE GAS [GHG] EMISSIONS, RELATED TO ENERGY GENERATION, TRANSMISSION OR DISTRIBUTION
- Y02E30/00—Energy generation of nuclear origin
- Y02E30/30—Nuclear fission reactors
Landscapes
- Physics & Mathematics (AREA)
- Life Sciences & Earth Sciences (AREA)
- General Life Sciences & Earth Sciences (AREA)
- General Physics & Mathematics (AREA)
- Geophysics (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本申请实施例提供一种氦气成因资源评价方法、装置、处理器及存储介质，属于氦气成因资源评价技术领域。所述方法包括：确定目标区域的氦源岩中已衰变的子体同位素的含量；基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；确定目标区域的氦源岩中氦气的运聚系数；基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。本申请通过上述方法能够实现氦气资源的有效评价，对于开展氦气成藏和富集规律研究具有重要指导作用。
Description
技术领域
本申请涉及氦气资源评价技术领域，具体涉及一种氦气成因资源评价方法、一种氦气成因资源评价装置、一种机器可读存储介质及一种处理器。
背景技术
氦气被称为“黄金气体”。发明人发现，由于氦源岩的厚度、运聚系数以及幔源岩的难以确定，目前国内外尚没有成熟的生氦法来开展氦气资源评价方法。发明人发现，氦气百分含量法虽然计算准确，但其依赖于氦气数据点的数量和质量、天然气储量的准确性。同时在成因法计算资源量时候，由于氦源岩的厚度、运聚系数以及幔源岩的难以确定，资源量的计算误差很大。发明人还发现，尽管氦气与天然气聚集在同一圈闭中，由于生氦和生烃机制完全不同，天然气与氦气的运聚系数不存在相关性，氦气藏的运聚有其自身的规律。
鉴于上述问题，本申请提出了克服上述问题或者至少部分地解决上述问题的一种氦气成因资源评价方法。
发明内容
本申请实施例的目的是提供一种氦气成因资源评价方法、一种氦气成因资源评价装置、一种机器可读存储介质及一种处理器。所述方法能够实现氦气资源的有效评价，对于开展氦气成藏和富集规律研究具有重要指导作用。
为了实现上述目的，本申请第一方面提供一种氦气成因资源评价方法，所述方法包括：
获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
确定目标区域的氦源岩的体积；
基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
确定目标区域的氦源岩中氦气的运聚系数；
基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
其中，所述确定目标区域的氦源岩的体积，包括：
基于预设温度确定目标区域的氦源岩的底界深度；
确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
在本申请实施例中，所述确定目标区域的氦源岩的底界深度，包括：基于预设温度确定目标区域的氦源岩的底界深度。
在本申请实施例中，所述基于预设温度确定目标区域的氦源岩的底界深度，包括：
将预设温度所对应的地层深度确定为氦源岩的底界深度，所述氦源岩包括：花岗岩和火山岩。
在本申请实施例中，所述预设温度的温度范围为：220℃～240℃。
在本申请实施例中，所述确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，包括：
基于地震资料和钻井资料确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，所述氦源岩包括：煤层、泥页岩。
在本申请实施例中，基于地震资料和钻井资料确定氦源岩的分布面积，包括：
对于厚度不稳定或者分布不连续的氦源岩，基于地震资料和钻井资料，通过小面元积分法求取氦源岩的分布面积。
在本申请实施例中，通过(1)式，将预设温度所对应的地层深度确定为氦源岩的底界深度：
H底＝(t预设温度-t地表温度)/K+50 (1)；
式中，H底表示氦源岩的底界深度，单位为m；
t地表温度表示地表温度，单位为℃；
K表示地温梯度，单位为℃/100m。
在本申请实施例中，通过(2)式，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度：
h＝H底-H顶 (2)；
式中，h表示氦源岩的厚度，单位为m；
H底表示氦源岩的底界深度，单位为m；
H顶表示氦源岩的顶面深度，单位为m。
在本申请实施例中，通过(3)式，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积：
V＝h×S (3)；
式中，V表示氦源岩的体积，单位为m3；
h表示氦源岩的厚度，单位为m；
S表示氦源岩分布的面积，单位为m2。
在本申请实施例中，通过(4)式，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量：
P＝v×ρ源岩×x (4)；
式中，P表示氦源岩中未衰变的母体同位素的含量，单位为mol；
v表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
在本申请实施例中，通过(5)式，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量：
D＝(eλt-1)×v×ρ源岩×x (5)；
式中，D表示氦源岩中已衰变的子体同位素的含量，单位为mol；
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
t表示氦源岩的绝对地质年龄，单位为百万年；
v表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
在本申请实施例中，氦源岩中的放射性元素包括：238U和232Th；
针对238U衰变，通过(6)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
n＝D×8/238 (6)；
针对232Th衰变，通过(7)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
n＝D×6/232 (7)；
式中，n表示氦源岩中氦气的物质的量，D表示氦源岩中已衰变的子体同位素的含量。
在本申请实施例中，通过(8)式，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量：
VHe生＝n×M/ρhe (8)；
式中，VHe生表示生氦总量，即生氦资源量，单位为m3；
n表示氦源岩中氦气的物质的量；
M表示生氦总量的摩尔质量；
ρhe表示氦气的密度，取值0.1786g/L。
在本申请实施例中，通过(9)式，确定目标区域的氦源岩中氦气的运聚系数：
Φ＝Φ1×Φ2 (9)；
式中，Φ1为氦源岩排氦率，单位为％；Φ2为水溶氦脱气率，单位为％。在本申请实施例中，通过(10)式，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量：
Vhe聚＝Φ×Vhe生 (10)；
式中，Vhe聚表示氦气聚集量，即氦气聚集资源量，单位为m3；
Φ表示运聚系数，单位为％；
Vhe生表示生氦总量，单位为m3。
在本申请实施例中，所述氦气包括3He和4He，所述方法还包括：
根据二元复合公式，通过3He/4He的比值，确定幔源氦气和壳源氦气的资源量。
本申请第二方面提供一种氦气成因资源评价装置，所述装置包括：
获取模块，用于获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
第一确定模块，用于确定目标区域的氦源岩的体积；
第二确定模块，用于基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
第三确定模块，用于基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
第四确定模块，用于基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
第五确定模块，用于基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
第六确定模块，用于确定目标区域的氦源岩中氦气的运聚系数；
第七确定模块，用于基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
其中，所述第一确定模块，包括：
底界深度确定模块，用于基于预设温度确定目标区域的氦源岩的底界深度；
顶面深度以及分布面积确定模块，用于确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
厚度确定模块，用于基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
体积确定模块，用于基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
本申请第三方面提供一种处理器，被配置成执行上述的氦气成因资源评价方法。
本申请第四方面提供一种机器可读存储介质，该机器可读存储介质上存储有指令，该指令在被处理器执行时使得所述处理器被配置成执行上述的氦气成因资源评价方法。
与现有技术相比，本发明的上述技术方案具有如下有益效果：
本申请提出的一种氦气成因资源评价方法、装置、存储介质及处理器，所述方法通过氦源岩中已衰变的子体同位素的含量确定出氦源岩中氦气的物质的量，然后基于氦源岩中氦气的物质的量确定出氦源岩的生氦总量，最后基于氦源岩中氦气的运聚系数，确定氦气聚集资源量。本申请通过上述方法能够实现氦气资源的有效评价，对于开展氦气成藏和富集规律研究具有重要指导作用。
本申请实施例的其它特征和优点将在随后的具体实施方式部分予以详细说明。
附图说明
附图是用来提供对本申请实施例的进一步理解，并且构成说明书的一部分，与下面的具体实施方式一起用于解释本申请实施例，但并不构成对本申请实施例的限制。在附图中：
图1示意性示出了根据本申请实施例的氦气成因资源评价方法的应用环境示意图；
图2示意性示出了根据本申请实施例的氦气成因资源评价方法的整体流程示意图；
图3示意性示出了根据本申请实施例的氦气成因资源评价方法中，计算目标区域的氦源岩的体积的流程示意图；
图4示意性示出了根据本申请实施例的He、N2和CH4亨利常数随温度的变化曲线图；
图5示意性示出了根据本申请实施例的氦源岩有效厚度示意图；
图6示意性示出了根据本申请实施例的氦气成藏模式及过程示意图；
图7示意性示出了根据本申请实施例的氦气成因资源评价装置的结构框图；
图8示意性示出了根据本申请实施例的计算机设备的内部结构图。
具体实施方式
为使本申请实施例的目的、技术方案和优点更加清楚，下面将结合本申请实施例中的附图，对本申请实施例中的技术方案进行清楚、完整地描述，应当理解的是，此处所描述的具体实施方式仅用于说明和解释本申请实施例，并不用于限制本申请实施例。基于本申请中的实施例，本领域普通技术人员在没有做出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。
需要说明，若本申请实施例中有涉及“第一”、“第二”等的描述，则该“第一”、“第二”等的描述仅用于描述目的，而不能理解为指示或暗示其相对重要性或者隐含指明所指示的技术特征的数量。由此，限定有“第一”、“第二”的特征可以明示或者隐含地包括至少一个该特征。另外，各个实施例之间的技术方案可以相互结合，但是必须是以本领域普通技术人员能够实现为基础，当技术方案的结合出现相互矛盾或无法实现时应当认为这种技术方案的结合不存在，也不在本申请要求的保护范围之内。
本申请提供的氦气成因资源评价方法，可以应用于如图1所示的应用环境中。其中，终端102通过网络与服务器104通过网络进行通信。其中，终端102可以但不限于是各种个人计算机、笔记本电脑、智能手机、平板电脑和便携式可穿戴设备，服务器104可以用独立的服务器或者是多个服务器组成的服务器集群来实现。
图2示意性示出了根据本申请实施例的氦气成因资源评价方法的流程示意图。如图2所示，在本申请一实施例中，提供了一种氦气成因资源评价方法，本实施例主要以该方法应用于上述图1中的终端102(或服务器104)
来举例说明，包括以下步骤：
步骤110，获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数。
在本实施例中，目标区域指的是氦气资源研究区域。在本实施例中，氦源岩包括花岗岩、火山岩和沉积岩等。花岗岩或者火山岩的年龄，可以根据测定岩石中锆石U-Pb的同位素年龄确定。沉积岩的年龄主要通过测定岩石中古生物年代确定。氦源岩的密度一般是常数，如花岗岩的密度为2.3g/cm3。氦源岩中放射性元素包括：铀(235U、238U)和钍(232Th)。可以通过采集不同层次、不同构造部位具有代表性的氦源岩样品，来测定氦源岩中铀或者钍等放射性元素的含量，可以多次测量铀或者钍的含量，同一样品误差小于5％，求取平均值。235U的衰变常数为9.8485×10-10/a，238U的衰变常数为1.55125 ×10-10/a，232Th的衰变常数为4.9745×10-11/a。
步骤120，确定目标区域的氦源岩的体积。
图3示意性示出了根据本申请实施例的氦气成因资源评价方法中，计算目标区域的氦源岩的体积的流程示意图。如图3所示，在本实施例中，所述目标区域的氦源岩的体积通过步骤121～步骤124确定：
步骤121，基于预设温度确定目标区域的氦源岩的底界深度。
示例性地，计算氦源岩的体积，首先，需要计算氦源岩的底界深度。
若氦源岩为花岗岩或者火山岩，则所述氦源岩的底界深度可以根据220℃～240℃所对应的地层深度确定。
图4示意性示出了根据本申请实施例的He、N2和CH4亨利常数随温度的变化曲线。如图4所示，依据亨利曲线，氦气(He)在180℃时与甲烷(CH4)、氮气(N2)溶解度、分压基本趋于一致，再进一步释放氦气的潜力基本消失， 240℃时氦源岩释放潜力趋于零，因此，本实施例将研究区220℃～240℃所对应的地层深度确定为花岗岩和火山岩等氦源岩的下限深度，即氦源岩的底界深度。
图5示意性示出了根据本申请实施例的氦源岩有效厚度示意图。
示例性地，若氦源岩为花岗岩或者火山岩，则氦源岩的底界深度可以通过(1)式计算：
H底＝(t预设温度-t地表温度)/K+50 (1)；
式中，H底表示氦源岩的底界深度，单位为m；
t预设温度表示预设温度，t预设温度的取值范围为220℃～240℃。
t地表温度表示地表温度，单位为℃；
K表示地温梯度，单位为℃/100m。
若氦源岩为沉积岩，例如，煤层或者泥页岩，则所述氦源岩的底界深度可以根据地震资料和钻井资料确定。
步骤122，确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积。
在本实施例中，可以通过地震资料和钻井资料确定氦源岩的顶面深度以及氦源岩的分布面积。对于厚度不稳定或者分布不连续的氦源岩，可以通过小面元积分法求取氦源岩的分布面积。
步骤123，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度。
示例性地，氦源岩的厚度可以通过(2)式计算：
h＝H底-H顶 (2)；
式中，h表示氦源岩的厚度，单位为m；
H底表示氦源岩的底界深度，单位为m；
H顶表示氦源岩的顶面深度，单位为m。
步骤124，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
示例性地，氦源岩的体积可以通过(3)式计算：
V＝h×S (3)；
式中，V表示氦源岩的体积，单位为m3；
h表示氦源岩的厚度，单位为m；
S表示氦源岩的分布面积，单位为m2。
本实施例提出基于亨利定律的生氦法，解决了氦源岩厚度参数确定的关键难点。常规方法氦源岩的下限采用主观推断法，根据氦源岩的露头厚度来推断地下氦源岩厚度，但地下氦源岩厚度变化大，导致氦源岩的体积计算不准确。本申请以氦气成藏理论为指导，将氦气“死亡线”220℃-240℃作为氦原岩深度的下限温度，将该温度对应的深度作为氦源岩的深度下限，由此确定氦源岩厚度，有效提高了氦源岩体积和生氦资源量准确性，误差下降10％。本实施例有效解决了现有技术中存在的氦源岩的厚度难以确定的问题。
步骤130，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量。
在本实施例中，可以通过(4)计算氦源岩中未衰变的母体同位素的含量：
P＝V×ρ源岩×x (4)；
式中，P表示氦源岩中未衰变的母体同位素的含量，单位为mol；
V表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
步骤140，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量。
在本实施例中，可以通过(5)计算氦源岩中已衰变的子体同位素的含量：
D＝(eλt-1)×V×ρ源岩×x (5)；
式中，D表示氦源岩中已衰变的子体同位素的含量，单位为mol；
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
t表示氦源岩的绝对地质年龄，单位为百万年；
v表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
步骤150，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量。
在本实施例中，氦源岩中的放射性元素包括：238U和232Th；
针对238U衰变，通过(6)式，计算氦源岩中氦气的物质的量：
n＝D×8/238 (6)；
针对232Th衰变，通过(7)式，计算氦源岩中氦气的物质的量：
n＝D×6/232 (7)；
式中，n表示氦源岩中氦气的物质的量，D表示氦源岩中已衰变的子体同位素的含量。
步骤160，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量。
在本实施例中，通过(8)式，计算氦源岩的生氦总量：
VHe生＝n×M/ρhe (8)；
式中，VHe生表示生氦总量，即生氦资源量，单位为m3；
n表示氦源岩中氦气的物质的量；
M表示生氦总量的摩尔质量；
ρhe表示氦气的密度，取值0.1786g/L。
示例性地，针对238U衰变，通过(8.1)式，计算氦源岩的生氦总量：
VHe生＝n×M/ρhe
＝D×M×8/238/ρhe
＝(eλt-1)×v×ρ源岩×x×M×8/238/ρhe (8.1)；
示例性地，针对232Th衰变，通过(8.2)式，计算氦源岩的生氦总量：
VHe生＝n×M/ρhe
＝D×M×6/232/ρhe
＝(eλt-1)×v×ρ源岩×x×M×6/232/ρhe (8.2)；
式(9)-(10)中，VHe生表示研究区生氦资源量，即氦源岩的生氦总量，单位为m3；
λ表示研究区铀或者钍同位素的衰变常数，无量纲；
t表示研究区氦源岩的绝对地质年龄，单位为百万年；
ρ源岩表示研究区氦源岩的密度，单位为t/m3；例如，花岗岩的密度取值2.5t/m3；
ρhe表示研究区氦气的密度，取值0.1786g/L；
v表示研究区氦源岩的体积，单位为m3；
x表示研究区氦源岩中铀或者钍的含量，单位为g/mol。
步骤170，确定目标区域的氦源岩中氦气的运聚系数。
图6示意性示出了根据本申请实施例的氦气成藏模式及过程示意图。如图6所示，本实施例根据研究区域构造活动期次和排氦特征，确定不同成藏模式下的氦源岩排氦率和水溶氦脱气率，两者相乘即可得到运聚系数。其中，不同成藏模式是指，不同来源的氦气在不同的圈闭中的聚集形式：按照气体成分，分为含氦氮气藏、含氦CO2气藏、含氦天然气藏等；按圈闭类型分为背斜、断层以及岩性；按照来源分为幔源、壳源或混源等等。
已经排氦的氦源岩与全部氦源岩体积之比即为氦源岩排氦率。通过测定刻度区已知水溶氦含量和气藏含氦量，计算已经聚集的氦气量与水溶氦的体积，两者之比即为水溶氦脱气率。
在本实施例中，通过(9)式，计算目标区域的氦源岩中氦气的运聚系数：
Φ＝Φ1×Φ2 (9)；
式中，Φ1为氦源岩排氦率，单位为％；Φ2为水溶氦脱气率，单位为％。
其中，氦源岩排氦率Φ1通过已经排氦的氦源岩与全部氦源岩体积之比来计算。测定刻度区已知水溶氦含量和气藏含氦量，计算已经聚集的氦气量与水溶氦的体积，两者之比即为水溶氦脱气率Φ2。
在不同的成藏模式下，氦气的初次运移和二次运移系数千差万别。已有的方法没有考虑不同成藏模式的差异，也没有给出具体的计算公式。本实施例将运聚系数与烃类成藏时间相匹配，解决运聚系数难以确定的问题。
步骤180，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
本实施例根据研究区花岗岩、火成岩、沉积岩等所有生氦总量，与运聚系数相乘，即可得到研究区氦气地质资源量。
在本实施例中，通过(10)式，计算氦气聚集资源量：
Vhe聚＝Φ×Vhe生 (10)；
式中，Vhe聚表示氦气聚集量，即氦气聚集资源量，单位为m3；
Φ表示运聚系数，单位为％；
Vhe生表示生氦总量，单位为m3。
在本实施例中，所述氦气包括3He和4He，所述方法还包括：
确定壳源氦和幔源氦中3He/4He的标准值，标准值如下：
壳源氦取值标准：3He/4He＝2.0×10-8；
幔源氦取值标准：3He/4He＝1.1×10-8；
壳幔混源氦取值标准：1.1×10-8＜3He/4He＜2.0×10-8。
在本实施例中，所述氦气包括3He和4He，所述方法还包括：根据二元复合公式，通过3He/4He比值计算幔源氦气和壳源氦气的资源量，壳源氦的比例计算公式如下：
幔源氦资源量比例的计算公式如下：
I4He＝1-I3He (12)；
式中，3He表示氦气中幔源氦气所占的比例，单位为％L；
4He表示氦气中壳源氦气所占的比例，单位为％L；
R表示实测氦气值体积，单位为m3；
Rc表示壳源氦气体积，单位为m3；
Rm表示幔源氦气体积，单位为m3。
图2为一个实施例中氦气成因资源评价方法的流程示意图。应该理解的是，虽然图2的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。
图7示意性示出了根据本申请实施例的氦气成因资源评价装置的结构框图。
在一个实施例中，如图7所示，提供了一种氦气成因资源评价装置200，包括获取模块210、第一确定模块220、第二确定模块230、第三确定模块 240、第四确定模块250、第五确定模块260、第六确定模块270以及第七确定模块280，其中：
获取模块210，用于获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
第一确定模块220，用于确定目标区域的氦源岩的体积；
第二确定模块230，用于基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
第三确定模块240，用于基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
第四确定模块250，用于基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
第五确定模块260，用于基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
第六确定模块270，用于确定目标区域的氦源岩中氦气的运聚系数；
第七确定模块280，用于基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
其中，所述第一确定模块220，包括：
底界深度确定模块，用于基于预设温度确定目标区域的氦源岩的底界深度；
顶面深度以及分布面积确定模块，用于确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
厚度确定模块，用于基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
体积确定模块，用于基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
所述氦气成因资源评价装置包括处理器和存储器，上述获取模块210、第一确定模块220、第二确定模块230、第三确定模块240、第四确定模块 250、第五确定模块260、第六确定模块270以及第七确定模块280等均作为程序单元存储在存储器中，由处理器执行存储在存储器中的上述程序模块中实现相应的功能。
处理器中包含内核，由内核去存储器中调取相应的程序单元。内核可以设置一个或以上，通过调整内核参数来实现氦气成因资源评价方法。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器 (RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)，存储器包括至少一个存储芯片。
本申请实施例提供了一种存储介质，其上存储有程序，该程序被处理器执行时实现上述氦气成因资源评价方法。
图8示意性示出了根据本申请实施例的计算机设备的内部结构图。
在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图8所示。该计算机设备包括通过系统总线连接的处理器A01、网络接口A02、显示屏A04、输入装置A05和存储器(图中未示出)。其中，该计算机设备的处理器A01用于提供计算和控制能力。该计算机设备的存储器包括内存储器A03和非易失性存储介质A06。该非易失性存储介质 A06存储有操作系统B01和计算机程序B02。该内存储器A03为非易失性存储介质A06中的操作系统B01和计算机程序B02的运行提供环境。该计算机设备的网络接口A02用于与外部的终端通过网络连接通信。该计算机程序被处理器A01执行时以实现一种氦气成因资源评价方法。该计算机设备的显示屏A04可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置 A05可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。
本领域技术人员可以理解，图8中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。
在一个实施例中，本申请提供的氦气成因资源评价装置可以实现为一种计算机程序的形式，计算机程序可在如图8所示的计算机设备上运行。计算机设备的存储器中可存储组成该氦气成因资源评价装置的各个程序模块，比如，图7所示的获取模块210、第一确定模块220、第二确定模块230、第三确定模块240、第四确定模块250、第五确定模块260、第六确定模块270 以及第七确定模块280。各个程序模块构成的计算机程序使得处理器执行本说明书中描述的本申请各个实施例的氦气成因资源评价方法中的步骤。
图8所示的计算机设备可以通过如图7所示的氦气成因资源评价装置中的获取模块210执行步骤110，计算机设备可通过第一确定模块220执行步骤120，计算机设备可通过第二确定模块230执行步骤130，计算机设备可通过第三确定模块240执行步骤140，计算机设备可通过第四确定模块250 执行步骤150，计算机设备可通过第五确定模块260执行步骤160，计算机设备可通过第六确定模块270执行步骤170，计算机设备可通过第七确定模块280执行步骤180。
本申请实施例提供了一种设备，设备包括处理器、存储器及存储在存储器上并可在处理器上运行的程序，处理器执行程序时实现以下步骤：
步骤110，获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
步骤120，确定目标区域的氦源岩的体积；
步骤130，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
步骤140，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
步骤150，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
步骤160，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
步骤170，确定目标区域的氦源岩中氦气的运聚系数；
步骤180，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量。
在一个实施例中，所述步骤120，包括：
步骤121，基于预设温度确定目标区域的氦源岩的底界深度；
步骤122，确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
步骤123，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
步骤124，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
在一个实施例中，所述确定目标区域的氦源岩的底界深度，包括：
基于预设温度确定目标区域的氦源岩的底界深度。
在一个实施例中，所述基于预设温度确定目标区域的氦源岩的底界深度，包括：
将预设温度所对应的地层深度确定为氦源岩的底界深度，所述氦源岩包括：花岗岩和火山岩。
在一个实施例中，所述预设温度的温度范围为：220℃～240℃。
在一个实施例中，所述确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，包括：
基于地震资料和钻井资料确定目标区域的氦源岩的底界深度、氦源岩的顶面深度以及氦源岩的分布面积，所述氦源岩包括：煤层、泥页岩。
在一个实施例中，所述基于地震资料和钻井资料确定氦源岩的分布面积，包括：
对于厚度不稳定或者分布不连续的氦源岩，基于地震资料和钻井资料，通过小面元积分法求取氦源岩的分布面积。
在一个实施例中，通过(1)式，将预设温度所对应的地层深度确定为氦源岩的底界深度：
H底＝(t预设温度-t地表温度)/K+50 (1)；
式中，H底表示氦源岩的底界深度，单位为m；
t地表温度表示地表温度，单位为℃；
K表示地温梯度，单位为℃/100m。
在一个实施例中，通过(2)式，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度：
h＝H底-H顶 (2)；
式中，h表示氦源岩的厚度，单位为m；
H底表示氦源岩的底界深度，单位为m；
H顶表示氦源岩的顶面深度，单位为m。
在一个实施例中，通过(3)式，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积：
V＝h×S (3)；
式中，V表示氦源岩的体积，单位为m3；
h表示氦源岩的厚度，单位为m；
S表示氦源岩的分布面积，单位为m2。
在一个实施例中，通过(4)式，基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量：
P＝v×ρ源岩×x (4)；
式中，P表示氦源岩中未衰变的母体同位素的含量，单位为mol；
v表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
在一个实施例中，通过(5)式，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量：
D＝(eλt-1)×v×ρ源岩×x (5)；
式中，D表示氦源岩中已衰变的子体同位素的含量，单位为mol；
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
t表示氦源岩的绝对地质年龄，单位为百万年；
v表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
在一个实施例中，氦源岩中的放射性元素包括：238U和232Th；
针对238U衰变，通过(6)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
n＝D×8/238 (6)；
针对232Th衰变，通过(7)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
n＝D×6/232 (7)；
式中，n表示氦源岩中氦气的物质的量，D表示氦源岩中已衰变的子体同位素的含量。
在一个实施例中，通过(8)式，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量：
VHe生＝n×M/ρhe (8)；
式中，VHe生表示生氦总量，即生氦资源量，单位为m3；
n表示氦源岩中氦气的物质的量；
M表示生氦总量的摩尔质量；
ρhe表示氦气的密度，取值0.1786g/L。
在一个实施例中，通过(9)式，确定目标区域的氦源岩中氦气的运聚系数：
Φ＝Φ1×Φ2 (9)；
式中，Φ1为氦源岩排氦率，单位为％；Φ2为水溶氦脱气率，单位为％。
在一个实施例中，通过(10)式，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量：
Vhe聚＝Φ×Vhe生 (10)；
式中，Vhe聚表示氦气聚集量，即氦气聚集资源量，单位为m3；
Φ表示运聚系数，单位为％；
Vhe生表示生氦总量，单位为m3。
在一个实施例中，所述氦气包括3He和4He，所述方法还包括：根据二元复合公式，通过3He/4He的比值，确定幔源氦气和壳源氦气的资源量。
本领域内的技术人员应明白，本申请的实施例可提供为方法、系统、或计算机程序产品。因此，本申请可采用完全硬件实施例、完全软件实施例、或结合软件和硬件方面的实施例的形式。而且，本申请可采用在一个或多个其中包含有计算机可用程序代码的计算机可用存储介质(包括但不限于磁盘存储器、CD-ROM、光学存储器等)上实施的计算机程序产品的形式。
本申请是参照根据本申请实施例的方法、设备(系统)、和计算机程序产品的流程图和/或方框图来描述的。应理解可由计算机程序指令实现流程图和/或方框图中的每一流程和/或方框、以及流程图和/或方框图中的流程和/ 或方框的结合。可提供这些计算机程序指令到通用计算机、专用计算机、嵌入式处理机或其他可编程数据处理设备的处理器以产生一个机器，使得通过计算机或其他可编程数据处理设备的处理器执行的指令产生用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的装置。
这些计算机程序指令也可存储在能引导计算机或其他可编程数据处理设备以特定方式工作的计算机可读存储器中，使得存储在该计算机可读存储器中的指令产生包括指令装置的制造品，该指令装置实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能。
这些计算机程序指令也可装载到计算机或其他可编程数据处理设备上，使得在计算机或其他可编程设备上执行一系列操作步骤以产生计算机实现的处理，从而在计算机或其他可编程设备上执行的指令提供用于实现在流程图一个流程或多个流程和/或方框图一个方框或多个方框中指定的功能的步骤。
在一个典型的配置中，计算设备包括一个或多个处理器(CPU)、输入/ 输出接口、网络接口和内存。
存储器可能包括计算机可读介质中的非永久性存储器，随机存取存储器 (RAM)和/或非易失性内存等形式，如只读存储器(ROM)或闪存(flash RAM)。存储器是计算机可读介质的示例。
计算机可读介质包括永久性和非永久性、可移动和非可移动媒体，可以由任何方法或技术来实现信息存储。信息可以是计算机可读指令、数据结构、程序的模块或其他数据。计算机的存储介质的例子包括，但不限于相变内存(PRAM)、静态随机存取存储器(SRAM)、动态随机存取存储器(DRAM)、其他类型的随机存取存储器(RAM)、只读存储器(ROM)、电可擦除可编程只读存储器(EEPROM)、快闪记忆体或其他内存技术、只读光盘只读存储器 (CD-ROM)、数字多功能光盘(DVD)或其他光学存储、磁盒式磁带，磁带磁磁盘存储或其他磁性存储设备或任何其他非传输介质，可用于存储可以被计算设备访问的信息。按照本文中的界定，计算机可读介质不包括暂存电脑可读媒体(transitory media)，如调制的数据信号和载波。
还需要说明的是，术语“包括”、“包含”或者其任何其他变体意在涵盖非排他性的包含，从而使得包括一系列要素的过程、方法、商品或者设备不仅包括那些要素，而且还包括没有明确列出的其他要素，或者是还包括为这种过程、方法、商品或者设备所固有的要素。在没有更多限制的情况下，由语句“包括一个……”限定的要素，并不排除在包括要素的过程、方法、商品或者设备中还存在另外的相同要素。
以上仅为本申请的实施例而已，并不用于限制本申请。对于本领域技术人员来说，本申请可以有各种更改和变化。凡在本申请的精神和原理之内所作的任何修改、等同替换、改进等，均应包含在本申请的权利要求范围之内。
Claims (15)
1.一种氦气成因资源评价方法，其特征在于，所述方法包括：
获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
确定目标区域的氦源岩的体积；
基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
确定目标区域的氦源岩中氦气的运聚系数；
基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量；
其中，所述确定目标区域的氦源岩的体积，包括：
基于预设温度确定目标区域的氦源岩的底界深度；
确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
2.根据权利要求1所述的方法，其特征在于，所述基于预设温度确定目标区域的氦源岩的底界深度，包括：
将预设温度所对应的地层深度确定为氦源岩的底界深度，所述氦源岩包括：花岗岩和火山岩。
3.根据权利要求2所述的方法，其特征在于，所述预设温度的温度范围为：220℃～240℃。
4.根据权利要求1所述的方法，其特征在于，所述确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积，包括：
基于地震资料和钻井资料确定氦源岩的顶面深度以及氦源岩的分布面积，所述氦源岩包括：煤层、泥页岩；
其中，基于地震资料和钻井资料确定氦源岩的分布面积，包括：对于厚度不稳定或者分布不连续的氦源岩，基于地震资料和钻井资料，通过小面元积分法求取氦源岩的分布面积。
5.根据权利要求3所述的方法，其特征在于，通过(1)式，将预设温度所对应的地层深度确定为氦源岩的底界深度：
H底＝(t预设温度-t地表温度)/K+50 (1)；
式中，H底表示氦源岩的底界深度，单位为m；
t地表温度表示地表温度，单位为℃；
K表示地温梯度，单位为℃/100m。
6.根据权利要求5所述的方法，其特征在于，通过(2)式，基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度：
h＝H底-H顶 (2)；
式中，h表示氦源岩的厚度，单位为m；
H底表示氦源岩的底界深度，单位为m；
H顶表示氦源岩的顶面深度，单位为m。
7.根据权利要求6所述的方法，其特征在于，通过(3)式，基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积：
V＝h×S (3)；
式中，V表示氦源岩的体积，单位为m3；
h表示氦源岩的厚度，单位为m；
S表示氦源岩的分布面积，单位为m2。
8.根据权利要求7所述的方法，其特征在于，通过(5)式，基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量：
式中，D表示氦源岩中已衰变的子体同位素的含量，单位为mol；
v×ρ源岩×x表示氦源岩中未衰变的母体同位素的含量；
λ表示氦源岩中放射性元素的同位素的衰变常数，无量纲；
t表示氦源岩的绝对地质年龄，单位为百万年；
v表示氦源岩的体积，单位为m3；
ρ源岩表示氦源岩的密度，单位为t/m3；
x表示氦源岩中放射性元素的含量，单位为mol/g。
9.根据权利要求8所述的方法，其特征在于，氦源岩中的放射性元素包括：238U和232Th；
针对238U衰变，通过(6)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
n＝D×8/238 (6)；
针对232Th衰变，通过(7)式，基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量：
n＝D×6/232 (7)；
式中，n表示氦源岩中氦气的物质的量，D表示氦源岩中已衰变的子体同位素的含量。
10.根据权利要求9所述的方法，其特征在于，通过(8)式，基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量：
VHe生＝n×M/ρhe (8)；
式中，VHe生表示生氦总量，即生氦资源量，单位为m3；
n表示氦源岩中氦气的物质的量；
M表示生氦总量的摩尔质量；
ρhe表示氦气的密度，取值0.1786g/L。
11.根据权利要求10所述的方法，其特征在于，通过(9)式，确定目标区域的氦源岩中氦气的运聚系数：
Φ＝Φ1×Φ2 (9)；
式中，Φ1为氦源岩排氦率，单位为％；Φ2为水溶氦脱气率，单位为％；
通过(10)式，基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量：
Vhe聚＝Φ×Vhe生 (10)；
式中，Vhe聚表示氦气聚集量，即氦气聚集资源量，单位为m3；
Φ表示运聚系数，单位为％；
Vhe生表示生氦总量，单位为m3。
12.根据权利要求1所述的方法，其特征在于，所述氦气包括3He和4He，所述方法还包括：
根据二元复合公式，通过3He/4He的比值，确定幔源氦气和壳源氦气的资源量。
13.一种氦气成因资源评价装置，其特征在于，所述装置包括：
获取模块，用于获取目标区域的氦源岩的绝对年龄、氦源岩的密度、氦源岩中放射性元素的含量以及氦源岩中放射性元素的衰变常数；
第一确定模块，用于确定目标区域的氦源岩的体积；
第二确定模块，用于基于所述氦源岩的密度、所述氦源岩中放射性元素的含量以及所述氦源岩的体积，确定氦源岩中未衰变的母体同位素的含量；
第三确定模块，用于基于所述氦源岩的绝对年龄、所述氦源岩中放射性元素的衰变常数以及所述氦源岩中未衰变的母体同位素的含量，确定氦源岩中已衰变的子体同位素的含量；
第四确定模块，用于基于所述氦源岩中已衰变的子体同位素的含量确定氦源岩中氦气的物质的量；
第五确定模块，用于基于所述氦源岩中氦气的物质的量确定氦源岩的生氦总量；
第六确定模块，用于确定目标区域的氦源岩中氦气的运聚系数；
第七确定模块，用于基于所述运聚系数和所述氦源岩的生氦总量确定氦气聚集资源量；
其中，所述第一确定模块，包括：
底界深度确定模块，用于基于预设温度确定目标区域的氦源岩的底界深度；
顶面深度以及分布面积确定模块，用于确定目标区域的氦源岩的顶面深度以及氦源岩的分布面积；
厚度确定模块，用于基于所述氦源岩的底界深度和氦源岩的顶面深度确定氦源岩的厚度；
体积确定模块，用于基于所述氦源岩的厚度以及所述氦源岩的分布面积确定氦源岩的体积。
14.一种处理器，其特征在于，被配置成执行根据权利要求1至12中任一项所述的氦气成因资源评价方法。
15.一种机器可读存储介质，该机器可读存储介质上存储有指令，其特征在于，该指令在被处理器执行时使得所述处理器被配置成执行根据权利要求1至12中任一项所述的氦气成因资源评价方法。
Priority Applications (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211271154.5A CN115685377B (zh) | 2022-10-17 | 2022-10-17 | 氦气成因资源评价方法及装置 |
| PCT/CN2022/138085 WO2024082401A1 (zh) | 2022-10-17 | 2022-12-09 | 氦气成因资源评价方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202211271154.5A CN115685377B (zh) | 2022-10-17 | 2022-10-17 | 氦气成因资源评价方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN115685377A true CN115685377A (zh) | 2023-02-03 |
| CN115685377B CN115685377B (zh) | 2023-05-19 |
Family
ID=85066765
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202211271154.5A Active CN115685377B (zh) | 2022-10-17 | 2022-10-17 | 氦气成因资源评价方法及装置 |
Country Status (2)
| Country | Link |
|---|---|
| CN (1) | CN115685377B (zh) |
| WO (1) | WO2024082401A1 (zh) |
Cited By (4)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117250661A (zh) * | 2023-09-19 | 2023-12-19 | 中国石油天然气股份有限公司 | 一种富氦气藏区带的评价方法、装置、设备及介质 |
| CN117332553A (zh) * | 2023-08-18 | 2024-01-02 | 中国地质大学（北京） | 页岩储层氦气同位素参数获取方法、终端及介质 |
| CN118114008A (zh) * | 2024-03-07 | 2024-05-31 | 中国地质大学（北京） | 自生自储型氦气有利区确定方法、终端及介质 |
| CN118330145A (zh) * | 2024-04-25 | 2024-07-12 | 中国石油大学(华东) | 一种基于瑞利分馏模型评价氦气二次运移过程中脱气程度的方法 |
Families Citing this family (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN118583907A (zh) * | 2024-06-04 | 2024-09-03 | 中国矿业大学(北京) | 富氦天然气成藏效率评价方法、装置、系统及存储介质 |
Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4378055A (en) * | 1977-08-05 | 1983-03-29 | Phillips Petroleum Company | Analyzing for helium in drilling muds to locate geothermal reservoirs |
| US20180202264A1 (en) * | 2016-02-08 | 2018-07-19 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
| US20190339407A1 (en) * | 2018-05-03 | 2019-11-07 | Chevron U.S.A. Inc. | System and method for mapping hydrocarbon source rock using seismic attributes |
| CN111308577A (zh) * | 2020-03-17 | 2020-06-19 | 中国石油化工股份有限公司 | 一种针对氦气气藏定量勘探中地区参数的确定方法 |
| CN111338001A (zh) * | 2020-03-17 | 2020-06-26 | 中国石油化工股份有限公司 | 一种氦气气藏的勘探方法 |
| CN114910976A (zh) * | 2022-04-18 | 2022-08-16 | 中国科学院西北生态环境资源研究院 | 低勘探程度区氦气资源潜力地质评价方法 |
-
2022
- 2022-10-17 CN CN202211271154.5A patent/CN115685377B/zh active Active
- 2022-12-09 WO PCT/CN2022/138085 patent/WO2024082401A1/zh unknown
Patent Citations (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US4378055A (en) * | 1977-08-05 | 1983-03-29 | Phillips Petroleum Company | Analyzing for helium in drilling muds to locate geothermal reservoirs |
| US20180202264A1 (en) * | 2016-02-08 | 2018-07-19 | Rs Energy Group Topco, Inc. | Method for estimating oil/gas production using statistical learning models |
| US20190339407A1 (en) * | 2018-05-03 | 2019-11-07 | Chevron U.S.A. Inc. | System and method for mapping hydrocarbon source rock using seismic attributes |
| CN111308577A (zh) * | 2020-03-17 | 2020-06-19 | 中国石油化工股份有限公司 | 一种针对氦气气藏定量勘探中地区参数的确定方法 |
| CN111338001A (zh) * | 2020-03-17 | 2020-06-26 | 中国石油化工股份有限公司 | 一种氦气气藏的勘探方法 |
| CN114910976A (zh) * | 2022-04-18 | 2022-08-16 | 中国科学院西北生态环境资源研究院 | 低勘探程度区氦气资源潜力地质评价方法 |
Cited By (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117332553A (zh) * | 2023-08-18 | 2024-01-02 | 中国地质大学（北京） | 页岩储层氦气同位素参数获取方法、终端及介质 |
| CN117332553B (zh) * | 2023-08-18 | 2024-03-29 | 中国地质大学（北京） | 页岩储层氦气同位素参数获取方法、终端及介质 |
| CN117250661A (zh) * | 2023-09-19 | 2023-12-19 | 中国石油天然气股份有限公司 | 一种富氦气藏区带的评价方法、装置、设备及介质 |
| CN117250661B (zh) * | 2023-09-19 | 2024-05-10 | 中国石油天然气股份有限公司 | 一种富氦气藏区带的评价方法、装置、设备及介质 |
| CN118114008A (zh) * | 2024-03-07 | 2024-05-31 | 中国地质大学（北京） | 自生自储型氦气有利区确定方法、终端及介质 |
| CN118330145A (zh) * | 2024-04-25 | 2024-07-12 | 中国石油大学(华东) | 一种基于瑞利分馏模型评价氦气二次运移过程中脱气程度的方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| WO2024082401A1 (zh) | 2024-04-25 |
| CN115685377B (zh) | 2023-05-19 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN115685377A (zh) | 氦气成因资源评价方法及装置 | |
| US20200211126A1 (en) | Prediction method for shale oil and gas sweet spot region, computer device and computer readable storage medium | |
| US11037084B2 (en) | Method and apparatus for evaluating exploitation value of a geological resource | |
| CN101806215B (zh) | 用束缚水饱和度数据判别储层流体类型的方法 | |
| EP3181804A1 (fr) | Procede pour caracteriser et exploiter une formation souterraine comprenant un reseau de fractures | |
| US20130282286A1 (en) | System and method for calibrating permeability for use in reservoir modeling | |
| CN110470581B (zh) | 确定储层应力敏感程度的方法、装置及存储介质 | |
| CN108005644B (zh) | 一种倾斜煤层动态渗透率预测方法及装置 | |
| CN116430474A (zh) | 一种地热资源储量计算方法、系统及电子设备 | |
| CN103345002A (zh) | 一种测井资料获取砂岩卤水层钾离子含量的方法 | |
| CN112862139A (zh) | 流体类型预测模型的构建方法、流体类型预测方法和装置 | |
| CN111413741A (zh) | 一种砂岩型铀矿资源量计算方法和装置 | |
| CN114428365A (zh) | 一种基于含气特征的致密砂岩气层测井识别方法及装置 | |
| CN106960264A (zh) | 低渗透砂砾岩地层压力预测方法及装置 | |
| CN108492014B (zh) | 一种确定地质资源量的数据处理方法及装置 | |
| CN110879920B (zh) | 一种获取沉积盆地坳陷区古地表热流的方法、装置及系统 | |
| Siegenthaler et al. | Model comparison of 14C and 10Be isotope records | |
| CN114580143A (zh) | 矿井涌水量动态预测方法及相关装置 | |
| Davis et al. | A dual porosity CNL logging system | |
| G. Ghedan et al. | An improved approach for generating saturation functions for simulation models using dynamic rock types | |
| CN115146889A (zh) | 气藏开发价值评估方法、系统、存储介质及电子设备 | |
| CN105528732B (zh) | 一种试气井产能预测方法 | |
| Laronga et al. | Time-Lapse Pulsed-Neutron Logs for CCS: What Have We Learned from All These Monitoring Runs? | |
| CN112861061A (zh) | 地质分层方法及系统、处理设备 | |
| CN117332553B (zh) | 页岩储层氦气同位素参数获取方法、终端及介质 |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |