CN118471358A - In-situ helium resource quantity evaluation method considering porosity and helium generation age correction - Google Patents
In-situ helium resource quantity evaluation method considering porosity and helium generation age correction Download PDFInfo
- Publication number
- CN118471358A CN118471358A CN202410634435.5A CN202410634435A CN118471358A CN 118471358 A CN118471358 A CN 118471358A CN 202410634435 A CN202410634435 A CN 202410634435A CN 118471358 A CN118471358 A CN 118471358A
- Authority
- CN
- China
- Prior art keywords
- helium
- situ
- pressure
- rock
- porosity
- 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.)
- Pending
Links
- 239000001307 helium Substances 0.000 title claims abstract description 130
- 229910052734 helium Inorganic materials 0.000 title claims abstract description 130
- SWQJXJOGLNCZEY-UHFFFAOYSA-N helium atom Chemical compound [He] SWQJXJOGLNCZEY-UHFFFAOYSA-N 0.000 title claims abstract description 130
- 238000011065 in-situ storage Methods 0.000 title claims abstract description 56
- 238000012937 correction Methods 0.000 title claims abstract description 18
- 238000011156 evaluation Methods 0.000 title claims abstract description 15
- 239000011435 rock Substances 0.000 claims abstract description 83
- 238000000034 method Methods 0.000 claims abstract description 30
- XLYOFNOQVPJJNP-UHFFFAOYSA-N water Substances O XLYOFNOQVPJJNP-UHFFFAOYSA-N 0.000 claims abstract description 25
- 229910052770 Uranium Inorganic materials 0.000 claims abstract description 20
- VNWKTOKETHGBQD-UHFFFAOYSA-N methane Chemical compound C VNWKTOKETHGBQD-UHFFFAOYSA-N 0.000 claims abstract description 14
- 238000012360 testing method Methods 0.000 claims abstract description 11
- 239000003345 natural gas Substances 0.000 claims abstract description 7
- 238000004458 analytical method Methods 0.000 claims abstract description 6
- 239000007789 gas Substances 0.000 claims description 34
- 238000004364 calculation method Methods 0.000 claims description 14
- 238000011160 research Methods 0.000 claims description 12
- 229910052776 Thorium Inorganic materials 0.000 claims description 10
- 230000015572 biosynthetic process Effects 0.000 claims description 9
- 230000006835 compression Effects 0.000 claims description 9
- 238000007906 compression Methods 0.000 claims description 9
- 239000011148 porous material Substances 0.000 claims description 8
- 238000005070 sampling Methods 0.000 claims description 6
- 238000005553 drilling Methods 0.000 claims description 5
- 239000011159 matrix material Substances 0.000 claims description 5
- 238000012546 transfer Methods 0.000 claims description 5
- 238000000691 measurement method Methods 0.000 claims description 4
- 238000006073 displacement reaction Methods 0.000 claims description 3
- 238000011067 equilibration Methods 0.000 claims description 3
- 238000007654 immersion Methods 0.000 claims description 3
- QSHDDOUJBYECFT-UHFFFAOYSA-N mercury Chemical compound [Hg] QSHDDOUJBYECFT-UHFFFAOYSA-N 0.000 claims description 3
- 229910052753 mercury Inorganic materials 0.000 claims description 3
- 238000007789 sealing Methods 0.000 claims description 3
- 238000001612 separation test Methods 0.000 claims description 3
- 230000000977 initiatory effect Effects 0.000 claims 1
- 238000003860 storage Methods 0.000 abstract description 2
- 239000004215 Carbon black (E152) Substances 0.000 description 6
- 229930195733 hydrocarbon Natural products 0.000 description 6
- 150000002430 hydrocarbons Chemical group 0.000 description 6
- 229910052500 inorganic mineral Inorganic materials 0.000 description 4
- 239000011707 mineral Substances 0.000 description 4
- 238000004519 manufacturing process Methods 0.000 description 3
- 230000006978 adaptation Effects 0.000 description 2
- 229910052586 apatite Inorganic materials 0.000 description 2
- 238000002474 experimental method Methods 0.000 description 2
- 239000010438 granite Substances 0.000 description 2
- 239000003673 groundwater Substances 0.000 description 2
- 229910052595 hematite Inorganic materials 0.000 description 2
- 239000011019 hematite Substances 0.000 description 2
- LIKBJVNGSGBSGK-UHFFFAOYSA-N iron(3+);oxygen(2-) Chemical compound [O-2].[O-2].[O-2].[Fe+3].[Fe+3] LIKBJVNGSGBSGK-UHFFFAOYSA-N 0.000 description 2
- 238000013178 mathematical model Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- VSIIXMUUUJUKCM-UHFFFAOYSA-D pentacalcium;fluoride;triphosphate Chemical compound [F-].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[Ca+2].[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O.[O-]P([O-])([O-])=O VSIIXMUUUJUKCM-UHFFFAOYSA-D 0.000 description 2
- 230000002285 radioactive effect Effects 0.000 description 2
- ZSLUVFAKFWKJRC-IGMARMGPSA-N 232Th Chemical compound [232Th] ZSLUVFAKFWKJRC-IGMARMGPSA-N 0.000 description 1
- UFHFLCQGNIYNRP-UHFFFAOYSA-N Hydrogen Chemical compound [H][H] UFHFLCQGNIYNRP-UHFFFAOYSA-N 0.000 description 1
- RTAQQCXQSZGOHL-UHFFFAOYSA-N Titanium Chemical compound [Ti] RTAQQCXQSZGOHL-UHFFFAOYSA-N 0.000 description 1
- 229910001570 bauxite Inorganic materials 0.000 description 1
- 230000009286 beneficial effect Effects 0.000 description 1
- 230000000903 blocking effect Effects 0.000 description 1
- 239000003245 coal Substances 0.000 description 1
- 239000010779 crude oil Substances 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000000151 deposition Methods 0.000 description 1
- 238000011161 development Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 239000001257 hydrogen Substances 0.000 description 1
- 229910052739 hydrogen Inorganic materials 0.000 description 1
- 239000007788 liquid Substances 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 239000003921 oil Substances 0.000 description 1
- 239000002245 particle Substances 0.000 description 1
- 230000000704 physical effect Effects 0.000 description 1
- 230000005258 radioactive decay Effects 0.000 description 1
- 238000001179 sorption measurement Methods 0.000 description 1
- 229910052719 titanium Inorganic materials 0.000 description 1
- 239000010936 titanium Substances 0.000 description 1
- JFALSRSLKYAFGM-UHFFFAOYSA-N uranium(0) Chemical compound [U] JFALSRSLKYAFGM-UHFFFAOYSA-N 0.000 description 1
- 238000012795 verification Methods 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
Classifications
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/10—Analysis or design of chemical reactions, syntheses or processes
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C10/00—Computational theoretical chemistry, i.e. ICT specially adapted for theoretical aspects of quantum chemistry, molecular mechanics, molecular dynamics or the like
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/20—Identification of molecular entities, parts thereof or of chemical compositions
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/30—Prediction of properties of chemical compounds, compositions or mixtures
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16C—COMPUTATIONAL CHEMISTRY; CHEMOINFORMATICS; COMPUTATIONAL MATERIALS SCIENCE
- G16C20/00—Chemoinformatics, i.e. ICT specially adapted for the handling of physicochemical or structural data of chemical particles, elements, compounds or mixtures
- G16C20/70—Machine learning, data mining or chemometrics
Landscapes
- Engineering & Computer Science (AREA)
- Computing Systems (AREA)
- Theoretical Computer Science (AREA)
- Chemical & Material Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Bioinformatics & Cheminformatics (AREA)
- Bioinformatics & Computational Biology (AREA)
- Crystallography & Structural Chemistry (AREA)
- Health & Medical Sciences (AREA)
- General Health & Medical Sciences (AREA)
- Spectroscopy & Molecular Physics (AREA)
- Physics & Mathematics (AREA)
- Artificial Intelligence (AREA)
- Databases & Information Systems (AREA)
- Evolutionary Computation (AREA)
- Medical Informatics (AREA)
- Software Systems (AREA)
- Data Mining & Analysis (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Analytical Chemistry (AREA)
- Chemical Kinetics & Catalysis (AREA)
- Investigating Strength Of Materials By Application Of Mechanical Stress (AREA)
Abstract
The invention relates to an in-situ helium resource quantity evaluation method taking both porosity and helium generation 'age' correction into consideration, belonging to the technical field of helium source identification and storage process analysis, comprising the following steps: defining a target producing interval, determining the helium rock type of the producing interval, and determining the volume and the quality of source rock; basic data are acquired, and the solubility coefficient of reservoir temperature, reservoir pressure and helium in water is defined; testing and analyzing to obtain the He content c in the natural gas of samples with different depths; obtaining an average abundance of initial U, th; the plunger sample helium method with highest accuracy and reliability is selected to represent the porosity of the shale, and the 4 He age is corrected, so that compared with the prior art, the calculated in-situ helium resource amount is more accurate and scientific, and better persuasion is achieved.
Description
Technical Field
The invention relates to an in-situ helium resource quantity evaluation method taking both porosity and helium generation 'age' correction into consideration, and further relates to a method for accurately evaluating the in-situ yield of shell source 4 He and the in-situ helium resource quantity, belonging to the technical field of helium source identification and reservoir forming process analysis.
Background
Helium is an element in the universe that is next to hydrogen and accounts for about 23% of the mass of the entire star system. Helium is widely available in the universe, but helium resources within the earth are scarce. Geochemical studies of helium in rare gases have mainly used two stable isotopes of 3 He and 4 He, which differ in their origin and production on earth. Among them, shell source 4 He is a main source of industrial helium, and is generated by decay of radioactive elements such as uranium (U) and thorium (Th) in rock. In view of the helium-rich gas fields that have been found to date, a sufficient and good quality helium source is the basis for the formation of helium-rich fields. Older sedimentary rocks (shale, bauxite, coal rock, etc.), magma rocks (granite, etc.), metamorphic rocks (gneiss, etc.), which are rich in radioactive U, th elements, can all be effective helium source rocks. In the prior in-situ helium resource amount calculating method, the in-situ helium yield is calculated only according to the stratum age, development characteristics, stratum porosity and U, th average content of helium source rock. However, in the calculation of the amount of helium resources in situ in a typical sedimentary helium source rock (shale), the following two problems are often ignored:
(1) Shale acts as an effective helium source rock, while in hydrocarbon-bearing reservoirs, it can act as a hydrocarbon source rock for hydrocarbon production. The hydrocarbon-based rock hydrocarbon-producing process can be divided into different stages according to temperature and maturity, namely: immature stage (temperature typically between 60-100 ℃, maturity between 0.5% -0.7%), low mature stage (temperature typically between 80-150 ℃, maturity between 0.7% -1.3%), high mature stage (temperature typically between 150-200 ℃, maturity between 1.3% -2.0%) and overmature stage (temperature typically between 200-300 ℃, maturity greater than 2.0%). While initial helium generation often requires breakthrough of the closure temperature of the different minerals bounding the helium element. For example, the blocking temperature of apatite is 55-110 ℃ (average 70 ℃), hematite is 90-250 ℃, zircon is 180-200 ℃, and titanium ore is about 200 ℃. Wherein, the old sedimentary rock is rich in apatite and hematite, namely, the old sedimentary rock can really enter into the helium generation stage when the temperature reaches 70 ℃, namely, the immature stage of the crude oil gas. This shows that helium will start to be released from minerals 2 km or more at a normal crust temperature gradient of 30 ℃/km and a surface temperature of 10 ℃. Therefore, unlike the conventional studies, the formation age cannot be regarded as the "age" of helium in general, and errors in helium resource evaluation are caused.
(2) The measurement of porosity is important for in situ helium production calculations, which represent the amount of helium that can be stored in groundwater as well as dissolved in groundwater. The shale minerals have complex composition, compact rock structure and strong heterogeneity, and the accurate test of the porosity is more difficult than the conventional experiment (CO 2 adsorption method). The rock porosity applied in the original in-situ helium yield calculation is not completely aimed at depositing helium source rock (hydrocarbon source rock), because the rock porosity which belongs to magma rock (granite) and metamorphic rock (gneiss) in helium source rock is high in brittleness index, large in mineral particles, strong in reservoir physical properties and other characteristics, and large in difference from dense and heterogeneous shale mainly developing micro-nano pores, and the porosity obtained by the conventional experiment is not completely applicable. That is, the porosity of different types of helium source rock cannot be approximated. In addition, in the previous calculation of in-situ helium resource amount, the porosity of the one-sided helium source rock of different types is calculated by 10%, and helium resource evaluation cannot be accurately performed.
Thus, there is a need for more accurate correction of the 4 He age and porosity parameters of shale to better service the calculation of in situ helium resource amounts.
Disclosure of Invention
Aiming at the defects of the existing evaluation analysis method, namely that the original shell source 4 He in-situ yield calculation formula is inaccurate in obtaining age and porosity parameters, the invention provides the in-situ helium resource amount evaluation method which takes account of the correction of the porosity and the raw helium 'age', corrects the age and the porosity, and calculates the obtained in-situ helium resource amount more accurately and scientifically.
The invention adopts the following technical scheme:
An in-situ helium resource amount evaluation method taking both porosity and helium generation 'age' correction into consideration comprises the following steps:
s1: defining a target producing interval, selecting a typical or key natural gas (rare gas) exploration horizon, determining the type of helium rock (deposited helium source rock) of the producing interval, and determining the volume and quality of the source rock;
S2: basic data are acquired, and the reservoir temperature, the reservoir pressure and the solubility coefficient of helium in water of a target interval sample are clearly researched;
S3: sampling at the wellhead of a typical well in a research area or sampling while drilling during logging, sealing to prevent helium leakage, and performing separation test to obtain the He content c in natural gas samples with different depths;
S4: u, th the half-life period of the main natural isotope is extremely long, so that the U, th concentration in the stratum nowadays can be used for approximately replacing the average abundance of the initial U, th deposit, and the U and Th contents can be obtained through geochemical analysis of the drilling core of a classical well;
s5: the in situ 4 He yield was calculated from equation (1):
Wherein ρ r is rock density in g/cm -3; t is the hiding time, i.e. 4 He age, per unit year; Is porosity; Λ=1, which is the efficiency of He transfer from the rock matrix to water; cm 3STPg-1 H2O is 4 He yield per gram of water at standard conditions; p (4 He) is the in situ 4 He yield obtained from equation (2):
P(4He)＝1.207×10-13[U]+2.867×10-13[Th] cm3STPg-1 rockyr-1 (2)
Wherein [ U ] and [ Th ] respectively represent average abundance of U and Th in a research area, and unit ppm is obtained in the step S4; substituting the average abundance of U and Th to obtain the in situ 4 He yield P (4He),cm3STPg-1 rockyr-1 is 4 He yield per gram of rock per year in standard state;
S6: correcting 4 He age T in the formula (1) to obtain corrected 4 He age T
In combination with the partial pressure in the study area, the solubility coefficient of helium in water S He and the annual helium content α in water entering the 1cm 3 pore from rock, the corrected 4 He "age" T is calculated by equation (2):
wherein c is the He content obtained in the step S3,%; p is the reservoir pressure obtained in step S2, atm; s He is the solubility coefficient of helium in water, cm 3/cm3. Atm; alpha is a fixed value, and 3.4X10 -12a-1 is taken; cp is the partial pressure of helium in the gas layer;
s7: porosity calculation using plunger-like helium method Substituting the porosity, corrected 4 He 'age' T, in-situ 4 He yield P (4 He), he transfer efficiency lambda (lambda=1) from the rock matrix to water and rock density into a formula (1) to calculate in-situ 4 He yield in the research area, and multiplying the obtained in-situ 4 He yield by source rock mass to obtain in-situ helium resource quantity in the research area.
Preferably, in step S1, the source rock volume is determined by the source rock thickness and area;
and collecting a deposited helium source rock sample, performing high-pressure mercury testing to obtain rock density, and obtaining the quality of the source rock according to the volume of the source rock and the rock density.
Preferably, in step S2, the reservoir temperature may be replaced with a measured downhole temperature or a drill bit temperature; reservoir pressure can be obtained by logging, formation pressure testing or mathematical model and the like, and can be replaced by formation pressure; the solubility coefficient of helium in water was calculated by PHREEQC software in combination with the reservoir temperature, reservoir pressure of the sample.
Preferably, in step S7, the process of calculating the porosity by using the plunger-like helium method is as follows:
Firstly, placing a sample into a sample chamber, vacuumizing the sample chamber to a fixed value p 1, injecting helium into a reference chamber to a certain pressure, recording the gas pressure p 2 in the reference chamber when the helium pressure in the reference chamber is balanced, then communicating the reference chamber with the sample chamber to fully saturate the sample pore with helium, recording the gas pressure p 3 in the sample chamber when the helium pressure in the sample chamber is balanced, and obtaining a formula (4) according to Boyle's law; under the conditions that the testing process is at constant temperature and the valve displacement volume is ignored, the shale sample skeleton volume Vg is calculated by the simultaneous formula (4) and the formula (5):
wherein, p 1 represents the absolute pressure of the sample chamber after vacuumizing, and MPa; v s represents the sample chamber volume, cm 3;Vg represents the sample skeleton volume, cm 3;Z1 represents the gas compression factor (i.e. the ratio of the molar volume of gas to the molar volume of ideal gas) at a pressure of p 1, and p 2 represents the initial absolute pressure of the reference chamber, MPa; v r denotes the reference chamber volume, cm 3;Z2 denotes the gas compression factor at p 2 pressure; p 3 represents the absolute pressure after equilibration, mpa; z 3 represents the gas compression factor at p 3 pressure;
obtaining total volume and porosity of the sample according to a caliper measurement method or an Archimedes immersion method The calculation is as follows:
In the middle of Representing the porosity of the rock sample; v z represents the total volume of the sample, cm 3.
Compared with the conventional porosity measurement method (plunger-type helium method, particle-type helium method and liquid saturation method) of the shale, the method provided by the invention has the advantages that the porosity of the shale is represented by adopting the plunger-type helium method which aims at the shale rich in organic matters and has the highest accuracy and reliability, and the age of 4 He is corrected by utilizing the formula (3), so that the in-situ yield of the shell source 4 He can be accurately calculated.
The invention is not exhaustive and can be seen in the prior art.
The beneficial effects of the invention are as follows:
According to the invention, on the basis that the former calculates the in-situ release of 4 He through a U, th radioactive decay formula, the 4 He age and the porosity parameter in the formula are corrected, so that compared with the former calculated in-situ helium resource amount, the method is more accurate and scientific, and has better convincing ability.
Drawings
The accompanying drawings, which are included to provide a further understanding of the application and are incorporated in and constitute a part of this specification, illustrate embodiments of the application and together with the description serve to explain the application.
FIG. 1 is a flow chart of an in situ helium resource quantity evaluation method of the present invention that takes into account both porosity and helium generation "age" correction.
Detailed Description
In order to better understand the technical solutions in the present specification, the following description will clearly and completely describe the technical solutions in the embodiments of the present invention in conjunction with the drawings in the implementation of the present specification, but not limited thereto, and the present invention is not fully described and is according to the conventional technology in the art.
Example 1
An in-situ helium resource quantity evaluation method which takes into account both porosity and raw helium 'age' correction, as shown in fig. 1, comprises the following steps:
S1: and (3) defining a target producing interval, selecting a typical or key natural gas (rare gas) exploration horizon, determining the volume of source rock through the thickness and the area of the source rock, collecting a deposited helium source rock sample, performing high-pressure mercury testing to obtain the rock density, and obtaining the quality of the source rock according to the volume of the source rock and the rock density.
S2: basic data are acquired, and the reservoir temperature, the reservoir pressure and the solubility coefficient of helium in water of a target interval sample are clearly researched;
reservoir temperature may be replaced with measured downhole temperature or bit temperature; reservoir pressure can be obtained by logging, formation pressure testing or mathematical model and the like, and can be replaced by formation pressure; the solubility coefficient of helium in water was calculated by PHREEQC software in combination with the reservoir temperature, reservoir pressure of the sample.
S3: sampling at the wellhead of a typical well in a research area or sampling while drilling during logging, sealing to prevent helium leakage, and performing separation test to obtain the He content c in natural gas samples with different depths;
S4: u, th the half-life period of the main natural isotope is extremely long, so that the U, th concentration in the stratum nowadays can be used for approximately replacing the average abundance of the initial U, th deposit, and the U and Th contents can be obtained through geochemical analysis of the drilling core of a classical well;
s5: the in situ 4 He yield was calculated from equation (1):
Wherein ρ r is rock density in g/cm -3; t is the hiding time, i.e. 4 He age, per unit year; Is porosity; Λ=1, which is the efficiency of He transfer from the rock matrix to water; cm 3STPg-1 H2O is 4 He yield per gram of water at standard conditions; p (4 He) is the in situ 4 He yield obtained from equation (2):
P(4He)＝1.207×10-13[U]+2.867×10-13[Th] cm3STPg-1 rockyr-1 (2)
Wherein [ U ] and [ Th ] respectively represent average abundance of U and Th in a research area, and unit ppm is obtained in the step S4; substituting the average abundance of U and Th to obtain the in situ 4 He yield P (4He),cm3STPg-1 rockyr-1 is 4 He yield per gram of rock per year in standard state;
S6: correcting 4 He age T in the formula (1) to obtain corrected 4 He age T
In combination with the partial pressure in the study area, the solubility coefficient of helium in water S He and the annual helium content α in water entering the 1cm 3 pore from rock, the corrected 4 He "age" T is calculated by equation (2):
wherein c is the He content obtained in the step S3,%; p is the reservoir pressure obtained in step S2, atm; s He is the solubility coefficient of helium in water, cm 3/cm3. Atm; alpha is a fixed value, and 3.4X10 -12a-1 is taken; cp is the partial pressure of helium in the gas layer;
s7: porosity calculation using plunger-like helium method Firstly, placing a sample into a sample chamber, vacuumizing the sample chamber to a fixed value p 1, taking p 1 to 0.1Pa in the embodiment, injecting helium into a reference chamber to 1.378MPa, recording the gas pressure p 2 in the reference chamber when the helium pressure in the reference chamber is balanced, then communicating the reference chamber with the sample chamber to fully saturate the sample pore with the helium, recording the gas pressure p 3 in the sample chamber after the helium pressure in the sample chamber is balanced, and obtaining a formula (4) according to Boyle's law; under the conditions that the testing process is at constant temperature and the valve displacement volume is ignored, the shale sample skeleton volume Vg is calculated by the simultaneous formula (4) and the formula (5):
wherein, p 1 represents the absolute pressure of the sample chamber after vacuumizing, and MPa; v s represents the sample chamber volume, cm 3;Vg represents the sample skeleton volume, cm 3;Z1 represents the gas compression factor (i.e. the ratio of the molar volume of gas to the molar volume of ideal gas) at a pressure of p 1, and p 2 represents the initial absolute pressure of the reference chamber, MPa; v r denotes the reference chamber volume, cm 3;Z2 denotes the gas compression factor at p 2 pressure; p 3 represents the absolute pressure after equilibration, mpa; z 3 represents the gas compression factor at p 3 pressure;
obtaining total volume and porosity of the sample according to a caliper measurement method or an Archimedes immersion method The calculation is as follows:
In the middle of Representing the porosity of the rock sample; v z represents the total volume of the sample, cm 3.
Substituting the porosity, corrected 4 He 'age' T, in-situ 4 He yield P (4 He), transfer efficiency of He from rock matrix to water Λ (Λ=1) and rock density into formula (1) to calculate in-situ 4 He yield in the research area, and multiplying the obtained in-situ 4 He yield by source rock mass to obtain in-situ helium resource amount in the research area.
Example 2
An in-situ helium resource amount evaluation method taking both porosity and helium generation 'age' correction into consideration selects a helium-containing/helium-rich natural gas field A of an oil-containing basin in a certain area for verification:
(1) Because the helium resource quantity is calculated on the basis of no data correction of the A gas field, the original data of the A gas field are obtained as follows: the rock density was 2.55g/cm 3, the mass was 149.71 ×10 12 g, the average abundance of U and Th was 25ppm and 20ppm, the formation age was 442Ma, and the porosity was 4%. Substituting the average abundance of U and Th into equation (2) can calculate the in situ 4 He yield P (4 He) to be 8.75×10 -12cm3 STP grock -1yr-1. Using the above data, substituting into equation (1) yields in situ 4 He of 2.37X10 -1cm3STP g-1 H2O were obtained. Finally, the yield of in-situ 4 He multiplied by the mass can obtain the in-situ helium resource amount of 3.54 multiplied by 10 10m3.
(2) Correction of 4 He age by the raw helium "age" formula, the He content c of 0.05%, the helium pressure in the gas layer of 400atm (reservoir pressure 40 Mpa), and the solubility coefficient of helium in water of 0.0063cm 3/cm3 atm, the helium content in water entering 1cm 3 pores from rock per year being constant, i.e. 3.4x10 -12a-1, were obtained by the above step S3. Substituting into equation (3) yields helium "age" of 371Ma.
(3) And correcting the porosity of the rock by a plunger-type helium method, taking a plurality of samples in the A gas field, measuring the porosity of the rock by the plunger-type helium method in the step S7, and taking the average value of the porosity of the plurality of samples as 6%.
(4) The corrected 4 He age and porosity parameters were re-substituted into equation (1) to give an in situ 4 He yield of 1.40 x 10 -1cm3 STP g-1 H2O. Also multiplying the in situ 4 He yield by the source rock mass gives an in situ helium source quantity of 1.94 x 10 10m3.
TABLE 1 results of in situ 4 He yield calculations and tables of parameters
It can be seen that the in-situ helium gas resource amount calculated from the corrected data is smaller than that before the correction because a great amount of escape of helium gas occurs in the process of generation and storage, and the result of calculation before the correction is only the maximum generation amount of helium gas and is not the resource amount of helium gas.
While the foregoing is directed to the preferred embodiments of the present invention, it will be appreciated by those skilled in the art that various modifications and adaptations can be made without departing from the principles of the present invention, and such modifications and adaptations are intended to be comprehended within the scope of the present invention.
Claims (4)
1. An in-situ helium resource amount evaluation method taking both porosity and helium generation 'age' correction into consideration is characterized by comprising the following steps:
S1: defining a target producing interval, determining the helium rock type of the producing interval, and determining the volume and the quality of source rock;
S2: basic data are acquired, and the solubility coefficient of reservoir temperature, reservoir pressure and helium in water is defined;
S3: sampling at the wellhead of a typical well in a research area or sampling while drilling during logging, sealing, and performing separation test to obtain the He content c in the natural gas of samples with different depths through test analysis;
s4: the average abundance of the deposit initiation U, th was approximated with U, th concentrations in the formation today;
s5: the in situ 4 He yield was calculated from equation (1):
Wherein ρ r is rock density in g/cm -3; t is the hiding time, i.e. 4 He age, per unit year; is porosity; Λ=1, which is the efficiency of He transfer from the rock matrix to water;
P (4 He) is the in situ 4 He yield obtained from equation (2):
P(4He)＝1.207×10-13[U]+2.867×10-13[Th] cm3STPg-1 rockyr-1 (2)
Wherein [ U ] and [ Th ] respectively represent average abundance of U and Th in a research area, and unit ppm is obtained in the step S4;
S6: correcting 4 He age T in the formula (1) to obtain corrected 4 He age T
In combination with the partial pressure in the study area, the solubility coefficient of helium in water S He and the annual helium content α in water entering the 1cm 3 pore from rock, the corrected 4 He "age" T is calculated by equation (2):
wherein c is the He content obtained in the step S3,%; p is the reservoir pressure obtained in step S2, atm; s He is the solubility coefficient of helium in water, cm 3/cm3. Atm; alpha is a fixed value, and 3.4X10 -12a-1 is taken;
s7: porosity calculation using plunger-like helium method Substituting the porosity and corrected 4 He 'age' T into a formula (1) to calculate the yield of in-situ 4 He in the research area, and multiplying the obtained yield of in-situ 4 He by the quality of source rock to obtain the in-situ helium resource amount in the research area.
2. The in-situ helium resource quantity evaluation method taking into account porosity and helium generation "age" correction according to claim 1, wherein in step S1, a source rock volume is determined by a source rock thickness and an area;
and collecting a deposited helium source rock sample, performing high-pressure mercury testing to obtain rock density, and obtaining the quality of the source rock according to the volume of the source rock and the rock density.
3. The in situ helium resource quantity evaluation method taking into account both porosity and helium generation "age" correction according to claim 1, wherein in step S2, reservoir temperature is replaced with measured downhole temperature or bit temperature; reservoir pressure is replaced with formation pressure; the solubility coefficient of helium in water was calculated by PHREEQC software in combination with the reservoir temperature, reservoir pressure of the sample.
4. The method for evaluating the in-situ helium resource amount taking into account both porosity and helium generation "age" correction according to claim 1, wherein in step S7, the process of calculating the porosity by using a plunger-like helium method is as follows:
Firstly, placing a sample into a sample chamber, vacuumizing the sample chamber to a fixed value p 1, injecting helium into a reference chamber to a certain pressure, recording the gas pressure p 2 in the reference chamber when the helium pressure in the reference chamber is balanced, then communicating the reference chamber with the sample chamber to fully saturate the sample pore with helium, recording the gas pressure p 3 in the sample chamber when the helium pressure in the sample chamber is balanced, and obtaining a formula (4) according to Boyle's law; under the conditions that the testing process is at constant temperature and the valve displacement volume is ignored, the shale sample skeleton volume Vg is calculated by the simultaneous formula (4) and the formula (5):
Wherein, p 1 represents the absolute pressure of the sample chamber after vacuumizing, and MPa; v s represents the sample chamber volume, cm 3;Vg represents the sample skeleton volume, cm 3;Z1 represents the gas compression factor under p 1 pressure, and p 2 represents the initial absolute pressure of the reference chamber, MPa; v r denotes the reference chamber volume, cm 3;Z2 denotes the gas compression factor at p 2 pressure; p 3 represents the absolute pressure after equilibration, mpa; z 3 represents the gas compression factor at p 3 pressure;
obtaining total volume and porosity of the sample according to a caliper measurement method or an Archimedes immersion method The calculation is as follows:
In the middle of Representing the porosity of the rock sample; v z represents the total volume of the sample, cm 3.
Applications Claiming Priority (2)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202410529689 | 2024-04-29 | ||
| CN2024105296890 | 2024-04-29 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN118471358A true CN118471358A (en) | 2024-08-09 |
Family
ID=92168451
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202410634435.5A Pending CN118471358A (en) | 2024-04-29 | 2024-05-22 | In-situ helium resource quantity evaluation method considering porosity and helium generation age correction |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN118471358A (en) |
Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130091925A1 (en) * | 2011-10-12 | 2013-04-18 | Thomas Henry Darrah | Method for Determining the Genetic Fingerpring of Hydrocarbons and Other Geological Fluids using Noble Gas Geochemistry |
| CN110687006A (en) * | 2019-09-30 | 2020-01-14 | 苏州冠德能源科技有限公司 | Rock gas content calculation method based on well site analytic experiment |
| CN116047624A (en) * | 2022-12-29 | 2023-05-02 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | Helium resource amount evaluation method |
-
2024
- 2024-05-22 CN CN202410634435.5A patent/CN118471358A/en active Pending
Patent Citations (3)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| US20130091925A1 (en) * | 2011-10-12 | 2013-04-18 | Thomas Henry Darrah | Method for Determining the Genetic Fingerpring of Hydrocarbons and Other Geological Fluids using Noble Gas Geochemistry |
| CN110687006A (en) * | 2019-09-30 | 2020-01-14 | 苏州冠德能源科技有限公司 | Rock gas content calculation method based on well site analytic experiment |
| CN116047624A (en) * | 2022-12-29 | 2023-05-02 | 中国地质调查局西安地质调查中心(西北地质科技创新中心) | Helium resource amount evaluation method |
Non-Patent Citations (1)
| Title |
|---|
| 王杰;刘文汇;陶成;腾格尔;席斌斌;王萍;杨华敏;: "海相油气成藏定年技术及其对元坝气田长兴组天然气成藏年代的反演", 地球科学, no. 06, 24 May 2018 (2018-05-24) * |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| Etheridge et al. | Natural and anthropogenic changes in atmospheric CO2 over the last 1000 years from air in Antarctic ice and firn | |
| Aeschbach-Hertig et al. | A paleotemperature record derived from dissolved noble gases in groundwater of the Aquia Aquifer (Maryland, USA) | |
| Lippmann et al. | Dating ultra-deep mine waters with noble gases and 36Cl, Witwatersrand Basin, South Africa | |
| Ballentine et al. | A Magnus opus: Helium, neon, and argon isotopes in a North Sea oilfield | |
| Hu et al. | Hydrocarbon generation and expulsion characteristics of Lower Permian P1f source rocks in the Fengcheng area, northwest margin, Junggar Basin, NW China: implications for tight oil accumulation potential assessment | |
| Moore et al. | Differentiating between biogenic and thermogenic sources of natural gas in coalbed methane reservoirs from the Illinois Basin using noble gas and hydrocarbon geochemistry | |
| EP2116690A1 (en) | Geochemical surveillance of gas production from tight gas fields | |
| Zhang et al. | Quantifying the helium and hydrocarbon accumulation processes using noble gases in the North Qaidam Basin, China | |
| Wen et al. | Assessing compositional variability and migration of natural gas in the Antrim Shale in the Michigan Basin using noble gas geochemistry | |
| Palcsu et al. | In-reservoir mixing of mantle-derived CO2 and metasedimentary CH4–N2 fluids–Noble gas and stable isotope study of two multistacked fields (Pannonian Basin System, W-Hungary) | |
| CN110687006B (en) | Rock gas content calculation method based on well site analytic experiment | |
| CN111060672B (en) | Method for reproducing whole history process of forming high-temperature overpressure natural gas reservoir | |
| Mayer et al. | A multi-tracer study of groundwater origin and transit-time in the aquifers of the Venice region (Italy) | |
| Li et al. | Tracing interaction between hydrocarbon and groundwater systems with isotope signatures preserved in the Anyue gas field, central Sichuan Basin, China | |
| Stalker et al. | Application of tracers to measure, monitor and verify breakthrough of sequestered CO2 at the CO2CRC Otway Project, Victoria, Australia | |
| CN113431537B (en) | Unsteady variable-flow-rate large-scale rock core water flooding gas relative permeability testing method | |
| Pang et al. | Hydrocarbon accumulation depth limit and implications for potential resources prediction | |
| Huo et al. | Geochemical characteristics and hydrocarbon expulsion of source rocks in the first member of the Qingshankou Formation in the Qijia‐Gulong Sag, Songliao Basin, Northeast China: Evaluation of shale oil resource potential | |
| CN110579810B (en) | Rare gas4Method for estimating shale gas reservoir age by He | |
| Serno et al. | Using oxygen isotopes to quantitatively assess residual CO2 saturation during the CO2CRC Otway Stage 2B Extension residual saturation test | |
| CN115356463B (en) | Method for judging whether substrate granite supplies helium to helium-rich natural gas reservoir | |
| Botor et al. | Petroleum generation in the Carboniferous strata of the Lublin Trough (Poland): an integrated geochemical and numerical modelling approach | |
| CN118471358A (en) | In-situ helium resource quantity evaluation method considering porosity and helium generation age correction | |
| Smith et al. | Testing helium equilibrium between quartz and pore water as a method to determine pore water helium concentrations | |
| Stock et al. | The Posidonia Shale of northern Germany: unconventional oil and gas potential from high-resolution 3D numerical basin modelling of the cross-junction between the eastern Lower Saxony Basin, Pompeckj Block and Gifhorn Trough |
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 |