Method for calculating natural hydrogen resource in sedimentary basin and application thereof

Abstract

The application provides a kind of sedimentary basin natural hydrogen resource quantity calculation method and its application, it is related to the technical field of resource exploration, including: by determining the total hydrogen production rate of target area, capture efficiency, hydrogen loss rate and hydrogen residence time, the geological resource quantity of natural hydrogen in target area is calculated by using material balance method;Wherein, the total hydrogen production rate includes the hydrogen generation rate from at least one source of serpentine reaction, radioactive radiation effect and deep mantle degassing;Capture efficiency is the proportion of effectively stored hydrogen gas in generated hydrogen gas;Hydrogen loss rate is the annual loss amount of hydrogen gas in underground system due to physical, chemical and biological loss mechanism;Hydrogen residence time is the time for hydrogen to stay and be saved in geological body.The application considers the key parameters such as total hydrogen production rate, capture efficiency, hydrogen loss rate and residence time, and first constructs a natural hydrogen geological resource quantity calculation model suitable for sedimentary basin scale.

Description

TECHNICAL FIELD

[0001] The present application relates to the technical field of resource exploration, and in particular to a method for calculating the amount of natural hydrogen resources in a sedimentary basin and application thereof. BACKGROUND

[0002] The calculation of resource amount is the basic basis for resource potential and evaluation of selected areas, development deployment and resource management in resource exploration and development, and is crucial to resource exploration and development. As a highly efficient and clean non-fossil energy, the development and utilization of natural hydrogen has attracted widespread attention.

[0003] However, the current exploration and research of natural hydrogen resources are still in the initial stage and face many challenges and risks. Among them, the lack of a method for calculating the amount of hydrogen resources is one of the important challenges restricting the discovery of natural hydrogen exploration, which is aimed at the main exploration and development target, i.e. a natural hydrogen system (an organic whole consisting of an effective hydrogen source body as the core, a reservoir cap rock and a migration system) and a sedimentary basin (composed of one or several natural hydrogen systems). Due to the particularity of the genetic type and physical and chemical properties of natural hydrogen, a hydrogen system is not equal to an oil and gas system, so the oil and gas resource amount calculation method is not applicable to natural hydrogen.

[0004] At present, foreign scholars have established a numerical model for predicting global natural hydrogen resource potential based on material balance method, which only provides a preliminary framework for evaluating global natural hydrogen resource potential, and the related parameters are all assumed, resulting in high uncertainty in resource potential evaluation, and the resource potential evaluation span is up to seven orders of magnitude. The root cause lies in the regional differences of natural hydrogen systems, and there are many natural hydrogen systems in different sedimentary basins and their periphery, and their formation and distribution differ significantly. This model cannot be applied to the calculation of resource amount and potential evaluation of regional or basin scale and natural hydrogen systems. Therefore, establishing a method for calculating the amount of natural hydrogen resources suitable for the exploration and development of sedimentary basins and their hydrogen systems is a key problem to be solved at present.

[0005] In view of this, the present application is proposed. SUMMARY

[0006] One of the purposes of the present application is to provide a method for calculating the amount of natural hydrogen resources in a sedimentary basin to at least solve one of the technical problems existing in the prior art.

[0007] The second purpose of the present application is to provide a method for calculating the amount of natural hydrogen resources in a sedimentary basin in the calculation of the amount of natural hydrogen resources in a sedimentary basin.

[0008] In order to achieve the above purposes of the present application, the following technical solutions are adopted:

[0009] In the first aspect, the present application provides a method for calculating the amount of natural hydrogen resources in a sedimentary basin, comprising:

[0010] Based on the mass dynamic balance relationship of the generation, migration, accumulation and loss of the underground abiogenic hydrogen, the geological resource amount of the natural hydrogen in the target region is calculated by determining the total hydrogen production rate, capture efficiency, hydrogen loss rate and hydrogen residence time of the target region through the material balance method;

[0011] The total hydrogen production rate includes the hydrogen production rate from at least one source of serpentinization reaction, radioactive radiation and deep mantle degassing; the capture efficiency is the proportion of the generated hydrogen effectively stored by geological trapping; the hydrogen loss rate is the annual loss amount of hydrogen in the underground system due to physical, chemical and biological loss mechanisms; and the hydrogen residence time is the time for hydrogen to stay and be preserved in the geological body.

[0012] Further, the method comprises the following steps:

[0013] (a) Based on the geological conditions of the sedimentary basin and the surrounding area, the type of the natural hydrogen system and its accumulation geological elements are identified;

[0014] (b) According to the generation mechanism of the underground abiogenic hydrogen, the hydrogen production rate of at least one source of serpentinization, radioactive radiation and deep mantle degassing is determined to obtain the total hydrogen production rate;

[0015] (c) Based on the accumulation ability and preservation conditions of hydrogen during migration, the capture efficiency of hydrogen captured into effective traps is determined;

[0016] (d) The diffusion, secondary conversion and microbial consumption process of hydrogen in the near-surface environment are analyzed to determine the loss rate of the hydrogen system;

[0017] (e) The residence time is determined by at least one of the following methods: geological events, inclusion analysis and isotope dating;

[0018] (f) The hydrogen resource amount of a single hydrogen system is calculated, and the calculation formula is: Q R ( G × η trap - L )× t res ;

[0019] wherein, Q R is the hydrogen resource amount; G is the total hydrogen production rate; η trap is the capture efficiency; L is the hydrogen loss rate; t res is the hydrogen residence time;

[0020] (g) adding the hydrogen resource amounts of the several hydrogen systems in the target region to obtain the total natural hydrogen resource amount of the target region.

[0021] Further, the hydrogen production rate from serpentinization is the serpentinization hydrogen production rate, and the calculation formula of the serpentinization hydrogen production rate is: G serp =S serp ×H serp ×F H2 ;

[0022] wherein, G serp is the serpentinization hydrogen production rate; S serp is the area of the iron-rich hydrogen source body; H serp is the thickness of the iron-rich hydrogen source body; F H2 is the hydrogen production flux.

[0023] Further, the hydrogen production rate from radioactive decomposition is the radioactive decomposition hydrogen production rate, and the calculation formula of the radioactive decomposition hydrogen production rate is: G rad =S rad ×H rad ×Φ rad ×ω U / Th ×F rad ;

[0024] wherein, G rad is the radioactive decomposition hydrogen production rate; S rad is the area of the uranium-rich hydrogen source body; H rad is the thickness of the uranium-rich hydrogen source body; Φ rad is the water-containing porosity; ω U / Th is the uranium / thorium content; F rad is the hydrogen production flux.

[0025] Further, the hydrogen production rate from deep mantle degassing is the mantle-derived gas annual contribution amount, and the mantle-derived gas annual contribution amount is obtained by an isotope mixing model; the isotope mixing model includes at least one of a hydrogen isotope mixing model and a helium isotope mixing model.

[0026] Further, the calculation formula of the hydrogen isotope mixing model is: G mant = (δD crust - δD mixe ) / (δD crust - δD mant ) × 100% × f surf ;

[0027] wherein, G mant is the mantle-derived gas annual contribution amount; δD crust is a crust source δD endmember value; δD mant is a mantle source δD endmember value; δD mixeThe measured value of the mixed gas is δD; f surf The natural hydrogen surface flux.

[0028] Further, the calculation formula of the helium isotope mixing model is G mant = (R crust -R mixe ) / (R crust -R mant ) x r H ;

[0029] Wherein, G mant is the annual contribution of the mantle source gas; R mant is the mantle source 3 He 4 end member value; R crust is the crust source 3 He 4 end member value; R mixe is the measured value of the mixed gas 3 He 4 ; r H is the mantle hydrogen H2 3 / He ratio.

[0030] Further, the annual contribution of the mantle source gas is constrained and corrected according to the tectonic environment and the associated mantle source gas ratio.

[0031] Further, the tectonic environment includes at least one of a rift, a subduction zone and a craton; and the associated mantle source gas ratio includes CO2 3 / He ratio.

[0032] In a second aspect, the present application provides an application of a sedimentary basin natural hydrogen resource quantity calculation method in calculating the sedimentary basin natural hydrogen resource quantity.

[0033] Compared with the prior art, the present application has the following beneficial effects:

[0034] The sedimentary basin natural hydrogen resource quantity calculation method provided by the present application is based on the mass dynamic balance relationship of abiogenic hydrogen generation-migration-aggregation-loss, and by comprehensively considering key parameters such as total hydrogen production rate, capture efficiency, hydrogen loss rate and residence time, a natural hydrogen geological resource quantity calculation model suitable for the scale of a sedimentary basin is first constructed. The method breaks through the limitation that the traditional oil and gas resource quantity estimation method has poor applicability in the hydrogen system, can accurately reflect the characteristics of multi-source genesis, high diffusivity and dynamic preservation of hydrogen, significantly improves the scientificity and reliability of resource quantity calculation, and provides effective technical support for natural hydrogen resource potential evaluation and exploration selection. BRIEF DESCRIPTION OF DRAWINGS

[0035] In order to more clearly illustrate the technical solutions in the specific embodiments or prior art of the present application, the drawings needed to be used in the description of the specific embodiments or prior art will be briefly introduced below. Obviously, the drawings in the following description are some embodiments of the present application, and other drawings can be obtained by those skilled in the art without any creative effort on the basis of these drawings.

[0036] Figure 1 The sedimentary basin natural hydrogen resource quantity calculation method and system provided by the present application are shown in the schematic diagram.

[0037] Figure 2 The natural hydrogen content distribution map of the northern Songliao Basin is shown in the figure.

[0038] Figure 3 The natural gas component characteristic map of different oil and gas layers in the northern Songliao Basin is shown in the figure.

[0039] Figure 4 The typical characteristic table of natural gas components with hydrogen content greater than 10% in the deep and lower oil and gas containing combination of the Songliao Basin is shown in the table. DETAILED DESCRIPTION

[0040] Unless otherwise defined herein, scientific and technical terms used in connection with the present application shall have the meanings that are commonly understood by those of ordinary skill in the art. Clear indications to the contrary are necessary to overcome the common meaning of the term when it appears alone or can be defined in accordance with its conventional meaning. In this application, the use of "or" means "and / or" unless otherwise stated. In addition, the use of the term "including" as well as other forms such as "include", "includes" for claiming consists of, or "comprise", "comprises" and / or "comprising", are non-limiting.

[0041] The technical solutions of the present application will be described below in conjunction with the embodiments. Obviously, the described embodiments are only a part of the embodiments of the present application, rather than all the embodiments. Based on the embodiments in the present application, all other embodiments obtained by those skilled in the art without any creative effort belong to the scope of protection of the present application.

[0042] The first aspect of the present application provides a sedimentary basin natural hydrogen resource quantity calculation method, comprising:

[0043] Based on the mass dynamic balance relationship of the generation, migration, accumulation and loss of underground non-biogenic hydrogen, the geological resource quantity of natural hydrogen in the target area is calculated by determining the total hydrogen production rate, capture efficiency, hydrogen loss rate and hydrogen residence time of the target area using material balance method.

[0044] The total hydrogen production rate includes hydrogen production rates from at least one source of serpentinization reaction, radioactive radiation and deep mantle degassing; the capture efficiency is the proportion of hydrogen effectively stored by geological trapping; the hydrogen loss rate is the annual loss amount of hydrogen in the underground system due to physical, chemical and biological loss mechanisms; and the hydrogen residence time is the time for hydrogen to stay and be preserved in the geological body.

[0045] The application discloses a natural hydrogen resource quantity calculation method applied to a sedimentary basin and a hydrogen system, and the method comprises the following steps: based on the revelation of the natural hydrogen accumulation mechanism of the sedimentary basin and the periphery, the hydrogen system and its accumulation geological elements are determined; based on the natural hydrogen generation-migration-accumulation-diffusion mass dynamic balance system of the underground non-biological causes (serpentinization, radiation, deep degassing), the hydrogen geological resource quantity calculation method of the natural hydrogen system is established, with "the hydrogen generation quantity in a given volume and time being equal to the sum of the accumulation quantity, the surface flux and the loss quantity" as the core; the hydrogen production rate of the hydrogen system of different sources is calculated through hydrogen source evaluation; the consumption rate and the capture efficiency are determined by using biological and non-biological consumption quantitative characterization and surface flux measurement; the residence time is determined by using the geological event and the isotopic dating method; the natural hydrogen resource quantity of different hydrogen systems is calculated according to the hydrogen generation rate of different sources; the total natural hydrogen resource quantity of the sedimentary basin is calculated based on the hydrogen resource quantity of different sources; and the natural hydrogen resource quantity calculation method of the sedimentary basin is constructed based on the calculation method and the parameter system of the generation quantity, the accumulation quantity and the consumption quantity of different source hydrogen systems. The application creates the hydrogen resource quantity calculation method of the natural hydrogen system, and provides a basic method and reliable basis for quantitatively evaluating the natural hydrogen resource potential and strategically selecting a region on the basin or regional scale.

[0046] In some preferred embodiments, the sedimentary basin natural hydrogen resource quantity calculation method comprises the following steps:

[0047] (1) Based on the revelation of the natural hydrogen accumulation mechanism of the sedimentary basin and the periphery, the hydrogen system and its accumulation geological elements are determined.

[0048] (2) Based on the natural hydrogen generation-migration-accumulation-diffusion mass dynamic balance system of the underground non-biological causes (serpentinization, radiation, deep degassing), the hydrogen geological resource quantity (In-situ geological resources) calculation method of the natural hydrogen system is established, with "the hydrogen generation quantity in a given volume and time being equal to the sum of the accumulation quantity, the surface flux and the loss quantity" as the core.

[0049] Core formula (mass balance method): Q R ( G × η trap - L ) ×t res ;

[0050] wherein, Q R is the hydrogen resource amount (also referred to as the amount of hydrogen resources in traps), in kg or m 3 ; G is the total hydrogen production rate (also referred to as the annual hydrogen production), in kg / yr or m 3 / yr, which is divided into G serp , G radi , and G mant ; η trap is the capture efficiency of traps (0-1, dimensionless); L is the hydrogen loss rate (kg / yr or m 3 / yr, loss rate); t res is the hydrogen residence time (yr, residence time).

[0051] (3) The hydrogen production rate of different source hydrogen systems is calculated through hydrogen source evaluation; the natural hydrogen genesis type, hydrogen system, and Figure 1 related parameters in the sedimentary basin are identified based on the content of (1) above, and the quality, hydrogen production rate, and scale of each type of hydrogen source body are evaluated.

[0052] wherein, the hydrogen production rate of the serpentinization source is the serpentinization hydrogen production rate (also referred to as the serpentinization hydrogen production amount), and the calculation formula of the serpentinization hydrogen production rate is: serp =G serp ×H serp ×F H2 ;

[0053] G serp is the serpentinization hydrogen production rate; S serp is the area of the iron-rich hydrogen source body (m 2 ), which is obtained through geological survey; H serp is the thickness of the iron-rich hydrogen source body (m), which is obtained through geological survey; and F H2 is the hydrogen production flux (kg H2 / m 3 ·yr), which is obtained through hydrogen production simulation experiment, reaction kinetics calculation, and / or analogy method.

[0054] wherein the hydrogen production rate from radioactive decay is a radioactive decay hydrogen production rate, and the radioactive decay hydrogen production rate is calculated according to the formula: G rad = S rad × H rad × Φ rad × ω U / Th × F rad ;

[0055] G rad is the radioactive decay hydrogen production rate; S rad is the area of the uranium-rich hydrogen source body (m 2 ), which is obtained by geological surveying; H rad is the thickness of the uranium-rich hydrogen source body (m), which is obtained by geological surveying; Φ rad is the water-bearing porosity (%), which is obtained by experimental testing; ω U / Th is the uranium / thorium content (ppm), which is obtained by experimental testing; and F rad is the hydrogen production flux (kg H2 / m 3 ·yr), which is obtained by radioactive hydrogen production simulation experiments and / or analogy methods (theoretical calculations).

[0056] wherein the hydrogen production rate from deep mantle degassing is a mantle-derived gas annual contribution (which can also be referred to as a mantle-derived hydrogen gas amount), and the mantle-derived gas annual contribution is obtained by an isotope mixing model; the mantle-derived hydrogen gas amount G mant The calculation is divided into two models, including at least one of a hydrogen isotope mixing model and a helium isotope mixing model.

[0057] In some preferred embodiments, the calculation formula of the hydrogen isotope (δD) mixing model is: G mant = (δD crust - δD mixe ) / (δD crust - δD mant ) × 100% × f surf ;

[0058] wherein G mant is the mantle-derived gas annual contribution (kg / yr or m 3 / yr); δD crust is a crust-derived δD end-member value (‰), which is obtained by a theoretical or empirical value; δD mant is a mantle-derived δD end-member value (‰), which is obtained by a theoretical or empirical value; δD mixe is a measured δD value of the mixed gas (‰), which is obtained by actual measurement; and f surf is a natural hydrogen gas surface flux (kg / yr or m 3 / yr), which is obtained by actual measurement (long-term continuous monitoring).

[0059] In some preferred embodiments, the helium isotope ( 3 He / 4 The calculation formula for the He R) mixture model is: G mant =(R crust -R mixe ) / (R crust -R mant )×r H ;

[0060] Among them, G mant Annual contribution of mantle-derived gas (kg / yr or m³) 3 / yr); R mant mantle source 3 He / 4 He endmember value (R / Ra=6-9, atmospheric) 3 He / 4 He ratio Ra = 1.4 × 10 -6 ), obtained through theory or analogy; R crust shell source 3 He / 4 He endmember values ​​(R / Ra = 0.01-0.05) are obtained through theory or analogy; R mixe For mixed gas 3 He / 4 He is the measured value, obtained through actual measurement; r H Mantle-derived hydrogen (H2) 3 The He ratio is obtained through theory or analogy.

[0061] Based on the tectonic environment (rift valley / subduction zone / craton) and associated mantle-derived gases such as CO2 / 3 The He ratio constrains and corrects the calculated amount of mantle-sourced hydrogen.

[0062] (4) Hydrogen loss rate L The hydrogen loss rate in the upper part of the hydrogen system within the sedimentary basin was obtained by comprehensively analyzing loss mechanisms such as diffusion, secondary transformation, and microbial consumption, combined with migration, diffusion, and leakage simulation studies under actual geological conditions, and by comparing and analyzing known hydrogen systems.

[0063] (5) Capture efficiency η trap The capture efficiency of hydrogen systems in potential reservoirs within sedimentary basins is obtained through a comprehensive analysis of the generation, storage and preservation conditions of natural hydrogen systems, their occurrence state and migration processes, and / or by analogy with similar known natural hydrogen systems.

[0064] (6) Duration of stay t resThe determination of the retention time utilizes the combined method of geological events, inclusion analysis and isotopic dating. The determination of the serpentinization-derived hydrogen gas focuses on the determination of the geological time (years) of the start and end of the serpentinization reaction and the duration. The determination of the radioactive decay hydrogen gas focuses on the determination of the age of the uranium-rich hydrogen source body and the formation time of the potential trap. The determination of the mantle-derived hydrogen gas focuses on the determination of the activity time of the deep fault zone.

[0065] (7) According to the hydrogen generation rates of different sources, the natural hydrogen gas resource amounts of different hydrogen systems in the sedimentary basin are calculated, and the calculation formula is:

[0066] Q R1 (G serp G radi G mant ) η trap L ) × t res

[0067] Q R2 (G serp G radi G mant ) η trap L ) × t res

[0068] According to the same method, the resource amounts Q R3 to Q R of the remaining hydrogen systems can be calculated in sequence; wherein, Q Rn (G serp G radi G mant ) η trap L ) × t res

[0069] (8) The resource amounts of multiple hydrogen systems in the target area are added to obtain the total natural hydrogen gas resource amount of the target area.

[0070] ​​​​​​​​​​​​​​​​​​Specifically, the total resource amount of natural hydrogen in the sedimentary basin is obtained by accumulation of the hydrogen system resource amount, and the calculation formula is:

[0071] Q R Q R1 Q R2 Q Rn

[0072] (9) Based on the calculation method and parameter system of the hydrogen system generation amount, accumulation amount and consumption amount of different sources, the calculation method of the natural hydrogen resource amount in the sedimentary basin is formed.

[0073] The third aspect of the present application provides an application of the calculation method of the natural hydrogen resource amount in the sedimentary basin in calculating the natural hydrogen resource amount in the sedimentary basin.

[0074] The present application will be further described by examples.

[0075] The present embodiment provides a calculation method of the natural hydrogen resource in the sedimentary basin, and calculates the natural hydrogen resource amount in the Songliao Basin, comprising the following steps:

[0076] (1) First, the basic geological parameters of the Songliao Basin are obtained by data collection. The total area of the Songliao Basin is about 26×104km 2 , which has a double-layer structure of lower subsidence and upper depression, and the base of the basin is mainly shallow metamorphic granite of Carboniferous-Permian system. Figure 2 The deep gas-bearing combination mainly includes four sets of gas layers of basement, Shahezi Formation, Yingcheng Formation and Denglouku Formation; the lower oil and gas-bearing combination includes Fuyu oil layer and Yangdachengzi oil layer; the middle oil and gas-bearing combination includes Saertu oil layer and Putaohua oil layer. Figure 3 Among them, the hydrogen-rich region of the deep gas-bearing combination is mainly distributed near the trap-controlling fault, with a depth of 2300-5300m, such as the Shengping-Songzhan uplift area near the Songxi fault zone (the hydrogen content of natural gas in wells Wangshen 1, Songshen 1 and Shengshen 6 is more than 10%, and the highest is 71.64%); the hydrogen content of natural gas in well Fangshen 4 on the Xuxi fault zone is 39.55%, and the hydrogen content of natural gas in well Zhaozhou subsidence 2 is as high as 85.54%; the hydrogen content of natural gas in well Putaohua 1 near the trap-controlling fault zone of Gulong fault subsidence is 38.71%. The hydrogen-rich region of the lower oil and gas-bearing combination is mainly located on the east side of the Sanjiaozai depression, with a depth of about 18000-2700m. The hydrogen content of natural gas in the middle oil and gas-bearing combination is generally low, and only in some wells in the Sanjiaozai depression, Aonan and west slope area, the hydrogen content of natural gas samples is more than 5%, with a depth of 480-1400m. Figure 4

[0077] ​​​​​According to the hydrogen isotope data of the hydrogen obtained from the deep gas-bearing combination of the Songliao Basin, the hydrogen formation mechanism in the deep layer of the Songliao Basin is mainly serpentinization, followed by radioactive decomposition and hydrogen generation and mantle degassing.

[0078] (2) The hydrogen production rate of different source hydrogen systems is calculated through hydrogen source evaluation; the hydrogen generation type and hydrogen system in the Songliao Basin are identified based on the content information of (1) above, and the quality, hydrogen generation rate and scale of each type of hydrogen source are evaluated according to the related parameters in (3) below, and the hydrogen generation rate (G) parameter is deduced, wherein: Figure 1

[0079] The serpentine hydrogen generation amount G serp Calculation: G serp =S serp ×H serp ×F H2 ;

[0080] S serp : The area of iron-rich hydrogen source (m 2 ), about 5×107 m 2 ;

[0081] H serp : The thickness of iron-rich hydrogen source (m), average 200m;

[0082] F serp : Hydrogen generation flux (kg H2 / m 3 ·yr), 2.0×10 -7 kg H2 / m 3 ·yr.

[0083] G serp =S serp ×H serp ×F H2 =2000 m 3 / yr (standard state).

[0084] The radioactive decomposition hydrogen generation amount G rad =S rad ×H rad ×Φ rad ×ω U / Th ×F rad ;

[0085] S rad : The area of uranium-rich hydrogen source (m 2 ), about 1×107 m 2 ;

[0086] H rad : The thickness of uranium-rich hydrogen source (m), average 100m;

[0087] Φ​rad Water content porosity (%), 5%;

[0088] ω U / Th Uranium / thorium content (ppm): 5 ppm;

[0089] F rad Hydrogen generation flux (kg H2 / m³) 3 ·yr), obtained through radioactive hydrogen production simulation experiments and / or analogy (theoretical calculations), approximately 5×10 -8 kg H2 / m 3 ·yr;

[0090] G rad = S rad ×H rad ×Φ rad ×ω U / Th ×F rad = 2.5 m 3 / yr.

[0091] Mantle-source hydrogen quantity G mant The calculation is divided into two models:

[0092] One is the hydrogen isotope (δD) mixing model, calculated using the formula: G mant =(δD crust -δD mixe ) / (δD crust -δD mant )×100%×f surf ;

[0093] However, based on the tectonic setting of the Songliao Basin (rift-depression basin), this embodiment uses a helium isotope mixing model for calculation: helium isotopes ( 3 He / 4 He R) hybrid model, calculation formula: G mant =(R crust -R mix e) / (R) crust -R mant )×r H ;

[0094] R mant Mantle Source 3 He / 4 He endmember value (R / Ra=6-9, atmospheric) 3 He / 4 He ratio Ra = 1.4 × 10 -6 ), 7 Ra;

[0095] R crust Shell Source 3 He / 4He end-member value (R / Ra = 0.01-0.05), 0.02 Ra (shell source end-member);

[0096] R mixe : mixed gas 3 He 4 He measured value, 0.8 Ra (measured mixed value);

[0097] r H : mantle source hydrogen gas (H2 3 He) ratio, about 4.5 x 10 4 m 3 / yr.

[0098] According to the tectonic environment (rift / subduction zone / craton) and the associated other mantle source gas such as CO2 3 He ratio to the calculated mantle source hydrogen amount, the flux correction factor is 0.10.

[0099] G mant = 502 m 3 / yr.

[0100] (3) other key parameters calculation and acquisition. First of all, the acquisition of the loss rate L, according to the closure of Songliao basin, considering the diffusion, microbial consumption, leakage and other factors, L = 2 m 3 / yr. Capture efficiency η trap acquisition, considering the development characteristics of multiple sets of reservoir and cap rock combination in Songliao basin, and hydrogen gas is mainly enriched in deep gas containing combination, controlled by deep fault, η trap =10%. The determination of residence time t res , using the method of geological events, inclusion analysis and isotope dating to determine the residence time of natural hydrogen resources in Songliao basin. For the hydrogen gas from serpentinization, the focus is to determine the geological time (years) of the start and end of serpentinization reaction. Previous studies have shown that serpentinization reaction began in the late Mesozoic, with a duration of about 70~30Ma; radioactive decay has been continuing since the formation of the basement, about 200Ma; the main active period of mantle source hydrogen gas is 140~120Ma. Comprehensive residence time tres = 25 Ma= 2.5 x 10 8 yr.

[0101] (4) according to the different source hydrogen generation rate, the amount of natural hydrogen resources in different hydrogen system of Songliao basin is calculated, and the calculation result is:

[0102] Total hydrogen generation rate: G = G serp + G radi + G mant= 2500 m 3 / yr.

[0103] The total resource amount of natural hydrogen in the sedimentary basin is obtained by accumulation of the resource amount of hydrogen system, and the calculation formula is:

[0104] Q R = Q R1 + Q R2 +…+ Q Rn ≈ 6.2×109 m 3 ;

[0105] These parameters are reasonably deduced and constrained based on the actual geological conditions, tectonic background and hydrogen system characteristics of Songliao Basin.

[0106] Finally, it should be noted that: the above examples are only used to illustrate the technical solutions of the present application, but not to limit it; although the present application has been described in detail with reference to the foregoing examples, those skilled in the art should understand that: it can still modify the technical solutions recorded in the foregoing examples, or make equivalent replacement for part or all of the technical features; and these modifications or replacements do not make the essence of the corresponding technical solutions deviate from the scope of the technical solutions of the embodiments of the present application.

Claims

1. A method for calculating the amount of natural hydrogen resources in a sedimentary basin, characterized by, Comprise: Based on the mass dynamic balance relationship of the generation, migration, accumulation and loss of underground abiogenic hydrogen, the geological resource quantity of natural hydrogen in the target area is calculated by determining the total hydrogen production rate, capture efficiency, hydrogen loss rate and hydrogen residence time of the target area using material balance method; Wherein, the total hydrogen production rate includes hydrogen production rate from at least one source of serpentinization reaction, radioactive radiation and deep mantle degassing; the capture efficiency is the proportion of effectively stored hydrogen among generated hydrogen; the hydrogen loss rate is the annual loss amount of hydrogen in the underground system due to physical, chemical and biological loss mechanisms; the hydrogen residence time is the time for hydrogen to stay and be preserved in the geological body; The sedimentary basin natural hydrogen resource quantity calculation method comprises the following steps: (a) Based on the geological conditions of the sedimentary basin and its surrounding area, the type of natural hydrogen system and its accumulation geological elements are identified; (b) According to the generation mechanism of underground abiogenic hydrogen, the hydrogen production rate of at least one source of serpentinization, radioactive radiation and deep mantle degassing is determined to obtain the total hydrogen production rate; (c) Based on the accumulation ability and preservation conditions of hydrogen during migration, the capture efficiency of hydrogen captured into effective traps is determined; (d) The diffusion, secondary conversion and microbial consumption process of hydrogen in the near-surface environment are analyzed to determine the loss rate of hydrogen system; (e) The residence time is determined by at least one of the following methods: geological events, inclusion analysis and isotope dating; (f) calculating the hydrogen resource amount of the single hydrogen system, and the calculation formula is: Q R ( G × η trap L )× t res ; wherein, Q R is the hydrogen resource amount; G is the total hydrogen production rate; η trap is the capture efficiency; L is the hydrogen loss rate; t res is the hydrogen residence time; (g) The hydrogen resource quantities of several hydrogen systems in the target area are added to obtain the total resource quantity of natural hydrogen in the target area.

2. The method of claim 1, wherein, The hydrogen production rate from the serpentinization source is the hydrogen production rate from serpentinization, and the formula for calculating the hydrogen production rate from serpentinization is: G serp = S serp × H serp × F H2 ; wherein G serp is the serpentinization hydrogen production rate; S serp is the iron-rich hydrogen source body area; H serp is the iron-rich hydrogen source body thickness; F H2 is the hydrogen production flux.

3. The method of claim 1, wherein, The hydrogen production rate from radioactive radiation source is the radioactive radiation hydrogen production rate, and the formula for calculating the radioactive radiation hydrogen production rate is: G rad = S rad × H rad × Φ rad × ω U / Th × F rad ; where G rad is the hydrogen production rate from radiolysis; S rad is the area of the uranium-rich hydrogen source; H rad is the thickness of the uranium-rich hydrogen source; Φ rad is the water-filled porosity; ω U / Th is the uranium / thorium content; F rad is the hydrogen production flux.

4. The method of claim 1, wherein, The hydrogen production rate from deep mantle degassing is the annual contribution of mantle source gas, and the annual contribution of mantle source gas is obtained by an isotope mixing model; the isotope mixing model includes at least one of hydrogen isotope mixing model and helium isotope mixing model.

5. The method of claim 4, wherein, The calculation formula of the hydrogen isotope mixing model is: G mant = (δD crust - δD mixe ) / (δD crust - δD mant ) x 100% x f surf ; where G mant is the mantle-derived gas contribution; δD crust is the crustal δD end-member value; δD mant is the mantle δD end-member value; δD mixe is the measured δD value of the mixed gas; f surf is the natural hydrogen gas surface flux.

6. The method of claim 4, wherein, The calculation formula of the helium isotope mixing model is: G mant = (R crust -R mixe ) / (R crust -R mant ) x r H ; where G mant is the mantle contribution; R mant is the mantle 3 He 4 He end-member value; R crust is the crustal 3 He 4 He end-member value; R mixe is the mixed gas 3 He 4 He measured value; r H is the mantle hydrogen H2 3 He ratio.

7. The method of claim 1, wherein, The annual contribution of mantle source gas is constrained and corrected according to the tectonic environment and the ratio of associated mantle source gas.

8. The method of claim 7, wherein, The tectonic environment includes at least one of a rift, a subduction zone, and a craton; and the associated mantle-derived gas ratio includes CO2 / He. 3 He ratio.

Images